Python cusignal.spectral_analysis.spectral.csd用法及代码示例

用法: cusignal.spectral_analysis.spectral.csd(x, y, fs=1.0, window='hann', nperseg=None, noverlap=None, nfft=None, detrend='constant', return_onesided=True, scaling='density', axis=- 1, average='mean')

使用 Welch 方法估计交叉功率谱密度 Pxy。

参数：

x：array_like: 测量值的时间序列
y：array_like: 测量值的时间序列
fs：浮点数，可选: x 和 y 时间序列的采样频率。默认为 1.0。
window：str 或 tuple 或数组，可选: 想要使用的窗口。如果window是字符串或元组，则传递给get_window生成窗口值，默认为DFT-even。有关窗口和所需参数的列表，请参阅get_window。如果window 是数组，它将直接用作窗口，其长度必须为nperseg。默认为 Hann 窗口。
nperseg：整数，可选: 每个段的长度。默认为 None，但如果 window 是 str 或 tuple，则设置为 256，如果 window 是数组，则设置为窗口的长度。
noverlap: int, optional：: 段之间重叠的点数。如果 None , noverlap = nperseg // 2 。默认为 None 。
nfft：整数，可选: 如果需要零填充 FFT，则使用的 FFT 的长度。如果 None ，则 FFT 长度为 nperseg 。默认为 None 。
detrend：str 或 function 或 False ，可选: 指定如何去除每个段的趋势。如果detrend 是一个字符串，它作为type 参数传递给detrend 函数。如果它是一个函数，它接受一个段并返回一个去趋势的段。如果 detrend 是 False ，则不进行去趋势。默认为‘constant’。
return_onesided：布尔型，可选: 如果 True ，返回真实数据的 one-sided 频谱。如果 False 返回 two-sided 频谱。默认为 True ，但对于复杂数据，始终返回 two-sided 频谱。
scaling：{ ‘density’, ‘spectrum’ }，可选: 选择计算交叉谱密度 (‘density’)，其中 Pxy 的单位为 V**2/Hz，计算交叉谱 (‘spectrum’)，其中 Pxy 的单位为 V**2，如果 x 和 y 以 V 为单位，fs 以 Hz 为单位。默认为‘density’
axis：整数，可选: 为两个输入计算 CSD 的轴；默认值在最后一个轴上(即 axis=-1 )。
average：{ ‘mean’, ‘median’ }，可选: 平均周期图时使用的方法。默认为‘mean’。

f：ndarray: 采样频率数组。
Pxy：ndarray: x,y 的交叉谱密度或交叉功率谱。

注意：

按照惯例，Pxy 的计算方法是 X 的共轭 FFT 乘以 Y 的 FFT。

如果输入序列的长度不同，则较短的序列将被零填充以匹配。

适当的重叠量取决于窗口的选择和您的要求。对于默认的 Hann 窗口，50% 的重叠是在准确估计信号功率和不过度计算任何数据之间的合理权衡。较窄的窗口可能需要较大的重叠。

参考：

1: P. Welch, “The use of the fast Fourier transform for the estimation of power spectra: A method based on time averaging over short, modified periodograms”, IEEE Trans. Audio Electroacoust. vol. 15, pp. 70-73, 1967.
2: Rabiner, Lawrence R., and B. Gold. “Theory and Application of Digital Signal Processing” Prentice-Hall, pp. 414-419, 1975

例子：

>>> import cusignal
>>> import cupy as cp
>>> import matplotlib.pyplot as plt

生成两个具有一些共同特征的测试信号。

>>> fs = 10e3
>>> N = 1e5
>>> amp = 20
>>> freq = 1234.0
>>> noise_power = 0.001 * fs / 2
>>> time = cp.arange(N) / fs
>>> b, a = cusignal.butter(2, 0.25, 'low')
>>> x = cp.random.normal(scale=cp.sqrt(noise_power), size=time.shape)
>>> # lfilter not currently implemented in cuSignal
>>> y = signal.lfilter(b, a, x)
>>> x += amp*cp.sin(2*cp.pi*freq*time)
>>> y += cp.random.normal(scale=0.1*cp.sqrt(noise_power), size=time.shape)

计算并绘制交叉谱密度的大小。

>>> f, Pxy = cusignal.csd(x, y, fs, nperseg=1024)
>>> plt.semilogy(cp.asnumpy(f), cp.asnumpy(cp.abs(Pxy)))
>>> plt.xlabel('frequency [Hz]')
>>> plt.ylabel('CSD [V**2/Hz]')
>>> plt.show()

相关用法

注：本文由纯净天空筛选整理自rapids.ai大神的英文原创作品 cusignal.spectral_analysis.spectral.csd。非经特殊声明，原始代码版权归原作者所有，本译文未经允许或授权，请勿转载或复制。

用法:

参数：

返回：

注意：

参考：

例子：