资源预览内容
第1页 / 共3页
第2页 / 共3页
第3页 / 共3页
亲,该文档总共3页全部预览完了,如果喜欢就下载吧!
资源描述
解:由x=1;2;3,算出b=-3; 7;6知:X1=(2*X3-3)/3 X2=(-X3+7)/2 X3=(2X1-X2+6)/21. 雅克比法:X1(1)=0;X2(1)=0;X3(1)=0;for i=1:50X1(i+1)= (2*X3(i)-3)/3;X2(i+1)=(-X3(i)+7)/2;X3(i+1)=(2*X1(i)-X2(i)+6)/2;End迭代50次结果如下:X1 = Columns 1 through 11 0 -1.0000 1.0000 -0.8333 1.0000 -0.6806 1.0000 -0.5405 1.0000 -0.4121 1.0000 Columns 12 through 22 -0.2945 1.0000 -0.1866 1.0000 -0.0877 1.0000 0.0029 1.0000 0.0860 1.0000 0.1622 Columns 23 through 33 1.0000 0.2320 1.0000 0.2960 1.0000 0.3547 1.0000 0.4085 1.0000 0.4577 1.0000 Columns 34 through 44 0.5029 1.0000 0.5444 1.0000 0.5823 1.0000 0.6171 1.0000 0.6490 1.0000 0.6783 Columns 45 through 51 1.0000 0.7051 1.0000 0.7297 1.0000 0.7522 1.0000X2 = Columns 1 through 11 0 3.5000 2.0000 3.3750 2.0000 3.2604 2.0000 3.1554 2.0000 3.0591 2.0000 Columns 12 through 22 2.9708 2.0000 2.8899 2.0000 2.8158 2.0000 2.7478 2.0000 2.6855 2.0000 2.6284 Columns 23 through 33 2.0000 2.5760 2.0000 2.5280 2.0000 2.4840 2.0000 2.4437 2.0000 2.4067 2.0000 Columns 34 through 44 2.3728 2.0000 2.3417 2.0000 2.3133 2.0000 2.2871 2.0000 2.2632 2.0000 2.2413 Columns 45 through 51 2.0000 2.2212 2.0000 2.2027 2.0000 2.1859 2.0000X3 = Columns 1 through 11 0 3.0000 0.2500 3.0000 0.4792 3.0000 0.6892 3.0000 0.8818 3.0000 1.0583 Columns 12 through 22 3.0000 1.2201 3.0000 1.3684 3.0000 1.5044 3.0000 1.6290 3.0000 1.7433 3.0000 Columns 23 through 33 1.8480 3.0000 1.9440 3.0000 2.0320 3.0000 2.1127 3.0000 2.1866 3.0000 2.2544 Columns 34 through 44 3.0000 2.3165 3.0000 2.3735 3.0000 2.4257 3.0000 2.4736 3.0000 2.5174 3.0000 Columns 45 through 51 2.5576 3.0000 2.5945 3.0000 2.6283 3.0000 2.65932. 高斯-赛德法:XX1(1)=0;XX2(1)=0;XX3(1)=0;for i=1:50XX1(i+1)= (2*X3(i)-3)/3;XX2(i+1)=(-X3(i)+7)/2;XX3(i+1)=(2*X1(i+1)-X2(i+1)+6)/2;End迭代50次结果如下:X1 = Columns 1 through 11 0 -1.0000 1.0000 -0.8333 1.0000 -0.6806 1.0000 -0.5405 1.0000 -0.4121 1.0000 Columns 12 through 22 -0.2945 1.0000 -0.1866 1.0000 -0.0877 1.0000 0.0029 1.0000 0.0860 1.0000 0.1622 Columns 23 through 33 1.0000 0.2320 1.0000 0.2960 1.0000 0.3547 1.0000 0.4085 1.0000 0.4577 1.0000 Columns 34 through 44 0.5029 1.0000 0.5444 1.0000 0.5823 1.0000 0.6171 1.0000 0.6490 1.0000 0.6783 Columns 45 through 51 1.0000 0.7051 1.0000 0.7297 1.0000 0.7522 1.0000X2 = Columns 1 through 11 0 3.5000 2.0000 3.3750 2.0000 3.2604 2.0000 3.1554 2.0000 3.0591 2.0000 Columns 12 through 22 2.9708 2.0000 2.8899 2.0000 2.8158 2.0000 2.7478 2.0000 2.6855 2.0000 2.6284 Columns 23 through 33 2.0000 2.5760 2.0000 2.5280 2.0000 2.4840 2.0000 2.4437 2.0000 2.4067 2.0000 Columns 34 through 44 2.3728 2.0000 2.3417 2.0000 2.3133 2.0000 2.2871 2.0000 2.2632 2.0000 2.2413 Columns 45 through 51 2.0000 2.2212 2.0000 2.2027 2.0000 2.1859 2.0000X3 = Columns 1 through 11 0 3.0000 0.2500 3.0000 0.4792 3.0000 0.6892 3.0000 0.8818 3.0000 1.0583 Columns 12 through 22 3.0000 1.2201 3.0000 1.3684 3.0000 1.5044 3.0000 1.6290 3.0000 1.7433 3.0000 Columns 23 through 33 1.8480 3.0000 1.9440 3.0000 2.0320 3.0000 2.1127 3.0000 2.1866 3.0000 2.2544 Columns 34 through 44 3.0000 2.3165 3.0000 2.3735 3.0000 2.4257 3.0000 2.4736 3.0000 2.5174 3.0000 Columns 45 through 51 2.5576 3.0000 2.5945 3.0000 2.6283 3.0000 2.6593结果分析:从上面的雅可比迭代法和高斯赛德尔迭代法这两种方法所得的实验结果可知,对于同样的矩阵,对于同样的精度,雅可比迭代法要比高斯-赛德尔迭代法慢。从这个例子可以得出结论,用高斯-赛德尔迭代法比雅可比迭代法收敛快,具体地说,在收敛的前提下,及时的更新迭代方程的数据可以获得更好的收敛速度,效率更好。这个结论在多数情况下是成立的,但也有相反的情况,即高斯-赛德尔迭代法比雅可比迭代法收敛慢,甚至还有雅可比迭代法收敛,高斯-赛德尔迭代法发散的情形。
收藏 下载该资源
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号