c51.m

来自「matlab例程MATLAB 6.x版的内建数据类型（Built-in data」· M 代码 · 共 11 行

T=0.015625;T0=4;N0=T0/T;
t=0:T:T*(N0-1);t=t';
f=T*exp(-2*t);
f(1)=T*0.5;
F=fft(f);
[Fp,Fm]=cart2pol(real(F),imag(F));
k=0:N0-1; k=k';
w=2*pi*k/T0;
subplot(211),plot(w(1:128),Fm(1:128))
subplot(212),plot(w(1:128),Fp(1:128))