雅可比迭代法

学完编程去卖煎饼果子

于 2023-03-25 16:34:50 发布

阅读量189

点赞数

分类专栏：学习作业文章标签： matlab

本文链接：https://blog.youkuaiyun.com/crucis0719/article/details/129768849

版权

学习作业专栏收录该内容

14 篇文章

订阅专栏

该代码实现了一个使用雅可比迭代法解决线性方程组的MATLAB函数。函数接受矩阵A、向量b、初始解X0、目标精度detx和最大迭代次数maxN作为参数。在每次迭代中，它更新解向量X，并检查是否达到预设的精度要求。如果超过最大迭代次数且未达到精度，函数会显示错误信息。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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)