exm0523_1.m

来自「这是本人收集的一些程序源代码」· M 代码 · 共 23 行

%exm0523_1.m
clear,clc,
N=[6 8 10 12 14];
for k=1:length(N)
	n=N(k);
	H=hilb(n);
	Hi=invhilb(n);
	b=ones(n,1);
	x_approx=H\b;
	x_exact=Hi*b;
	ndb=norm(H*x_approx-b);nb=norm(b);
	ndx=norm(x_approx - x_exact);nx=norm(x_approx);
	er_actual(k)=ndx/nx;
	K=cond(H);
	er_approx(k)=K*eps;
	er_max(k)=K*ndb/nb;
end  
disp('Hilbert矩阵阶数'),disp(N)
format short e
disp('实际误差  er_actual'),disp(er_actual),disp('')
disp('近似的最大可能误差  er_approx'),disp(er_approx),disp('')
disp('最大可能误差  er_max'),disp(er_max),disp('')