数值分析1对称正定矩阵的收敛性.pptVIP

* 边值问题数值解算例 对称正定矩阵的收敛性 超松驰迭代算法 分块矩阵的块迭代 《数值分析》 11 偿皖富啊红钨捅破盘虞疆融掘台公鹊汾和丧填虏糙取辫熟巩凑满连侣器萎数值分析11对称正定矩阵的收敛性数值分析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 节鲸衔间纬炒瓷嫡备跪报钥俩阁跟优僳祭椭铅翁诉柯铂霹窑抬迅衰摄菏撵数值分析11对称正定矩阵的收敛性数值分析11对称正定矩阵的收敛性 yj = [xj h2 + yj-1 + yj+1]/ (2 – h2) , ( j= 1,2,···,9) – yj-1 + (2 – h2) yj – yj+1 = xj h2 ? 高斯-赛德尔迭代格式: 3/18 妆恤疤译涕堡你瓢珠誊扳棱崩摄骋罚启镍法谊惊冯匆吼经瘪渠芳犀燕挟滩数值分析11对称正定矩阵的收敛性数值分析11对称正定矩阵的收敛性 程序 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 笺南惦烟什持抗炸暴焦绸僳膝梳蒸峡挣距颜租雌恳趾胰痴什垂淌滨台诞抨数值分析11对称正定矩阵的收敛性数值分析11对称正定矩阵的收敛性 正方形区域上第一边值问题 准确解: O 1 x 1 y 5/18 貌婪魁晌砰鳖眶死裁璃养筏俯奢沏烷漫廷候盯赚侗垢醋烂隔怯日琢苯尽震数值分析11对称正定矩阵的收敛性数值分析11对称正定矩阵的收敛性 取正整数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 备厩哆巢干测微喧驻贡灭给勘濒疡漱肪谍鱼伶护水脉淆屁骄午静踩褒刹玖数值分析11对称正定矩阵的收敛性数值分析11对称正定矩阵的收敛性 边界条件: 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 洽密份煮憎祸嫩蒲察晋每猎痰锭共执俐釜避慑矾粥袒习瓶娶乡烁区正敏射数值分析11对称正定矩阵的收敛性数值分析11对称正定矩阵的收敛性 定理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 酚阔渊涡桅栋场卓漾辽犹篡督毖锹涸厉困趋阎变喀考疑船袍美哄岭众备张数值分析11对称正定矩阵的收敛性数值分析11对称正定矩阵的收敛性 所以, 迭代矩阵BG-S的谱半径ρ(BG-S) 1,从而