从信号处理角度进行分析

  简单的时间序列直接做各种谱分析频谱,包络谱,平方包络谱,功率谱,倒谱等等

  比如一些简单的旋转机械振动时间序列信号

  如果频谱不好分析,那可以分析如下图所示的时间序列的时频谱

  给个简单的模拟信号的例子

  其连续小波变换时频谱如下

  一个模拟的轴承内圈故障振动信号,带有明显的周期性

  相应的频谱如下,红色虚线代表故障特征频率及相应的倍频

  包络谱如下

  看一下相应的CWT时频谱,很明显能看出冲击性

  还可以试试小波相干与交叉小波分析

  小波相干交叉小波,可以很好地反映两个不同时间序列变化之间的相关性。小波相干分析,一般反映序列间周期性变化趋势的一致性,但不直接反映变化周期的强度关系;交叉小波分析,一般反映序列间共有周期的强度。

  此外,如果时频谱线能量发散,时频脊线模糊,还可以试试同步压缩之类的算法

  时间序列信号处理系列基于Python的同步压缩变换 哥廷根数学学派的文章 知乎 httpszhuanlan.zhihu.comp

单向轴承工作原理(时间序列如何分析周期性?)

  当时间序列信号中噪声较大时,为了有利于周期性分析,不可避免地要进行降噪前处理

  比如KSVD降噪

  样条框架降噪

  Morlet小波降噪

  当待分析的时间序列过于复杂时,那可能要引入模态分解多分辨分析,比如小波分解,经验模态分解及其变体,变分模态分解,经验小波变换,局部均值分解,辛几何模态分解,各种各样的自适应分解算法

  基于小波脊线的时间序列分解

  好多同学都对各种模态分解方法的时间序列处理感兴趣,那就随便说一下

  实际上,时间序列通常由多个具有物理意义的分量组成,在很多时候,为了更容易的研究信号,我们希望在与原始数据相同的时间尺度上单独研究这些分量中的一个或多个,理想情况下,我们希望这些经MRA分解到的多个分量在物理上是有意义的,可容易解释的。多分辨率分析MRA通常与小波或小波包相关联,但诸如经验模态分解EMD,变分模态分解VMD等模态分解方法也可以构成MRA。

  先给一个简单的合成信号,信号以Hz的频率采样秒钟。

单向轴承工作原理(时间序列如何分析周期性?)

  该信号由个主要分量组成频率为 Hz的时间局部振荡分量频率为 Hz的时间局部振荡分量和趋势项分量。趋势项分量为正弦曲线,频率为.Hz。Hz的振荡分量发生在 .到 . 秒之间,而 Hz的振荡分量发生在 . 到 秒之间。

  但这些分量从时域波形中无法分辨,因此进行频域变换。

  从频率中可以更容易地辨别振荡分量的频率,但时间局部性信号却丢失。为了同时定位时间和频率信息,使用连续小波变换进行分析。

  从CWT时频谱图中可以看出Hz和Hz分量的时间范围,但没有发现趋势项分量。为了分离出信号的分量并单独进行分析,接下来使用多分辨分析,直接在时域中进行相关操作。

  多分辨分析通过将信号分成不同分辨率的分量进而缩小分析范围,而提取不同分辨率的信号分量相当于分解数据在不同时间尺度上的变化,或等效地在不同频带上进行分析。首先,采用离散小波变换的变体最大重叠离散小波变换对信号进行多分辨分析,分解层数为。关于最大重叠离散小波变换的相关内容,请查看如下文献。

  最大重叠离散小波变换的层多分辨分析分解如下

  如果从上向下看,会看到所分解的分量变得越来越平滑,即分量频率越来越低。回想一下,原始信号包含个主要分量,一个 Hz 的高频振荡成分一个 Hz 的低频振荡成分和一个趋势成分,它们都被加性噪声破坏了。

  从D 图中可以看出时间局部化的高频分量被分解出来,而下面的两个图包含较低频率的振荡分量,这是多分辨率分析的一个重要方面,最后S子图中包含了趋势项分量。

  除了小波多分辨分析,经验模态分解 EMD 是一种所谓的数据自适应多分辨技术。 EMD 在不使用固定基函数的情况下递归地从数据中提取不同的分辨率成分,关于EMD相关文献浩如烟海,不做赘述了。EMD的多分辨分析分解如下所示

  虽然MRA分解分量的数目不同,但 EMD MRA和小波 MRA会产生相似的信号波形,在 EMD MRA分解中,高频振荡成分位于第个本征模态函数中 IMF,低频振荡成分主要位于IMF和IMF中,IMF 中的趋势项分量与小波技术提取的趋势分量非常相似。

  自适应多分辨分析的另一种技术是变分模态分解 VMD,VMD 从信号中提取固有模式函数或振荡模式,并不使用固定基函数进行分析。EMD在时域上递归,以逐步提取低频IMF分量,而VMD 首先识别频域中的信号峰值并同时提取所有模式,相关文献如下

  Dragomiretskiy, Konstantin, and Dominique Zosso. Variational Mode Decomposition. IEEE Transactions on Signal Processing , no. February –. httpsdoi.org.TSP...

  VMD的多分辨分析分解如下所示

  由上图可知,与小波和EMD类似,VMD将个分量基本分离了出来。

  还有一种数据自适应多分辨分析技术经验小波变换 EWT ,EWT根据分析信号的频率构造 Meye小波进而进行自适应小波,之前写过EWT相关的内容

  经验小波变换在信号处理及轴承故障诊断中的应用 哥廷根数学学派的文章 知乎httpszhuanlan.zhihu.comp

  EWT的多分辨分析分解如下所示

  与之前的EMD和小波MRA类似,EWT分解出了相关的振荡分量,用于执行分析的滤波器及其通带信息如下

  下面考虑一段神户地震信号,源于 年 月 日在澳大利亚霍巴特的塔斯马尼亚大学记录,从 GMT 开始,以 秒的间隔持续 分钟。

  以最大重叠离散小波变换为例,其层MRA分解如下

  从D和D子图中可以看出初级与延迟次级波分量,地震波中的分量以不同的速度传播,初级波比次级剪切波传播的更快。

  将信号分解为若干分量的目的通常是去除某些分量以减轻对信号分析的影响,MRA技术的关键是重建原始信号的能力,如下

  每种方法的最大重建误差约为 或更小,表明它们可以完美地对信号进行重建。在很多研究中我们对趋势项不感兴趣,由于趋势项一般位于最后一个 MRA 分解分量中,因此只需将该分量去除,然后进行重建。

  此外,再删除第个MRA分解分量看起来主要是噪声

  在前面我们将趋势项删除,然而在许多应用中,趋势项可能是我们的主要研究部分,因此可视化几种MRA方法所提取的趋势项分量。

  根据