地球物理计算方法第一章教材.ppt

*/47 §6.8 解线性方程组的迭代法 可见，迭代序列逐次收敛于方程组的解，而且迭代7次得到精确到小数点后两位的近似解。 k x1(k) x2(k) x3(k) ‖x(k)-x*‖? 0 1 2 3 4 5 6 7 0 1.4 1.11 0.929 0.9906 1.01159 1.000251 0.9982364 0 0.5 1.20 1.055 0.9645 0.9953 1.005795 1.0001255 0 1.4 1.11 0.929 0.9906 1.01159 1.000251 0.9982364 1 0.5 0.2 0.071 0.0355 0.01159 0.005795 0.0017636 计算结果列表如下: */47 G-S迭代法的计算公式为: §1.8 解线性方程组的迭代法 */47 同样取初始向量x(0)=(0,0,0)T，计算结果为： 由计算结果可见，G-S迭代法收敛较快.取精确到小数点后两位的近似解，G-S迭代法只需迭代 3 次，而J迭代法需要迭代 7 次。 k x1(k) x2(k) x3(k) ‖x(k)-x*‖? 0 1 2 3 0 1.4 1.0634 0.9951044 0 0.78 1.02048 00 1.026 0.987516 11 0.4 0.0634 0.0048956 §1.8 解线性方程组的迭代法 */47 若使‖x(k) –x*‖? ,只需, 即 可以事先估计达到某一精度需要迭代多少步。 §1.8 解线性方程组的迭代法 */47 例2 用J迭代法求例1中方程组的解，取x(0)=(0,0,0)T，若使误差??x(k)-x*???10-5，问需要迭代多少次？ 解 由例1知，x(1)=(1.4,0.5,1.4)T， 于是有，??x(1)-x(0)???=1.4，??B???=0.5 。 k应满足 故取k=19，即需要迭代19次。 §1.8 解线性方程组的迭代法 */47 §6.8 解线性方程组的迭代法 试建立一个收敛的迭代格式，并说明收敛性。 解 按如下方法建立迭代格式 由于迭代矩阵的行范数小于1，故此迭代法收敛。 例3 已知解线性方程组 */47 §6.8 解线性方程组的迭代法 例4 用SOR方法解线性方程组 解 SOR方法迭代公式为 方程组的精确解是 x*=(2,1,-1)T。 */47 §6.8 解线性方程组的迭代法 k x1(k) x2(k) x3(k) 0 1 2 3 … 20 0 3.65 22.5661399 …… 1.9999987 0 0.8845882 0.4230939 0.6948261 …… 1.0000013 0 -0.2021098 -0-0.4952594 …… -1.0000034 从结果可见，迭代 20 次时已获得精确到小数点后五位的近似解。如果取 ?=1.25，则需要迭代 56 次才能得到具有同样精度的近似解；如果取 ?=1，则需迭代 110 次以上。 取x(0)=(0,0,0)T，?=1.46,计算结果如下: */47 §6.8 解线性方程组的迭代法 */47 §6.8 解线性方程组的迭代法 注意：未必Seidel方法一定比Jacobi方法好。 */47 §1.9 病态方程组和迭代改善法 本节内容 一. 矩阵的条件数 二. 误差分析 返回章节目录 */47 §1.9 病态方程组和迭代改善法 一.矩阵的条件数 */47 §6.10 病态方程组和迭代改善法 */47 §6.10 病态方程组和迭代改善法 右端项b的扰动对解的影响 */47 §6.10 病态方程组和迭代改善法 系数矩阵A的扰动对解的影响 */47 §6.10 病态方程组和迭代改善法 条件数的定义 */47 §6.10 病态方程组和迭代改善法 条件数的性质 */47 §1.9 病态方程组和迭代改善法 */47 §1.9 病态方程组和迭代改善法 */47 §6.10 病态方程组和迭代改善法 */47 §1.9 病态方程组和迭代改善法 */47 “病态”方程的经验判断 */47 §1.9 病态方程组和迭代改善法 */47 §1.9 病态方程组和迭代改善法 二. 误差分析 */47 §6.10 病态方程组和迭代改善法 * * * * */47 §1.7 向量和矩阵的范数 */47 §1.7 向量和矩阵的范数 */47 §1.7 向量和矩阵的范数 */47 §1.7 向量和矩阵的范数 */47 §1.7 向量和矩阵的范数 */47 §1.7