matlab用雅可比(Jacobi)迭代法求解方程组
·
function x=myjacob(A,b,x0)
myeps=0.0001;
n=length(b);
err0=1;
while err0>myeps
x1=zeros(n,1);
for i=1:n
s1=A(i,i:i-1)*x0(1:i-1);
s2=A(i,i+1:n)*x0(i+1:n);
x1(i)=(b(i)-s1-s2)/A(i,i);
end
err0=norm(x1-x0);x0=x1;
end
x=x1;
end
魔乐社区(Modelers.cn) 是一个中立、公益的人工智能社区,提供人工智能工具、模型、数据的托管、展示与应用协同服务,为人工智能开发及爱好者搭建开放的学习交流平台。社区通过理事会方式运作,由全产业链共同建设、共同运营、共同享有,推动国产AI生态繁荣发展。
更多推荐


所有评论(0)