Question

怎么用matlab求线性方程组的雅可比迭代矩阵

Answer 1

function [x,n]=Jacobi_Solve(A,b,x0,dalt)
% Jacobi 跌代法解线性方程组
%[x,n]=Jacobi_Solve(A,b,x0,dalt)
% A 方程组系数
% b 常数项(列向量)
% x0 初始值,默认为 0
% dalt 精度,默认为 10
% x 返回跌代结果
% n 返回跌代次数
e=1; i=0;
r=size(b);%将矩阵b的行数及列数赋值给r
a=b;
if nargin<4 %输入参数个数＜4
dalt=1e-8;
end
if nargin<3
x=zeros(r);%创建一个r行全0的矩阵
else
x=x0;
end
r=r(1);
for t=1:r
a(t)=A(t,t);%选出主对角线上的元素
A(t,t)=0;
A(t,:)=A(t,:)/a(t);
end
b=b./a;
while e>=dalt
Y=b-A*x;
e=max(abs(Y-x));
x=Y; i=i+1;
end
if nargout>1 %函数输出变量数的个数>1
n=i;
end

望采纳！