求：雅可比矩陣叠代法MATLAB編程

x = zeros(size(b)); %初始解設置為與b同型的零向量

k = 0; %叠代次數的記數變量，初始量設為0

r = 1; %前後項之差的無窮範數

% % % % % % % % % % % % % % % %

D = diag(diag(A));

B = inv(D)*(D-A);

f = inv(D)*b;

% % % % % % % % % % % % % % % %

p = max(abs(eig(B))); %譜半徑大於等於1就不收斂

if p >= 1

'叠代法不收斂'

return

end

while r >e

x0 = x;

x = B*x0 + f;

k = k + 1;

r = norm (x-x0,inf);

end

'所求解為'

x

'叠代次數為'

k

自己以前編的。。。。

上一篇:寶馬7系底盤升降怎麽調

下一篇:i310105f參數