GNU倍频程极短信号的频率分析

Question

我有一些非常短的信号，从示波器(50k-200 k样本)注册超过2ms的时间长度。这些是带有静电放电(静电放电)火花信号的声信号。我想得到该信号的一些频率数据，在近声频率范围内(高达30 that )，具有尽可能高的时间分辨率。

我试过绘制一个光谱图(八度图)来查看信号，但是输出并没有真正的用处。使用specgram( x, N, fs );，其中x是我的fs采样率信号，当N的低值时，我接收从大约500 the的高频开始的图形，对于大N值(如2^12-13)，获得更好的频率分辨率，但是窗口太宽，在整个信号长度上只接收到2个频谱值。

据我所知，这可能是傅立叶变换的局限性，它很可能是由散斑函数使用的(实际上，我对信号分析不太了解)。

还有其他方法来获取这种信号的频率(随时间的变化)信息吗？我读过一些关于小波的内容，但是当我尝试使用信号包的dwt函数时，我收到了这样的错误：

error: 'wfilters' undefined near line 51 column 14
error: called from
    dwt at line 51 column 12

即使这样可行，我也不太确定我是否知道如何使用这些小波函数的输出.

Answer 1

要从如此高的采样率中获取音频信息，您将需要获得一个样本向量，其长度足以包含音频频率上的至少几个完整的周期，例如，许多10毫秒的连续样本，这可能超出了您的范围，也可能不超过您的范围。为了合理地处理这一数量的数据，您可能会对样本数据进行低通滤波，使其只包含音频，然后将其重采样到较低的采样率，但超过滤波器截止频率的两倍。然后，你将得到一个更短的样本向量来为你的音频频谱分析提供一个FFT。