func2p3.m

来自「数值分析最常用的四十种算法」· M 代码 · 共 12 行

%func2p3.m 验证病态性
H=[1,1/2,1/3;1/2,1/3,1/4;1/3,1/4,1/5]
b=[11/6,13/12,47/60]'
b2=[1.83,1.08,0.783]'
x=[1,1,1]'
x=H\b
x=H\b2
norm(inv(a),inf)*norm(a,inf)
%b轻微变化后（各分量减少0.003333），解已变为(1.0800,0.5400,1.4400)
%原矩阵有严重的病态性;cond(infinite)(a)=3.4077,可见右边的扰动将对解造成重大影响