如何使用光谱图和/或周期图绘制PSD？

Question

我有一个来自加速度计的m/s^2数据的<1,000,000x1 double>矢量。我想创建此数据的功率谱密度图(功率/赫兹与频率)。我可以通过频谱图运行向量来找到频率，但我可能不正确地使用窗口。我得到一个<131073x8复杂的double>返回，但我不知道它的内容是什么。我可以简单地使用矢量上的函数周期图来绘制功率/hz与频率的关系图吗？我试过了，matlab返回一个空白图。

对此的支持将非常感谢..谢谢你，亚历克斯

Answer 1

如果您只需要功率谱密度图，您可以运行

periodogram(x)

你会得到归一化频率上的功率谱密度(单位为dB)的曲线图。如果希望密度高于实际频率(以Hz为单位)，则需要指定数据的采样率。语法是

periodogram(x, [], [], Fs)

其中Fs是采样频率(以Hz为单位)。然后，您会得到一个横轴从0到Fs / 2 (奈奎斯特频率)的曲线图。

作为[]给出的两个参数可用于指定周期图法的窗口长度和用于底层快速傅立叶变换的数据点数量。您可以使用这些参数，特别是窗口长度，以调节频谱的频率分辨率和估计功率谱密度的精度之间的权衡。

使用这些语法时，periodogram不会返回任何输出。但是，您可以使用以下语法

[P, f] = periodogram(x, [], [], fs);

它以P为单位返回功率谱密度的估计值，并以f为单位返回相应的频率(以Hz为单位)。在此基础上，您可以使用以下命令生成一个类似于periodogram生成的无输出的图

semilogy(f, P)

即用垂直轴的对数刻度在f上绘制P。