scipy.signal.csd¶
- scipy.signal.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。
- 参数
- xarray_like
测量值的时间序列
- yarray_like
测量值的时间序列
- fs浮动，可选
的采样频率 x 和 y 时间序列。默认为1.0。
- window字符串或元组或array_like，可选
要使用的所需窗口。如果 window 是字符串或元组，则将其传递给
get_window以生成窗口值，该窗口值是DFT-即使在默认情况下也是如此。看见get_window有关窗口和所需参数的列表，请执行以下操作。如果 window is array_like它将直接用作窗口，其长度必须为nperseg。默认为Hann窗口。- nperseg整型，可选
每个线段的长度。默认值为None，但如果Window为字符串或元组，则设置为256，如果Window为array_like，则设置为窗口长度。
- noverlap：int，可选
线段之间要重叠的点数。如果 None ，
noverlap = nperseg // 2。默认为 None 。- nfft整型，可选
如果需要填充零的FFT，则使用的FFT长度。如果 None ，FFT长度为 nperseg 。默认为 None 。
- 下降趋势 ：字符串或函数或 False ，可选字符串或函数或
指定如何对每段数据段进行趋势调整。如果
detrend是字符串，则将其作为 type 参数设置为detrend功能。如果它是一个函数，它接受一个段并返回一个去势的段。如果detrend是 False ，没有进行去趋势性操作。默认值为“Constant”。- return_onesided布尔值，可选
如果 True ，返回真实数据的单边频谱。如果 False 返回双面光谱。默认为 True ，但是对于复杂的数据，总是返回双面频谱。
- scaling{‘密度’，‘光谱’}，可选
在计算交叉谱密度(‘密度’)之间进行选择，其中 Pxy 单位为V 2/Hz and computing the cross spectrum ('spectrum') where `Pxy` has units of V 2，如果 x 和 y 是以V和V为单位测量的 fs 是以赫兹为单位测量的。默认为“密度”
- axis整型，可选
沿其计算两个输入的CSD的轴；默认值在最后一个轴上(即
axis=-1)。- average{‘均值’，‘中位数’}，可选
平均周期图时使用的方法。如果谱是复数，则分别计算实部和虚部的平均值。默认为“Mean”。
1.2.0 新版功能.
- 退货
- fndarray
采样频率数组。
- Pxyndarray
x，y的交叉谱密度或交叉功率谱。
参见
periodogram简单的、可选修改的周期图
lombscargle非均匀采样数据的Lomb-Scarger周期图
welch功率谱密度按Welch方法计算。 [等价于CSD(x，x)]
coherence用韦尔奇的方法进行幅度平方相干性计算。
注意事项
按照惯例，PXY是用X的共轭FFT乘以Y的FFT来计算的。
如果输入序列的长度不同，将对较短的序列进行补零以匹配。
适当的重叠量将取决于窗口的选择和您的要求。对于默认的汉恩窗口，50%的重叠是准确估计信号功率，同时不过度计算任何数据之间的合理折衷。较窄的窗口可能需要较大的重叠。
0.16.0 新版功能.
参考文献
- 1
P.Welch，“快速傅立叶变换在功率谱估计中的应用：一种基于短周期图时间平均的方法”，IEEE Transans。音频电声。卷。15年，第70-73页，1967年。
- 2
题名/责任者：Realthcare/by J.“数字信号处理理论与应用”，Prentice-Hall，第414-419页，1975
示例
>>> from scipy import signal >>> import matplotlib.pyplot as plt >>> rng = np.random.default_rng()
生成具有一些共同特征的两个测试信号。
>>> fs = 10e3 >>> N = 1e5 >>> amp = 20 >>> freq = 1234.0 >>> noise_power = 0.001 * fs / 2 >>> time = np.arange(N) / fs >>> b, a = signal.butter(2, 0.25, 'low') >>> x = rng.normal(scale=np.sqrt(noise_power), size=time.shape) >>> y = signal.lfilter(b, a, x) >>> x += amp*np.sin(2*np.pi*freq*time) >>> y += rng.normal(scale=0.1*np.sqrt(noise_power), size=time.shape)
计算并绘制交叉谱密度的大小。
>>> f, Pxy = signal.csd(x, y, fs, nperseg=1024) >>> plt.semilogy(f, np.abs(Pxy)) >>> plt.xlabel('frequency [Hz]') >>> plt.ylabel('CSD [V**2/Hz]') >>> plt.show()
