《高斯列主元消元法》由会员分享，可在线阅读，更多相关《高斯列主元消元法（5页珍藏版）》请在装配图网上搜索。
1、1. 用高斯列主元消去法解线性方程组 Ax=b。0.3 x 101559.143159.17、A =5.291-6.130-12,b =46.7811.29521 1211丿 2丿编程思路：在消元之前，在剩余的元素中选取绝对值最大的为主元，并将其换到主元的 位置，以防止出现小数被大数吃掉的现象。然后再将矩阵化为三角矩阵，按照高斯消去 法进行计算，最后进行回代求解。程序代码：#include #include void main()/ 定义变量/float Ab45,x;int i,j,m,k;/输入增广矩阵 Ab/printf(please input coefficient matrix A
2、b by linen);for(i=0;i4;i+)for(j=0;j5;j+)scanf(%f,&Abij);/ 输出已输入矩阵，验证是否有误/for(i=0;i4;i+)for(j=0;j5;j+)printf(%8.4f ,Abij);printf(n);printf(nnnnnn);/消元过程/ for(k=0;k=3;k+) /寻找列主元行号/ m=k;for(i=k+1;ix)m=i;printf(%dnnnn,m);if(Abmk=0)printf(error);return;else/交换两行/ for(j=k;j=4;j+) x=Abkj;Abkj=Abmj;Abmj=x;f
3、or(i=0;i4;i+)for(j=0;j5;j+)printf(%8.4f ,Abij); printf(n);printf(nnnnnn);for(i=k+1;i=3;i+)Abik=Abik/Abkk;for(j=k+1;j=4;j+)Abij=Abij-Abik*Abkj;/输出变换后的矩阵/ for(i=0;i4;i+)for(j=0;j=0;i-)x=0;for(j=i+1;j4;j+) x=x+Abij*Abj4;Abi4=(Abi4-x)/Abii;/输出结果/for(k=0;k=3;k+)printf(%f ,Abk4);数据输入运行结果2用追赶法解正定矩阵 Ax=b。-1-130-20、0r 6iA =，b =0-24-3-2 00-35丿 1丿编程思路：将各阶顺序主子式都不为0的矩阵A,分解为一个非奇异下三角阵和一个单位上三角阵, 将元方程Ax=b化为Ly=d和Ux=y,先求出y,在求出X。运行结果
高斯列主元消元法
