数值分析实验二(列主元Gauss消去法).docVIP

《数值分析》实验报告 实验编号：实验二 课题名称：列主元Gauss消去法 一、算法介绍 1、输入矩阵的阶数n，方程组的增广矩阵A； 2、对k=0,1,…,n-2，循环：选取列中绝对值最大的元素，将主元所在的行的元素保存在数组temp[n+1]中。若主元为零，则系数矩阵奇异，计算停止；否则，顺序进行。如果绝对值最大的元素就在矩阵的对角线上，则进行普通高斯消元法的第一大步，否则将方程组系数换行之后再进行普通高斯消元法的第一大步； 3、然后利用回代法求解线性方程组。 二、程序代码 #includeiostream #includecmath #includeiomanip using namespace std; int main() { int n=0,k=0,i=0,j=0,h=0,g=0,flag=0,i1,j1; double max=0,m=0; cout***利用列主元Gauss消元法求解线性方程组***endl; cout请输入矩阵的阶数：endl; cinn; double a[n][n+1]; double t[n+1]; double x[n]; memset(a,0,sizeof(a)); memset(x,0,sizeof(x)); cout请输入方程组的增广矩阵：endl; for(i=0;in;i++) { for(j=0;jn+1;j++) { cina[i][j]; } } for(k=0;kn-1;k++) { max=0;j1=0; for(i=k;in;i++) { if(fabs(a[i][k])max) { max=fabs(a[i][k]); i1=i; j1=k; } } if(max==0) { cout该系数矩阵为奇异矩阵，计算停止endl; flag=1; break; } else { cout第j1+1列中绝对值最大的元素是a[i1][j1] ，在线性方程组的第i1+1行endl; if(i1!=k) { for(j=0;j=n;j++) { t[j]=a[i1][j]; a[i1][j]=a[k][j]; a[k][j]=t[j]; } } for(i=k+1;i=n-1;i++) { m=a[i][k]/a[k][k]; for(j=k;j=n;j++) a[i][j]=a[i][j]-m*a[k][j]; for(g=0;gn;g++) { for(h=0;hn+1;h++) coutsetiosflags(ios::fixed)setprecision(2)a[g][h] ; coutendl; } coutendl; } } } if(flag==0) { x[n-1]=a[n-1][n]/a[n-1][n-1] ; double sum=0; for(k=n-2;k=0;k--) { sum=0; for(i=n-1;i=k;i--) sum+=a[k][i]*x[i]; x[k]=(a[k