使用scipy.signal.csd进行跨谱密度估计

如何解决使用scipy.signal.csd进行跨谱密度估计

我第一次使用scipy.signal.csd，发现有些令人困惑。文档指出返回了交叉功率谱密度（两个输入x和y的Pxy）。这种描述使我认为应该返回实数而不是复数（因此，如果这是不正确的，那么问题就差不多了！）

看看源代码，csd调用_spectral_helper，其中包括以下内容：

if not same_data:
    # All the same operations on the y data
    result_y = _fft_helper(y,win,detrend_func,nperseg,noverlap,nfft,sides)
    result = np.conjugate(result) * result_y
elif mode == 'psd':
    result = np.conjugate(result) * result

如果将相同的数据传递给csd（Welch函数执行此操作），则“结果”（最后一行）显然是真实的且符合预期。如果数据不同，则使用行结果= np.conjugate（result）* result_y，它是一组复数。

因此，如上所述，如果这是众所周知的和显而易见的，那么可以。否则，np.abs（result）应该在_spectral_helper中，还是应该更改文档以使不熟练的用户明白这一点？