ft.m

来自「在matlab环境下的向前和向后euler算法价偏微分方程问题」· M 代码 · 共 14 行

function ft(N)
%x(n)为输入信号;N为所取点的个数
y=zeros(1,N);
for n=1:N
    x(n)=cos(0.48*pi*n)+cos(0.52*pi*n);
end
for k=1:N
    for m=1:N
        y(k)=y(k)+x(m)*exp(-j*2*pi*(k-1)*(m-1)/N);
        z(k)=norm(y(k));
    end
end
plot(z)