频率提取（FFT）的 Matlab 实现，代码如下clear %编写骆遥 fs=1000 t=0:1/fs:0.6; f1=100; f2=300; x=sin(2*pi*f1*t)+sin(2*pi*f2*t); subplot(711) plot(x); title(f1（100Hz）f2(300Hz)的正弦信号，初相 0) xlabel(序列（n）) grid onnumber=512y=fft(x,number); n=0:length(y)-1; f=fs*n/length(y); subplot(713) plot(f,abs(y); title(f1f2 的正弦信号的 FFT（512 点）) xlabel(频率 Hz) grid onx=x+randn(1,length(x); subplot(715) plot(x); title(原 f1f2 的正弦信号（含随机噪声）) xlabel(序列（n）) grid ony=fft(x,number); n=0:length(y)-1;f=fs*n/length(y); subplot(717) plot(f,abs(y); title(原 f1f2 的正弦信号（含随机噪声）的 FFT（512 点）) xlabel(频率 Hz) grid on