列主元消去法解方程组.docVIP

数值分析上机实验报告（二） 问题描述：利用列主元消去法求解下列方程组 2X1+5X2+3X3 - 2X4 7 2X1- 2X2+3X3+5X4 -1 X1+3X2+2X3+3X4 0 X1+2X2+ X3 - 2X4 4 算法原理： 由高斯消去法知道，在消去过程中可能出现akk k 0的情况，这时候消去法将无法进行，所以最好选取系数矩阵（或消元后的低阶矩阵）中绝对值最大的元素作为主元，以使高斯消去法具有较好的数值稳定性。 实验步骤： det 1; 对于k 1,2,···,n-1 按列选主元 |aik.k| max|aik| （2）如果ai.k 0,则计算停止（det（A） 0） （3）如果ik k则转（4） 换行：akj aik·j j k,k+1,···,n bk bik det -det （4）消元计算 对于i k+1,···，n amik aik/akk 对于j k+1,···,n aij aij—mik*akj bi bi-mik*bik 5 det akk*det 如果则计算停止（det A 0） 回代求解 （1）bn bn/ann （2）对于i n-1···，2,1 bi bi-∑aij*bj /aii 5.det ann*det 四、实验框图 五、源程序 # include # include # d