弟哥素来。。

N=128;n=0:N-1;t=0.01*n;
x=2*sin(4*pi*t)+5*cos(16*pi*t);
subplot(2,1,1)
>> N=128;n=0:N-1;t=0.01*n;
x=2*sin(4*pi*t)+5*cos(16*pi*t);
subplot(2,1,1);stem(n,x,'.');title('signal x(n)');grid;
X=fft(x,N);
magx=abs(X(1:1:N/2+1));k=[0:1:N/2];w=2*pi/N*k;
subplot(2,1,2);plot(w/pi,magx);title('FFT N=128');
xlabel('频率(单位：pi)');ylabel('|X|');grid