mat7labs的个人空间 http://blog.eetop.cn/1631436 [收藏] [复制] [分享] [RSS]

日志

雅克比迭代法解方程

已有 534 次阅读2019-11-26 01:53 |系统分类:芯片设计

function [X,k,err]=jacobi(A,B,P,delta,max1)
N=length(B);
D=zeros(N);
for i=1:N
    D(i,i)=A(i,i);
end
%求出迭代式中的两个参数
Bj=inv(D)*(D-A);
Fj=inv(D)*B;
%迭代过程
for k=1:max1
    X=Bj*P+Fj;
    if (norm(abs(X-P))<delta)|((norm(abs(X-P))/(norm(abs(X-P))+eps))<delta)
        err=norm(abs(X-P));
        break
    end
        P=X;
    end

评论 (0 个评论)

facelist

您需要登录后才可以评论 登录 | 注册

关闭

站长推荐上一条 /4 下一条

小黑屋| 关于我们| 联系我们| 在线咨询 |  EETOP 创芯网 ( 京ICP备:10050787号 京公网安备:11010502037710 )

GMT+8, 2020-1-28 19:49 , Processed in 0.028116 second(s), 6 queries , Gzip On, Redis On.

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

返回顶部