现代科学与工程计算-孟大志-第二篇、第6章.pptVIP

6.3.2 列主元的LU分解 在Doolittle分解法中，若 时，计算无法进行，或者 绝对值很小时，按分解公式计算时可能会引起舍入误差的增大， 因此，与Gauss列主元消去法一样，为了保证运算能顺利进行 以及方法的稳定性，矩阵三角分解一般也采用选列主元的技术。 在计算机进行LU分解时，通常采用紧凑格式进行存贮， 即将 存放在 位置上，将 存放在 的位 置上，并将 存放在 的位置上。 紧凑格式的列主元Doolittle分解 第1步：令 ，选主元 ，则将 的第1行与第 行交换，然后执行紧凑格式Doolittle分解的 的第1行和 L的第1列； 第r步： 计算 ，选主元 , 则交换 的第r行与 行元素，然后执行紧凑格式 Doolittle分解的U的第r行和L的第列； 第n步：不需选主元，直接分解即可。 6.3.3 Cholesky分解 设A为对称正定阵，则A的各阶顺序主子式均大于零，因此存在 唯一的Doolittle分解 ，为了利用 A的对称性，将U1再分解为 对角阵D与单位上三角阵U0的乘积，即 ，其中 ， 于是 ，由 A的对称性可知： 为单位上三角阵， 为单位下三角阵，由矩阵 三角分解的唯一性知 故对称正定矩阵A 有唯一的分解式： (6.3.2) 称为 分解，利用矩阵的乘法可得 中元素 的计算公式，对 有 由A的正定性知 ，则对角阵D可以继续分解为 因此(6.3.2)式可改写为 称为 分解（或称Cholesky分解），其中 ，为对角线元素为正的下三角阵， 利用矩阵的乘 法可得矩阵 L中元素 的计算公式 因此，可以对j从1到n逐列计算L的元素，可以求出矩阵L。 分解较 分解的优点在于不必进行开方运算，既适合于 对称正定矩阵，也适合对称且各阶顺序主子式不为零的矩阵。 而 分解只适合于对称正定矩阵。 6.3.4 方程组的矩阵三角分解法 1. 分解法 对于待解的方程组 ，当系数矩阵 实现了LU分解后， 原方程组可写成 ，设 ，则将方程组 变成相继 求解两个三角形方程组 先由下三角方程组 可以顺次求出向量 的所有分量 再由上三角方程组 ，回代求出解向量x的所有分量 LU分解法与高斯消去法的计算量完全相同。 2. 平方根法 在科学与工程计算中，许多实际问题应用有限元法解结构 力学问题可归结为求解线性方程组，其系数矩阵大多是对称正 定的。平方根法就是利用对称正定矩阵的Cholesky分解而得到 的一种求解对称正定方程组的有效方法，特点是分解过程无需 选主元，且具有良好的数值稳定性。 (1) 类似于LU分解法，当系数矩阵 A实现了 分解后,可将方程组 改写成两个三角形方程组 先利用 顺次求出向量y的所有分量，再利用 回代 求出解向量x的所有分量，这种求解方程组的方法称为 分解法（或平方根法）。平方根法和 分解法的计算量都约 为高斯消去法计算量的一半 ，但平方根法的运算量要比 分解法多一点，而且要开平方，因此 分解法也称为改进的 平方根法。 例3 分别用 分解法和 分解法求解下列方程组，计算 过程中保留2位小数。 解：注意到该方程组的系数矩阵是对称的且所有顺序主子式都 大于0，因此可以用 分解法和 分解法求解，记方程组的 系数矩阵为 ，右端项为 。 (1) 分解法：首先对系数矩阵进行 分解 Step3(i=3):得 Step1(i=1):得 Step2(i=2):得 即系数矩阵的 分解为 利用 ，求出y 然后 求出x (2) 分解法：首先对系数矩阵进行 分解 Step1(i=1):得 Step2(i=2):得 Step3(i=3):得 即系数矩阵的 分解为 * * * 第6章 解线性方程组的直接方法 6.1 问题驱动：投入产出分析 投入产出分析是20世纪30年代由美国经济学家首先提出的， 它是研究整个经济系统各部门之间“投入”与“产出”关系的线性模 型，一般称为投入产出模型。国民经济各个部门之间存在着相 互依存的关系，每个部门在运转中将其它部门的成品或半成品 经过加工（称为投入）变为自己的产品（称为产出），如何根 据各部门之间的投入-产出关系，确定各部门的产出水平，以满 足社会的需求，是投入产出综合平衡模型研究的问题，试讨论 如下简化问题。 设国民经济仅由农业、制造业和服务业三个部门构成， 已知某年它们之间的投入和产出关系、外部需求、 初始投入等如表6.1.1所示（数字表示产值，单位 为亿元）。 表6.1.1 国民经济各个部门间的关系 表中