原创力文档- 1、本文档共18页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
* 边值问题数值解算例 对称正定矩阵的收敛性 超松驰迭代算法 分块矩阵的块迭代 《数值分析》 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 = xj h2 ( j= 1,2,···,9) 2/18 yj = [xj h2 + yj-1 + yj+1]/ (2 – h2) , ( j= 1,2,···,9) – yj-1 + (2 – h2) yj – yj+1 = xj h2 ? 高斯-赛德尔迭代格式: 3/18 程序 h=0.1;x=0:h:1; y=zeros(size(x)); r1=h*h;r=2-r1; er=1;k=0; while e0.0001 er=0; for j=2:10 s=(y(j-1)+y(j+1) + r1*x(j))/r; er=max(er,abs(s-y(j)));y(j)=s; end k=k+1 end 准确解: y(x)=sin x/sin 1 - x ----- y(x) o yj 4/18 正方形区域上第一边值问题 准确解: O 1 x 1 y 5/18 取正整数n,记 对区域做网格剖分: xi = i h ( i = 0，1，……，n+1 ) yj = j h ( j = 0，1，……，n+1) 在点(xi，yj ) 处记 uij = u(xi ,yj) ,五点差分格式 整理 6/18 边界条件: u0, j = 0 (j = 1, ···, n); un, j = ( j = 1, ···, n); ui,0 = 0 ( i = 1, ···, n); ui,n = 0 ( i = 1, ···, n) 结点数n2 102 202 402 迭代次数 91 303 978 CPU时间(s) 0.14 1.843 24.6720 误差 0.0028 5.6995e-4 6.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)-1LT x =λx ? LT x =λ(D – L)x ? A正定,故 p = xTDx0, 记 xTLTx = a , 则有 xTLTx =λxT(D – L)x xTAx=xT(D – L – LT)x=p – a – a =p – 2a 0 8/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.3 用SOR方法解方程组(ω=1.4) w=input(input: w:=); A=[2 -1 0;-1 2 -1;0 -1 2]; b=[1;0;1.8]; x=[1;1;1]; er=1;k=0; while er0.0005 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)=(1-w)*t+w*(b(i)-s)/A(i,i);
文档评论（0）