matlab雅克比迭代法算例.doc

《matlab雅克比迭代法算例.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、雅克比迭代法算例子函数：m文件%雅克比Jacobi.mfunctionx=Jacobi(A,b,x0,eps,N)%用雅克比迭代法解方程n=length(b);x=ones(n,1);k=0;whilek<=Nfori=1:nx(i)=(b(i)-A(i,[1:i-1,i+1:n])*x0([1:i-1,i+1:n]))/A(i,i);endk=k+1;ifnorm(x-x0,inf)Nwarning('算法超出迭代的次数');elsedisp(['迭代次数=',num2s

2、tr(k)])xend主程序：（调用以上程序）。M文件clcformatlongA=[10,-1,2,0;-1,11,-1,3;2,-1,10,-1;0,3,-1,8];b=[6;25;-11;15];x0=[0;0;0;0];eps=1e-6;N=30;x=Jacobi(A,b,x0,eps,N);