列主元高斯消去法和列主元三角分解法解线性方程计算方法实验报告1 【课题名称】用列主元高斯消去法和列主元三角分解法解线性方程【目的和意义】高斯消去法是一个古老的求解线性方程组的方法,但由它改进得到的选主元的高斯消去法则是目前计算机上常用的解低阶稠密矩阵方程组的有效方法。用高斯消去法解线性方程组的基本思想时用矩阵行的初等变换将系数矩阵A约化为具有简单形式的矩阵(上三角矩阵、单位矩阵等),而三角形方程组则可以直接回带求解用高斯消去法解线性方程组bAx(其中 A∈Rn×n)的计算量为:乘除法运算步骤为32(1)(1)(21)(1)(1)262233nnnnnnnnnnnMDn,加减运算步骤为(1)(21)(1 )(1 )(1 )(25)6226nnnn nnnnnnAS。相比之下,传统的克莱姆法则则较为繁琐,如求解20阶线性方程组,克莱姆法则大约要195 10 次乘法,而用高斯消去法只需要3060 次乘除法。在高斯消去法运算的过程中,如果出现abs(A(i,i))等于零或过小的情况,则会导致矩阵元素数量级严重增长和舍入误差的扩散,使得最后的计算结果不可靠, 所以目前计算机上常用的解低阶稠密矩阵方程的快速有效的方法时列主元高斯消去法,从而使计算结果更加精确。2、列主元三角分解法高斯消去法的消去过程, 实质上是将 A分解为两个三角矩阵的乘积A=LU,并求解 Ly=b 的过程。回带过程就是求解上三角方程组Ux=y。所以在实际的运算中, 矩阵 L和 U可以直接计算出,而不需要任何中间步骤, 从而在计算过程中将高斯消去法的步骤进行了进一步的简略,大大提高了运算速度,这就是三角分解法采用选主元的方式与列主元高斯消去法一样,也是为了避免除数过小, 从而保证了计算的精确度【计算公式】1、 列主元高斯消去法设有线性方程组 Ax=b,其中设 A为非奇异矩阵。方程组的增广矩阵为1112112122221[ ,]nnlaaaabaaaba bLLMLLMMMMMMMMM第 1 步(k=1):首先在 A的第一列中选取绝对值最大的元素1la ,作为第一步的主元素 : 然后交换( A,b)的第 1 行与第 l 行元素,再进行消元计算。设列主元素消去法已经完成第1 步到第 k-1 步的按列选主元, 交换两行, 消元计算得到与原方程组等价的方程组 A(k)x=b(k) 111max0lii naa(1)(1)(1)(1)(1)1112111(2)(2)(2)(2)22222( )(( )1)()( )( )( )()1,1()(,)( )[, ][,]kkkknkknknkkkkkkkknkkkknkkknnnaaaabaaabaababbaaaA bAbLMLMOMMMLMMM MMMML第 k 步计算如下:对于 k=1,2,⋯, n-1 (1)按列选主元...