数值分析9解线性方程组迭代法.pptVIP

课件 * 课件 解线性方程组的迭代法 ( Iterative Methods ) 雅可比迭代法 赛德尔迭代 迭代法的矩阵表示 《数值分析》 9 课件 ( i,j = 1,···,n ) 记 X=[u11··· u1,n u21··· u2,n······ un,1··· un,n]T AX = b A=LU ? (L + U – I ) 边值问题: 2/18 课件 ( i,j = 1,···,n ) 3/18 课件 例4.1 特点:系数矩阵主对角元均不为零 计算格式 X(1)=B X(0) + f 取 X(0) = 4/18 课件 X* 1.0000 1.0000 1.0000 X(0) 0 0 0 X(1) 0.7778 0.8000 0.8667 X(2) 0.9630 0.9644 0.9778 X(3) 0.9929 0.9935 0.9952 计算格式: X(k+1)=BX(k)+f 准确解 X(4) ········ 0.9987 0.9988 0.9991 5/18 课件 || X(4) – X(3) ||1 = 0.0150, 反例(不收敛) -78.8667 -79.8667 -0.1333 -806.5333 -716.6667 -9.7156 ······ ······ ······ -8.0000 -7.0000 0.8667 6/18 课件 雅可比迭代法 (i = 1,2,…n; k=1,2,……) 取初始向量X(0)=[x1(0) x2(0) ··· xn(0)]T, 迭代计算 (i = 1,2,…,n) 7/18 课件 分裂: AX = b 任取 x(0), 迭代计算产生向量序列: 若 则X* 是方程组 AX = b 的解 迭代矩阵 迭代计算格式 ( k=0,1,2,······ ) X(k+1)= BX(k) + f X(1),X(2),···,X(k),······ splitting D + ( A – D ) X = BX + f DX = (D – A )X + b ? 8/18 课件 迭代法适用于解大型稀疏方程组 (万阶以上的方程组,系数矩阵中零元素占很大比例,而非零元按某种模式分布) 背景: 电路分析、边值问题的数值解和数学物理方程 问题: (1)如何构造迭代格式？ (2)迭代格式是否收敛？ (3)收敛速度如何？ (4)如何进行误差估计？ 9/18 课件 (i = 1,2,…,n) 高斯-赛德尔迭代法 (i = 1,2,…n; k =1,2,……) 取初始向量x(0)=[x1(0) x2(0) ··· xn(0)]T, 迭代计算 10/18 课件 例 11/18 课件 雅可比迭代算法 A=[9 -1 -1;-1 10 -1;-1 -1 15]; b=[7;8;13];x=[0;0;0]; er=1;k=0; while er0.00005 er=0;k=k+1; for i=1:3 s=0;t=x(i);x(i)=0; for j=1:3 s=s+A(i,j)*x(j); end x(i)=t; y(i)=(b(i)-s)/A(i,i); er=max(abs(x(i)-y(i)),er); end x=y;x end 0.7778 0.8000 0.8667 0.9630 0.9644 0.9719 0.9929 0.9935 0.9952 0.9987 0.9988 0.9991 0.9998 0.9998 0.9998 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 12/18 课件 高斯-赛德尔迭代算法 A=[9 -1 -1;-1 10 -1;-1 -1 15]; b=[7;8;13];x=[0;0;0]; er=1;k=0; while er0.00005 er=0;k=k+1; for i=1:3 s=0;t=x(i);x(i)=0; for j=1:3 s=s+A(i,j)*x(j); end x(i)=(b(i)-s)/A(i,i); er=max(abs(x(i)-t),er); end x end 0.7778 0.8778