matlab的gccphat转C语言,Python中的GCCPHAT互相关

我试图用python实现GCC-PHAT。在

该方法类似于以下两个环节：

link1和link2

GCC-PHAT和使用FFT的正常互相关之间的唯一区别似乎是除以幅度。在

这是我的代码：import numpy as np

import matplotlib.pyplot as plt

from scipy.fftpack import rfft, irfft, fftfreq, fft, ifft

def xcorr_freq(s1,s2):

pad1 = np.zeros(len(s1))

pad2 = np.zeros(len(s2))

s1 = np.hstack([s1,pad1])

s2 = np.hstack([pad2,s2])

f_s1 = fft(s1)

f_s2 = fft(s2)

f_s2c = np.conj(f_s2)

f_s = f_s1 * f_s2c

denom = abs(f_s)

denom[denom < 1e-6] = 1e-6

f_s = f_s / denom # This line is the only difference between GCC-PHAT and normal cross correlation

return np.abs(ifft(f_s))[1:]

我通过注释fs = fs / denom检查了这个函数产生的结果与宽带信号的正常互相关相同。在

下面是一个示例测试代码，显示上面的GCC-PHAT代码的性能比正常的互相关差：

^{pr2}$

以下是GCC-PHAT的结果：

以下是正常互相关的结果：

由于GCC-PHAT应该能为宽带信号提供更好的互相关性能，我知道我的代码有问题。非常感谢任何帮助！在