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

实验五 矩阵的lu分解法，雅可比迭代法 班级： 学号： 姓名： 实验五 矩阵的LU分解法，雅可比迭代 一、目的与要求： 熟悉求解线性方程组的有关理论和方法； 会编制列主元消去法、LU 分解法、雅可比及高斯—塞德尔迭代法德程序； 通过实际计算，进一步了解各种方法的优缺点，选择合适的数值方法。 二、实验内容： 会编制列主元消去法、LU 分解法、雅可比及高斯—塞德尔迭代法德程序，进一步了解各种方法的优缺点。 三、程序与实例 列主元高斯消去法 算法：将方程用增广矩阵[A∣b]=(表示 消元过程 对k=1,2,…,n-1 ①选主元，找使得 = ②如果，则矩阵A奇异，程序结束；否则执行③。 ③如果，则交换第k行与第行对应元素位置， j=k,┅,n+1 ④消元，对i=k+1, ┅,n计算 对j=l+1, ┅,n+1计算 回代过程 ①若，则矩阵A奇异，程序结束；否则执行②。 ②；对i=n-1, ┅,2,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; for(int i=start;iend;i++){ if(abs(p[i][start])abs(p[max][start])) max=i; } return max; } void swapRow(double** p,int one,int other,int col){ double temp=0; for(int i=0;icol;i++){ temp=p[one][i]; p[one][i]=p[other][i]; p[other][i]=temp; } } bool dispel(double** p,int row,int col){ for(int i=0;irow;i++){ int flag=findMax(p,i,row); //找列主元行号 if(p[flag][i]==0) return false; swapRow(p,i,flag,col); //交换行 for(int j=i+1;jrow;j++){ double elem=p[j][i]/p[i][i]; //消元因子 p[j][i]=0; for(int k=i+1;kcol;k++){ p[j][k]-=(elem*p[i][k]); } } } return true; } double sumRow(double** p,double* q,int row,int col){ double sum=0; for(int i=0;icol-1;i++){ if(i==row) continue; sum+=(q[i]*p[row][i]); } ret