预览加载中，请您耐心等待几秒...

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

预览
1/3
2/3
3/3

在线预览结束，喜欢就下载吧，查找使用更方便

下载提示 文本预览

如果您无法下载资料，请参考说明：

1、部分资料下载需要金币，请确保您的账户上有足够的金币

2、已购买过的文档，再次下载不重复扣费

3、资料包下载后请先用软件解压，在使用对应软件打开

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=D\b;x=B*x0+f;n=1;whilenorm(x-x0)>=epsx0=x;x=B*x0+f;n=n+1;if(n>=m)disp('可能不收敛');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=B*x0+f;n=n+1;if(n>=m)disp('迭代次数过多，可能不收敛');return;end;endPS:有兴趣小朋友试一下，结果与x0无关