作者TsaoNiMa (絕對不酸郭嚴文)
看板MATLAB
標題[問題] 使用xcorr與fft計算PSD
時間Wed May 15 22:46:29 2013
我用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)
1F：→ coolfox0719:加上絕對值再算呢? 05/16 10:14
我在網路找到參考資料了
計算出來應該是對的
不過結果跟我預期的有落差Q___Q
※ 編輯: TsaoNiMa 來自: 140.113.144.35 (05/16 10:37)
※ 編輯: TsaoNiMa 來自: 140.113.144.35 (05/16 10:52)
2F：→ labbat:離散偶函數不可能有中點 -T~T的0或1~2T的0.5T+0.5都不存在 05/16 13:03