日志

雅克比迭代法解方程

已有 2857 次阅读| 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

GMT+8, 2025-11-4 15:33 , Processed in 1.024797 second(s), 18 queries , Gzip On, Redis On.

eetop公众号

创芯大讲堂

创芯人才网