fs=100; x=-2:1/fs:2; y=sin(3*pi*x); z=rectpuls(x); figure;plot(x,y,x,z,':r');
my=abs(fft(y)); mz=abs(fft(z)); my=my/max(my); %归一化 mz=mz/max(mz); %归一化
f=(0:1/length(x):1)*fs; figure;plot(f(1:fs/2),my(1:fs/2),f(1:fs/2),mz(1:fs/2),':r'); %为什么这么做? %还记得Nyquist采样定理吗?哎~~你对了,就是这个原因。 %频谱图同理。
my=20*log10(my+eps);mz=20*log10(mz+eps); figure;plot(f(1:fs/2),my(1:fs/2),f(1:fs/2),mz(1:fs/2),':r');
%图形见全文
信号时域波形:
傅立叶变换的幅度图:
频谱图:
(责任编辑:泉水) |