- 149
- 0
- 约2.35万字
- 约 99页
- 2017-07-21 发布于湖北
- 举报
* */97 五、特征值与特征向量 为使消元后的对角元为正,用x1构造镜面反射阵H1 ,使得H1 x1 = y1 = ( , 0, 0)T。令u1 = x1 ? y1 = (2 ? , ?1, 0)T, 有 A1 的第2列对角线以下(含对角线)为x2 = (–3/ , –1)T, * */97 五、特征值与特征向量 ||x2|| = 。用x2构造镜面反射阵 ,使得 = ( , 0)T, 令u2 = x2 ? y2 = ( , ?1)T, 易得 * */97 五、特征值与特征向量 (3)用平面旋转变换 * */97 五、特征值与特征向量 返回 人有了知识,就会具备各种分析能力, 明辨是非的能力。 所以我们要勤恳读书,广泛阅读, 古人说“书中自有黄金屋。 ”通过阅读科技书籍,我们能丰富知识, 培养逻辑思维能力; 通过阅读文学作品,我们能提高文学鉴赏水平, 培养文学情趣; 通过阅读报刊,我们能增长见识,扩大自己的知识面。 有许多书籍还能培养我们的道德情操, 给我们巨大的精神力量, 鼓舞我们前进。 * * * */97 4. 带参数的迭代法 逐次(超)松弛法中的参数? 的选取不易。另一种方法是Richardson同步迭代法(简称RF法)。 当A对称正定时,易验证当0 ? 2/?1时,RF法收敛, 而当 时收敛最快。其中?1和?n分别为A 的最大和最小特征值。 四、方程组迭代法 * */97 解题指导: 1 . 本部分主要内容是用迭代法求线性方程组的解,重点是Jacobi和Gauss-Seidel迭代法,掌握相关的计算公式、迭代矩阵以及收敛的充要条件和充分条件,并会用相关定理判别迭代法的收敛性。 (1) 对所构造的迭代法能写出具体的迭代矩阵B,并利用||B|| 1判别收敛性; (2) 对不满足充分条件的方程组或A带有参数的方程组,判别收敛性通常需求迭代矩阵B的特征值和谱半径? (B),并由? (B) 1判别是否收敛。 (3) 掌握与迭代法相关的向量序列、矩阵序列的收敛性 四、方程组迭代法 * */97 引理3: 设A为n阶方阵,则 (4) 利用渐进收敛速度R(B) = – ln ? (B) 比较收敛快慢; 2 . J法和GS法 (1) 对给定的方程,写出迭代公式; (2) 写出对应的迭代矩阵,用范数或谱半径判别收敛性; (3) 从方程组的系数矩阵A直接判别收敛性: 注意:矩阵B的特征多项式两种形式: | A – ? I | 与 | ? I – A | 等价(仅差一个符号),用于求特征值没区别。 四、方程组迭代法 * */97 三、典型例题 1、设 ,证明 (1)当 时,矩阵A是正定的 。 (2)用Jacobi迭代法求解方程组Ax = b是收敛的。 提示:(1)顺序主子式大于0,得2a2 1,整理即得。 (2)||BJ|| 不一定小于1,但f (?) = ? (?2 – 2a2), 知? = 0,?2 1,故? (BJ) 1 四、方程组迭代法 * */97 2、设方程组 ,其中 (1) 用迭代法收敛的充要条件求出使Jacobi迭代法和G-S迭代法均收敛的a的取值范围; (2)当a ? 0时,给出这两种迭代法的收敛速度之比。 解:(1) 四、方程组迭代法 * */97 一般地,判别方程组Ax = b的迭代法的收敛性,可先从A本身的特征进行:如严格对角占优,对称正定;其次是利用迭代矩阵的范数||B||∞,注意到Jacobi法的严格对角占优与||B||∞ 1等价,最后使用谱半径。 四、方程组迭代法 (2) * */97 此外,对于J法和GS法,求迭代矩阵B的谱半径,可不必先求迭代矩阵B。事实上, 对于J法,B = – D -1(L+U),特征方程为 |? I – B| = |? I + D -1(L+U)| = | D -1[? D + (L+U)]| =| D -1|·|? D + (L+U)| = 0 即 |? D + (L+U)| = 0 特征多项式就是A的对角线元素乘? 后的行列式。 四、方程组迭代法 *
原创力文档

文档评论(0)