Jacobi迭代法求解线性方程组实验报告.docxVIP

仿真平台与工具应用实践 Jacobi迭代法求解线性方程组 A 实验报告 院 系： 专业班级： 姓 名： 学 #: 指导老师： 二实验目的 熟悉Jacobi迭代法原理； 学习使用Jacobi迭代法求解线性方程组; 编程实现该方法； 三、实验内容 应用Jacobi迭代法解如下线性方程组: 4 ?勺_心+兀3 =7 4xi-8x2 + x3=-21 ,要求计算精度为IO -2%! +x2+ 5xy = 15 四、 A 实验过程 ⑴、算法理论 儿“?。仞迭代格式的引出是依据迭代法的基本思想：构造一个向量系列 仅（）｝,使其收敛至某个极限X*,则X*就是要求的方程组的准确解。 Jacobi迭代 ⑷內+%勺+…+5耳=6 ⑷內+%勺+…+5耳=6 ⑷內+%勺+…+5耳=6 ⑷內+%勺+…+5耳=6 将方程组:勺內+如吃+…+色儿=% 将方程组: 勺內+如吃+…+色儿=% 〔°川旺+(52欠2+…+勺“丿舁=仇 “=勺2兀2+勺3习+?+九￡+8| 在假设?H0,改写成/从乜冲…+从“ ⑵ 亠=乞內+仇2勺+…+乞(,1)兀T + Sn 如果引用系数矩阵 bxg】 bx g】 ? A _ ? z — ? ? 9 O — ? ,g 一 ■ ￡ Sn 及向量X = 方程组(2)和(2)分别可写为：AX=h及X = 3X + g,这样就得到了 jacobi 迭代格式X^=BXk +g。用jacobi迭代解方程组AX=b时，就可任意取初值 X°带入迭代可知式X=BX*g,然后求limXS但是，〃比较大的时候, 写方程组⑴和⑵是很麻烦的，如果直接llM, b能直接得到B, g就是矩阵 与向量的运算了，那么如何得到B, g呢实际上，如果引进非奇异对角矩阵 ?!1 0 仏 H 0) D = 0 \ 一。…叽 将A分解成：A = A-D + D.要求AX=b的解，实质上就有 AX = {A-D)X + DXy而 D是非奇异的，所以 Z)T存在，DX = AX+ (D-A)X, 从而有X = D~lAX+D~lb,我们在这里不妨令B = I-D-lA,g = D~lb就得到 jacobi 迭代格式：X+l =BXk +g （2）算法框图 000x = 0 0 0 x = (3入算法程序 m文件： function x=jacobi(A,b,prdelta,n) N=length(b); %返回矩阵b的最大长度 for k=l:n for j=l:N x(j)=(b(j)-A(L[l:j-lj+l:N])*P([l:j-lJ+l:N]))/A(jj); end err=abs(norm(x-P)); %求(x，-P)模的绝对值 P=xf; ? ?? 讦(errdelta) %判断是否符合精度要求 break; end end E=eye(N,N); %产生N行N列矩阵 D=diag(diag(A)); f=A*inv(D); %彳是八乘D的逆矩阵 B=E-f; P x=x; k,err B MATLAB 代码： ? clear all A二[4,1;4,l;-2, 1, 5]; ! b=[7, -21, 15] p=[o,o,oy； x=jacobi(A,bFle?7,20) (4人算法实现 4xn-x12+xI3 =7 用jacohi迭代法求解方程组: 4x^| — 8x7? + X23 = —21 一 2x31 + x32 + 5^33 = 15 正常计算结果是2, 3, 4 ,下面是程序输出结果: P = k = 17 err = B = 实验体会 MATLAB是非常实用的软件，能够避免大量计算，简化我们的工作, 带来便捷。通过本次试验，我了解了 MATLAB软件，提高了解决实际问 题的能力。 参考文献 《科学计算与数学建模实验报告Jacobi迭代法求解线性方程组》