MATLAB实现雅克比迭代.docx

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

1、%雅克比迭代clearall;closeall;clc;tic;formatlongedisp('请输入参数');K=input('维数K=');ticA=100*rand(K);%A元素是0-100fori=1:KA(i,i)=sum(abs(A(i,:)))+25*rand(1);%对角占优的量为0~25endb=zeros(K,1);fori=1:K;x=0;forr=1:K;x=x+A(i,r);endb(i,1)=x;end%产生b矩阵，b中的元素为A中对应行的和，目的是使方程解全为1jd=input('控制精度jd=');I=e

2、ye(K);%单位阵L=-tril(A,-1);%下三角U=-triu(A,1);%上三角D=diag(diag(A));%对角矩阵B=D(L+U);%迭代矩阵或B=I-DA;f=Db;%fx0=zeros(K,1);%初始迭代矩阵y=B*x0+f;%迭代公式s=1;%迭代次数whilenorm(y-x0)>=jd&&s<3000;%迭代条件x0=y;y=B*x0+f;%迭代公式s=s+1;endystoct=1:K;yy=abs(y'-1)/1;plot(t,yy);title('绝对误差图')