资源描述：
《迭代法实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、一、课题名称Matlab解线性方程组的迭代法二.目的和意义1、通过上机计算体会迭代法求解线性方程组的特点，并能和消去法比较；2、运用所学的迭代法算法，解决各类线性方程组，编出算法程序；3、体会上机计算时，终止步骤：
2、Wy(TL<*或k>(予给的迭代次数)，对迭代法敛散性的意义；4、体会初始解x(,松弛因子的选取，对计算结果的影响三、实验要求1、体会迭代法求解线性方程组，并能与消去法做以比较；2、分别对不同精度要求，如由迭代次数体会该迭代法的收敛快慢；54310,10,10.8722.-s23=£3、对方程组2,3使用SOR方法时，选取松弛因子3二0.8,0.9,1,1.1,1.2等
3、，试看对算法收敛性的影响，并能找出你所选用的松弛因子的最佳者；4、给出各种算法的设计程序和计算结果。四、问题描述给出下列儿个不同类型的线性方程组，请用适当算法计算其解1、设线性方程组42_3-1r1s6-5一36542-1320-7Av15-13-476-16786—8571702-13—421610-11-91734462-713900-18_3-24・4—4022-1_21—4-1141-80Av16-121—8-12243-3-44025-3-100063-32、设对称正定阵系数阵线方程组0000•vi5010012-1031丫33-1194—-333326-35丫64653
4、01X7132-122•勺3820124丫919—863-1_Yio_-21■■400■■0■320-6-356丫320-4-33丫4234-10-3丫59111-4丫6-??1142-15-4719453、三对角形线性方程组_4-100000000_'■Yi'7_-14-10000000v250-14-1000000£-1300-14-100000□2000-14-1000060000-14-1000-1200000-14-100Y714000000-14-10*8-40000000-14-1*9500000000-14.V10_-5a=[4000000004・10000000;
5、04000000;00400000;00040000;0000-14-1000;00000-14-100;000000-14-10;0000000-14-1;00000000・14];b=[75-1326-1214-45-5];n=10;x=[0000000000];六、实验程序设计Jacobi迭代法functionz=chap8_7(a,b,n,x)fork=l:19fori二l:ns=0;forj=l:nifj'i,s二s+a(i,j)*x(j);endendy(i)=(b(i)-s)/a(i,i);endx=y;endz=x计算方程组一»a=[4,2,-3,-1,2,1,0,
6、0,0,0;8,6,-5,一3,6,5,0,1,0,0;4,2,-2,T,3,2,-1,0,3,1;0,-2,1,5,-1,3,-1,1,9,4;一4,2,6,-1,6,7,-3,3,2,3;&6,-&5,7,17,2,6,-3,5;0,2,-1,3,-4,2,5,3,0,1;16,10,-11,-9,17,34,2,-1,2,2;4,6,2,-7,13,9,2,0,12,4;0,0,-1,&-3,-24,-&6,3,-11；»b=[5;12;3;2;3;46;13;38;19;-21J;»n=10;»x二[0000000000]>>z=chap87(a,b,n,x);1.0e+
7、012*Colurnns1through9-0.0809-o.1452-0・39550.87640.55090.21610.17321.93530.2550Column10-1.8123计算方程组二»a=[42-402400;22-1・21320;・4・1141・8-356;0-216-1・4・33;21・8-1224-10_3；43-3-44111-4;025-3-101142;0063-3-4219];»b=[0-620239-22-1545J;»n=8;»x=[0()000000]';»z=chap8_7(a,b,n,x);z=1,0e+005*-2.6705-4.02851
8、.09571.5976-0.6513-1.7472-0.05560.5478计算方程组三»a=[40000000040000000;0・14・1000000;00・14・100000;000・14-10000;0000-14-1000;00000-14-100;000000-14-10;0000000-14-l;00000000-14J;»b二[75・1326・1214-45-5];»n=10;»x=[0000000000];»z=chap8_7(a,b,n,x);z
