您好,欢迎访问三七文档
当前位置:首页 > 行业资料 > 酒店餐饮 > LU分解法(C语言)
LU分解法求解线性方程:#includestdio.hvoidsolve(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;ji;j++){t=-l[i][j];s1=s1+t*y[j];}y[i]=(b[i]+s1)/l[i][i];}for(i=n;i=1;i--)/*第二次回代过程开始*/{s2=0;for(j=n;ji;j--){t=-u[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:\n);/*输入方程组的个数*/scanf(%d,&n);printf(inputarrayA:\n);/*读取原矩阵A*/for(i=1;i=n;i++)for(j=1;j=n;j++)scanf(%f,&a[i][j]);printf(inputarrayB:\n);/*读取列矩阵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=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:\n);/*输出矩阵L*/for(i=1;i=n;i++){for(j=1;j=n;j++)printf(%7.3f,l[i][j]);printf(\n);}printf(arrayU:\n);/*输出矩阵U*/for(i=1;i=n;i++){for(j=1;j=n;j++)printf(%7.3f,u[i][j]);printf(\n);}solve(l,u,b,x,n);printf(解为:\n);for(i=1;i=n;i++)printf(x%d=%f\n,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.000000Pressanykeytocontinue
本文标题:LU分解法(C语言)
链接地址:https://www.777doc.com/doc-4607324 .html