这是用户在 2025-4-12 24:27 为 https://app.immersivetranslate.com/pdf-pro/3280963b-6d89-4c39-927a-ae1415634196/ 保存的双语快照页面,由 沉浸式翻译 提供双语支持。了解如何保存?

基于深度学习的振动信号自适应模态分解与瞬时频率估计

郭一冰 a , b , c a , b , c ^(a,b,c){ }^{a, b, c} ,鲍跃全 a , b , c , a , b , c , ^(a,b,c,**){ }^{a, b, c, *} ,李惠 a,b,c a,b,c  ^("a,b,c "){ }^{\text {a,b,c }} ,张宇峰 d,e d,e  ^("d,e "){ }^{\text {d,e }} a ^("a "){ }^{\text {a }} 哈尔滨工业大学,工业和信息化部智能防灾减灾重点实验室,哈尔滨 150090,中国 b b ^(b){ }^{\mathrm{b}} 哈尔滨工业大学,结构工程灾变与控制教育部重点实验室,哈尔滨 150090,中国 c ^("c "){ }^{\text {c }} 哈尔滨工业大学土木工程学院,哈尔滨 150090,中国 d ^("d "){ }^{\text {d }} 在役长大桥梁安全与健康国家重点实验室,江苏 211112,中国 e ^("e "){ }^{\text {e }} 苏交科集团,江苏 211112,中国

  文章信息

  关键词:

  时频分析
非线性和非平稳
自适应模态分解
  单分量
  瞬时频率
  深度学习

  摘要

在时变环境和载荷下,结构的振动监测数据呈现多频率成分的非平稳特性,极端载荷下还可能表现出显著非线性。时频分析方法非常适合分析这类非线性和非平稳信号。然而,大型结构的振动信号通常包含多种频率成分,使得现有自适应模态分解方法难以有效实现模态分离,且易受虚假单分量干扰。虚假单分量会严重降低瞬时频率的估计精度。针对模态分解中的模态混叠和虚假单分量问题,本文提出了一种基于深度学习的信号自适应模态分解与瞬时频率估计数据驱动方法。通过将时频分析问题转化为深度神经网络的优化任务,设计了内置冗余自适应滤波器与稀疏正则化机制,从而能准确估计本征模态函数和时变频率。 所提方法通过数值非平稳信号及长缆索模型实验振动信号验证。结果表明,该方法能自适应地将信号分解为单分量且无虚假分量,同时实现对应瞬时频率的高精度估计。在频谱重叠的无噪声信号案例中,瞬时频率估计精度达 98.52 % 98.52 % 98.52%98.52 \% 。该方法对噪声更具鲁棒性,在信噪比为 5 dB 时,瞬时频率识别误差仅为 4.90 % 4.90 % 4.90%4.90 \% ,远低于 EEMD 与 VMD 方法。

  1. 引言

结构健康监测(SHM)是理解工程结构演化规律、保障结构安全最直接先进的方法,已在众多大型结构中得到广泛应用与研究[1-3]。监测数据的有效分析是结构参数识别、损伤检测与状态评估的基础[4]。在时...
  术语表
A   自适应滤波器
B   傅里叶变换矩阵
C   中心矩阵
D   带宽矩阵
e e ee 相对估计误差
ε ε epsi\varepsilon   少量
f f ff   合成信号
f ( t ) f ( t ) f(t)\mathbf{f}(t) 结构振动信号
f ^ r e c ( ω ) f ^ r e c ( ω ) widehat(f)_(rec)(omega)\widehat{\mathbf{f}}_{r e c}(\omega) 重构信号频谱
F F F\mathscr{F}   傅里叶变换
F   非线性变换
F   实测索力
F ^ F ^ widehat(F)\widehat{\mathrm{F}} 信号的频谱矩阵
F ^ IMF F ^ IMF  widehat(F)_("IMF ")\widehat{\mathbf{F}}_{\text {IMF }} 本征模态函数的频谱矩阵
F ^ red F ^ red  widehat(F)_("red ")\widehat{\mathbf{F}}_{\text {red }} 冗余单分量矩阵
H   滤波器组
IF   估计缆索力
I r , k I r , k I_(r,k)I_{r, k}   正交值
K   滤波器数量
L L LL   损失函数
L c L c L_(c)L_{c}   电缆长度
m m mm   质量密度
M 本征模态函数数量
N N NN   数据点数量
r ( t ) r ( t ) r(t)\mathbf{r}(t)   残差
u k ( t ) u k ( t ) u_(k)(t)\mathbf{u}_{k}(t)   本征模态函数
ω ~ ( t ) ω ~ ( t ) widetilde(omega)(t)\widetilde{\boldsymbol{\omega}}(t)   瞬时频率
W   稀疏权重矩阵
z k ( t ) z k ( t ) z_(k)(t)\mathrm{z}_{\mathrm{k}}(t)   解析信号
Z   解析信号矩阵
Z ^ Z ^ widehat(Z)\widehat{\mathbf{Z}} 解析信号的频谱矩阵
Nomenclature A adaptive filters B Fourier transform matrix C center matrix D bandwidth matrix e relative estimation error epsi a small number f synthetic signal f(t) structural vibration signal widehat(f)_(rec)(omega) spectrum of reconstructed signal F Fourier transform F nonlinear transformation F measured cable force widehat(F) spectrum matrix of signal widehat(F)_("IMF ") spectrum matrix of intrinsic mode functions widehat(F)_("red ") redundant mono-component matrix H filter set IF estimated cable force I_(r,k) orthogonal value K number of filters L loss function L_(c) cable length m mass density M number of intrinsic mode functions N number of data points r(t) residual u_(k)(t) intrinsic mode function widetilde(omega)(t) instantaneous frequency W sparse weights matrix z_(k)(t) analytic signal Z analytic signal matrix widehat(Z) spectrum matrix of analytic signal| Nomenclature | | | :---: | :---: | | A | adaptive filters | | B | Fourier transform matrix | | C | center matrix | | D | bandwidth matrix | | $e$ | relative estimation error | | $\varepsilon$ | a small number | | $f$ | synthetic signal | | $\mathbf{f}(t)$ | structural vibration signal | | $\widehat{\mathbf{f}}_{r e c}(\omega)$ | spectrum of reconstructed signal | | $\mathscr{F}$ | Fourier transform | | F | nonlinear transformation | | F | measured cable force | | $\widehat{\mathrm{F}}$ | spectrum matrix of signal | | $\widehat{\mathbf{F}}_{\text {IMF }}$ | spectrum matrix of intrinsic mode functions | | $\widehat{\mathbf{F}}_{\text {red }}$ | redundant mono-component matrix | | H | filter set | | IF | estimated cable force | | $I_{r, k}$ | orthogonal value | | K | number of filters | | $L$ | loss function | | $L_{c}$ | cable length | | $m$ | mass density | | M | number of intrinsic mode functions | | $N$ | number of data points | | $\mathbf{r}(t)$ | residual | | $\mathbf{u}_{k}(t)$ | intrinsic mode function | | $\widetilde{\boldsymbol{\omega}}(t)$ | instantaneous frequency | | W | sparse weights matrix | | $\mathrm{z}_{\mathrm{k}}(t)$ | analytic signal | | Z | analytic signal matrix | | $\widehat{\mathbf{Z}}$ | spectrum matrix of analytic signal |
在多变环境与荷载作用下,结构振动监测数据呈现多频率成分的非平稳特性。某些复杂环境效应和荷载条件下,振动数据往往具有非线性特征,例如重载车辆作用下的桥梁拉索振动数据及强震作用下受损结构的振动数据[5-7]。这些信号的瞬时频率为结构健康状态评估与损伤检测提供了丰富信息,因此精确的时频分析估计对结构健康监测(SHM)至关重要[8-10]。此类非线性非平稳 SHM 信号的时频分析关键在于精准的自适应模态分解过程,该过程可将振动信号分解为一组称为本征模态函数(IMFs)的单分量项[11]。每个 IMF 对应一个独特的时变频率成分,进而可通过基于希尔伯特变换的解析信号 [ 12 , 13 ] [ 12 , 13 ] [12,13][12,13] 计算瞬时频率,用于结构的进一步损伤检测或状态评估。 因此,将信号准确分解为固有模态函数(IMFs)是后续有意义时变参数识别的关键步骤。
传统的基于基展开的模式分解方法,如傅里叶变换,采用形式已知的正交基对信号进行分解,存在灵活性差和适应性不足的缺点[14,15],无法为非线性非平稳信号提供足够的精度。为此,Huang 等人[11]提出了经验模态分解(EMD),其具有自适应性和操作简便的优势,已被广泛应用于结构健康监测(SHM)中的模态参数识别 [ 16 , 17 ] [ 16 , 17 ] [16,17][16,17] 和损伤检测[18-20]。Smith 提出的局部均值分解(LMD)[21]摒弃了包络提取中的插值操作,在一定程度上避免了端点效应。Huang 等人提出的集合经验模态分解(EEMD)[22]提高了白噪声鲁棒性,解决了 EMD 的模态混叠问题[23-25]。随后,Gilles 提出的经验小波变换(EWT)[26]与 Dragomiretskiy 和 Zosso 提出的变分模态分解(VMD)[27]均基于频谱分割,相比 EMD 具有坚实的数学基础,并在衍生方法中针对模态数量进行了改进[28,29]。 侯和史[30,31]提出了基于非线性基追踪和非线性匹配追踪的数据驱动自适应稀疏时频分析(ASTFA)方法,这些方法在理论上比 EMD 更具基础性,其优势在于 IMF 具有稀疏性。然而,由包络线的三次样条插值引起的端点效应[11]、间歇性信号和噪声导致的模态混叠[32,33],以及受限于频率分辨率而无法分离的单分量[34],都可能导致自适应模态分解中出现虚假单分量,其中 EMD 得到的 IMF 包含不同频率的分量[35]。
这些虚假分解的单分量存在于 EEMD 和 LMD[35]中。预先选定的分解模式会导致 VMD[27,36]产生过分割或欠分割的 IMF,而不恰当的初始值会影响 ASTFA 的收敛性并导致虚假单分量。然而,对于大型结构而言,振动信号通常包含丰富的模态信息。在周围环境存在显著背景噪声的情况下,很难从振动信号中精确分解这些模态,现有方法容易获得虚假单分量,从而为有效的瞬时频率估计带来挑战。
自适应模态分解和瞬时频率估计的时频分析理论上是受非线性目标函数控制的优化问题。因此,深度学习网络(DNNs)等机器学习方法因其更强大的全局最优性和逼近任意非线性模型的特性而得以应用[37-39]。在结构健康监测(SHM)领域,已有研究将机器学习与结构模态识别、系统辨识、损伤检测等任务相结合[40-42]。Liu 等[43]提出了一种结合独立特性的机器学习方法用于结构模态参数识别;Guo 等[44]在损失函数中嵌入物理信息的长短期记忆神经网络,用于振动台系统辨识;Hesser 等[45]利用二维卷积神经网络进行特征提取,并通过迁移学习训练数据,实现了声发射源的精准识别。前期自适应稀疏时频分析的研究工作已表明机器学习在精确瞬时频率估计方面的良好潜力[46]。 然而,该方法中目标 IMF 的手动筛选及基于傅里叶的字典使用,对于非三角信号的自适应性存在明显不足。
为提升自适应性并克服模态混叠与虚假单分量问题,本文提出一种基于深度学习的信号自适应稀疏模态分解与瞬时频率估计新方法。通过将时频分析转化为深度学习中的优化问题,构建了内置自适应冗余滤波器及稀疏正则化的深度神经网络。借助 DNN 强大的优化能力,实现了 IMFs 的精确分解与高分辨率瞬时频率获取。

2. 面向时频分析的深度学习方法

2.1 理论基础

考虑一个非线性非平稳的结构振动信号 f ( t ) f ( t ) f(t)\mathbf{f}(t) ,包含 N N NN 个离散点。该原始信号可分解为 M M MM 个单分量和一个残差,其中每个单分量是一个 IMF,代表包含独特自然频率的独立模态。IMF 的数量应尽可能少,因此仅需最少数量的 IMF 组合即可最大程度重构信号。这一模态分解过程可表述为:
{ minimize M subject to : f ( t ) = k = 1 M u k ( t ) + r ( t ) , t R , M R , f ( t ) R N minimize M  subject to :  f ( t ) = k = 1 M u k ( t ) + r ( t ) , t R , M R , f ( t ) R N {[minimize M],[" subject to : "f(t)=sum_(k=1)^(M)u_(k)(t)+r(t)","t inR","M inR","f(t)inR^(N)]:}\left\{\begin{array}{l} \operatorname{minimize} M \\ \text { subject to : } \mathbf{f}(t)=\sum_{k=1}^{M} \mathbf{u}_{k}(t)+\mathbf{r}(t), t \in \mathbb{R}, M \in \mathbb{R}, \mathbf{f}(t) \in \mathbb{R}^{N} \end{array}\right.
其中 f ( t ) f ( t ) f(t)\mathbf{f}(t) 为待分析的非线性非平稳时间序列信号, M M MM 是 IMF 的数量, u k ( t ) u k ( t ) u_(k)(t)\mathbf{u}_{k}(t) r ( t ) r ( t ) r(t)\mathbf{r}(t) 分别为第 k k kk 个分解出的 IMF 及残差。
模态确定后,通过频谱方法获取每个 IMF 的解析形式,仅选取其傅里叶变换的正频率成分[8,9]。因此,每个 IMF 解析信号的频谱为:
z ^ k ( ω ) = { u ^ k ( ω ) , for ω > 0 0 , otherwise z ^ k ( ω ) = u ^ k ( ω ) ,  for  ω > 0 0 ,  otherwise  widehat(z)_(k)(omega)={[ widehat(u)_(k)(omega)","," for "omega > 0],[0","," otherwise "]:}\widehat{\mathbf{z}}_{k}(\omega)= \begin{cases}\widehat{\mathbf{u}}_{k}(\omega), & \text { for } \omega>0 \\ 0, & \text { otherwise }\end{cases}
其中 u ^ k ( ω ) = F [ u k ( t ) ] u ^ k ( ω ) = F u k ( t ) widehat(u)_(k)(omega)=F[u_(k)(t)]\widehat{\mathbf{u}}_{k}(\omega)=\mathscr{F}\left[\mathbf{u}_{k}(t)\right] F F F\mathscr{F} 表示傅里叶变换。随后如式(3)和式(4)[12,13]所示,从解析信号中识别各 IMF 的瞬时频率:
z k ( t ) = F 1 [ z ^ k ( ω ) ] = z k , r e ( t ) + j z k , i m ( t ) ω ~ k ( t ) = d d t ( arctan z k , i m ( t ) z k , r e ( t ) ) z k ( t ) = F 1 z ^ k ( ω ) = z k , r e ( t ) + j z k , i m ( t ) ω ~ k ( t ) = d d t arctan z k , i m ( t ) z k , r e ( t ) {:[z_(k)(t)=F^(-1)[ widehat(z)_(k)(omega)]=z_(k,re)(t)+j*z_(k,im)(t)],[ widetilde(omega)_(k)(t)=(d)/(dt)(arctan((z_(k,im)(t))/(z_(k,re)(t))))]:}\begin{aligned} & \mathbf{z}_{k}(t)=\mathscr{F}^{-1}\left[\widehat{\mathbf{z}}_{k}(\omega)\right]=\mathbf{z}_{k, r e}(t)+j \cdot \mathbf{z}_{k, i m}(t) \\ & \widetilde{\boldsymbol{\omega}}_{k}(t)=\frac{d}{d t}\left(\arctan \frac{\mathbf{z}_{k, i m}(t)}{\mathbf{z}_{k, r e}(t)}\right) \end{aligned}
其中 z k , r e ( t ) z k , r e ( t ) z_(k,re)(t)\mathbf{z}_{k, r e}(t) z k , i m ( t ) z k , i m ( t ) z_(k,im)(t)\mathbf{z}_{k, i m}(t) 分别是第 k k kk 个解析信号的实部与虚部, F 1 F 1 F^(-1)\mathscr{F}^{-1} 表示傅里叶逆变换, ω ~ k ( t ) ω ~ k ( t ) widetilde(omega)_(k)(t)\widetilde{\boldsymbol{\omega}}_{k}(t) 则是第 k k kk 个瞬时频率。

2.2. 机器学习任务形式化

上述问题的核心在于如何有效分解模态而不产生虚假单分量信号,因为分解结果直接影响瞬时频率的准确性。每个模态中应仅呈现单一特定频率成分。式(1)中的模态分解问题可转化为如式(5)所示的非凸最小二乘优化问题:
minimize M , subject to : f ( t ) k = 1 M u k ( t ) 2 2 ε minimize M ,  subject to :  f ( t ) k = 1 M u k ( t ) 2 2 ε minimize M," subject to : "||f(t)-sum_(k=1)^(M)u_(k)(t)||_(2)^(2) <= epsi\operatorname{minimize} M, \text { subject to : }\left\|\mathbf{f}(t)-\sum_{k=1}^{M} \mathbf{u}_{k}(t)\right\|_{2}^{2} \leqslant \varepsilon
其中 ε ε epsi\varepsilon 是一个小参数,表示信号与本征模态函数(IMFs)之和的差异。由于这一非凸最小二乘问题不存在唯一解,传统优化算法专注于从局部展开求解以降低计算成本和算法复杂度,容易陷入局部最优。在此情况下,深度学习凭借非线性模型强大的优化能力和无限逼近特性,能更有效地找到理想解。因此,式(1)-(4)中的时频分析可转化为机器学习任务,这可以视为对已知待分析信号施加非线性变换 F F FF ,并如式(6)所示输出瞬时频率:
ω ~ ( t ) = F ( f ( t ) ) , f ( t ) R N , t [ 0 , T ] ω ~ ( t ) = F ( f ( t ) ) , f ( t ) R N , t [ 0 , T ] widetilde(omega)(t)=F(f(t)),quadf(t)inR^(N),quad t in[0,T]\widetilde{\boldsymbol{\omega}}(t)=F(\mathbf{f}(t)), \quad \mathbf{f}(t) \in \mathbb{R}^{N}, \quad t \in[0, T]
上述内容构成了基于深度学习的数据驱动型信号自适应稀疏模态分解与瞬时频率估计技术的基石,其程序流程图如图 1 所示。首先,将结构健康监测系统(SHM)采集的振动信号 f ( t ) f ( t ) f(t)\mathbf{f}(t) 作为输入离散时域信号,通过傅里叶变换得到频谱 F ^ F ^ widehat(F)\widehat{\mathbf{F}} 以便进行频域操作。接着,应用一组冗余的自适应滤波器 A A A\mathbf{A} (其可变中心频率、带宽及稀疏系数 W W W\mathbf{W} 作为深度神经网络的权重)进行最稀疏单分量频谱提取,从而分解出本征模态函数(IMFs) F ^ I M F F ^ I M F widehat(F)_(IMF)\widehat{\mathbf{F}}_{I M F} 。随后提取 IMFs 的正频率成分作为解析信号频谱 Z ^ Z ^ widehat(Z)\widehat{\mathbf{Z}} ,最终通过解析信号的相位估计瞬时频率 ω ~ ( t ) ω ~ ( t ) widetilde(omega)(t)\widetilde{\boldsymbol{\omega}}(t) 。与传统方法不同,IMFs 提取与瞬时频率估计的流程均由设计的深度神经网络实现。
图1. 程序流程图

2.3 深度神经网络构建

式(5)中的优化问题类似于深度神经网络的目标函数。所设计的八层全连接 DNN 优化器细节如图 2 所示,其实现了时频分析中的两个关键步骤:从原始信号到本征模态函数(IMFs)的自适应稀疏模态分解,以及 IMFs 的瞬时频率估计。
自适应稀疏模态分解在设计的 DNN 前四层完成。时间序列信号首先输入至输入层,随后通过傅里叶基矩阵转换至频域。其次,在第二层将信号分离为实部与虚部后,引入具有可变中心频率和带宽的自适应滤波器组,并筛选单分量以实现自适应滤波。为避免虚假单分量,在第三层输出端施加稀疏正则化(亦称 1 1 ℓ_(1)\ell_{1} 正则化)以获取最稀疏的单分量。该过程通过剔除冗余分量,并利用最稀疏单分量——即本征模态函数(IMFs)重构信号。
固有模态函数(IMFs)的瞬时频率估计通过第 5、6、7 和 8 层的希尔伯特变换实现。首先,对解析信号的频谱进行正频率滤波。接着,乘以逆傅里叶矩阵并合并结果,从解析信号的相位中计算出瞬时频率。
在该设计的深度神经网络中,三类权重至关重要。前两类是不同自适应滤波器的峰值中心 C C C\mathbf{C} (余下部分简称中心)和带宽 D D D\mathbf{D} ,第三类则是稀疏正则化层生成的滤波器稀疏系数 W W W\mathbf{W}
具体细节如下。对于设计的深度神经网络的第二层,执行如式(7)所示的傅里叶变换,其中 F ^ F ^ widehat(F)\widehat{\mathbf{F}} f ( t ) f ( t ) f(t)\mathbf{f}(t) 的傅里叶变换矩阵形式, B B B\mathbf{B} 为傅里叶矩阵:
F ^ = f ( t ) B = [ f ^ 1 ( ω ) , f ^ 2 ( ω ) , , f ^ N ( ω ) ] R N × N B = ( b r , k ) r = 1 , , N , k = 1 , , N F ^ = f ( t ) B = f ^ 1 ( ω ) , f ^ 2 ( ω ) , , f ^ N ( ω ) R N × N B = b r , k r = 1 , , N , k = 1 , , N {:[ widehat(F)=f(t)*B=[ widehat(f)_(1)(omega), widehat(f)_(2)(omega),dots, widehat(f)_(N)(omega)]inR^(N xx N)],[B=(b_(r,k))_(r=1,dots,N,k=1,dots,N)]:}\begin{aligned} & \widehat{\mathbf{F}}=\mathbf{f}(t) \cdot \mathbf{B}=\left[\widehat{\mathbf{f}}_{1}(\omega), \widehat{\mathbf{f}}_{2}(\omega), \ldots, \widehat{\mathbf{f}}_{N}(\omega)\right] \in \mathbb{R}^{N \times N} \\ & \mathbf{B}=\left(b_{r, k}\right)_{r=1, \ldots, N, k=1, \ldots, N} \end{aligned}
在设计的 DNN 第三层中,引入了冗余自适应滤波器以提取频谱单分量。此步骤中,冗余滤波器提取的单分量数量多于实际固有模态函数(IMFs)。滤波是一种提取有意义信号成分并去除不理想部分的方法[47]。本方法的一个主要差异在于采用了一组具有不同中心频率和带宽的自适应滤波器,其设计基于波数域中的 Meyer 尺度函数,如式(8)所示:
Φ = { 1 , | l | 1 1 2 ( 1 cos π l ) , 1 < | l | < 2 0 , otherwise Φ = 1 , | l | 1 1 2 ( 1 cos π l ) , 1 < | l | < 2 0 ,  otherwise  Phi={[1",",|l| <= 1],[(1)/(2)(1-cos pi l)",",1 < |l| < 2],[0","," otherwise "]:}\Phi=\left\{\begin{array}{lr} 1, & |l| \leqslant 1 \\ \frac{1}{2}(1-\cos \pi l), & 1<|l|<2 \\ 0, & \text { otherwise } \end{array}\right.
其中 l l ll 表示波数。对于所提出的包含 K ( K > M ) K ( K > M ) K(K > M)K(K>M) 个元素的自适应滤波器组 A A A\mathbf{A} ,其矩阵形式表达式为:
图 2. 所提方法的 DNN 结构示意图。
A = m 1 m 2 2 ( 1 cos | l | C D D π ) + ( 1 m 1 ) m 2 A = [ a ^ 1 ( ω ) , a ^ 2 ( ω ) , , a ^ K ( ω ) ] R N × K , C = [ c 1 , c 2 , , c K ] R K , D = [ d 1 , d 2 , , d K ] R K A = m 1 m 2 2 1 cos | l | C D D π + 1 m 1 m 2 A = a ^ 1 ( ω ) , a ^ 2 ( ω ) , , a ^ K ( ω ) R N × K , C = c 1 , c 2 , , c K R K , D = d 1 , d 2 , , d K R K {:[A=(m_(1)m_(2))/(2)(1-cos((|l|-C∣-D)/(D)pi))+(1-m_(1))m_(2)],[A=[ widehat(a)_(1)(omega), widehat(a)_(2)(omega),dots, widehat(a)_(K)(omega)]inR^(N xx K)","C=[c_(1),c_(2),dots,c_(K)]inR^(K)","D=[d_(1),d_(2),dots,d_(K)]inR^(K)]:}\begin{aligned} & \mathbf{A}=\frac{\mathbf{m}_{1} \mathbf{m}_{2}}{2}\left(1-\cos \frac{|l|-\mathbf{C} \mid-\mathbf{D}}{\mathbf{D}} \pi\right)+\left(1-\mathbf{m}_{1}\right) \mathbf{m}_{2} \\ & \mathbf{A}=\left[\widehat{\mathbf{a}}_{1}(\omega), \widehat{\mathbf{a}}_{2}(\omega), \ldots, \widehat{\mathbf{a}}_{K}(\omega)\right] \in \mathbb{R}^{N \times K}, \mathbf{C}=\left[c_{1}, c_{2}, \ldots, c_{K}\right] \in \mathbb{R}^{K}, \mathbf{D}=\left[d_{1}, d_{2}, \ldots, d_{K}\right] \in \mathbb{R}^{K} \end{aligned}
式中 a ^ i ( ω ) , i = 1 , , K a ^ i ( ω ) , i = 1 , , K widehat(a)_(i)(omega),i=1,dots,K\widehat{\mathbf{a}}_{i}(\omega), i=1, \ldots, K 表示傅里叶域中的第 i i ii 个自适应滤波器; l l ll 为水平方向坐标; m 1 = sgn ( l l C ) + 1 2 m 1 = sgn ( l l C ) + 1 2 m_(1)=(sgn(ll-C∣)+1)/(2)\mathbf{m}_{1}=\frac{\operatorname{sgn}(l l-\mathbf{C} \mid)+1}{2} m 2 = sgn ( D | | l | C | ) + 1 2 m 2 = sgn ( D | | l | C | ) + 1 2 m_(2)=(sgn(D-||l|-C|)+1)/(2)\mathbf{m}_{2}=\frac{\operatorname{sgn}(\mathbf{D}-||l|-\mathbf{C}|)+1}{2} 是将余弦形函数转换为双峰形函数的两个参数; C C C\mathbf{C} D D D\mathbf{D} 分别是 K K KK 个滤波器的中心频率矩阵和带宽矩阵,这些参数控制着自适应滤波器的形状,作为设计 DNN 第三层的权重参数。
冗余单分量矩阵 F ^ red F ^ red  widehat(F)_("red ")\widehat{\mathbf{F}}_{\text {red }} 在设计的 DNN 第三层后按式(10)提取,其中 f ^ i ( ω ) f ^ i ( ω ) widehat(f)_(i)(omega)\widehat{\mathbf{f}}_{i}(\omega) 表示由第 i i ii 个自适应滤波器检测到的信号第 i i ii 个频谱分量:
F ^ red = A F ^ = [ f ^ 1 ( ω ) , f ^ 2 ( ω ) , , f ^ K ( ω ) ] C N × K F ^ red  = A F ^ = f ^ 1 ( ω ) , f ^ 2 ( ω ) , , f ^ K ( ω ) C N × K widehat(F)_("red ")=A* widehat(F)=[ widehat(f)_(1)(omega), widehat(f)_(2)(omega),dots, widehat(f)_(K)(omega)]inC^(N xx K)\widehat{\mathbf{F}}_{\text {red }}=\mathbf{A} \cdot \widehat{\mathbf{F}}=\left[\widehat{\mathbf{f}}_{1}(\omega), \widehat{\mathbf{f}}_{2}(\omega), \ldots, \widehat{\mathbf{f}}_{K}(\omega)\right] \in \mathbb{C}^{N \times K}
在第四层应用稀疏矩阵以确保信号重建所需 IMF 数量最小化,此步骤将剔除冗余不必要的单分量。重建信号 f ^ rec ( ω ) f ^ rec  ( ω ) widehat(f)_("rec ")(omega)\widehat{\mathbf{f}}_{\text {rec }}(\omega) 的频谱如式(11)和(12)所示:
F ^ I M F = F ^ red W f ^ r e c ( ω ) = k = 1 M F ^ I M F ( r , k ) , f ^ rec ( ω ) R N F ^ I M F = F ^ red  W f ^ r e c ( ω ) = k = 1 M F ^ I M F ( r , k ) , f ^ rec  ( ω ) R N {:[ widehat(F)_(IMF)= widehat(F)_("red ")*W],[ widehat(f)_(rec)(omega)=sum_(k=1)^(M) widehat(F)_(IMF)(r","k)"," widehat(f)_("rec ")(omega)inR^(N)]:}\begin{aligned} & \widehat{\mathbf{F}}_{I M F}=\widehat{\mathbf{F}}_{\text {red }} \cdot \mathbf{W} \\ & \widehat{\mathbf{f}}_{r e c}(\omega)=\sum_{k=1}^{M} \widehat{\mathbf{F}}_{I M F}(r, k), \widehat{\mathbf{f}}_{\text {rec }}(\omega) \in \mathbb{R}^{N} \end{aligned}
其中 W W W\mathbf{W} 为稀疏权重矩阵(同时也是待优化的 DNN 权重), F ^ IMF F ^ IMF  widehat(F)_("IMF ")\widehat{\mathbf{F}}_{\text {IMF }} 表示 IMF 频谱矩阵, r r rr k k kk 分别对应矩阵 F ^ IMF F ^ IMF  widehat(F)_("IMF ")\widehat{\mathbf{F}}_{\text {IMF }} 的行列索引。
本方法的核心在于从原始信号中获取自适应稀疏 IMF 分量,进而精确识别瞬时频率。该优化过程通过最小化重建时序信号(即 IMF 分量之和)与输入原始信号之间的差异来实现,其机器学习任务的最优化问题可表述为式(13),精确表达式见式(14):
min ( f ( t ) f rec ( t ) 2 2 + W 1 ) min A C N × K , W R N f A F ^ W B 1 2 2 + W 1 min f ( t ) f rec  ( t ) 2 2 + W 1 min A C N × K , W R N f A F ^ W B 1 2 2 + W 1 {:[ min(||f(t)-f_("rec ")(t)||_(2)^(2)+||W||_(1))],[min_(AinC^(N xx K),WinR^(N))||f-A*( widehat(F))*W*B^(-1)||_(2)^(2)+||W||_(1)]:}\begin{aligned} & \min \left(\left\|\mathbf{f}(t)-\mathbf{f}_{\text {rec }}(t)\right\|_{2}^{2}+\|\mathbf{W}\|_{1}\right) \\ & \min _{\mathbf{A} \in \mathbb{C}^{N \times K}, \mathbf{W} \in \mathbb{R}^{N}}\left\|\mathbf{f}-\mathbf{A} \cdot \widehat{\mathbf{F}} \cdot \mathbf{W} \cdot \mathbf{B}^{-1}\right\|_{2}^{2}+\|\mathbf{W}\|_{1} \end{aligned}
其中 f ( t ) f r e c ( t ) 2 2 = r = 1 N ( f ( t r ) f r e c ( t r ) ) 2 f ( t ) f r e c ( t ) 2 2 = r = 1 N f t r f r e c t r 2 ||f(t)-f_(rec)(t)||_(2)^(2)=sum_(r=1)^(N)(f(t_(r))-f_(rec)(t_(r)))^(2)\left\|\mathbf{f}(t)-\mathbf{f}_{r e c}(t)\right\|_{2}^{2}=\sum_{r=1}^{N}\left(f\left(t_{r}\right)-f_{r e c}\left(t_{r}\right)\right)^{2} W 1 = i = 1 K | w i | W 1 = i = 1 K w i ||W||_(1)=sum_(i=1)^(K)|w_(i)|\|\mathbf{W}\|_{1}=\sum_{i=1}^{K}\left|w_{i}\right|
为实现设计 DNN 对信号的精确模态分解,式(15)提出了一种新颖的双重损失函数:
L = 1 N r = 1 N ( f ( t r ) F 1 [ k = 1 M A F ^ W ( r , k ) ] ) 2 + i = 1 K | w i | L = 1 N r = 1 N f t r F 1 k = 1 M A F ^ W ( r , k ) 2 + i = 1 K w i L=(1)/(N)sum_(r=1)^(N)(f(t_(r))-F^(-1)[sum_(k=1)^(M)A*( widehat(F))*W(r,k)])^(2)+sum_(i=1)^(K)|w_(i)|L=\frac{1}{N} \sum_{r=1}^{N}\left(f\left(t_{r}\right)-\mathscr{F}^{-1}\left[\sum_{k=1}^{M} \mathbf{A} \cdot \widehat{\mathbf{F}} \cdot \mathbf{W}(r, k)\right]\right)^{2}+\sum_{i=1}^{K}\left|w_{i}\right|
其中第一项用于最小化输入信号与重构信号之间的重建误差,第二项则用于在冗余分量上施加最稀疏[48]的权重,从而在倒数第二层获得最稀疏的 IMF; N N NN 表示离散时间域信号的总长度; k k kk 为列数, Im [ ] Im [ ] Im[*]\operatorname{Im}[\cdot] 代表复矩阵的虚部。
所提出的方法为数据驱动型,依赖于数据训练,但其模型架构具有通用性,可适用于不同场景。对于该显式神经网络模型,训练过程同时也是获取解的优化过程。在优化过程中,我们采用随机梯度下降法(SGD)更新权重。由于损失函数是非凸函数,不仅存在一个全局最小值,还可能存在多个局部极小值,因此每一步参数更新仅使用单个样本。这种方式下,梯度与传统梯度下降法相比可能存在轻微偏差。
在完成本征模态函数(IMFs)的自适应稀疏分解后,瞬时频率估计流程被嵌入后续无需优化权重的深度神经网络层中。通过如式(16)所示在频域构建解析信号并提取正频率分量,在获得最优自适应滤波器组 A A A\mathbf{A} 及其关联稀疏权重 W W W\mathbf{W} 后,即可实现瞬时频率估计:
Z ^ = A F ^ W H = [ z ^ 1 ( ω ) , z ^ 2 ( ω ) , , z ^ M ( ω ) ] C N × M Z ^ = A F ^ W H = z ^ 1 ( ω ) , z ^ 2 ( ω ) , , z ^ M ( ω ) C N × M widehat(Z)=A* widehat(F)*W*H=[ widehat(z)_(1)(omega), widehat(z)_(2)(omega),dots, widehat(z)_(M)(omega)]inC^(N xx M)\widehat{\mathbf{Z}}=\mathbf{A} \cdot \widehat{\mathbf{F}} \cdot \mathbf{W} \cdot \mathbf{H}=\left[\widehat{\mathbf{z}}_{1}(\omega), \widehat{\mathbf{z}}_{2}(\omega), \ldots, \widehat{\mathbf{z}}_{M}(\omega)\right] \in \mathbb{C}^{N \times M}
其中 Z ^ Z ^ widehat(Z)\widehat{\mathbf{Z}} 表示本征模态函数解析信号的频谱矩阵, H H H\mathbf{H} 为滤波器组,其定义为:
H = { f ^ r e c ( ω ) , for ω = 0 2 f ^ r e c ( ω ) , for 1 ω N 2 1 f ^ r e c ( N 2 ) , for ω = N 2 0 , for N 2 + 1 ω N 1 H = f ^ r e c ( ω ) ,  for  ω = 0 2 f ^ r e c ( ω ) ,  for  1 ω N 2 1 f ^ r e c N 2 ,  for  ω = N 2 0 ,  for  N 2 + 1 ω N 1 H={[ widehat(f)_(rec)(omega)","," for "omega=0],[2 widehat(f)_(rec)(omega)","," for "1 <= omega <= (N)/(2)-1],[ widehat(f)_(rec)((N)/(2))","," for "omega=(N)/(2)],[0","," for "(N)/(2)+1 <= omega <= N-1]:}\mathbf{H}= \begin{cases}\widehat{\mathbf{f}}_{r e c}(\omega), & \text { for } \omega=0 \\ 2 \widehat{\mathbf{f}}_{r e c}(\omega), & \text { for } 1 \leqslant \omega \leqslant \frac{N}{2}-1 \\ \widehat{\mathbf{f}}_{r e c}\left(\frac{N}{2}\right), & \text { for } \omega=\frac{N}{2} \\ 0, & \text { for } \frac{N}{2}+1 \leqslant \omega \leqslant N-1\end{cases}
频谱 Z ^ Z ^ widehat(Z)\widehat{\mathbf{Z}} 通过逆傅里叶矩阵 B 1 B 1 B^(-1)\mathbf{B}^{-1} (式(18))转换回解析信号 Z Z Z\mathbf{Z} ,因此瞬时频率可直接从解析信号中计算得出[12,13],如式(19)所示:
Z = Z ^ B 1 ω ~ ( t ) = 1 2 π d d t ( arctan Z i m Z r e ) Z = Z ^ B 1 ω ~ ( t ) = 1 2 π d d t arctan Z i m Z r e {:[Z= widehat(Z)*B^(-1)],[ widetilde(omega)(t)=(1)/(2pi)(d)/(dt)(arctan((Z_(im))/(Z_(re))))]:}\begin{aligned} & \mathbf{Z}=\widehat{\mathbf{Z}} \cdot \mathbf{B}^{-1} \\ & \widetilde{\boldsymbol{\omega}}(t)=\frac{1}{2 \pi} \frac{d}{d t}\left(\arctan \frac{\mathbf{Z}_{i m}}{\mathbf{Z}_{r e}}\right) \end{aligned}
其中 Z i m Z i m Z_(im)\mathbf{Z}_{i m} Z r e Z r e Z_(re)\mathbf{Z}_{r e} 分别是 Z Z Z\mathbf{Z} 的虚部和实部。

  3. 示例

3.1. 示例1 具有幅频调制的合成信号

以一个包含五个调幅或调频本征模态函数(IMF)的非平稳信号为例。每个 IMF 的解析表达式均已给出,原始输入信号为 f f ff (即五个 IMF 之和):
图3. 示例1中用于合成信号的深度神经网络结构。
f 1 = sin ( 500 π t + cos ( 4 π t ) ) , f 2 = ( 0.5 + cos ( 2 π t ) ) cos ( 380 π t ) , f 3 = 1 1.5 + sin ( 2 π t ) cos ( 240 π t + sin ( 16 π t ) ) , f 4 = sin ( 150 π t + 2 sin ( 6 π t ) ) , f 5 = cos ( 30 π t + 10 sin ( 2 π t ) ) , f 1 = sin ( 500 π t + cos ( 4 π t ) ) , f 2 = ( 0.5 + cos ( 2 π t ) ) cos ( 380 π t ) , f 3 = 1 1.5 + sin ( 2 π t ) cos ( 240 π t + sin ( 16 π t ) ) , f 4 = sin ( 150 π t + 2 sin ( 6 π t ) ) , f 5 = cos ( 30 π t + 10 sin ( 2 π t ) ) , {:[f_(1)=sin(500 pi t+cos(4pi t))","],[f_(2)=(0.5+cos(2pi t))*cos(380 pi t)","],[f_(3)=(1)/(1.5+sin(2pi t))*cos(240 pi t+sin(16 pi t))","],[f_(4)=sin(150 pi t+2sin(6pi t))","],[f_(5)=cos(30 pi t+10 sin(2pi t))","]:}\begin{gathered} f_{1}=\sin (500 \pi t+\cos (4 \pi t)), \\ f_{2}=(0.5+\cos (2 \pi t)) \cdot \cos (380 \pi t), \\ f_{3}=\frac{1}{1.5+\sin (2 \pi t)} \cdot \cos (240 \pi t+\sin (16 \pi t)), \\ f_{4}=\sin (150 \pi t+2 \sin (6 \pi t)), \\ f_{5}=\cos (30 \pi t+10 \sin (2 \pi t)), \end{gathered}
f = f 1 + f 2 + f 3 + f 4 + f 5 f = f 1 + f 2 + f 3 + f 4 + f 5 f=f_(1)+f_(2)+f_(3)+f_(4)+f_(5)f=f_{1}+f_{2}+f_{3}+f_{4}+f_{5}
根据图2所示的网络结构,本例采用的深度神经网络如图3所示,其预训练网络参数列于表1。神经网络超参数基于经验选取。信号总长度为 N N NN = 1024 = 1024 =1024=1024 ,初始冗余滤波器及系数数量设为20。为加速收敛,初始中心频率值围绕解析中心频率设定而非随机设置。此外,自适应滤波器带宽可先设为合理数值,随后在优化过程中调整以匹配各分量的频谱宽度。
网络重构的时频信号如图 4 所示。根据设计的 DNN 框架,输出的重构信号是各 IMF 分量的总和。IMF 分量的提取结果如图 5 和图 6 所示,表明所提方法的分解结果更接近 IMF 的解析表达式。图 5(b)和图 6(b)展示了 EMD 的相应波形,其中提取的 IMF 分量含有额外的固有模态。对于 EMD 而言,虚假的 IMF 分量几乎不可避免,因为它们是通过三次样条插值从包络中提取的[11]。此外,EMD 的另一个局限是当某些 IMF 分量的频率过于接近时(即两个 IMF 分量的频率比高于 0.75[34]),它们难以被分离。如图 6(b)所示, f 1 f 1 f_(1)f_{1} f 2 ( f 1 / f 2 = 0.76 > 0.75 ) f 2 f 1 / f 2 = 0.76 > 0.75 f_(2)(f_(1)//f_(2)=0.76 > 0.75)f_{2}\left(f_{1} / f_{2}=0.76>0.75\right) 对于 EMD 是不可分离的,而所提方法(图 6(a))则能清晰分离,没有模态混叠或虚假单分量。图 5(c)和图 6(c)展示了 VMD 的结果,同样获得了准确的 IMF 分量。
在详细说明 DNN 中自适应稀疏滤波器如何提取 IMF 时,图 7 展示了频域中每个 IMF 对应的滤波器波形。当训练轮次为 0 时,中心频率和带宽的初始值设置如表 2 和表 3 所示。每个滤波器的初始带宽均为 20,其中五个滤波器的初始中心频率经过粗略调整接近解析值,其余则随机选择。图 7 表明自适应滤波器能同步进行滤波操作,每个滤波器独立提取一个 IMF,且滤波器的最终中心频率接近解析解,与初始值相差较大。此外,尽管初始赋予相同的滤波器带宽,但根据各滤波器提取的 IMF 带宽特性,其最终带宽值会动态调整。这些特性验证了所提方法具有良好的自适应能力。
合成信号 f f ff 的神经网络模型损失曲线如图 8 所示。为清晰展示形态变化, IMF 3 IMF 3 IMF_(3)\mathrm{IMF}_{3} IMF 1 IMF 1 IMF_(1)\mathrm{IMF}_{1} 的自适应滤波提取过程展示于图 9(a)。优化过程中,中心频率和带宽两个参数被调整以匹配固有模态函数(IMFs)。所有 IMFs 参数列于表 2 和表 3 中。由表可见, IMF 3 IMF 3 IMF_(3)\mathrm{IMF}_{3} IMF 1 IMF 1 IMF_(1)\mathrm{IMF}_{1} 的中心频率经优化后更接近解析值,最终分别达到 119.94 和 253.63,精度分别为 99.95 % 99.95 % 99.95%99.95 \% 98.55 % 98.55 % 98.55%98.55 \% 。最终带宽值分别为 34.36 和 29.15,完全覆盖了这两个 IMFs 在频域的波段范围。
冗余自适应滤波器的稀疏系数如图 9(b)所示。在 30 个训练周期时,几乎所有权重已完成优化且方差极小,表明优化接近完成。当周期达到 100 时,权重呈现最稀疏分布。仅五个滤波器的最终权重约为 1,说明 IMFs 的最稀疏数量为 5 个,同时...
  表1
合成信号 f f ff 的深度神经网络相关参数。
  参数   数量
  初始滤波器数量 20
初始滤波器系数数量 20
  学习率 1
  动量 0.9
  最大训练轮数 100
  批次大小 64
Parameter Number Initial filter numbers 20 Initial filter coefficient numbers 20 Learning rate 1 Momentum 0.9 Max epoch 100 Batch size 64| Parameter | Number | | :--- | :--- | | Initial filter numbers | 20 | | Initial filter coefficient numbers | 20 | | Learning rate | 1 | | Momentum | 0.9 | | Max epoch | 100 | | Batch size | 64 |
图 4. 输入合成信号 f f ff 及通过所提方法重构的信号:(a)波形与(b)傅里叶频谱。
  表2
合成信号 f f ff 的自适应滤波器中心频率(Hz)。
  迭代次数 固有模态函数
Epochs IMF| Epochs | | :--- | | IMF |
25 50 75
1 270.00 259.74 256.60 254.82
2 175.00 182.86 185.19 186.61
3 135.00 121.10 120.15 119.98
4 90.00 80.60 78.16 78.01
5 5.00 14.99 14.99 14.99
"Epochs IMF" 25 50 75 1 270.00 259.74 256.60 254.82 2 175.00 182.86 185.19 186.61 3 135.00 121.10 120.15 119.98 4 90.00 80.60 78.16 78.01 5 5.00 14.99 14.99 14.99| Epochs <br> IMF | | 25 | 50 | 75 | | :--- | ---: | ---: | ---: | ---: | | 1 | 270.00 | 259.74 | 256.60 | 254.82 | | 2 | 175.00 | 182.86 | 185.19 | 186.61 | | 3 | 135.00 | 121.10 | 120.15 | 119.98 | | 4 | 90.00 | 80.60 | 78.16 | 78.01 | | 5 | 5.00 | 14.99 | 14.99 | 14.99 |
  表3
合成信号 f f ff 的自适应滤波器带宽(Hz)。
  训练周期 0 25 50 75
IMF
1 20.00 27.46 28.52 28.94
2 20.00 24.25 24.88 29.14
3 20.00 29.35 31.43 33.03
4 20.00 27.31 28.58 29.23
5 20.00 29.59 32.32 34.27
Epochs 0 25 50 75 IMF 1 20.00 27.46 28.52 28.94 2 20.00 24.25 24.88 29.14 3 20.00 29.35 31.43 33.03 4 20.00 27.31 28.58 29.23 5 20.00 29.59 32.32 34.27| Epochs | 0 | 25 | 50 | 75 | | :--- | :--- | :--- | :--- | :--- | | IMF | | | | | | 1 | 20.00 | 27.46 | 28.52 | 28.94 | | 2 | 20.00 | 24.25 | 24.88 | 29.14 | | 3 | 20.00 | 29.35 | 31.43 | 33.03 | | 4 | 20.00 | 27.31 | 28.58 | 29.23 | | 5 | 20.00 | 29.59 | 32.32 | 34.27 |
图8. 合成信号 f f ff 的神经网络模型损失曲线。
其他值接近于0。
图 10 对比了不同方法下分解得到的 IMF 分量之间的正交性。每个方块代表通过如下公式[11]计算的正交性结果:
I O r , k = | IMF r , IMF k | IMF r 2 IMF k 2 I O r , k = IMF r , IMF k IMF r 2 IMF k 2 IO_(r,k)=(|(:IMF_(r),IMF_(k):)|)/(||IMF_(r)||_(2)*||IMF_(k)||_(2))I O_{r, k}=\frac{\left|\left\langle\mathrm{IMF}_{r}, \mathrm{IMF}_{k}\right\rangle\right|}{\left\|\mathrm{IMF}_{r}\right\|_{2} \cdot\left\|\mathrm{IMF}_{k}\right\|_{2}}
其中下标 r r rr k k kk 分别对应正交矩阵中的行与列。 I O r , k I O r , k IO_(r,k)I O_{r, k} 值越小表示正交性越好。图 10(b)、图 10(c)及图 10(d)表明,尽管 EMD 和 VMD 方法均获得了可接受的正交性,但所提方法的正交性显著优于这两种方法。
在获取本征模态函数(IMFs)后,相关瞬时频率通过深度神经网络(DNN)识别并输出(图 11(a))。作为对比方法,图 11(b)展示了结合经验模态分解(EMD)与希尔伯特变换的希尔伯特-黄变换(HHT)处理结果。图 11(c)则显示对变分模态分解(VMD)获得的 IMFs 进行希尔伯特变换实现的瞬时频率识别。通过计算各瞬时频率与其解析解之间的相对识别误差( e e ee ),结果表明所提方法与 VMD 均具有较高精度。由于不准确的 IMFs 导致 HHT(图 11(b))识别结果不可用,这印证了精确模态分解结果的重要性。
图 11. 合成信号 f f ff 的瞬时频率识别结果及误差( e e ee ):(a)本文方法,其中 a k ( k = 1 , , 5 a k ( k = 1 , , 5 a_(k)(k=1,dots,5\mathbf{a}_{k}(k=1, \ldots, 5 ω ~ k ( k = 1 , , 5 ) ω ~ k ( k = 1 , , 5 ) widetilde(omega)_(k)(k=1,dots,5)\widetilde{\boldsymbol{\omega}}_{k}(k=1, \ldots, 5) 分别代表解析频率和识别频率,(b)HHT 方法,(c)基于 HT 的 VMD 方法。

3.2 示例2 不同信噪比下频率重叠的合成信号分析

为验证所提方法的抗噪性能,分析了不同信噪比(SNR)下的非线性非平稳信号。相较于前一示例,该信号的单分量具有频谱重叠特性且频带更宽。原始输入信号 s s ss 及其单分量的解析表达式如式(22)所示:
s 1 = { 4 t sin ( ( 50 + 50 t ) 2 π t ) , 0 t 0.25 4 3 ( t 1 ) sin ( ( 50 + 50 t ) 2 π t ) , 0.25 t 1 s 2 = { 2 t sin ( ( 120 + 100 t ) 2 π t ) , 0 t 0.5 2 ( t 1 ) sin ( ( 120 + 100 t ) 2 π t ) , 0.5 t 1 s 3 = ( 2 + cos ( 8 π t ) ) cos ( 300 π ( t + 1 ) 2 ) , s 1 = 4 t sin ( ( 50 + 50 t ) 2 π t ) , 0 t 0.25 4 3 ( t 1 ) sin ( ( 50 + 50 t ) 2 π t ) , 0.25 t 1 s 2 = 2 t sin ( ( 120 + 100 t ) 2 π t ) , 0 t 0.5 2 ( t 1 ) sin ( ( 120 + 100 t ) 2 π t ) , 0.5 t 1 s 3 = ( 2 + cos ( 8 π t ) ) cos 300 π ( t + 1 ) 2 , {:[s_(1)={[4t*sin((50+50 t)*2pi t)",",0 <= t <= 0.25],[-(4)/(3)(t-1)*sin((50+50 t)*2pi t)",",0.25 <= t <= 1]:}],[s_(2)={[2t*sin((120+100 t)*2pi t)",",0 <= t <= 0.5],[-2(t-1)*sin((120+100 t)*2pi t)",",0.5 <= t <= 1]:}],[s_(3)=(2+cos(8pi t))*cos(300 pi(t+1)^(2))","]:}\begin{aligned} & s_{1}= \begin{cases}4 t \cdot \sin ((50+50 t) \cdot 2 \pi t), & 0 \leqslant t \leqslant 0.25 \\ -\frac{4}{3}(t-1) \cdot \sin ((50+50 t) \cdot 2 \pi t), & 0.25 \leqslant t \leqslant 1\end{cases} \\ & s_{2}= \begin{cases}2 t \cdot \sin ((120+100 t) \cdot 2 \pi t), & 0 \leqslant t \leqslant 0.5 \\ -2(t-1) \cdot \sin ((120+100 t) \cdot 2 \pi t), & 0.5 \leqslant t \leqslant 1\end{cases} \\ & s_{3}=(2+\cos (8 \pi t)) \cdot \cos \left(300 \pi(t+1)^{2}\right), \end{aligned}
图 12. 合成信号 s s ss 的波形:(a)信噪比为 20 dB , 10 dB , 5 dB 20 dB , 10 dB , 5 dB 20dB,10dB,5dB20 \mathrm{~dB}, 10 \mathrm{~dB}, 5 \mathrm{~dB} 时的时域波形,(b)对应的频谱图。
s = s 1 + s 2 + s 3 s = s 1 + s 2 + s 3 s=s_(1)+s_(2)+s_(3)s=s_{1}+s_{2}+s_{3}
三个单分量的瞬时频率均为线性单调递增函数,频率区间分别为[50 150]、[120 320]和[300 600]。原始时间序列信号及其在添加白噪声( SNR = 20 dB , 10 dB SNR = 20 dB , 10 dB SNR=20dB,10dB\mathrm{SNR}=20 \mathrm{~dB}, 10 \mathrm{~dB} 和 5 dB)后的频谱波形如图 12 所示。采样频率为 2000 Hz,持续时间为 1 秒。
预训练的 DNN 参数列于表 4。为验证噪声干扰下的识别精度,考虑了 SNR 20 dB , 10 dB SNR 20 dB , 10 dB SNR20dB,10dB\mathrm{SNR} 20 \mathrm{~dB}, 10 \mathrm{~dB} 和 5 dB 两种情况,并将所提方法、EEMD 和 VMD 的信号分解结果对比展示于图 13 至图 18。
如图 13 和图 14 所示,与 VMD(图 17 和图 18)相比,所提方法的单分量在时域和频域上分离效果更佳,而 VMD 在处理高度非平稳信号时失效[35]。然而,EEMD 产生了超出预期的单分量数量,且部分分量并非 IMF(图 16)。
初始中心频率和带宽值在 epoch 为 0 时按表 5 和表 6 设定,各滤波器初始带宽均为 50。最终 IMFs 的中心频率分别为 91.84 Hz , 217.29 Hz 91.84 Hz , 217.29 Hz 91.84Hz,217.29Hz91.84 \mathrm{~Hz}, 217.29 \mathrm{~Hz} 和 449.90Hz,对应带宽为 80.22 Hz , 104.70 Hz 80.22 Hz , 104.70 Hz 80.22Hz,104.70Hz80.22 \mathrm{~Hz}, 104.70 \mathrm{~Hz} 和 297.53Hz。由于 VMD 在此例中失效,无需比较分解后 IMF 的正交性。以噪声干扰最强的 5dB 信噪比情况为例,展示了 IMF 选择的自适应滤波器(图 19)、神经网络损失曲线(图 20)及自适应滤波器稀疏系数(图 21)。原始 10 个滤波器最终缩减至 3 个。
图 22 对比了所提方法、EEMD 结合 HT 及 VMD 结合 HT 对三个单分量瞬时频率的识别结果。时变频率识别考虑了两种典型场景:无噪声原始信号及 5dB 信噪比环境。各方法与解析解之间的识别误差 ( e ) ( e ) (e)(e) 通过逐点计算得出。
  表4
合成信号的 DNN 相关参数 s s ss
  参数   数量
  初始滤波器数量 10
初始滤波器系数数量 10
  学习率 1
  动量 0.9
  最大训练轮数 1000
  批次大小 256
Parameter Number Initial filter numbers 10 Initial filter coefficient numbers 10 Learning rate 1 Momentum 0.9 Max epoch 1000 Batch size 256| Parameter | Number | | :--- | :--- | | Initial filter numbers | 10 | | Initial filter coefficient numbers | 10 | | Learning rate | 1 | | Momentum | 0.9 | | Max epoch | 1000 | | Batch size | 256 |
图 13. 所提方法在不同噪声水平下对合成信号 s s ss 的 IMF 分量对比:(a)无噪声,(b) SNR = 20 dB SNR = 20 dB SNR=20dB\mathrm{SNR}=20 \mathrm{~dB} ,(c) SNR = 10 SNR = 10 SNR=10\mathrm{SNR}=10 dB 及(d) ( d ) SNR = 5 dB ( d ) SNR = 5 dB (d)SNR=5dB(\mathrm{d}) \mathrm{SNR}=5 \mathrm{~dB}
图 14. 所提方法在不同噪声水平下对合成信号 s s ss 的 IMF 频谱对比:(a)无噪声,(b)信噪比 = 20 dB = 20 dB =20dB=20 \mathrm{~dB} ,(c)信噪比 = 10 dB = 10 dB =10dB=10 \mathrm{~dB} ,(d)信噪比 = 5 dB = 5 dB =5dB=5 \mathrm{~dB}
图 15. EEMD 方法在不同噪声水平下对合成信号 s s ss 的 IMF 分量对比:(a)无噪声,(b) SNR = 20 dB SNR = 20 dB SNR=20dB\mathrm{SNR}=20 \mathrm{~dB} ,(c) SNR = 10 dB SNR = 10 dB SNR=10dB\mathrm{SNR}=10 \mathrm{~dB} ,(d) SNR = 5 dB SNR = 5 dB SNR=5dB\mathrm{SNR}=5 \mathrm{~dB}
图 16. EEMD 方法在不同噪声水平下对合成信号 s s ss 的 IMF 频谱对比:(a)无噪声,(b) SNR = 20 dB SNR = 20 dB SNR=20dB\operatorname{SNR}=20 \mathrm{~dB} ,(c) SNR = 10 dB SNR = 10 dB SNR=10dB\mathrm{SNR}=10 \mathrm{~dB} ,(d) SNR = 5 dB SNR = 5 dB SNR=5dB\mathrm{SNR}=5 \mathrm{~dB}
图 17. 合成信号 s s ss 经 VMD 在不同噪声水平下的 IMF 分量对比:(a)无噪声,(b)信噪比 = 20 dB = 20 dB =20dB=20 \mathrm{~dB} ,(c) SNR = 10 dB SNR = 10 dB SNR=10dB\mathrm{SNR}=10 \mathrm{~dB} ,(d)信噪比 = 5 dB = 5 dB =5dB=5 \mathrm{~dB}
实例中,结果显示无噪声和 5dB 信噪比时,所提方法的数值分别为 1.48 % 1.48 % 1.48%1.48 \% 4.90 % 4.90 % 4.90%4.90 \% ,而 EEMD 为 40.43 % 40.43 % 40.43%40.43 \% 53.12 % 53.12 % 53.12%53.12 \% ,VMD 为 30.37 % 30.37 % 30.37%30.37 \% 36.56 % 36.56 % 36.56%36.56 \% 。所提方法即使在 5dB 信噪比下也能更精确识别瞬时频率(图 22(b))。EEMD 由于冗余 IMF 产生了更多瞬时频率,且误差随噪声增大而增加。VMD 因图 17 中获取的 IMF 不准确而失效。
本例中,考虑 30 % 30 % 30%30 \% 频谱重叠以验证所提方法。土木工程结构的振动数据通常具有弱频谱重叠和窄频带特性,因此该方法适用于土木工程监测数据分析。
图 18. 合成信号 s s ss 经 VMD 在不同噪声水平下的 IMF 频谱对比:(a)无噪声,(b) SNR = 20 dB SNR = 20 dB SNR=20dB\mathrm{SNR}=20 \mathrm{~dB} ,(c) SNR = 10 dB SNR = 10 dB SNR=10dB\mathrm{SNR}=10 \mathrm{~dB} ,(d) SNR = 5 dB SNR = 5 dB SNR=5dB\mathrm{SNR}=5 \mathrm{~dB}
  表5
合成信号 s s ss 自适应滤波器的中心频率(Hz)。
训练周期 0 250 IMF
Epochs 0 250 IMF | Epochs | 0 | 250 | | :--- | :--- | :--- | | IMF | | |
500 750
1 300 369.65 439.90 448.69
2 200 213.51 216.64 216.01
3 100 90.94 90.91 91.41
"Epochs 0 250 IMF " 500 750 1 300 369.65 439.90 448.69 2 200 213.51 216.64 216.01 3 100 90.94 90.91 91.41| Epochs 0 250 <br> IMF | | 500 | 750 | | | :--- | :--- | :--- | ---: | ---: | | 1 | 300 | 369.65 | 439.90 | 448.69 | | 2 | 200 | 213.51 | 216.64 | 216.01 | | 3 | 100 | 90.94 | 90.91 | 91.41 |
  表6
合成信号 s s ss 自适应滤波器的带宽 ( Hz ) ( Hz ) (Hz)(\mathrm{Hz})
  本征模态函数(IMF)的迭代周期
Epochs IMF| Epochs | | :--- | | IMF |
0 250 500 750
1 50 130.61 228.24 297.38
2 50 71.68 85.14 104.53
3 50 62.84 69.93 79.50
"Epochs IMF" 0 250 500 750 1 50 130.61 228.24 297.38 2 50 71.68 85.14 104.53 3 50 62.84 69.93 79.50| Epochs <br> IMF | 0 | 250 | 500 | 750 | | :--- | :--- | ---: | ---: | ---: | | 1 | 50 | 130.61 | 228.24 | 297.38 | | 2 | 50 | 71.68 | 85.14 | 104.53 | | 3 | 50 | 62.84 | 69.93 | 79.50 |
图 19. 合成信号 s s ss IMF 分量提取的自适应滤波器。

3.3 实例3 电缆模型实验示例

如图 23 所示进行缆索实验以验证所提方法的瞬时频率估计精度,并分别采用所提方法、EEMD 和 VMD 对缆索振动信号的 IMF 分量进行分析。测试缆索置于缆索制造厂区坑道内,一侧安装液压千斤顶加载,另一侧设置紧固机构。缆索质量密度为 22 kg m 1 22 kg m 1 22kg*m^(-1)22 \mathrm{~kg} \cdot \mathrm{~m}^{-1} ,长度 191.952 米,两端锚固长度各 1.0 米。如图 23(a)所示,缆索上布置了五个加速度计。
图20. 合成信号 s s ss 的神经网络模型损失曲线。
图21. 合成信号 s s ss 的自适应滤波器稀疏系数。
试验中,缆索分三步张拉。首先进行预张拉阶段,施加 1250 kN 的力以测量缆索实际长度,随后释放。第二步,将缆索复位至初始状态后重新加载至 1250 kN,并逐步增加荷载至 1500 kN , 2000 kN 1500 kN , 2000 kN 1500kN,2000kN1500 \mathrm{kN}, 2000 \mathrm{kN} 及 2500 kN,同时保持荷载一段时间(见图 24)。最后对缆索卸载。由于高度限制及张力不足,缆索中部因重力下垂,影响了其振动模态。分析表明,当缆索力达到 2240 kN 时(1312 秒时刻),缆索脱离地面处于张紧状态。因此,仅选取 1312 秒至 1531 秒间的缆索力数据进行分析,如图 24 所示。
由 1 号加速度计监测的缆索振动数据波形与傅里叶谱如图 25 所示。所提方法重构信号与原输入信号在时域和频域上差异均较小。采用与图 3 相同的 DNN 结构,相关参数列于表 7。
在此实例中,15 个滤波器的初始带宽为 0.4,中心频率大致按 0.88 的倍数设置。冗余滤波器的系数初始值设为 1。信号长度为 4096 个监测离散点。图 26 和图 27 展示了所提方法、EEMD 和 VMD 分解得到的 IMF 分量,分别对应 9 个、4 个和 9 个 IMF。此处 VMD 的 IMF 数量人为设置为与所提方法相同,因为若预设模态数不准确,该方法会失效。值得注意的是,所提方法能良好分离 IMF 分量,频谱中未出现虚假单分量,而 EEMD 无法实现有效分离(图 27(b))。对于 VMD 而言,该方法未能捕获频率在 2.5 Hz , 5.5 Hz 2.5 Hz , 5.5 Hz 2.5Hz,5.5Hz2.5 \mathrm{~Hz}, 5.5 \mathrm{~Hz} 和 9.5Hz 附近的单分量。
提取出的 IMF 优化自适应滤波器如图 28 所示,神经网络模型的损失曲线展示于图 29。经过 1000 轮神经网络训练后,自适应滤波器从 15 个冗余滤波器缩减至 9 个,完全覆盖了 IMFs 的频谱范围,同时相关滤波器系数被稀疏优化,其中 9 个滤波器权重值接近 1(图 30)。
图 31. 缆索振动信号 IMFs 的正交矩阵:(a)本文方法,(b)EEMD 及(c)VMD。
图 31 检验了 IMFs 间的正交性,显示本文方法的 I O r , k I O r , k IO_(r,k)I O_{r, k} 值远小于 EEMD 和 VMD,这表明所提方法分解得到的 IMFs 具有更优的正交性。
接下来,在获得电缆信号的自适应稀疏 IMF 后,通过深度神经网络识别相关的瞬时频率。首先对电缆振动信号进行短时傅里叶变换的频谱分析,图 32 揭示了 9 个频率成分,其中能量最大的频率集中在 6 至 8 赫兹之间。选取了五个离散时间点(标记为 t 1 t 5 t 1 t 5 t_(1)∼t_(5)t_{1} \sim t_{5} ),均显示出在 6 8 Hz 6 8 Hz 6∼8Hz6 \sim 8 \mathrm{~Hz} 范围内较高的频率幅值。根据平坦张紧弦理论[49]计算得到时变电缆力的识别结果:
I F = 4 m L c 2 ( ω ~ k ( t ) 2 π k ) 2 I F = 4 m L c 2 ω ~ k ( t ) 2 π k 2 IF=4mL_(c)^(2)(( widetilde(omega)_(k)(t))/(2pi k))^(2)\mathbf{I F}=4 m L_{c}^{2}\left(\frac{\widetilde{\boldsymbol{\omega}}_{k}(t)}{2 \pi k}\right)^{2}
式中 L c L c L_(c)L_{c} 表示电缆长度, m m mm 为线质量密度, ω ~ k ( t ) ω ~ k ( t ) widetilde(omega)_(k)(t)\widetilde{\boldsymbol{\omega}}_{k}(t) 是电缆的第 k k kk 阶固有频率。为验证该方法对多信号处理的可行性,基于加速度计监测信号识别的电缆力——
图32. 电缆振动信号的频谱图。右图显示了五个时间点的绝对振幅。
图 33. 由 1 号、3 号和 5 号加速度计监测信号识别的时变电缆力及其识别误差(e)。
图 33 对比了加速度计 1 号、3 号和 5 号的数据。识别出的索力按 IMF 分量顺序排列,其中最低固有频率对应最大阶数。各索力估计值与实测真实索力之间的相对误差通过以下公式计算:
e = I F k F 2 F 2 × 100 % e = I F k F 2 F 2 × 100 % e=(||IF_(k)-F||_(2))/(||F||_(2))xx100%e=\frac{\left\|\mathbf{I} \mathbf{F}_{k}-\mathbf{F}\right\|_{2}}{\|\mathbf{F}\|_{2}} \times 100 \%
其中 I F k I F k IF_(k)\mathbf{I F}_{k} 表示第 k k kk 个识别出的索力, F F F\mathbf{F} 为实测索力。计算了1号加速度计的识别误差。从最低到最高固有频率对应的九个识别索力相对误差分别为 2.80 % , 1.18 % , 1.16 % , 0.98 % , 1.93 % , 1.52 % , 1.70 % , 2.43 % 2.80 % , 1.18 % , 1.16 % , 0.98 % , 1.93 % , 1.52 % , 1.70 % , 2.43 % 2.80%,1.18%,1.16%,0.98%,1.93%,1.52%,1.70%,2.43%2.80 \%, 1.18 \%, 1.16 \%, 0.98 \%, 1.93 \%, 1.52 \%, 1.70 \%, 2.43 \% 2.18 % 2.18 % 2.18%2.18 \% 。本例中,神经网络参数在分析三个传感器信号时保持完全一致。三个信号源识别出的索力与实测值均具有良好吻合度,进一步验证了所提方法在多源数据下的可行性。
如图所示, I F 3 I F 3 IF_(3)\mathbf{I F}_{3} I F 4 I F 4 IF_(4)\mathbf{I F}_{4} 具有最高的识别精度,这与图 32 和图 25(b)中 6 8 Hz 6 8 Hz 6∼8Hz6 \sim 8 \mathrm{~Hz} 之间频率的最大能量一致,其中较高的振幅描绘了峰值和振荡细节。这表明 IMF 的能量越大,瞬时频率的估计就越准确。
为了说明 IMF 能量对瞬时频率估计精度的影响,图 34 展示了不同时间点 IMF 的振幅及相对误差。识别出的索力相对误差与频谱图中五个不同时间点电缆振动信号的绝对振幅相关联。结果显示,对于 I F 3 I F 3 IF_(3)\mathbf{I F}_{3} I F 4 I F 4 IF_(4)\mathbf{I F}_{4} ,观察到振幅最大时相对误差最小,这表明索力估计精度与 IMF 的能量高度相关。

  4. 结论

本文提出了一种基于深度学习的自适应模态分解与瞬时频率估计方法,用于结构健康监测中的非线性和非平稳信号分析。为实现无模态混叠和虚假分量的精确本征模态函数提取以及高精度的瞬时频率估计,构建了一种新型深度神经网络框架。在该网络中,采用自适应滤波器进行单分量频谱提取,并将希尔伯特变换嵌入网络层中以实现

图34. 索缆振动信号频谱图幅值与相对误差的关系。
瞬时频率估计。通过考虑频谱重叠和不同噪声干扰水平的数值信号,以及具有变化张力工况的长索模型实验振动数据,验证了所提方法的有效性。
该方法能够自适应地分解具有紧密间隔单分量的合成信号。此外,在示例中实现了 1.48 % 1.48 % 1.48%1.48 \% 的信号分解误差,对于频谱重叠较弱的单分量,该误差低于 EEMD 和 VMD。此外,该方法对噪声具有鲁棒性。在信噪比为 5 dB 的示例中,瞬时频率的识别误差为 4.90 % 4.90 % 4.90%4.90 \% ,而 EEMD 和 VMD 分别为 53.12 % 53.12 % 53.12%53.12 \% 36.56 % 36.56 % 36.56%36.56 \% 。此外,与 EEMD 和 VMD 相比,获得了更多正交分解的 IMF,合成信号和实验信号的正交性值分别为 10 5 10 5 10^(-5)10^{-5} 10 6 10 6 10^(-6)10^{-6}
从缆索振动数据的结果来看,时变缆力估计值与实测值之间的识别误差小于 2.80 % 2.80 % 2.80%2.80 \% 。同时研究了本征模态函数(IMF)能量与缆力识别误差的关系,能量较高的 IMF 会带来更高的识别精度。对于能量最高的 IMF,识别误差可低至 0.98 % 0.98 % 0.98%0.98 \% 。此外,所提出的基于深度神经网络的时频方法对实际监测数据具有更好的适应性。该 DNN 模型可重新训练,并能很好地适应各种结构的新输入数据。

  基金资助

本研究得到国家自然科学基金资助[项目编号:51978216、52192664]。

CRediT 作者贡献声明

郭一兵:概念化、方法论、软件、验证、形式分析、初稿撰写。鲍跃全:概念化、方法论、文稿审阅与编辑、监督指导。李惠:概念化、方法论、文稿审阅与编辑、监督指导。张玉峰:资源提供。

利益冲突声明

作者声明,对于可能影响本研究报告的已知竞争性经济利益或个人关系不存在。

  数据可用性

作者无权限共享数据。

  参考文献

[1] 鲍勇,陈志,魏松,徐毅,唐震,李华,结构健康监测中的数据科学与工程研究现状,《工程》5 (2019) 234 242 234 242 234-242234-242
[2] 李华楠,任力,贾志国,易天昊,李大胜,大型复杂土木基础设施结构健康监测研究进展,《土木结构健康监测杂志》6 (2016) 3-16。
[3] 欧进萍, 李惠, 中国大陆结构健康监测研究综述与未来趋势, 结构健康监测 9 (2010) 219-231.
[4] 侯锐, 夏勇, 土木工程结构基于振动的损伤识别新进展综述: 2010-2019, 声与振动杂志 491 (2021) 115741.
[5] 杨勇, 李素娟, S. Nagarajaiah, 李惠, 周培, 基于复杂度追踪的加速度实时输出唯一识别时变索力, 结构工程杂志 142 (2016) 04015083.
[6] B. Bhowmik, T. Tripura, B. Hazra, V. Pakrashi, 基于递归典型相关分析的鲁棒线性和非线性结构损伤检测, 机械系统与信号处理 136 (2020) 106499.
[7] S. Panda, T. Tripura, B. Hazra, 基于近实时包络提取的地震激励结构在线损伤检测, 结构健康监测. 21 (2022) 298-319.
[8] L. Li, H. Zhou, H. Liu, C. Zhang, J. Liu, 一种耦合经验模态分解与长短期记忆网络的混合方法用于预测 SHM 系统缺失测量信号数据, 结构健康监测. 20 (2021) 1778-1793.
[9] L. Lian, W. Ren, S. Wang, 分数阶傅里叶变换:时频表征与结构瞬时频率识别, 机械系统与信号处理. 178 (2022) 109305.
[10] D. Goyal, B. Pabla, 结构健康监测中的振动监测方法与信号处理技术综述, 计算机方法工程档案. 23 (2016) 585-594.
[11] 黄锷、沈正、龙慎仪等,《基于经验模态分解与希尔伯特谱的非线性非平稳时间序列分析方法》,《皇家学会会刊 A 辑-数学物理》454 卷(1998 年)903-995 页。
[12] 丹尼斯·盖博,《通信理论》,《电气工程师学会杂志》93卷(1946年)42-57页。
[13] 让·阿莱克西·维尔,《解析信号理论及其应用》,《电缆与传输》2卷(1948年)61-74页。
[14] 斯特凡·马莱,《小波基》,载于《信号处理的小波之旅:稀疏方法》(第3版),爱思唯尔出版社,伯灵顿,2009年,263-376页。
[15] W.J. Fitzgerald,《傅里叶分析》,载于:《电信工程师参考手册》,Butterworth-Heinemann 出版社,牛津,1993 年,第 2-17 页。
[16] H. Aied, A. González, D. Cantero,《基于集合经验模态分解的桥梁加速度响应中移动荷载引起的刚度突变识别》,《机械系统与信号处理》66 卷(2016 年)314-338 页。
[17] M. Lazhari, AJJoS Sadhu,《采用自适应经验模态分解法的结构分散模态识别》,《声音与振动》447 卷(2019 年)20-41 页。
[18] P. Singh, M. Keyvanlou, A. Sadhu,《一种基于有限传感器的结构状态评估改进时变经验模态分解方法》,《工程结构》232 卷(2021 年)111882 页。
[19] C. Delprete, E. Brusa, C. Rosso, F. Bruzzone, 基于正交经验模态分解的轴承健康监测, 《振动与冲击》2020 年 (2020).

    • 通讯作者单位:哈尔滨工业大学土木工程学院,哈尔滨 150090,中国。
    电子邮箱地址:baoyuequan@hit.edu.cn (Y. Bao).
    https://doi.org/10.1016/j.ymssp.2023.110463
    收稿日期:2022年10月31日;修回日期:2023年4月6日;录用日期:2023年5月20日
    2023年5月31日在线可获取
    0888-3270/© 2023 爱思唯尔有限公司版权所有。