实验报告2：直接求解线性方程组的解.docVIP

实验报告2：线性方程组求解 Gauss列主元消去法 和 LU三角分解法 解线性方程组。 掌握用Gauss列主元消去法 和 LU三角分解法设计程序，从而实现线性组方程求根1.Gauss列主元消去法是通过Gauss消去，将线性方程组的系数矩阵化为简单的上三角阵，然后进行回代求出x； 2.LU三角分解法的实质是将方程组的系数矩阵化为一个下三角矩阵L和一个上三角阵U，利用 A=LU将Ax=b 化简为LUx=b ,进一步化简为 Ly=b 和 Ux=y;这是一次回代求出y和x； 1.Gauss选主元消去法的核心算法： 首先是进行选主元，每次消元前，通过比较得出每次主子式中的最大项，采用item_i和item_j[n]保存其行标和列标，如果行标和列标与正常的主元不一致，则通过列变换和行变换（被打乱的列顺序将在得到解最后重新恢复），采取主子式的最大项为主元，进行该列的消元。 消元方法的伪代码如下： for(k=0;kn-1;k++) { for(j=k+1;jn;j++) a[k][j] = a[k][j]/a[k][k]; bb[k] = bb[k]/a[k][k]; （此处直接将m[j]的值存入a[k][j]和b[k]中，不仅节省了内存，还方便后面的回代时的累加） for(i=k+1;in;i++) { for(j=k+1;jn;j++) a[i][j] = a[i][j]-a[i][k]*a[k][j]; bb[i] = bb[i] - a[i][k]*bb[k]; } } 然后根据消元所得的上三角阵，回代求值，伪代码大体如下： x[n-1] = bb[n-1]/a[n-1][n-1]; for(i=n-2;i-1;i--) { t = 0; for(j=i+1;jn;j++) t = t+a[i][j]*x[j]; x[i] = bb[i] - t; } （消元过程中的赋值已使a[i][j]实际等于原先的a[i][j]/a[i][i]） 在回代后，方程组的解向量就被自动保存到私有数据成员x[n]中了。 LU三角形分解法的核心算法： （1）实质是将方程组的系数矩阵化为一个下三角矩阵L和一个上三角阵U 中得倒y[n]，代入中得到x[n]暨求得其解。 （2）关于L[n][n]与U[n][n]的各项求法： 对于L的对角线上的元素，设定其为1，暨L[i][i]=1； 对于U的第一行和L的第一列，有： 3.对于U的第k行和L的第k列，有： （3）关于回代求解，基本上是重复Gauss消去法中的回代过程，故不做赘述。 C++语言程序及变量说明 程序建立了一个线性方程组类（Linear_Equations.h），类的数据成员与函数成员如下： class LineEqt { private: int n; //方程组的维数 double** A; //系数矩阵A double* b; //常数向量b double* x; //方程组的解x public: LineEqt(int number); ~LineEqt(); void Input(); //输入方程组 char file_name[]); void Show(); //显示方程组 void ShowRoot();//显示方程组的根 void GetRoot(); void Gauss(char file_name[]); //高斯求解 (将保存结果的文件名作为形参) void LU(char file_name[]); //LU分解求解 (将保存结果的文件名作为形参) void (char file_name[]); //Cholesky求解 (将保存结果的文件名作为形参) void (char file_name[]); //Jacobi迭代求解 (将保存结果的文件名作为形参) void (char file_name[]); //G-S迭代求解 (将保存结果的文件名作为形参) }; 类中的构造、析构函数、Input()、Show()、ShowRoot()与核心算法无关，故不做详述，仅从GetRoot()、Gauss(char [])、void LU(char [])详述： const double ep = 1e-8; const char* Error[] = { Error_0:\n Failed to open file for getting values, Check to make sure to Be correct !, Error_1:\n Failed to open file for outputing, Check to make sure to Be correct