线性代数方程组的直接法材料.ppt

1.消去法解方程组的基本思想： 逐次消去未知数 把Ax=b,化为与之等价的三角方程组 现考虑n阶线性方程组 7.2.2矩阵的三角分解 从矩阵变换的观点看： 从矩阵变换的观点看： 从矩阵变换的观点看：(一般有) 作业 2，4 LU 分解求解线性方程组 矩阵A的谱半径 2. 基本原理 原理1 (对称阵的三角分解定理Th7.7) 设A为n阶对称阵，且A的所有顺序主子式均 不为零，则A可唯一分解为 A=LDLT 其中L为单位下三角阵，D为对角阵. 证： A所有顺序主子式不为零,则A有唯一的LU分解。 将U再分解为 其中D为对角阵，U0为单位上三角阵，于是 又 由分解的唯一性得 即 原理2 (对称正定阵的三角分解定理Th7.8) 设A为n阶对称阵，且A的所有顺序主子式均 大于零，则存在一个非奇异下三角阵 L 使 A=LLT，当限定L的对角元素为正时，这种 分解是唯一的。 矩阵的这种分解称为Cholesky分解 乔里斯基 证明： 由原理1可知A可分解为 则 由于对于任意非零向量x， y=(LT)-1x也为非零向量，于是 由A的正定性 即D正定，所以D的对角元素均为正数， 令 则 原理2 (对称正定阵的三角分解定理Th7.8) 设A为n阶对称阵，且A的所有顺序主子式均 大于零，则存在一个非奇异下三角阵 L 使 A=LLT，当限定L的对角元素为正时，这种 分解是唯一的。 矩阵的这种分解称为Cholesky分解 乔里斯基 3 平方根法思路 (与直接三角分解法解线性方程组类似推导公式) 思路： 比较A与LLT的相应元素，可得计算公式 分解唯一Th7.8 设n阶对称正定矩阵A有分解 ,先用待定系数法求L的元素 。 4.用平方根法解线性代数方程组的算法 （1）对矩阵A进行Cholesky分解，即A=LLT，由矩阵乘法： j=1 j=2,3,…,n （2）求解Ly=b （3）求解LTx = y Ax = b (LLT)x=b Ly = b LTx = y 用平方根法求解对称正定方程组时 由 可知 因此 2.平方根法是数值稳定的 5. 平方根法稳定性和计算量 平方根法计算量:n3/6 L求出,LT也求出来了,所以是LU分解的计算量的一半。 1.不需选取主元 三、改进的平方根法 1. 引入 平方根法lii要用到开方，为了避免开根方，利用Th7.7的分解式A=LDLT,这就是改进的平方根法。 优点： 1、数值稳定。 2、计算量小，大约为 次乘除法，是一般矩阵 A的LU分解计算量的一半。 缺点： 计算lii时要开平方。 6. 平方根法优缺点： 设n阶对称正定矩阵A有分解 。其中L为单位下三角阵， D为对角阵。即 2．改进平方根法 （为避免开平方） =1 求解对称正定方程组Ax=b的改进平方根法(计算公式)： 1．分解计算 令为tij 优点： 1、计算量小，同平方根方法一样，是一般矩阵A的LU分解(消元法)计算量的一半，是目前解对称正定矩阵方程组的有效方法。 3、精度较高。 2、计算简单(没有开方运算)。 例 用改进的平方根法求解方程组 分析：对称很显然,由TH7.7可LDLT分解 解 1、分解 2、求解计算 到方程组的解： 作业： 习题 14 三、追赶法 1. 初步介绍 追赶法求解系数矩阵为三对角阵的方程组。 将系数矩阵A分解为两个三角阵的乘积，A=LU，其中L为下三角矩阵，U为单位上三角矩阵。 对角占优矩阵: 三对角线方程组,其形式为: 其中 2.追赶法 设A为满足对角占优条件的n阶三对角阵 ， 理论基础： 则有唯一三角分解A=LU 3. 追赶法的本质 4. 追赶法的计算量和数值稳定性 追赶法的计算量比较小 追赶法是数值稳定的 就是把Gauss消去法用到三对角方程组上 消元过程称为追。回代过程称为赶 一、向量的范数 §7.5 向量和矩阵的范数 二、矩阵的范数 数量积(内积) 称向量的数量积 Euclid范数 Euclid范数P52页 一、向量的范数 1. 向量范数的定义(定义7.2) 设对任意向量 x∈Rn，按一定的规则有一实数与之对应，记为‖x‖，若‖x‖满足 则称‖x‖为向量的范数 在Rn上的向量x =(x1,…,xn)T∈Rn，三种常用的范数为： 称为∞－范数或最大范数 称为1－范数 称为2－范数 称为p－范数 2. 常用的向量范数 举例：计算向量 x=(1, -2, 3)T的各种范数. 解： 设给定Rn中的向量序列{ }，即 其中 若对任何i (i = 1, 2,…, n )都有 则向量 称为向量序列{ }的极限， 或者