a24.m

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

m=101;
n=20;
t=linspace(-1,1,m);
e=10^(-12);
y=exp(t)+e*randn(1,m);
c=polyfit(t,y,n-1);
y=y';
A=zeros(m,n);
for i=1:n
    for j=1:m
            A(j,i)=y(j)^(i-1);

    end
end
[U,S,V]=svd(A);
b=U'*y;

plot(b)