数值分析11.ppt

数值分析11.ppt

ID:49395087

大小:278.00 KB

页数:18页

时间:2020-02-06

数值分析11.ppt_第1页
数值分析11.ppt_第2页
数值分析11.ppt_第3页
数值分析11.ppt_第4页
数值分析11.ppt_第5页
资源描述:

《数值分析11.ppt》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、边值问题数值解算例对称正定矩阵的收敛性超松驰迭代算法分块矩阵的块迭代《数值分析》11例1常微分方程边值问题在x1=0.1,x2=0.2,···,x9=0.9处的数值解解:令h=0.1,记yj=y(xj)(j=1,2,···,9),将代入微分方程,整理得–yj-1+(2–h2)yj–yj+1=xjh2(j=1,2,···,9)2/18yj=[xjh2+yj-1+yj+1]/(2–h2),(j=1,2,···,9)–yj-1+(2–h2)yj–yj+1=xjh2高斯-赛德尔迭代格式:3/18程序h=0.1;x=0:h:1;y=zeros(size(x));r1=h*h;r=2

2、-r1;er=1;k=0;whilee>0.0001er=0;forj=2:10s=(y(j-1)+y(j+1)+r1*x(j))/r;er=max(er,abs(s-y(j)));y(j)=s;endk=k+1end准确解:y(x)=sinx/sin1-x-----y(x)oyj4/18正方形区域上第一边值问题准确解:O1x1y5/18取正整数n,记对区域做网格剖分:xi=ih(i=0,1,……,n+1)yj=jh(j=0,1,……,n+1)在点(xi,yj)处记uij=u(xi,yj),五点差分格式整理6/18边界条件:u0,j=0(j=1,···,n);un,j=(j

3、=1,···,n);ui,0=0(i=1,···,n);ui,n=0(i=1,···,n)结点数n2102202402迭代次数91303978CPU时间(s)0.141.84324.6720误差0.00285.6995e-46.6671e-4高斯-赛德尔迭代法实验:7/18定理4.4方程组Ax=b中,若A是实对称正定矩阵,则Gauss-Seidel迭法收敛证明:由A=D–L–LTBG-S=(D–L)-1LT设λ为BG-S的任一特征值,x为其特征向量,则(D–L)-1LTx=λxLTx=λ(D–L)xA正定,故p=xTDx>0,记xTLTx=a,则有xTLTx=λxT(

4、D–L)xxTAx=xT(D–L–LT)x=p–a–a=p–2a>08/18所以,迭代矩阵BG-S的谱半径ρ(BG-S)<1,从而当方程组Ax=b的系数矩阵A是实对称正定矩阵时,Gauss-Seidel迭代法收敛称R=–lnρ(B)为迭代法的渐近收敛速度9/18(i=1,2,···,n;k=1,2,3,··········)超松驰(SOR)迭代法Gauss-Seidel迭代格式10/18定理4.8若A是对称正定矩阵,则当0<ω<2时SOR迭代法解方程组Ax=b是收敛的定理4.9若A是严格对角占优矩阵,则当0<ω<1时SOR迭代法解方程组Ax=b是收敛的迭代矩阵11/18例4

5、.3用SOR方法解方程组(ω=1.4)w=input('input:w:=');A=[2-10;-12-1;0-12];b=[1;0;1.8];x=[1;1;1];er=1;k=0;whileer>0.0005er=0;k=k+1;fori=1:3s=0;t=x(i);x(i)=0;forj=1:3s=s+A(i,j)*x(j);endx(i)=(1-w)*t+w*(b(i)-s)/A(i,i);er=max(abs(x(i)-t),er);endendkk=10x=1.19991.39991.5999ω=1.2,只需k=612/18块迭代法简介设A∈Rn×n,x∈Rn,b

6、∈Rn将方程组Ax=b中系数矩阵A分块其中,Aii∈Rni×ni,Aij∈Rni×nj,xi∈Rni,bi∈Rni13/18将A分解,A=DB–LB–UBJacobi块迭代DBX(k+1)=(LB+UB)X(k)+bi=1,2,···,r(2)Gauss-Seidel块迭代DBX(k+1)=LBX(k+1)+UBX(k)+bi=1,2,···,r14/18块迭代求解X1=[x1x2x3]TX2=[x4x5x6]Tb1=[050]Tb2=[6–26]TDX1–X2=b1–X1+DX2=b2DX1(k+1)=b1+X2(k)DX2(k+1)=b2+X1(k)15/18(i,j

7、=1,···,n)边值问题:(i,j=1,···,n;k=1,2,3,······)保留三对角块16/18取正整数n,h=1/(n+1)离散点xi=ihyj=jhzm=mh(i,j,m=0,1,···,n+1)用七点差分格式计算求解,用slice命令绘四维图17/18高斯-赛德尔迭代法18/18

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。