2.m

来自「学习使用各种数值方法计算数值代数问题。其中包括计算矩阵的特征值问题」· M 代码 · 共 23 行

clear all
m=101;
n=20;
t=linspace(-1,1,m);

for i=1:m
    for j=1:n
        y(i)=exp(t(i));
        A(i,j)=exp(t(i));
    end
end
y=y';
[U,S,V]=svd(A);
[Q,R]=qr(A);
for k=1:rank(A)
    a(k)=S(k,k);
end
semilogy(a)

%for i=1:101
%plot(A(:,i),t)
%hold on
%end