>> fs=250;N=128;
n=-N/2:1:N/2;t=n/fs;
f=sin(40*t).*cos(140*2*pi*t);
figure(3),stem(t,f);
xlabel('t');ylabel('f(n)');
title('采样信号波形图');
grid;
F=abs(fftshift(fft(f,N+1)));
figure(4),plot(n,F);
xlabel('n');
title('采样频谱波形图');
grid;