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

* 边值问题数值解算例 对称正定矩阵的收敛性 超松驰迭代算法 分块矩阵的块迭代 《数值分析》 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 正方形区域上第一边值问题 准确解: