1.m

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

clear all
m=101;
n=20;
t=linspace(-1,1,m);
alpha=rand(n,1);
lambda=rand(n,1);
for i=1:m
    for j=1:n
        y(i)=alpha(j)*exp(-lambda(j)*t(i));
        A(i,j)=exp(-lambda(j)*t(i));
    end
end
y=y';
r=y-A*alpha;
[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