数值分析雅克比高斯迭代法MATLAB程序.doc

《数值分析雅克比高斯迭代法MATLAB程序.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、1.雅可比迭代法：function[x,n]=jaccbi(A,b,x0,eps,t)ifnargin==3;eps=1e-6;m=200;elseifnargin<3error('输入的数有误');return;elseifnargin==5m=t;endD=diag(diag(A));L=-tril(A,-1);U=-triu(A,1);B=D(L+U);f=Db;x=B*x0+f;n=1;whilenorm(x-x0)>=epsx0=x;x=B*x0+f;n=n+1;if(n>=m)disp('

2、可能不收敛');return;end;end2.高斯—赛德尔迭代法：function[x,n]=gsdddy(A,b,x0,eps,t)ifnargin==3;eps=1e-6;m=200;elseifnargin<3error('输入有误');return;elseifnargin==5m=t;endD=diag(diag(A));L=-tril(A,-1);U=-triu(A,1);B=(D-L)U;f=(D-L)b;x=B*x0+f;n=1;whilenorm(x-x0)>=epsx0=x;x=

3、B*x0+f;n=n+1;if(n>=m)disp('迭代次数过多，可能不收敛');return;end;endPS:有兴趣小朋友试一下，结果与x0无关