t2_7.m

来自「用matlab编写的追赶法程序」· M 代码 · 共 23 行

function[K]=T2_7(n,m)
A=hilb(n);
B=[];
if(m<1)
    B=eye(n);
end
if(m>=1)
    for j=1:n
         B(j,j)=A(j,j);
    end   
end
for i=1:(m-1)/2
    for j=1:n
        if(j-i>=1)
            B(j-i,j)=A(j-i,j);
        end
        if(j+i<=n)
             B(j+i,j)=A(j+i,j);
        end
    end
end
C=inv(B)*A;
K=norm(C)*norm(inv(C));