作者gggould (evanescent)
看板MATLAB
標題怎麼畫出Fs/2以上的頻譜??
時間Wed Mar 27 11:55:28 2013
比如以下的code
Fs = 1000;
t = 0:1/Fs:1-1/Fs;
x = cos(2*pi*100*t)+randn(size(t));
plot(psd(spectrum.periodogram,x,'Fs',Fs,'NFFT',length(x)));
有什麼方法可以畫出Fs/2以上的頻譜嗎
主要想看一下信號通過升頻或降頻digital filter 後 剩下的aliasing effect
GOOGLE後找不到
請問有這方面的example嗎
thanks a lot..
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 67.169.80.46
1F:→ labbat:時間的解析度提高則頻率的範圍拉長 03/27 12:18
2F:→ jatj:如果你修過信號與系統你會知道fft是循環的 03/27 13:57
3F:→ jatj:而且對real signal, fft is even 03/27 13:57
4F:→ jatj:so, 0-Fs/2 is enough, though you can see something > Fs/2 03/27 13:59
5F:→ labbat:循環的前提是稀疏 整倍數刻度之外的值皆為零 03/27 18:50