第四讲线性方程组的求解.ppt

引言 举例 2.1.2 高斯消去法的运算量 高斯消去法比起克莱姆法则和约当消去法，主要长处就是运算量较少。 用克莱姆法则求解n阶线性方程组，需计算n+1个行列式；每个行列式是n!个乘积之和，而每个乘积是n个元素相乘，因此共需乘法次数 ( n + 1）· n!( n － 1)= ( n + 1) ! ( n － 1) 当n=20时( n + 1) ! ( n － 1)≈9.7×1020。要完成这么多次乘法，在每秒做一亿次乘法运算的计算机上，也需30.8万年。因此克莱姆法则在实际计算中不适用。 高斯消去法消去过程分n-1步，第k步变换n－k行：对这些行先求倍数，再从n+1－k元素减去第k行对应列的倍数；因此共需乘、除次数 回代过程求xn需1次除法，求xn-1需1次乘法、1次除法，……，求需n –1次乘法1次除法，因此共需乘除次数 两过程共需乘除次数 。 当n=20时，N=3 060，显然比克莱姆法则少得多。 线性代数教材讲授线性方程组解法时，往往提倡用行初等变换把增广矩阵化为最简形。如上例方程组，多用如下变换： ［例］ 求下面矩阵的LU分解 解： 先求二阶主子矩阵的LU分解 u11=a11=2，u12=a12=-1，l21=a21／u11=0， u22=a22-l21u12=-4-0=-4 本题i=3，由公式 得 1.追赶法 追赶法与稀疏线性方程组 追赶法仍然保持LU分解特性,它是一种特殊的LU分解。充分利用了系数矩阵的特点，而且使之分解更简单,得到对三对角线性方程组的快速解法。 因三对角矩阵的非零元素呈“带状”，我们也因此将它叫做带状矩阵。 定理3：若对称矩阵A各阶顺序主子式不为零时, 则A可以唯一分解为A= LDLT ，这里 （2）求解下三角形方程组 （3）求解LTX = y 3．改进平方根法 其中 改进平方根法解对称正定方程组的算法 令LTX = y，先解下三角形方程组LDY = b得 解上三角形方程组LTX = Y得 ? 特殊方程组的解法 1.追赶法 2.LDLT分解法 三对角线性方程组： 设有方程组Ax=d，其中A为三对角矩阵。 假设系数矩阵A满足条件：对A作Crout分解形式为： 第i个分量 第j个分量 追赶法计算公式 “追” “赶” 定理 如果上带宽为q，下带宽为p的n阶带状矩阵A有Doolittle分解。A=LU，则L是下带宽为p的单位下三角矩阵，U是上带宽为q的上三角矩阵。 下面举实例用追赶法来解三对角方程组。 实际问题中，当求解方程组的系数矩阵是对称矩阵时，则用下面介绍的LDLT 分解法可以简化程序设计并减少计算量. 从定理可知,当矩阵A的各阶顺序主子式不为零时,A有唯一的Doolittle分解A=LU.此时,当然有,所以矩阵U的对角线元素uii ?0,(i=1,2,?,n),将矩阵U的每行依此提出uii 2. LDLT分解法 由A=AT，得 由分解的唯一性有， 即，于是可得下面的结论。 LT为L的转置矩阵。 当A有LDLT分解时，利用矩阵运算法则及相等原理易得计算ljk及dk的公式为 例4 用列主元消去法解方程组 解 第一次消元对 因列主元素为a31(1),故先作行变换r2__r3,然后进行消元计算可得 第二次消元对[A(2) |b(2) ],因列主元素为a32(2) ,故先作行变换r2__r3,然后进行消元计算可得 由此回代,得x=(1.9272,-0.69841,0.90038)T与精确解 x=(1.9273,-0.69850,0.90042)T相比较是比较准确的. -0.002x1+2x2+2x3 =0.4 x1+0.78125x2 =1.3816 3.996x1+5.5625x2+4x3=7.4178 -0.002 2 2 0.4 [A(1) |b(1) ]= 1 0.78125 0 1.3816 3.996 5.5625 4 7.4178 3.9