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

数值分析对称正定矩阵的收敛性第一页，共十八页，2022年，8月28日 例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第二页，共十八页，2022年，8月28日 yj = [xj h2 + yj-1 + yj+1]/ (2 – h2) , ( j= 1,2,···,9)– yj-1 + (2 – h2) yj – yj+1 = xj h2 ?高斯-赛德尔迭代格式:3/18第三页，共十八页，2022年，8月28日 程序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+1end准确解: y(x)=sin x/sin 1 - x----- y(x)o yj4/18第四页，共十八页，2022年，8月28