数值分析08解线性方程组的迭代法.ppt

* 下面描述这一过程。 SOR迭代法的分量形式 推导SOR迭代格式的矩阵形式（以三阶方程为例） 推导SOR迭代格式的矩阵形式 SOR法收敛性的结论： （1）SOR方法收敛的必要条件为0 ?2 （2）若系数阵A对称正定，则当0 ?2时， SOR方法收敛 （3）若系数阵A严格对角占优，则当0 ??1时， SOR方法收敛。 在计算机上采用动态计算形式 (1) x(i)=0 (i=1,2,…,n) (2)对k =1,…,Kmax,循环计算到第(7)步 (3)置ER=0 4、对称超松弛迭代法（SSOR法） SSOR迭代法的矩阵形式： 注： （1）关于 的收敛条件和准则与SOR方法相同； （2）收敛快慢对 的选取不敏感。 5、块超松弛迭代法（BSOR法） 例 3-17 (热传导问题) 设有一维热传导方程的初边值问题 试用数值方法求出t=0.2时刻金属杆的温度分布. 三．应用实例 解:(1)对空间进行离散. (2)对微分算子进行离散. t 0.02 0.01 0 0 0.1 0.2 …. 0.9 1 x 采用无条件稳定的Crank-Nicholson格式,则有 或 加上边界条件后有 加上边界条件后有 其矩阵形式为 ATn+1=dn 解:(1)剖分求解域. 0 1 2 …. N N+1 X Y N+1 N ： 2 1 0 例 3-18 数值求解正方形域上的Poisson方程边值问题 解:(1)剖分求解域. (2)对微分算子进行离散. 0 1 2 …. N N+1 X Y N+1 N ： 2 1 0 在每个点(xi,yj)上的有限差分方程为 在边界上 又称为五点差分格式 (i,j) (i,j+1) (i,j-1) (i-1,j) (i+1,j) 对非边界点进行编号: 顺序为-----从下往上,从左往右 相应的解向量和右端向量分别为 Gauss-Seidel迭代法参考程序: n=9; b(2:n+1,2:n+1)=0.02; U=zeros(n+2,n+2); e=0.000000001; for k=1:1000 %迭代求解 er=0; for j=2:n+1 for i=2:n+1 Ub=U(i,j); U(i,j)=(U(i-1,j)+U(i+1,j)+U(i,j-1)+U(i,j+1)+b(i,j))/4; 　　er=er+abs(Ub-U(i,j)); %估计当前误差 　　end end if er/n^2e,break;end %判断是否达到计算精度，如果达到则退出循环 end f(x,y)=2, h=0.1 差分方程组的矩阵形式为 Au=f 如果把每一条线上的网点看作一个组,如 其中 可用块迭代法 (即线迭代法) 求解Au=f function [u,k]=kgs2(n) f=2*1/(n+1)^2*ones(n+2,n+2); a=-1*ones(1,n); b=4*ones(1,n);c=-1*ones(1,n); u=zeros(n+2,n+2);e=0.00001; for k=1:2000 er=0; for j=2:n+1 Ub=u(:,j); d(1:n)=f(2:n+1,j)+u(2:n+1,j-1)+u(2:n+1,j+1) ; x=zg(a,b,c,d);u(2:n+1,j)=x; er=er+norm(Ub-u(:,j),1); end if er/n^2e ,break;end end 习题三　P122-----17，18（1），21，25 上机作业： 数值实验题三--P124--3、4（任选一题） 习题　 P122-----21，24，25，26 上机作业： 数值实验题三--P124--3、4（任选一题） * 让M=D,N=A-D, 即可得Jacobi 迭代。 Jacobi迭代矩阵 推导其分量形式 第i个方程除以aii(i =1,2,…,n),得 Jacobi迭代法的分量形式 则 x(k+1)=BJx(k)+g ， 这里 BJ=D-1(L+U) , g=D-1b Jacobi迭代法的分量形式