ex1903.m

来自「在信号中傅里叶的变换以及其逆变换的相互转换」· M 代码 · 共 15 行

%例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)')