计算方法与2线性方程组直接法 .ppt

方程组 此方程组的准确解为x1=0, x2=-1。现将其右端加以微小的扰动使之变为： 经计算可得准确解为x1=2, x2=-3. 这两个方程组的解相差很大，说明方程组的解对常数项b的扰动很敏感。 病态方程组： 如果方程组AX=b由于A或b的小扰动而导致解严重失真，则此方程组称为病态方程组，否则称为良态方程组。 判定一个病态方程组的简单方法； 病态方程组一般不能用解方程组的常有方法求解，而采用“迭代求精法”来计算。（列主元消元法的应用） 2.6 LU分解 当A的所有顺序主子式均不为零时，矩阵A可唯一地分解为两个三角矩阵的乘积 A＝LU 其中，L是单位下三角矩阵，U是上三角矩阵 (2.13) 定义 叫 的三角（因子）分解，其中 是 是上三角。 下三角, 为单位下三角阵（对角元全为1）， 为上三角阵，则称 为Doolittle分解; 若 是下三角， 是单位上三角，则称 定理 n阶阵 有唯一Doolittle分解(Crout） 的前n-1个顺序主子式不为0.（证略） 三角分解不唯一,为此引入 定义 若 为Crout分解。 为什么要讨论三角分解？ 若在消元法进行前能实现三角分解A＝LU，则 容易回代求解 (2.14) L是单位下三角矩阵，因此 (2.15) 其中，yi是向量Y的分量，Y＝(y1,y2,…,yn)T,再从UX＝Y中解出X (2.16) 1.直接三角分解法（以Doolittle分解为例） 设 　　 　 ＝ 由矩阵乘法 ……………………… (2.18) (2.17) [解]先LU分解系数矩阵，由2.17式得 再由2.18式得 对第一个方程组，由2.15式得 再由2.16式得 对第二个方程组，由2.15式得 再由2.16式得 2．平方根法 定理 设A对称正定，则有非奇异下三角阵L，使 ---- 理论基础 (证略） 分解方法：设 ( choleskey分解) 3 追赶法 追赶法仍然保持LU分解特性，它是一种特殊的LU分解。充分利用了系数矩阵的特点，而且使之分解更简单，得到对三对角线性方程组的快速解法。 因三对角矩阵的非零元素呈“带状”，我们也因此将它叫做带状矩阵。 三对角线性方程组： 设有方程组Ax=d，其中A为三对角矩阵。 假设系数矩阵A满足条件：对A作Crout分解形式为： 第i个分量 第j个分量 追赶法计算公式 定理 如果上带宽为q，下带宽为p的n阶带状矩阵A有Doolittle分解。A=LU，则L是下带宽为p的单位下三角矩阵，U是上带宽为q的上三角矩阵。 下面举实例用追赶法来解三对角方程组。 第2章 解线性代数方程组的直接法 本章研究的对象是 n 阶线性代数方程组 对象 (2.1) 线性系统广泛存在于工程、科学以及社会科学、商业和经济问题的定量分析等领域中 用矩阵和向量的记法来表示，(2.1)式可写成 (2.2) 其中A＝(aij)是方程组(2.1)的系数aij构成的n×n阶矩阵，称为系数矩阵。B＝{bi},X ＝{xi}是n维向量，X是未知量，B称为右端项。 使方程组(2.1)中每一个方程都成立的一组数x1*,x2*, …,xn* 称为式，(2.1)的解，把它记为向量的形式，称为解向量。 克莱姆(cramer)法则 如果方程组(2.1)的系数矩阵A的行列式不等于零,那么,这个方程组有唯一解,而且它们可以表示为 按上面的等式求解, 就要做 N=(n2-1)n!+n 次乘除法运算,这个计算量是大得惊人的. 例如,当n=10时，乘除法的运算次数共 当n=100时， 1033次/秒的计算机要算10120年; 解线性方程组的方法可以分为2类： ①直接法：在没有舍入误差的情况下，用有限步的四则运算得出精确解的方法。目前常用的是列主元消去法和矩阵三角分解法 ②迭代法：先给一个初始值，按一定法则逐步求解出各个更准确的近似值的方法。 本章讲解直接法 准确，可靠，理论上得到的解是精确的 速度快，但有误差 2.1 消元法 我们知道，下面有3种方程的解我们可以直接求出： ① n次运算 ② （n＋1）n/2次运算 ③ （n＋1）n/2次运算 对方程组（2.1)，作如下的变换，解不变 ①交换两个方程的次序 ②一个方程的两边同时乘以一个非0的数 ③一个方程的两边同时乘以一个非0数，加到另一个方程上 因此，对应的对增广矩阵(A，B)，作如下的变换，解不变 ①交换矩阵的两行 ②某一行乘以一个非0的数 ③某一个乘以一个非0数，加