计算方法分解及列主元高斯消去解线性方程组.doc

北京信息科技大学实验报告 课程名称： 数值分析 实验项目： LU分解及列主元高斯消去解线性方程组 学 院： 光电信息与通信工程学院 专 业： 通信工程 姓名/学号： 赖裕东 2009010476 班 级： 通信0904 实验时间： 2010.12.28 成 绩： 指导教师： 孟 祥 花 实验报告内容 实验目的与要求（实验题目） 用LU分解及列主元高斯消去法解线性方程组 {10 -7 0 1 {X1 {8 -3 2.099999 6 2 X2 5.900001 5 -1 5 -1 X3 = 5 2 1 0 2} X4} 1} 输出AX=b中系数A=LU分解的矩阵L及U，解向量X及detA；列主元法的行交换次序，解向量X及detA；比较两种方法所得的结果。 模型建立（相关主要计算公式） 高斯消去法 三、实验过程、步骤（程序） #includeiostream.h #includemath.h #includeconio.h #includeiomanip.h class glz { public: int row,col,i,j,k,t; double mat [100][100]; double tmp; public: void initialize(); void solve(); }; void glz :: initialize() { cout数值分析课 高斯列主元消去法 endl; cout输入方程组的个数endl; cinrow; cout请输入 未知数个数+1 endl; cincol; cout按行输入各未知数前系数 和 等号右边常数 (按行输入哦~~~)endl; for(int i=0; irow; i++) for(int j=0; jcol; j++) cinmat[i][j]; } void glz::solve() { double x[100]; for(i=0; icol-2; i++) { k=i; double max=mat[k][i]; for(j=i+1; jcol-1; j++) { if(abs(mat[j][i])abs(mat[k][i])) { max=mat[j][i]; k=j; } } if (k!=i) { for(j=i; jcol; j++) { double tmp=mat[k][j]; mat[k][j]=mat[i][j]; mat[i][j]=tmp; } } for (j=i+1; jcol-1; j++) { double tmp=mat[j][i]/mat[i][i]; for(t=i+1; tcol; t++) { mat[j][t]-=tmp*mat[i][t]; } } } x[col-2]=mat[col-2][col-1]/mat[col-2][col-2]; for (i=col-3; i=0; i--) { for(j=col-2; j=i+1; j--) mat[i][col-1]-=x[j]*mat[i][j]; x[i]=