列主元Gauss消去法.docx

列主元Gauss消去法：程序：#include stdio.h#include stdlib.h#include math.h#define M 9#define N (M+1)float num[M][N];//增广矩阵数组float x[M];void Read_Xishu(){ FILE *fp; int i,j; if((fp = fopen(../shuzu.txt,r))==NULL) { printf(打开错误!); exit(1); }for(i=0;iM;i++)//从文件中读取数据for(j=0;jN;j++)fscanf(fp,%f,num[i][j]);printf(该方程组的增广矩阵num[%d][%d]为：\n,M,N); for(i=0;iM;i++)//将数据输出到显示屏{for(j=0;jN;j++)printf(%3.0f ,num[i][j]);printf(\n);}fclose(fp);}void Lie_Zhu_Yuan(){int i,j,k;int Lmax;//列主元所在的行float temp;for(i=0;iM;i++){Lmax=i;temp=(float)fabs(num[i][i]);for(j=i+1;jM;j++){if(temp(float)fabs(num[j][i])){temp=(float)fabs(num[j][i]);Lmax=j;}}if(Lmax!=i){for(k=i;kN;k++){temp=num[Lmax][k];num[Lmax][k]=num[i][k];num[i][k]=temp;}}for(j=i+1;jM;j++)//将系数矩阵化成上三角阵{temp=num[j][i]/num[i][i];for(k=i;kN;k++)num[j][k]=num[j][k]-num[i][k]*temp;}}}void Result(){int i,j;x[M-1]=num[M-1][N-1]/num[M-1][N-2];for(i=M-2;i=0;i--){x[i]=num[i][N-1];for(j=N-2;ji;j--)x[i]=x[i]-num[i][j]*x[j];x[i]=x[i]/num[i][i];}printf(该线性方程组的解为：\n);for(i=0;iM;i++)printf(x[%d]=%f\n,i+1,x[i]);printf(\n);}int main(void){Read_Xishu();//从文件中读取系数Lie_Zhu_Yuan();//列主元，转化成等价的三角方程组Result();//求解 return 0;}运行结果：上机体会：通过此次上机编程我学会了如何从文本中读取数据，嵌套语句的使用，用列主元Gauss消去法解线性方程的方法