fft_change.m

来自「快速傅立叶变换的程序」· M 代码 · 共 11 行

 N=128;
 T=1;
 t=linspace(0,T,N);
 x=12*sin(2*pi*10*t+pi/4)+5*cos(2*pi*40*t);
 dt=t(2)-t(1);
 f=1/dt;
 X=fft(x);
 F=X(1:N);
  f=f*(0:N-1)/N;
  
 plot(f,abs(F),'-*');