a23.m

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

m=101;
n=20;
t=linspace(-1,1,m);
for k=(-16):4:0
    e=10^k;
    y=exp(t)+e*randn(1,m);
    y=y';
    plot(y)
    hold on
end