2022年高斯列主消元法法和分解法 .docx

《2022年高斯列主消元法法和分解法 .docx》由会员分享，可在线阅读，更多相关《2022年高斯列主消元法法和分解法 .docx（7页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精品_精品资料_高斯列主消元法： #include #include #define N 10int mainint i=0,j=0,r=0,n=0,k=1;float aNN,xN,max= -1,sum1=0,sum2=0;/ 输入矩阵维数printfPlease enter n 0n11:n; scanf%d,&n;/ 输入矩阵元素,最终一列为bfori=1;i=n;i+forj=1;j=n+1;j+printfPlease enter a%d%d:n,i,j;scanf%f,&aij;dofori=k;imaxmax=aik; r=i;ifark=0可编辑资料 - - - 欢迎下载精品

2、_精品资料_/ 换行printfCannot solve.n; return 0;可编辑资料 - - - 欢迎下载精品_精品资料_forj=1;j=n+1;j+sum1=akj; akj=arj; arj=sum1;/ 将矩阵化为上三角矩阵fori=k+1;i=n;i+可编辑资料 - - - 欢迎下载精品_精品资料_k=k+1;forj=k+1;j=n+1;j+aij=aij -aik*akj/akk;可编辑资料 - - - 欢迎下载精品_精品资料_max=-1;whilek=1;k -forj=k+1;j=n;j+ sum2=sum2+akj*xj;xk=akn+1 -sum2/akk; s

3、um2=0;/ 输出结果fori=1;i=n;i+printfx%d=%fn,i,xi;while1; return 0;可编辑资料 - - - 欢迎下载精品_精品资料_LU 分解法： #include #include #define N 10 / 矩阵大小范畴/* 使用已经求出的x,向前运算 x供 getx 调用* float a矩阵 U* float x方程组解* inti解的序号数组X 元素序号* intn矩阵大小* return公式中需要的和*/float getmxfloat aNN, float xN, int i, int nfloat mx = 0; int r;forr=i

4、+1; rn; r+mx += air * xr;return mx;/* 使用已经求出的 y,向前运算 y供 gety 调用* float a矩阵 L* float y数组 Y* inti数组 Y 元素序号* intn矩阵大小* return公式中需要的和*/float getmyfloat aNN, float yN, int i, int nfloat my = 0; int r;forr=0; rNprintf 矩阵过大； n; continue;flag=0;printf 请输入系数矩阵值： n; fori=0; in; i+forj=0; jn; j+printfa%d%d: ,

5、i, j;scanf%f, &aij;printf 请输入右端项数组： n; fori=0; in; i+printfb%d: , i;scanf%f, &bi;/* 显示原始矩阵 */ printfn原始矩阵： n; fori=0; in; i+forj=0; jn; j+printf%0.3f ,aij; printfn;printfnn;/* 初始化矩阵 l*/可编辑资料 - - - 欢迎下载精品_精品资料_fori=0; in; i+forj=0; jn; j+ifi=j lij = 1;/* 开头 LU 分解 */* 第一步：对矩阵 U 的首行进行运算 */ fori=0; in;

6、i+u0i = floata0i/l00;/* 其次步：逐步进行LU 分解 */ fori=0; in -1; i+/* 对“ L 列”进行运算 */ forj=i+1; jn; j+fork=0,sum=0; kn; k+ifk .= i sum += ljk*uki;lji = floataji-sum/uii;/* 对“ U 行”进行运算 */ forj=i+1; jn; j+fork=0,sum=0; kn; k+ifk .= i+1 sum += li+1k*ukj;ui+1j = floatai+1j-sum;/* 输出矩阵 l*/ printf 矩阵 L： n; fori=0;

7、in; i+forj=0; jn; j+printf%0.3f , lij;可编辑资料 - - - 欢迎下载精品_精品资料_printfn;/* 输出矩阵 u*/ printfn矩阵 U：n; fori=0; in; i+forj=0; jn; j+printf%0.3f , uij;printfn;/* 回代方式运算数组Y*/ fori=0; in; i+yi = getyl,b,y,i,n;/* 显示数组 Y*/ printfnn数组 Y：n; fori=0; i=0; i-xi = getxu,y,x,i,n;/* 显示数组 X*/ printfnn数组 X： n; fori=0; in; i+printfx%d = %0.3fn, i+1,xi;return 0;可编辑资料 - - - 欢迎下载