[問題] 使用xcorr與fft計算PSD
我用xcorr分別計算出了x1與x2的autocrrelation function
(因為某些需求 只能從ACF下手)
若R_x1 的lag從-n ~ +n
如果原本的x1屬於實數，R_x1會為even function，S_x1會屬於實數
但直接用fft計算PSD時
我發現PSD會有複數出現!!!
因為我需要計算s_x1.*s_x2
所以不希望受到虛部的影響
所以希望一開始產生PSD時就能夠準確計算出PSD且沒有虛部
請問我要怎麼修改程式讓虛部消失
謝謝
如果考慮fft(R_x1 lag從0~+N) 這樣可行嗎?
--
或者這樣問好了
如果我現有2個ACF R_x1 R_x2 與一個CCR R_x1x2
我要如何算出它們的coherence function??
因為一直有虛部出現 有時候還會超出1....
所以上來請教是否能給我一些建議
謝謝
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.144.35
※ 編輯: TsaoNiMa 來自: 1.169.150.94 (05/15 23:58)
→
05/16 10:14, , 1F
05/16 10:14, 1F
我在網路找到參考資料了
計算出來應該是對的
不過結果跟我預期的有落差Q___Q
※ 編輯: TsaoNiMa 來自: 140.113.144.35 (05/16 10:37)
※ 編輯: TsaoNiMa 來自: 140.113.144.35 (05/16 10:52)
→
05/16 13:03, , 2F
05/16 13:03, 2F
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章