雅可比迭代法

最新推荐文章于 2024-06-14 19:00:00 发布

学完编程去卖煎饼果子

最新推荐文章于 2024-06-14 19:00:00 发布

阅读量137

点赞数

分类专栏：学习作业文章标签： matlab Powered by 金山文档

本文链接：https://blog.csdn.net/crucis0719/article/details/129768849

版权

学习作业专栏收录该内容

14 篇文章 0 订阅

订阅专栏

function X=jacdd(A,b,X0,detx,maxN)
[n m]=size(A);
for k=1:maxN
for j=1:m
X(j)=(b(j)-A(j,[1:j-1,j+1:m])*X0([1: j-1,j+1:m]))/A(j,j);
end
djwcX=norm(X'-X0); 
xdwcX=djwcX/(norm(X')+eps);
X0=X'; 
if (djwcX<detx)&(xdwcX<detx)
return
end
end
if (djwcX>detx)&(xdwcX>detx)
disp('雅可比迭代次数已经超过最大迭代次数maxN')
end

举例应用：

A=[27 6 -1;6 15 2;1 1 54];

b=[85;5;110];

X0=[0;0;0];

detx=0.00001;

maxN=20;

X=jacdd(A,b,X0,detx,maxN)