原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
《数值分析》第七章实验报告
2010/5/18
EXERCISE SET 7.3 P451
1、3、Jacobi 法
e)
4x x x x 6
1 2 3 5
x 3x x x 6
1 2 3 4
2x x 5x x x 6
1 2 3 4 5
x x x 4x 6
1 2 3 4
2x x x 4x 6
2 3 4 5
编写MATLAB 程序如下
function [x,k,index] = Jacobi(A,b,ep,it_max)
%求解线性方程组的Jacobi迭代法，其中
% A 方程组的系数矩阵
% b 方程组的右端项
% ep 精度要求。省缺为1e-5
% it_max 最大迭代次数，省缺为100
% x 方程组的解
% k 迭代次数
% index index=1表示迭代收敛到指定要求；
% index=0表示迭代失败
if nargin 4 it_max = 100 ;end
if nargin 3 ep = 1e-5;end
n = length(A);k = 0;
x = zeros (n,1);y = zeros(n,1);index = 1;
while 1
for i = 1:n
y(i) = b(i);
for j = 1:n
if j ~= i
y(i)=y(i)-A(i,j)*x(j);
end
end
if abs(A(i,i))1e-10 | k == it_max
index = 0;return;
end
1
y(i) = y(i)/A(i,i);
end
if norm(y-x,inf)ep
break;
end
x = y;k = k+1;
end
在命令行窗口运行
A=[4 1 1 0 1;-1 -3 1 1 0;2 1 5 -1 -1;-1 -1 -1 4 0;0 2 -1 1 4];b=[6 6 6 6 6];ep = 1e-3;
[x,k,index] = Jacobi (A,b,ep)
x =
0.7873
-1.0029
1.8654
1.9128
1.9891
k =
11
index =
1
所以方程的近似解为
x1=0.7873
x2=-1.0029
x3=1.8654
x4=1.9128
x5=1.9891
f)
4x x x 0
1 2 4
x 4x x x 5
1 2 3 5
x 4x x 0
2 3 6
x 4x x 6
1 4 5
x x 4x x 2
2 4 5 6
2
x x 4x 6
3 5 6
仍是运用上面的Jacobi.m 程序，在命令行窗口运行
A=[4 -1 0 -1 0 0;-1 4 -1 0 -1 0;0 -1 4 0 0 -1;-1 0 0 4 -1 0;0 -1 0 -1 4 -1;0 0 -1 0 -1 4];b=[0 5 0 6 -2
6];ep=1e-3;
[x,k,i
文档评论（0）