实验四 线性方程组数值解法.docxVIP

实验四 线性方程组数值解法 一、目的与要求 熟悉求解线性方程组的有关理论和方法； 会编制列主元消去法、LU分解法、雅可比及高斯-塞德尔迭代法的程序； 通过实际计算，进一步了解各种方法的优缺点，选择合适的数值方法。 二、算法描述 1、列主元高斯消去法 算法 将方程用增广矩阵[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,计算 2、矩阵直接三角分解法 将方程组中的A分解为，其中L为单位下三角矩阵，U为上三角矩阵，则方程组化解为2个方程组，，具体算法如下 3、雅可比迭代法 设方程组的系数矩阵的对角元素，M为迭代次数容许的最大值， 3、 4、高斯-塞尔德迭代法 设方程组Ax=b的系数矩阵的对角线元素,,M为迭代次数容许的最大值,ε为容许误差 ①取初始向量,令k=0。 ②对i=1,2,…,n计算 ③如果，则输出,结束；否则执行④。 ④如果则不收敛，终止程序；否则,转②。 三、源代码 #includestdio.h #includemath.h #includestdlib.h #define e 0.001 //列主元高斯消去法 void Gaussian(); //列主元消去法的排序函数 int MAX(int n,int m,int l,double b[10][10]); //矩阵直接三角分解法 void Matrix(); //雅可比迭代法 void Jacobi(); //高斯-塞德尔迭代法 void Serdel(); int main() { char ch1,ch2; do{ printf(***********************************************************\n); printf(********************请选择矩阵求解算法:********************\n); printf(********************1.列主元高斯消去法*********************\n); printf(********************2.直接矩阵三角分解法*******************\n); printf(********************3.雅可比迭代法*************************\n); printf(********************4.高斯-塞德尔迭代法********************\n); printf(***********************************************************\n); printf(请选择:); fflush(stdin); ch1 = getchar(); if(ch14||ch11){ printf(选择错误，请重新选择!!!!\n); system(Pause); } else{ switch(ch1){ case 1: Gaussian(); break; case 2: Matrix(); break; case 3: Jacobi(); break; case 4: Serdel(); break; } } printf(***********************************************************\n); printf(请选择是否继续该程序(Y/N):); fflush(stdin); ch2 = getchar(); }while(ch2 == Y); return 0; } //寻找数组每第m列最大值元素的位置，并返回 int MAX(int n,int m,int l,double b[10][10]){ double t = fabs(b[l][m]); int i; for(i = l;in;i++){ if(tfabs(b[i][m])) t = fabs(b[i][m]); } for(i = l;in;i++){ if(fabs(b[i][m]) == t){ break; } } return i; } //列主元高斯消去法 void Gaussian