ft1.m

来自「我们的作业」· M 代码 · 共 11 行

   t=0:1/99:1;
   x=sin(2*pi*15*t)+sin(2*pi*40*t);
   y=fft(x);
   m=abs(y);
   p=unwrap(angle(y));
   f=(0:length(y)-1)*99/length(y);
   subplot(1,2,1),plot(f,m);
   set(gca,'Xtick',[15,40,60,85]);
   subplot(1,2,2),plot(f,p*180/pi);
   set(gca,'Xtick',[15,40,60,85]);