第6章 解线性方程组的直接方法.ppt

　　由于在计算机实现时当uri计算好后ari就不用了，因此计算好L,U的元素后就存放在A的相应位置.例如 最后在存放A的数组中得到L,U的元素. 　　由直接三角分解计算公式，需要计算形如∑aibi的式子，可采用“双精度累加”，以提高精度. 　　如果已经实现了A=LU的分解计算，且L,U保存在A的相应位置，则用直接三角分解法解具有相同系数的方程组Ax=(b1,b2,…,bm)是相当方便的，每解一个方程组Ax=bj仅需要增加次乘除法运算. 　　矩阵A的分解公式(3.2),(3.3)又称为杜里特尔(Doolittle)分解. 　　直接分解法大约需要n3/3次乘除法，和高斯消去法计算量基本相同. 2. 选主元的三角分解法 　　从直接三角分解公式了看出当urr=0时计算将中断，或者当urr绝对值很小时，按分解公式直接计算可能引起舍入误差的积累. 但如果A非奇异，我们可通过交换A的行实现矩阵PA的LU分解，因此可采用与列主元消去法类似的方法(可以证明下述方法与列主元消去法等价)，将直接三角分解法修改为(部分)选主元的三角分解法. 　　设第r-1步分解已完成，这时有 　　设r步分解需用到(3.2)式及(3.3)式，为了避免用小的数urr作除数，引进量 于是有 　　取 交换A的r行与ir行元素，将　调到(r, r)位置(将(i, j)位置的新元素仍记为lij及aij ), 于是有