a13.m

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

m=101;
n=20;
t=linspace(-1,1,m);
y=rand(1,m);
c=polyfit(t,y,n-1);
A=zeros(m,n);
y=y';
for i=1:n
    for j=1:m
         A(j,i)=y(j)^(i-1);

    end
end
[U,S,V]=svd(A);
plot(U(:,2))
hold on
plot(U(:,3))