线性方程组地迭代法.docVIP

习题题三：求解Ax=b的迭代法 第 PAGE 9页 共 NUMPAGES 19页 习题题三：求解的迭代法 核心内容： 了解为什么已经有了直接法后，还要研究迭代法； 基本迭代法与近代迭代法； 基本迭代法——不动点迭代法 基本迭代法的构造：分裂 迭代 迭代 迭代 基本迭代法的收敛性与收敛率 引理： 定理：收敛 常用迭代法的收敛性 迭代：，由此判断 迭代：，由此判断 至少有一个特征值为0 特殊情况下，迭代、迭代、迭代的收敛性 1、对角占优矩阵； 2、对称正定矩阵。 6、加速迭代法 7、共轭梯度法 一、用迭代法、迭代法和超松弛迭代法求解，并判断收敛性。 1、，讨论用迭代法、迭代法的收敛性。 解：（1）迭代法 迭代法发散 特征方程可用 原因是， 故。 即 （2）迭代法 迭代法收敛 更方便地 原因是， 故。 迭代法收敛 2、 解：（1）迭代法 特征值，故迭代法收敛。 取 迭代到此为止 能解释这种规律吗？ 因（3重根），故存在非奇异矩阵使得 故 这是因为迭代收敛，故必为常数，此为问题的解。 一般地，，则，即迭代法第步为收敛！（本例迭代3步） （2）迭代法 特征值，于是 迭代发散。 迭代法至少有一个零特征根。 3、 （1）用迭代法和迭代法收敛； （2）用迭代法到为止，取； （3）用迭代法，取同（2）一样的终止准则与初值； 解：（1）为对角占优矩阵，故迭代法和迭代法收敛； （2）迭代格式 ，取 迭代结果为 1 0.6000 2.2727 -1.1000 1.875 2.2727 2 1.0473 1.7159 -0.8052 0.8852 0.9898 3 0.9326 2.0860 -1.1214 1.1309 0.3701 4 … … … … … 4、是否可约？如果是，请求出置换矩阵。 解：欲使，须找出左下角的零子块。已知中的零元素在 把的第2行与第3行对换， 第2列与第3列对换，即可实现 得 ，故是可约的。 6、H为n阶实对称阵，A为n阶实对称正定阵，研究迭代格式（*）如果正定，证明：迭代格式（*）从任一初始向量均收敛. 证明：如欲证（*）收敛，只需证明 设H的任一特征向量为，y为对应特征向量，即. 又正定，故对于y必有 另一方面由的正定性知 故 解，从而. 这样就证明了，从而迭代必收敛 7、矩阵A对称正定，，证明：若2D-A正定，解Jacobi迭代求解Ax=b必收敛. 证明：设求解Ax=b的Jacobi迭代阵，特征值为，相应特征向量为y，即 两边乘D：. 两边再加上，即 故 于是 (*) 由于A正定，故. 从而D也正定; 又2D-A正定，故（*）右端为正. . 这样 则 8、设矩阵A对称正定，证明迭代对任意初始向量都收敛到的解. 证明：改写迭代式为显式 因A对称正定，，故必非奇异，有逆，故 迭代阵 设B的特征值为，对应特征向量为u，故，从而有 为了求，在上式两边与u做内积： 因A正定，，故 于是，即. 所以该迭代收敛，且一定收敛于方程的解. 9、设为对称正定矩阵，对角元素为1. 构作的高斯-塞德尔迭代（SG-S迭代）法如下：分解 证明该迭代收敛. 证明：把SG-S迭代写成不动点迭代的标准形式 把(1)改写为 再把(3)代入(2): 故 欲证SG-S迭代收敛，须证明. L为严格下三角矩阵 又 也即 故 于是 在B的左右两端各乘和，得，即， 注意到是正定对称矩阵，故与B的特征值皆为实数，相应特征向量也u也为实数，故 B的特征向量u可写为： 把B的表达式代入： 即 两边做关于u的内积： 左边： 右边： 故 所以迭代收敛. 10、设矩阵A非奇异，试证用Gauss-Seidel迭代求解必收敛。 证明：设，因为A非奇异，故 即正定 正定阵为系数阵的线性方程组，用G-S迭代必收敛 11、设A为正交阵，B=2I-A, 证明：线性方程组用G-S迭代必收敛. 证明：[分析]如果证明B非奇异，即正定，则G-S迭代收敛. 而证明B非奇异的方法之一，是证明B的特征值均非零. 设为A的特征值，对应特征向量为u. 即： 两边与做内积： 故 如果