作者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