数值分析__线性方程组直接解法.pptVIP

第五章 线性方程组的直接解法 第五章 线性方程组的直接解法 §1 Gauss消去法 1.1 顺序Gauss消去法 1.2 列主元Gauss消去法 §2 直接三角分解方法? 2.1 Gauss消去法的矩阵运算 2.2 Doolittle分解法 2.3 平方根法 2.4 追赶法 §1 Gauss消去法 一、顺序Gauss消去法 二 列主元Gauss消去法 §2 直接三角分解方法 一 、 Gauss消去法的矩阵运算 二、Doolittle分解法 例7-3 利用Doolittle三角分解法分解矩阵 例7-4 利用Doolittle三角分解方法解线性方程 例7-5 利用Doolittle三角分解方法解线性方程组 三、平方根法 四、追赶法 五、三角分解方法的优点 第七章 总 结 宝山壁画 宝山壁画是引人注目的昂贵文物。此壁画发现于阿鲁科尔沁旗东沙布乡境内。1994年列为“全国十大考古新发现”之一。宝山壁画中最引人注目的是《杨贵妃教鹦鹉图》。该画高0.7米、宽2.3米，用于笔重彩绘制，最突出的表现了 晚唐风格。唐代擅长绘贵妇仕女的大师周昉绘制了《杨贵妃教鹦鹉图》，不仅享誉中原，而且还影响全国各地。发现于阿旗宝山古墓里的这幅画，就是契丹人聘请中原画家按照周氏风格绘制的， 技法深得周氏画风的真传。在唐人真迹稀如星风的今天，能够从中完整了解唐代人物画的杰出成就，堪称美术史研究的辛事。这幅壁画现今保存在阿鲁科尔沁旗博物馆，历经千年，恍如新绘，是该馆的镇馆之宝。 欢迎大家观看！ 关于方程组 Ax=b , 如果对系数矩阵进行了平方根分解 A＝LLT，则将方程组化为： Ly＝b , LTx＝y 解得 于是，关于系数矩阵是对称正定矩阵的线性方程组Ax=b的求解，分两步进行： 第一步：系数矩阵的平方根分解 第二步：解等价方程组 例7-6 用平方根法求解对称正定方程组 解: 首先进行A 的Cholesky 分解 A＝LLT 2 -0.5 0.5 2 1.5 1 2 -0.5 0.5 2 1.5 1 得 y1＝2，y2＝3.5 ，y3＝1 得 x1＝1，x2＝1，x3＝1 求解Ly＝b: 再求解 LTx＝y: 关于对称正定方程组 也可以用一般的三角分解法求解，这时由 4 -1 1 4 -0.25 0.25 4 3 7 0.75 1 1 求得 x1＝1，x2＝1，x3＝1 追赶法是专门用于求解三对角方程组的。这类方程组经常出现于用差分方法或有限元方法求解二阶常微分方程边值问题、热传导问题及三次样条函数插值等问题，三对角方程组Ax＝b的系数矩阵具有如下形式： 设A为一个三对角矩阵，那么它的顺序主子式均不为零的一个充分条件是： 在此条件下，可对A进行三角分解，设 比较矩阵的对应元素，根据矩阵乘法规则，可得到 i-1列 i 行 第3列 i-1列 i -1行 第3列 i-1列 第3列 i -1行 i-1列 于是，由以上结果： 分别得到： 也就是说,用这一组公式可以对三对角矩阵进行三角分解： 对于三对角方程组Ax＝b，设A的三角分解为A＝LU，则原方程组等价于 Ly＝b，Ux＝y 由Ly＝b，即 解得 解得 由Ux＝y，即 于是，对于三对角矩阵方程组 Ax=b，如下的两组公式便构成了构成了解三对角方程组的追赶法： 例7-7 用追赶法求解三对角方程组 解 : 首先进行系数矩阵的三角分解 求解方程组Ly＝y，即 得 到 y1＝1/2，y2＝1/3，y3＝1/4，y4＝1 再求解方程组 Ux＝y，即 得 到 x1＝1，x2＝1，x3＝1，x4＝1 当三对角矩阵A满足对角占优条件时，追赶法是数值稳定的。 追赶法具有计算程序简单，存贮少，计算量小的优点。 解三对角方程组Ax＝b的追赶法，用四个一维数组存放方程组数据 { ai，bi， ci， di } , di——常数项。 输入 方程组数据 {ai， bi， ci，di，n} 输出 计算解x＝（x1，x2，…，xn）T 根据以下算法编写程序 1． 2 对于i＝1，2，…，n－1，计算 3 4 5 输出 x1，x2，…xn 此时，在完成并存贮矩阵L和U后，右端项第改变一次仅需增加 n2 次运算。 1. 当需要求解具有同系数矩阵的一系列方程组 Ax＝bi，i＝1，2，…，m 时，可大节省计算量。