请看代码：

1 clf;fs=100; %采样频率 2 Ndata=128; %数据长度 3 N=256; %FFT的数据长度 4 n=0:Ndata-1;t=n/fs;   %数据对应的时间序列 5 %x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t);   %时间域信号 6 x=0.5*sin(2*pi*20*t); 7 y=fft(x,N);   %信号的Fourier变换 8 mag=abs(y);    %求取振幅 9 f=(0:N-1)*fs/N; %真实频率10 % subplot(2,2,1),11 plot(f(1:N/2),mag(1:N/2)*2/N); %绘出Nyquist频率之前的振幅12 xlabel('频率/Hz');ylabel('振幅');13 title('Ndata=128 Nfft=256');grid on;

请看图：