三对角阵的LU分解和三对角方程组的求解（C语言）

C语言实现三对角阵LU分解与解三对角方程组

最新推荐文章于 2021-05-20 08:30:22 发布

zhangchao3322218

最新推荐文章于 2021-05-20 08:30:22 发布

阅读量7.3k

点赞数 1

CC 4.0 BY-SA版权

分类专栏： C/C++ 数值分析文章标签： c 语言 float ini 出版算法

本文链接：https://blog.csdn.net/zhangchao3322218/article/details/7412676

本文展示了如何使用C语言实现三对角矩阵的LU分解以及求解三对角方程组的过程。通过输入矩阵的大小N，程序将生成一个三对角矩阵，并进行LU分解，然后利用分解结果求解方程组。程序中包含了计算L和U矩阵、打印L和U矩阵以及求解向量的步骤。

/*三对角阵的LU分解和三对角方程组的求解

-------------A=LU的分解算法-------
参考教材：《数值分析》李乃成，梅立泉，科学出版社
《计算方法教程》第二版凌永祥，陈明逵
*/
#include<stdio.h>
#include<stdlib.h>
#include<math.h>

int main(void)
{
int i,j,n;
int N;

printf("请输入 N(10,20,30或任意值): ");
scanf("%d",&N);
float *a=(float *)malloc(sizeof(float)*N);
float *b=(float *)malloc(sizeof(float)*N);
float *c=(float *)malloc(sizeof(float)*(N-1));
float *d=(float *)malloc(sizeof(float)*N);
float *x=(float *)malloc(sizeof(float)*N);
float *y=(float *)malloc(sizeof(float)*N);
float *u=(float *)malloc(sizeof(float)*N);
float *l=(float *)malloc(sizeof(float)*N);
a[0]=0;
for(n=0;n<N-1;n++)
{
b[n]=2;a[n+1]=1;c[n]=1;
}
b[N-1]=2;
d[0]=1;d[N-1]=-1;
for(i