第2章 线性方程组的直接法.pptVIP

定理2.2.2:设矩阵A对称正定，则存在唯一的对角元全为正的下三角阵G使得A＝GGT计算格式为:平方根法的优点：(1)乘除法运算量比一般LU分解要小得多；不选主元的平方根法是数值稳定的。缺点：有n个开方运算。第30页，共54页，星期日，2025年，2月5日为避免过多的开方运算，在更多情况下是将A分解为：A=LDLT，其中L为下三角，D为对角阵。事实上，由前面的讨论可知：A=LU，其中L为单位下三角阵，U为上三角阵。记U的元素为uij，D=diag(uii)，由于A对称正定，必有uii≠0(i=1,2,…,n)，所以U=DU*，其中注：将矩阵A作Doolittle分解或Crout分解,由矩阵乘法可得A的LDLT分解。第31页，共54页，星期日，2025年，2月5日?解三对角方程组的追赶法定理2.2.3：若A为对角占优的三对角阵，且满足则方程组有唯一的LU分解。第32页，共54页，星期日，2025年，2月5日第一步:对A作Doolittle分解追赶法公式的推导：（以四阶为例）第33页，共54页，星期日，2025年，2月5日该过程称为“追”的过程。第34页，共54页，星期日，2025年，2月5日该过程称为“赶”的过程。第35页，共54页，星期日，2025年，2月5日一般情形的三对角方程组计算公式：计算次序为：最好牢记第36页，共54页，星期日，2025年，2月5日第37页，共54页，星期日，2025年，2月5日直接比较等式两边的元素，可得到计算公式：第一步:对A作Crout分解：注:也可通过对A作Crout分解进行求解第38页，共54页，星期日，2025年，2月5日第1页，共54页，星期日，2025年，2月5日解线性方程组的两类方法:直接法:经过有限次运算后可求得方程组精确解的方法(不计舍入误差!)迭代法：从解的某个近似值出发，通过构造一个无穷序列去逼近精确解的方法(一般有限步内得不到精确解)第2页，共54页，星期日，2025年，2月5日一、高斯消去法思路首先将方程组Ax=b化为上三角方程组，此过程称为消去过程，再求解上三角方程组，此过程称为回代过程.§2.1高斯消去法和选主元高斯消去法第3页，共54页，星期日，2025年，2月5日将增广矩阵的第i行+li1?第1行，得到：消去过程：第一步：设，计算因子其中第4页，共54页，星期日，2025年，2月5日第k步：设，计算因子且计算共进行n?1步，得到第5页，共54页，星期日，2025年，2月5日定理2.1：若A的所有顺序主子式均不为0，则高斯消去法能顺序进行消元，得到唯一解。回代过程：第6页，共54页，星期日，2025年，2月5日例2.1.1高斯消元法x1=-13x2=8x3=2第7页，共54页，星期日，2025年，2月5日二、选主元消去法为避免这种情况的发生,可通过交换方程的次序，选取绝对值大的元素作主元.基于这种思想导出了主元素法在高斯消去法消去过程中可能出现的情况，这时高斯消去法将无法进行；即使主元素但很小，其作除数，也会导致其它元素数量级的严重增长和舍入误差的扩散第8页，共54页，星期日，2025年，2月5日列主元消去法在第k步消元前，在系数矩阵第k列的对角线以下的元素中找出绝对值最大的元。若p≠k,交换第k个与第p个方程后,再继续消去计算.这种方法称为列主元Gauss消去法。列主元Gauss消去法保证了｜lik｜≤1(i=k+1,k+2,…，n).第9页，共54页，星期日，2025年，2月5日例2.1.2列主元法第一列中绝对值最大是10，取10为主元n阶方程组第k轮消元时，选第k列的后(n–k+1)个元素中绝对值最大作主元。第10页，共54页，星期日，2025年，2月5日x3=6.2/6.2=1x2=(2.5-5x3)/2.5=-1x1=(7+7x2-0x3)/10=0x1=0x2=-1x3=1第二列的后两个数中选出主元2.5第11页，共54页，星期日，2025年，2月5日全主元消去法在第k步消去前，在系数矩阵右下角的n-k+1阶主子阵中，选绝对值最大的元素作为主元素。(1)If