列主元素高斯消去法.docVIP

实验报告 题 目 ： 列主元素高斯消去法 学生姓名 陈玉霞 学 号 20111325015 学 院 信息与控制学院 专 业 信息工程（系统工程方向） 指导教师 殷传洋 2013 年 11 月 18日 一、实验目的 （1）掌握高斯消去法的基本思路和迭代步骤； （2）培养编程与上机调试能力。 二、实现功能 本程序采用GAUSS列主元消去法求解线性方程组。 AX=b. 其中A为N阶矩阵，X，b均为N围列向量。 三、算法描述 高斯消去法基本思路设有方程组Ax=b，设A是可逆矩阵。高斯消去法的基本思想就是将矩阵的初等行变换作用于方程组的增广矩阵B?[A?b]，将其中的A变换成一个上三角矩阵，然后求解这个三角形方程组。 实验内容 X1-0.5*X2+1.5*X3=0.5 X2-0.25*X3=0.5 X3=-6 实验步骤 (1) 写出增广矩阵A; (2) 选主元; (3) 判断是否奇异; (4) 交换对应行元素; (5) 消元; (6) 回代,计算出X[1], X[2], X[3]. 六、代码 # includestdio.h # includemath.h # define delta 1e-6 #define N 100 void main() { int i,j,t,r,n,u,c=0; float p,L,max,s; float X[N]; float a[N][N+1]; printf(请输入方程的阶数\n); scanf(%d,n); printf(输入的原方程系数,中间用空格隔开\n); for(i=0;in;i++) for(j=0;jn+1;j++) scanf(%f,a[i][j]); printf(方程系数为\n); for(i=0;in;i++) { for(j=0;jn+1;j++) { printf(%.2f ,a[i][j]); if(j==n) {printf(\n);} } } for(j=0;jn-1;j++) { { max=fabs(a[j][j]); r=j; } for(i=j+1;in;i++) if(fabs(a[i][j])max) { max=a[i][j]; r=i; } if(fabs(a[i][j])delta) printf(矩阵奇异); for(t=j;tn+1;t++) { } p=a[j][t]; a[j][t]=a[r][t]; a[r][t]=p; } for(i=j+1;iN;i++) { L=a[i][j]/a[j][j]; for(t=j;tn+1;t++) a[i][t]=a[i][t]-L*a[j][t]; } printf(输出原方程的解为：\n); X[n-1]=a[n-1][n]/a[n-1][n-1]; for(i=n-2;i=0;i--) { s=a[i][n]; for(j=i+1;jn;j++) s=s-a[i][j]*X[j]; X[i]=s/a[i][i]; } for(u=0;un;u++) for(j=0;jn+1;j++) { printf(%12f,a[u][j]); c++; if(c%(n+1)==0) printf(\n); } for(i=0;in;i++) { printf(x(%d)=%.4f\n,i+1,X[i]); if(i==n-1) printf(\n); } } 实验结果 八、实验体会 从实验中，我学会了很多，首先要弄懂高斯消去法的原理和计算步骤,使我更加了解Gauss列主元素消去法的实质，掌握了新知识的同时,也为以后的学习打下一个好的基础。