LU分解法(C语言)

《LU分解法(C语言)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、LU分解法求解线性方程：#includevoidsolve(floatl[][100],floatu[][100],floatb[],floatx[],intn){inti,j;floatt,s1,s2;floaty[100];for(i=1;i<=n;i++)/*第一次回代过程开始*/{s1=0;for(j=1;j=1;i--)/*第二次回代过程开始*/{s2=0;for(j=n;j>i;j--){t=-u

2、[i][j];s2=s2+t*x[j];}x[i]=(y[i]+s2)/u[i][i];}}voidmain(){floata[100][100],l[100][100],u[100][100],x[100],b[100];inti,j,n,r,k;floats1,s2;for(i=1;i<=99;i++)/*将所有的数组置零，同时将L矩阵的对角值设为1*/for(j=1;j<=99;j++){l[i][j]=0,u[i][j]=0;if(j==i)l[i][j]=1;}printf("inputn:");/*输入方程组的个数*/scanf("%d",&

3、n);printf("inputarrayA:");/*读取原矩阵A*/for(i=1;i<=n;i++)for(j=1;j<=n;j++)scanf("%f",&a[i][j]);printf("inputarrayB:");/*读取列矩阵B*/for(i=1;i<=n;i++)scanf("%f",&b[i]);for(r=1;r<=n;r++)/*求解矩阵L和U*/{for(i=r;i<=n;i++){s1=0;for(k=1;k<=r-1;k++)s1=s1+l[r][k]*u[k][i];u[r][i]=a[r][i]-s1;}for(i=

4、r+1;i<=n;i++){s2=0;for(k=1;k<=r-1;k++)s2=s2+l[i][k]*u[k][r];l[i][r]=(a[i][r]-s2)/u[r][r];}}printf("arrayL:");/*输出矩阵L*/for(i=1;i<=n;i++){for(j=1;j<=n;j++)printf("%7.3f",l[i][j]);printf("");}printf("arrayU:");/*输出矩阵U*/for(i=1;i<=n;i++){for(j=1;j<=n;j++)printf("%7.3f",u[i][j]);p

5、rintf("");}solve(l,u,b,x,n);printf("解为:");for(i=1;i<=n;i++)printf("x%d=%f",i,x[i]);}运行结果：inputn:3inputarrayA:223477-245inputarrayB:31-7arrayL:1.0000.0000.0002.0001.0000.000-1.0002.0001.000arrayU:2.0002.0003.0000.0003.0001.0000.0000.0006.000解为:x1=2.000000x2=-2.000000x3=1.00000

6、0Pressanykeytocontinue