第3章 线性代数方程组的数值解法.ppt

3.1 引言 在这一章，我们将建立解线性方程组的一般方法，然后分析与计算机解题有关的误差并研究控制和降低这些误差的方法，最后介绍解线性方程组的迭代法。 本章讨论形如式(3.1.1) 的线性方程组的数值解法。 这是一个有n个方程， n个未知数x1, x2, ?, xn的系统，元素aij和bi是实数。可表示为矩阵形式 记 x=(x1， x2， ?， xn)T b=(b1， b2， ?， bn)T 为n 阶方阵 方程成为简单形式 Ax=b (3.1.2) 解线性方程组的方法，在线性代数中有根据Gram法则的方法。若A为非奇异矩阵，即|A|≠0，则xi=?i/?，方程组有唯一解。其中，?为A的行列式，?i为?中第i列元素被b替代的行列式。 根据行列式的定义计算一个n阶行列式需要n!(n?1)次乘法，因此，就目前而言，当n稍大一些，在电脑上用Gram法则解线性方程组是不可行的。 回顾中学课本中解一次方程组的方法：加减消元法。 例 3.2.1 r3?2r1：-3x1=-3 ?x1=1； r2?r3：x2=1 行初等变换 x1=1， x2=1 代入r1： x3=4?2?1=1 这里由于人可以方便的判断出r3?2r1和r2?r3可以同时消去2个未知数，很容易的得到了解。 3.2 Gauss消去法 计算机不能这样灵活的加减消元，但能按照既定的规则按部就班地进行加减消元。首先，分别以第1行乘以适当的数加到2~n行消去第1列除a11之外的所有元素，然后，第2行分别乘以适当的数加到3~n行消去第2列a22以下的所有元素，??。这就是顺序Gauss消去法的消元过程。 本节重点讨论Gauss消去法的基本原理、与矩阵三角分解的关系、计算量和可行性、稳定性。 3.2.1 顺序消去过程和LU三角分解 将变换之前的A、b记为A(1)、b(1)，则方程组(3.1.2)可写为A(1)x = b(1)，增广矩阵为[A，b]= [A(1)，b(1)]。 第1步消元： 丛r2,r3,? ,rn中消去x1项，条件a11(1)≠0，使得 A(2)x = b(2)?A(1)x = b(1) 其中 具体方法： -(r1(1)/a11(1))a21(1)加到第2行， -(r1(1)/a11(1))a31(1)加到第3行，?， -(r1(1)/a11(1))ai1(1)加到第i行，?， -(r1(1)/a11(1)) an1(1) 加到第n行。这就使第1列a11(1)以下元素皆为0。 记 l21 = a21(1)/a11(1)；?； li1 = ai1(1)/a11(1)，i=2,3,?,n 则上述过程可描述为：以矩阵[A(1)，b(1)]中的第1行分别乘以 -li1加到2~n行，即 ri(1)?li1r1(1) ?ri(2)，i=2,3, ?,n 其中第i行 aij(2) = aij(1)? li1a1j(1)，i, j=2,3,?,n bi(2) = bi (1)? li1b1 (1)，i=2,3,?,n 第2步消元： 丛r3,r4,? ,rn中消去x2项，条件a22(2)≠0，使得 A(3)x = b(3)?A(2)x = b(2) 具体方法： li2 = ai2(2)/a22(2)， ri(2)?li2r2(2) ?ri(3)， i=3,4,?,n ?? 第k步消元： 丛rk+1,rk+2,? ,rn中消去xk项，条件akk(k)≠0，使得 A(k+1)x = b(k+1)?A(k)x = b(k) 其中 (3.2.1) 具体方法： 记 lik = aik(k)/akk(k)，i=k+1,k+2,?,n (3.2.2) 以矩阵[A(k)，b(k)]中的第k行乘以-lik加到i行，即 ri(k)?likrk(k) ?ri(k+1)，i= k+1,k+2, ?,n 其中第i行 aij(k+1) = aij(k)? likakj(k)，i, j= k+1,k+2,?,n (3.2.3) bi(k+1) = bi (k)? likbk (k)，i= k+1,k+2,?,n 当完成第k=n?1步时， A(1)变为上三角阵A(n) ，Gauss消元过程结束，得到与原方程组等价的方程组 A(n) x= b(n) (3.2.4) 其中 (3.2.5) 这个过程称为Gauss消去法的回代过程。写成简洁的形式 (3.2.6) 消元过程的矩阵描述 分别以-li1乘以矩阵[A(1)，b