数值计算实验二课案.docVIP

附：报告封面要求 计算方法C （2014-2015-2） 【 解线性方程组 】实验报告 教学班级：__软件4班_____ _2014112218_________ 姓名：______何华均_____ ______________ 实验二 LU分解法的优点 实验目的： 掌握LU分解法的思想和方法。 实验题目： 给定矩阵A与向量b 1、求A的LU分解。 2、利用A的LU分解解下列方程组： ① Ax=b ② A2x=b ③ A3x=b 对第③题分析一下，如果先求M=A3，在解Mx=b，有和缺点？ 3、利用A的LU分解法求A-1，其中n由自己选择，例如取n=5、10。 提示：A-1的第k列为方程组之解，A-1的准确解为： 算法流程图如下： 源程序清单如下： #includecstdio #includecmath #includecassert class LUDec { public: void LU(double **A, int N, int M, double **L, double **U) { double **B = new double *[N]; assert(B != NULL); for(int i = 0; i N; i ++) { B[i] = new double [M]; assert(B[i] != NULL); } // init L U and B for(int i = 0; i N; i ++) { for(int j = 0; j N; j ++) { L[i][j] = 0.0; } for(int j = 0; j M; j ++) { U[i][j] = 0.0; B[i][j] = A[i][j]; } } // please avoid abnormal data here for(int i = 0; i N; i ++) { if( fabs(B[i][i]) 1e-9 ) continue; for(int j = i + 1; j N; j ++) { B[j][i] /= B[i][i]; for(int k = i + 1; k M; k ++) B[j][k] -= B[j][i] * B[i][k]; } } for(int i = 0; i N; i ++) { L[i][i] = 1.0; for(int j = 0; j M; j ++) { if( j i ) L[i][j] = B[i][j]; else U[i][j] = B[i][j]; } } if( B ) { for(int i = 0; i N; i ++) delete[] B[i]; delete[] B; B = NULL; } } }; int main() { LUDec cls; double **A = NULL, **L = NULL, **U = NULL; int N, M; scanf(%d %d, N, M); assert((N 0 N = M)); // A: N x M A = new double *[N]; assert(A != NULL); for(int i = 0; i N; i ++) { A[i] = new double [M]; assert(A[i] != NULL); } // make sure that A can be LU-ed for(i