雅可比迭代法的MATLAB程序.doc

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

1、雅可比迭代法的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表示计算成功;ifnargin<4it_max=100;endifnargin<3ep=le-5;endn=length(A);k=0;x=zeros(n,1);y=zeros(n,1);index=1;whilek<=it

2、_maxfori=1:nifabs(A(i,i))

3、代次数;%index为指标变量,index=0表示计算失败,index=1表示计算成功;ifnargin<4it_max=100;endifnargin<3ep=le-5;endn=length(A);k=0;x=zeros(n,1);y=zeros(n,1);index=1;whilek<=it_maxfori=1:nifabs(A(i,i))

4、i)=(b(i)-A(i,1:i-1)*y(1:i-1)-A(i,i+1:n)*x(i+1:n)/A(i,i);endendifnorm(y-x,inf)