清华第五版数值分析第5章ppt课件费下载.pptVIP

三角分解的存在唯一性定理 定理 引理：若 则 存在唯一性证明 证明：存在性：高斯消去法能进行下去，则三角分解存在。 唯一性：设 则 下三角阵=上三角阵，所以 所以 二 矩阵的直接三角分解 直接计算 A 的 LU 分解(例) （续） 1) 2) 3) 一般计算公式 计算量与 Gauss 消去法同. LU 分解求解线性方程组 a11 a12 a13 … a1n u11 u12 u13 … u1n a21 a22 a23 … a2n l21 u22 u23 … u2n a31 a32 a33 … a3n l31 l32 u33 … u3n an1 an2 an3 … ann ln1 ln2 ln3 … unn … … … ... … … … … ... ... (1) (3) (5) (2n-1) (2) (4) (6) (2n) 紧凑格式 ①旧元素减去左边行与 顶上列向量的点积 ②计算行不用除法 ③计算列要除主对角元 注： 例 求矩阵的Doolittle分解 ? ? ? ? 例:用矩阵的直接三角分解法解方程组 或 用 Doolittle 分解法 三 解三对角线性方程组的 追赶法 设三对角方程组 为Ax=d,即: 其中|i-j|1时,aij=0,且满足如下的条件: (1)|b1||c1|0,|bn||an|0 (2)|bi|≥|ai|+|ci|, aici≠0, i=2,3,…,n-1. 可以证明A的Crout分解形式为： L U 由矩阵乘法 用矩阵乘法比较和Crout 分解可推导总结算法步骤如下: 1.分解计算 2.解Ly=d 3.解Ux=y 追的过程 赶的过程 定理 设A为对称正定矩阵，则存在唯一分解A=LDLT，其中L为单位下三角阵，D=diag(d1,d2,…,dn)且di0(i=1,…,n) 四 系数矩阵为对称正定矩阵方程组的平方根法 A对称：AT=A A正定：A的各阶顺序主子式均大于零。即 证明： 本题考查学生综合知识解决问题的能力。最早研制成功的火车、飞机分别是在第一、二次工业革命中。火车动力来源于蒸汽机，飞机动力来源于内燃机 本题考查学生综合知识解决问题的能力。最早研制成功的火车、飞机分别是在第一、二次工业革命中。火车动力来源于蒸汽机，飞机动力来源于内燃机 第五章 线性代数方程组的数值解法 引言 关于线性方程组的数值解法一般有两类。 直接法：经过有限步算术运算，可求得方程组的精确解的方法（若在计算过程中没有舍入误差）. 代表性的算法是高斯(Gauss)消去法。计算代价高. 迭代法：用某种极限过程去逐步逼近线性方程组精确解的方法. 简单实用。 思路 首先将A通过初等行变换化为上三角阵，此过程称为消去过程，再求解如下形状的方程组，此过程称为回代求解。 5.1 高斯消去法 = 上三角方程组的回代求解： 二 高斯消去法步骤如下： 第一步： 则 设 Ax=b. 记A(1)=A b(1)=b 第二步： 则 重复上述过程, 最后得 原方程组的同解方程组 消元过程 Gauss 消去法 上三角方程组的回代求解： 定理：若A的各阶顺序主子式不等于0，则高斯消去法能顺序进行消元，得到唯一解。 高斯消去法条件(高斯消去法适用范围) 证明：高斯消去法能进行下去的充要条件是 而消元过程不会改变顺序主子式 的值。所以 解的唯一性又要求 高斯消去法存在的问题： 如果某个 ,但很小，会引入较大的误差。 1.在高斯消去法消去过程中可能出现 的情况，这时高斯消去法将无法进行 三 高斯列主元消去法 具体步骤为： 基本思想：在每轮消元之前，选列主元素（绝对值最大的元素） 以第二步为例: 挑选从第二行开始的第二列中的最大元素,交换行将其变为主元 如此至多经过n-1步,就得到与之同解的上三角形方 程组的增广矩阵,再用回代过程即可得方程组的解. 高斯列主元消去法的适用范围 定理：系数矩阵非奇异，这高斯列主元消去法可行。 高斯消元法每一消元步骤相当于对系数矩阵左乘一个行初等矩阵，具体为 第一行乘 –m21 加到第二行 ··········· 第一行乘 –