matplotlib.axes.Axes.csd

轴。csd ( x , y , NFFT = None , Fs = None , Fc = None , detrend = None , window = None , noverlap = None , pad_to = None , side = None , scale_by_freq = None , return_line = None , * , data =无，** kwargs )

绘制交叉光谱密度。

交叉谱密度\(P_{xy}\)Welch 的平均周期图法。向量x和y被分成 NFFT长度段。每个段由函数 detrend 去趋势并由函数window加窗。noverlap给出段之间重叠的长度。x和y的直接 FFT 的乘积在每个段上平均计算\(P_{xy}\)，通过缩放来校正由于开窗引起的功率损耗。

如果 len( x ) < NFFT或 len( y ) < NFFT，它们将被零填充到NFFT。

参数：

x, y一维数组或序列

包含数据的数组或序列。

Fs浮点数，默认：2

采样频率（每个时间单位的样本数）。它用于计算傅立叶频率freqs，以每时间单位的周期为单位。

窗口可调用或 ndarray，默认值：window_hanning

长度为NFFT的函数或向量。要创建窗口向量，请参见 window_hanning、window_none、numpy.blackman、numpy.hamming、 numpy.bartlett、scipy.signal、scipy.signal.get_window等。如果将函数作为参数传递，它必须将数据段作为参数并返回该段的窗口版本。

边{'默认'，'单面'，'双面'}，可选

返回频谱的哪一边。'default' 对于真实数据是单面的，对于复杂数据是双面的。'oneside' 迫使单方面的频谱回归，而 'twolateral' 则迫使双方回归。

pad_to int，可选

执行 FFT 时数据段填充到的点数。这可能与NFFT不同，后者指定使用的数据点数。虽然不会增加光谱的实际分辨率（可分辨峰之间的最小距离），但这可以在图中提供更多点，从而获得更多细节。这对应于调用中的nfft参数。默认为 None，将pad_to设置为NFFT

NFFT int，默认值：256

FFT 的每个块中使用的数据点数。幂 2 是最有效的。这不应该用于获得零填充，否则结果的缩放将不正确；使用pad_to代替。

detrend {'none', 'mean', 'linear'} 或可调用的，默认值：'none'

在 fft-ing 之前应用于每个段的函数，旨在消除均值或线性趋势。与 MATLAB 中的detrend参数是一个向量不同，在 Matplotlib 中它是一个函数。该mlab 模块定义detrend_none、detrend_mean和detrend_linear，但您也可以使用自定义函数。您还可以使用字符串来选择功能之一： 'none' calls detrend_none。“卑鄙”的电话 detrend_mean。'线性'调用detrend_linear。

scale_by_freq布尔值，默认值：True

生成的密度值是否应按缩放频率缩放，以 1/Hz 为单位给出密度。这允许对返回的频率值进行积分。对于 MATLAB 兼容性，默认值为 True。

noverlap int，默认值：0（无重叠）

段之间的重叠点数。

Fc整数，默认值：0

x的中心频率，它偏移绘图的 x 范围，以反映在采集信号、然后过滤和下采样到基带时使用的频率范围。

return_line布尔值，默认值：False

是否包括在返回值中绘制的线对象。

返回：

Pxy一维数组

交叉谱的值\(P_{xy}\)缩放之前（复值）。

频率一维数组

与Pxy中的元素对应的频率。

线Line2D

此函数创建的行。仅在return_line为 True 时返回。

其他参数：

数据可索引对象，可选

如果给定，以下参数也接受一个字符串s，它被解释为data[s]（除非这引发异常）：

x , y

**kwargs

关键字参数控制Line2D属性：

财产	描述
agg_filter	一个过滤器函数，它接受一个 (m, n, 3) 浮点数组和一个 dpi 值，并返回一个 (m, n, 3) 数组和距图像左下角的两个偏移量
alpha	标量或无
animated	布尔
antialiased或aa	布尔
clip_box	Bbox
clip_on	布尔
clip_path	补丁或（路径，变换）或无
color或 c	颜色
dash_capstyle	CapStyle或 {'butt', 'projecting', 'round'}
dash_joinstyle	JoinStyle或 {'miter', 'round', 'bevel'}
dashes	浮动序列（以点为单位的开/关墨水）或（无，无）
data	(2, N) 数组或两个一维数组
drawstyle或 ds	{'default', 'steps', 'steps-pre', 'steps-mid', 'steps-post'}，默认值：'default'
figure	Figure
fillstyle	{'full', 'left', 'right', 'bottom', 'top', 'none'}
gapcolor	颜色或无
gid	字符串
in_layout	布尔
label	目的
linestyle或 ls	{'-', '--', '-.', ':', '', (offset, on-off-seq), ...}
linewidth或 lw	漂浮
marker	标记样式字符串，Path或MarkerStyle
markeredgecolor或机械	颜色
markeredgewidth或者喵喵	漂浮
markerfacecolor或 mfc	颜色
markerfacecoloralt或 mfcalt	颜色
markersize或毫秒	漂浮
markevery	无或 int 或 (int, int) 或 slice 或 list[int] 或 float 或 (float, float) 或 list[bool]
mouseover	布尔
path_effects	AbstractPathEffect
picker	浮动或可调用[[艺术家，事件]，元组[布尔，字典]]
pickradius	未知
rasterized	布尔
sketch_params	（比例：浮动，长度：浮动，随机性：浮动）
snap	布尔或无
solid_capstyle	CapStyle或 {'butt', 'projecting', 'round'}
solid_joinstyle	JoinStyle或 {'miter', 'round', 'bevel'}
transform	未知
url	字符串
visible	布尔
xdata	一维数组
ydata	一维数组
zorder	漂浮

也可以看看

psd

相当于设置。y = x

笔记

对于绘图，功率被绘制为 \(10 \log_{10}(P_{xy})\)分贝，虽然\(P_{xy}\)本身被退回。

参考

Bendat & Piersol -- 随机数据：分析和测量程序，John Wiley & Sons (1986)

的示例matplotlib.axes.Axes.csd

惩教署演示

惩教署演示