原创力文档- 1、本文档共2页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
雅可比迭代法的MATLAB 程序：
Function[x，k,index ］=Jacobi （A,b ，ep,it—max ）
％ 求线性方程组的雅可比法；
% A 为方程组的系数矩阵;
％ b 为方程组的右端项；
％ x 为方程组的解；
％ ep 为精度要求，缺省值为le-5 ；
％ it_max 为最大迭代次数,缺省值为100；
% k 为迭代次数；
% index 为指标变量，index=0 表示计算失败,index=1 表示计算成功;
if nargin 〈4
it_max=100;
end
if nargin<3
ep=le—5 ；
end
n=length(A ）;k=0;x=zeros （n ，1);y=zeros （n，1)；index=1;
while k<=it_max
for i=1 ：n
if abs (A(i,i ））<le—10
index=0 ；
return;
end
y(i)=(b(i ）—A （i ，1:n)*x(1:n)+A(i,i）*x （i))/A （i,i ）；
end
if norm （y-x，inf ）<ep
break ；
end
k=k+1 ；
x=y ；
end
高斯—赛德尔迭代的MATLAB 程序
Function ［x,k ，index]=Gau-seidel(A，b,ep,it—max)
％ 求线性方程组的高斯-赛德尔迭代法；
％ A 为方程组的系数矩阵；
% b 为方程组的右端项;
％ x 为方程组的解；
% ep 为精度要求，缺省值为le—5；
％ it_max 为最大迭代次数,缺省值为100;
% k 为迭代次数;
％ index 为指标变量,index=0 表示计算失败,index=1 表示计算成功；
if nargin<4
it_max=100；
end
if nargin<3
ep=le-5；
end
n=length(A);k=0;x=zeros(n,1）；y=zeros （n ，1）；index=1;
while k 〈=it_max
for i=1 ：n
if abs （A （i,i ））<le—10
index=0;
return;
End
if i==1
y （i)= （b(i)-A(i，i+1 ：n) ＊x （i+1:n)/A （i ，i ）；
elseif i==n
y （i)=(b （i ）-A(i，1：i-1 ）*y(1：i—1）/A(i ，i ）;
else
y （i)= （b(i)—A(i,1 ：i—1) ＊y(1:i—1）-A （i ，i+1:n)*x （i+1:n ）/A(i ，i) ；
end
end
if nor （y-x,inf) 〈ep
break ；
end
k=k+1;
x=y ；
end
文档评论（0）