实验线性方程方程组迭代法收敛速度.docxVIP

武汉理工大学理学院数学系课程实验报告 课程名称： 数值分析 班级 信计1202 H期 20144)6-()6 成纸评定 姓名 陈细军 实验室 老师签名 学号 0121214410219 实验六：线性方程方程组迭代法收敘速度 实验 名称 Matlab 软件 1、实验内容 研究耕线性方程方程组Av =乃迭代法收敛速度■给定人w R”z为五对角矩阵 实 脸 冃 的 及 内 容3-1/2-1/4-1/2 实 脸 冃 的 及 内 容 3 -1/2 -1/4 -1/2 3 -1/2 -1/4 -1/4 -1/2 ■ ? 3 ? ■ -1/2 ? ? -1/4 ■ ■ ? -1/4 ■ -1/2 ■ 3 ? -1/2 -1/4 -1/2 3 -1/4 -1/2 A = -1/4 -1/2 3 (1)选取不冋的初始向凤X?及右端项向Zb ,给定迭代谋斧要求?用雅可比迭代和尚斯 赛德尔迭代法求解,观察础到的序列足否收敛？若收敛，记录迭代次数，分析计算纠果并 得出你的结论. (2)用SOR迭代法求上述方程组的解，松弛系数e以102的不同值，在 卜⑷_兀如)性10」时停止迭代.记录迭代次数，分析计算结果并得出你的结论. 实验要求：用c涪言或者Mathb涪言编程?但不能克接调用求做命令。 实验6实验步骤： 程序设计 计算实例 摄写实验报告 雅可比迭代程序代码如卜?： function [xf kJ =jacobi(afb) n=length (b); m-^200; e=le-6; xO=zeros(nr1); k=0; x=xO; x0=x+2*e; d-diag(diag(a)); lw-tril(a,-1); u=-triu (a,1); 实 while norm(xO-xrinf)eikm k-k+1;xO^x; 验 x=inv(d)*(1+u)*x+inv(d)*b; k, 原 disp(x), end; 理 if k==m error (?失败、? 「次数卜线?); 步 end 主程序： 骡 A=1 inspace (3, 3y 20); 、 B=1inspace(-1/2, -1/2, 19); C=1inspace(-1/4,-1/4,18); M=diag(A)+diag(B, 1) +diag(B,-l) -diag (C, 2) +diag (C, -2); b=[l:20] b = Columns 1 through 17 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 Columns 18 through 20 18 19 20 [x, k] -jacobi (M, b) 0. 724659702143361 344383768859657 007149458117254 668987888979304 334391909533303 1.000390017119131 666799906930891 333311335950194 999807268587987 666131976683678 331998265054314 996701088892924 658503676668257 313278141752161 9? 950091274B27077 545450436392965 11.025137935197158 281671786490090 697322157117078 389692934985428 k 二 23 高斯赛德尔迭代法程序代码如下： function [xr k]-^Gau^Sed(a,b,x0# e#m) n=length (b); if nargin5 m=200; end; if nargin4 e-le-5; end; if nargin3 xO-zeros(nr1); end; k?0; x-xO; x0=x+2*e; al^tril(a); ial=inv (al); while norm(x0-xr inf)ekm XH—ial斗(a—al) *xo+iml*b; disp (X二 end i f kHHm error enQ H前氓 AJLinspacep3? 20)- ￥1 inspace ( 1、2;1、2?19二 culinspace(ll、4?ll、4?18二 M=diag(A)+diag?l)+diag(w—1Tdiagp2)+diagr—2)- b=〔=20〕 Columns 1 Through 17 1 2 3 4 5 6 7 8 9S11 12 13 M 15 16 17 Columns 18 through 20 一8 】9 20 xo=r2ol; 〉〉〔x?kl=Gaulseds5g??pm) 0- 72亠66】M0895