[工学]第4章 解线性方程组的直接法.pptVIP

例6 用Doolittle分解求解例5的方程组 解 对增广矩阵作统一处理，得 即 使用回代公式求得 平方根法——设A为n阶对称正定矩阵，则A 可分解为 其中 用直接分解法，可得到 例7 解 改进的平方根法——为了避免开方运算 若A0,则有唯一的分解式 其中 有矩阵乘法规则可得 求得矩阵L、D的元素后，线性方程组就转化为 这是两个三角形方程组，可以用逐步递推法求得其解， 具体计算公式为 解实三对角线方程组的追赶法 其中空白部分均为零元素，并且|b1| |c1|0 方程组的矩阵形式为 利用矩阵的直接三角分解法来推导其计算公式， 将A分解为两个三角阵的乘积 其中 其中 按照矩阵的乘法规则有： 由此可推出计算公式 这样，求解方程组就转化为求解两个三角形方程组 求解公式为 例8 解 §5 误差分析 由于系数矩阵A或者右端常数项b的微小变化引起解的很大变化的方程组，称为病态方程组，系数矩阵称为病态矩阵。反之，称为良态。 矩阵A准确，右端常数项b有误差 ，相应解的改变量为 ， 原方程变为： 如果b准确，而A有误差 ，相应解的误差为 ，则有 因此， 能用来刻画矩阵“病态”特性，称作非 奇异矩阵A的条件数，记作 cond(A)，即： §6 本章小结 高斯消去法 主元高斯消去法 三角分解法（LU） 误差分析 第4章 解线性代数方程组的直接方法 §0 引 言 §1 向量和矩阵的范数 §2 Gauss 消去法 §3 高斯主元素消去法 §4 高斯消去法的变形－三角分解法 §5 误差分析 §6 本章小结 §0 引 言 解线性方程组的两类方法： 直接法: 经过有限次运算后可求得方程组精确解的方法(不计舍入误差) 迭代法：从解的某个近似值出发，通过构造一个无穷序列去逼近精确解的方法。（一般有限步内得不到精确解） §1 向量和矩阵的范数 4.1 向量范数 定义：n维向量 ，对应非负实数 满足：（1）非负性 ，当且仅当 x＝0时，（2）齐次性（3）三角不等式 常用的向量范数， （1）向量的1-范数 （2）向量的2-范数 （3）向量的 -范数 4.2 矩阵范数 与向量范数类似，定义n阶方阵A的范数。 定义： 设A为n阶方阵，对应的非负实数 满足： 常用的矩阵范数： （1）方阵A的行范数（无穷范数） （2）方阵A的列范数（1-范数） （3）方阵A的谱范数（2-范数） （4）算子范数 §2 Gauss 消去法 转化为等价（同解）的三角形方程组 2.1 Gauss 消去法计算过程 从而得到其等价方程组 从而得到其等价方程组 如此计算下去，进行n-1步消元过程后得到的等价方程组 系数矩阵与常数项为： 2.2 消去过程算法 2.3 回代过程算法 消去第一列的 n-1 个系数要计算 n×(n-1) 个乘法。 2.4 Gauss消去法乘法计算量 4.3三角分解法 LU分解次序 以三角矩阵为例给出的记忆公式 1.先解三角形方程组Ly=b,即 2，再解三角形方程组Ux=y,即 例5 Doolittle分解法—把矩阵A分解成单位下三角矩阵 与上三角矩阵的乘积，即 方程组 其计算公式为 其回代过程公式为