gusucode.com > 《精通matlab 7》包括配套光盘 > 配套光盘里的程序附件/Ex-19/Ex1903.m

    %例19-3：傅立叶变换
t = 0:0.001:0.6;                         %采样频率为1000Hz
x = sin(2*pi*50*t)+sin(2*pi*120*t);    %产生正弦波
y = x + 2*randn(size(t));             %叠加随机噪声
plot(1000*t(1:50),y(1:50))            %绘制y的曲线
title('信号受零均值的随机噪声干扰')
xlabel('时间(毫秒)')
%...
Y = fft(y,512);         %对y做傅立叶变换，取512个点
Pyy = Y.* conj(Y) / 512;   
f = 1000*(0:256)/512;   %设置频率轴（横轴）坐标，1000是采样频率
plot(f,Pyy(1:257))       %绘制频域图
title('y的频率成分')
xlabel('频率 (Hz)')