func12.asv

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

m=[5,-3,2;6,-4,4;4,-4,5]
maxtimes=1000
a=m
start=now
for i=1:maxtimes
   [Q,R]=qr(a)
   b=a
   a=R*Q
   if (abs(norm(a,inf)-norm(b,inf))<1e-4)
       disp('iteration times:')
       disp(i)
       break
   end
end
finish=now
disp('Elapse time:')
disp(finish
  
%a =

%    2.9999   -1.0001  -12.0002
%   -0.0001    2.0001   -1.9989
%    0.0000   -0.0000    1.0000

%iteration times:
%    22