雅可比迭代法与矩阵的特征值.pdfVIP

实验五 矩阵的lu分解法，雅可比迭代法 班级： 学号： 姓名： 实验五 矩阵的 LU 分解法，雅可比迭代 一、目的与要求： 熟悉求解线性方程组的有关理论和方法； 会编制列主元消去法、LU 分解法、雅可比及高斯—塞德尔迭代法德程 序； 通过实际计算，进一步了解各种方法的优缺点，选择合适的数值方法。 二、实验内容： 会编制列主元消去法、LU 分解法、雅可比及高斯—塞德尔迭代法德程 序，进一步了解各种方法的优缺点。 三、程序与实例 列主元高斯消去法 算法：将方程用增广矩阵 [A ∣b]=(a ) 表示 ij n(n1) 1) 消元过程 对 k=1,2,…,n-1 ①选主元，找  使得   i k,k 1, ,n k max = a a i ,k kin ik k ②如果a 0，则矩阵A 奇异，程序结束；否则执行③。 i ,k k ③如果i k，则交换第 k 行与第 行对应元素位置，i k k a a j=k,┅,n+1 kj i j k ④消元，对 i=k+1, ┅,n 计算 l a /a ik ik kk 对 j=l+1, ┅,n+1 计算 a a l a ij ij ik kj 2) 回代过程 ①若a 0，则矩阵A 奇异，程序结束；否则执行②。 nn ②x a /a ；对 i=n-1, ┅,2,1,计算 n n,n1 nn  n     x  a  a x /a  i  i,n1 ij j  ii ji1 程序与实例 程序设计如下： #include iostream #include cmath using namespace std; void disp(double** p,int row,int col){ for(int i=0;irow;i++){ for(int j=0;jcol;j++) coutp[i][j] ; coutendl; } } void disp(double* q,int n){ cout=====================================endl; for(int i=0;in;i++) coutX[i+1]=q[i]endl; cout=====================================endl; } void input(double** p,int row,int col){ for(int i=0;irow;i++){ cout输入第i+1行:; for(int j=0;jcol;j++) cinp[i][j]; } } int findMax(double** p,int start,int end){ int max=start; f