MATLAB实现雅克比迭代(共2页).docx

精选优质文档-倾情为你奉上%雅克比迭代clear all;close all;clc;tic;format longe disp(请输入参数);K=input(维数K=);ticA=100*rand(K);% A元素是0-100for i=1:K A(i,i)=sum(abs(A(i,:)+25*rand(1); %对角占优的量为025endb=zeros(K,1);for i=1:K; x=0; for r=1:K; x=x+A(i,r); end b(i,1)=x;end %产生b矩阵，b中的元素为A中对应行的和，目的是使方程解全为 1jd=input(控制精度jd=);I=eye(K); %单位阵L=-tril(A,-1);%下三角U=-triu(A,1);%上三角D=diag(diag(A);%对角矩阵B=D(L+U); %迭代矩阵 或B=I-DA;f=Db;%fx0=zeros(