范数是应用
在学习雅可比迭代法中,实现雅可比迭代法时:
%fileName:jacobi.m:
function [y,n]=jacobi(A,b,x0,ep)
if nargin==3 %当调用函数时没有输入精度要求,自动设置为1e-8
ep=1e-8;
end
D=diag(diag(A));
L=-tril(A,-1);
U=-triu(A,1);
B=D\(L+U);
f=D\b;
y=B*x0+f;
n=1;
while norm(y-x0)>=ep %这里不能用y-x0>=ep
x0=y;
y=B*x0+f;
n=n+1;
end
第13行的norm是求y-x0范数,其实就是两者之间的距离,不能直接用两者相减,主要原因是,两者相减时,只要两者中有一个方向值相同,则条件就不能成立了,就无法循环了。
人就像是被蒙着眼推磨的驴子,生活就像一条鞭子;当鞭子抽到你背上时,你就只能一直往前走,虽然连你也不知道要走到什么时候为止,便一直这么坚持着。
浙公网安备 33010602011771号