怎麼畫出Fs/2以上的頻譜??

看板MATLAB作者 (evanescent)時間13年前 (2013/03/27 11:55), 編輯推噓0(005)
留言5則, 2人參與, 最新討論串1/1
比如以下的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

03/27 12:18, , 1F
時間的解析度提高則頻率的範圍拉長
03/27 12:18, 1F

03/27 13:57, , 2F
如果你修過信號與系統你會知道fft是循環的
03/27 13:57, 2F

03/27 13:57, , 3F
而且對real signal, fft is even
03/27 13:57, 3F

03/27 13:59, , 4F
so, 0-Fs/2 is enough, though you can see something > Fs/2
03/27 13:59, 4F

03/27 18:50, , 5F
循環的前提是稀疏 整倍數刻度之外的值皆為零
03/27 18:50, 5F
文章代碼(AID): #1HKcspKR (MATLAB)
文章代碼(AID): #1HKcspKR (MATLAB)