列主元高斯消去法和列主元三角分解法解线性方程(20200709204233).pdfVIP

实用文案 计算方法实验报告 1 【课题名称】 用列主元高斯消去法和列主元三角分解法解线性方程 【目的和意义】 高斯消去法是一个古老的求解线性方程组的方法， 但由它改进得到的选主元的高斯消去 法则是目前计算机上常用的解低阶稠密矩阵方程组的有效方法。 用高斯消去法解线性方程组的基本思想时用矩阵行的初等变换将系数矩阵 A 约化为具 有简单形式的矩阵（上三角矩阵、单位矩阵等） ，而三角形方程组则可以直接回带求解 用高斯消去法解线性方程组 Ax b （其中A ∈Rn×n ）的计算量为：乘除法运算步骤为 3 n ( n 1) ( n 1) n (2 n 1) n ( n 1) n ( n 1) n 2 n MD n 2 6 2 2 3 3 ，加减运算步骤为 ( n 1) n (2 n 1 ) n ( n 1 ) n ( n 1) ( n 1) n (2 n 5) AS 6 2 2 6 。相比之下，传统的克莱姆 19 法则则较为繁琐， 如求解 20 阶线性方程组， 克莱姆法则大约要 5 10 次乘法， 而用高斯消 去法只需要 3060 次乘除法。 在高斯消去法运算的过程中，如果出现 abs(A(i,i)) 等于零或过小的情况，则会导致矩 阵元素数量级严重增长和舍入误差的扩散， 使得最后的计算结果不可靠， 所以目前计算机上 常用的解低阶稠密矩阵方程的快速有效的方法时列主元高斯消去法， 从而使计算结果更加精 确。 2、列主元三角分解法 高斯消去法的消去过程，实质上是将 A 分解为两个三角矩阵的乘积 A=LU，并求解 Ly=b 的过程。回带过程就是求解上三角方程组 Ux=y。所以在实际的运算中，矩阵 L 和 U 可以直 接计算出， 而不需要任何中间步骤， 从而在计算过程中将高斯消去法的步骤进行了进一步的 简略，大大提高了运算速度，这就是三角分解法 采用选主元的方式与列主元高斯消去法一样， 也是为了避免除数过小， 从而保证了计算的精 确度 【计算公式】 1、 列主元高斯消去法 设有线性方程组 Ax=b，其中设 A 为非奇异矩阵。方程组的增广矩阵为 a a a b 11 12 1n 1 a a a b 21 22 2n 2 [a, b] al 1 a a a b n1 n 2 nn n 第 1 步（ k=1 ）：首先在 A 的第一列中选取绝对值最大的