实验一Gauss消元法解方程组.docVIP

实验一 : Gauss消元法解方程组 学号: 姓名: 指导老师：马季骕 专业班级:计算机科学与技术（非师范） 1、算法说明： Gauss消元法的基本做法就是把方程组转化成为一个如下图的等价的三角方程组，这个过程叫做消元。得到三角方程组后，就可以逐个求出Xn，Xn-1，…，X1，这个过程叫回代。 　　　,对计算 (2) 回代过程 建立两个数组a和b，通过循环语句将ｎ阶增广矩阵输入进去，通过对列的循环对每一列进行消去未知数，通过ｎ小步ｎ大步把矩阵化简成上三角形矩阵，最后通过迭代法解得方程组得解。 2、源代码： #includeiostream #includecmath #includeiomanip #include fstream int const N=100; using namespace std; int main() { int n,i,j,k,m,h,t,q; double a[N][N],b[N],x[N],w,v,p; cout请输入矩阵的阶数：endl; while(cinn) //以下为数据输入，并显示所求方程 { cout请输入增广矩阵：endl; for(i=1;i=n;i++) { for(j=1;j=n;j++) { cina[i][j]; } cinb[i]; } k=1; for(i=1;i=n;i++) //消元的第一重循环 { p=0; q=0; for(m=i;mn+1;m++) { if(pabs(a[m][i])) { p=abs(a[m][i]); q=m; } } cout第k个系数矩阵为：endl; //开始矩阵消元的过程 for(m=1;m=n;m++) { for(h=1;h=n;h++) { coutsetiosflags(ios::left)setw(10)a[m][h]; } coutb[m]endl; } for(t=i;t=n-1;t++) //从此处开始时具体的消元算法 { w=a[t+1][i]; for(j=i;j=n;j++) { a[t+1][j]=a[t+1][j]-a[i][j]*(w/a[i][i]); } b[t+1]=b[t+1]-b[i]*(w/a[i][i]); } coutendl; k++; } for(i=1;i100;i++) //从此处开始为回代过程 x[i]=0; for(i=n;i=1;i--) {