MATLAB设计Jacobi迭代法运行错误,设计函数：function [x,k]=jacobi_f(A,b,x0,tol,max)[n,m]=size(A);xold=x0;C=-A;for i=1:nC(i:i)=0;C(i,:)=C(i,:)/A(i,i);d(i,:)=b(i)/A(i,i);endi=1;while i

Question

MATLAB设计Jacobi迭代法运行错误,
设计函数：
function [x,k]=jacobi_f(A,b,x0,tol,max)
[n,m]=size(A);xold=x0;C=-A;
for i=1:n
C(i:i)=0;
C(i,:)=C(i,:)/A(i,i);
d(i,:)=b(i)/A(i,i);
end
i=1;
while i

Answer 1

你的x没有初始化,导致当你迭代不收敛时x没有赋值,肯定x作为返回值报错.
函数开始加一句
x=[];
还有我觉得你的算法有问题,
C(i:i)=0;//i:弄错了吧
d(i,:)=b(i)/A(i,i);//冒号是什么意思
具体请自行检查.