-
舰船机械设备运转产生的振动会在海水中形成复杂的低频线谱噪声[1],该类噪声是反潜设备探测的主要目标特征。主动隔振可以有效抑制低频振动,目前国内外已展开诸多研究[2-4]。
在算法领域,传统滤波最小均方(filter-X least mean square, FXLMS)算法因结构简单且易于实现而被广泛应用[5],但对多线谱复杂激励等情况难以实现有效控制。为此,赵洪亮等[6]提出一种选频有源控制算法,通过构造二阶正交滤波结构实现了对多个线谱的控制。李彦等[7]提出一种多通道窄带Fx-Newton控制算法,用带通滤波器提取多个线谱分别予以控制。高伟鹏[8]提出了一种变步长的小波包自适应算法,将采集信号分解到不同带频上实现了多频信号控制。然而,上述算法都依赖具有精确频率的参考信号,在工程中存在传感器不易安装和通道耦合等问题,使得参考信号难以采集或者存在误差,即参考信号失配,从而影响控制效果。为此,张志谊[9]提出了通过测量信号自相关序列获得信号频率,翟晓军等[10]提出了利用快速傅里叶变换(FFT)方法对信号进行频谱估计。但是,上述算法均存在估计误差较大的问题,难以合成信噪比较高的参考信号。研究表明,基于自适应陷波器频率估计方法能够准确实现,且其陷波特性适合进行线谱控制[11-13]。
基于上述研究,本文将提出一种基于自适应陷波器反馈式多线谱控制算法,即首先将误差信号通过级联自适应陷波器,根据自适应算法更新陷波器参数,以估计多个信号频率,进行相位补偿并合成各参考信号。然后,将信号进行Hilbert变换得到另一路参考信号。最后,进入并行控制器完成幅值的更新,实现多线谱振动抑制,以解决参考信号失配影响。
-
陷波器是一种特殊带阻滤波器,可以实现对特定频率分量的剥离,但对其他频率分量的影响较小,图1(a)所示为其幅频特性。陷波器结构通常可用二阶IIR(infinite impulse response)滤波器表示,其传递函数为
$$H\left( {\textit{z}} \right) = \frac{{1 + a \cdot {{\textit{z}}^{ - 1}} + {{\textit{z}}^{ - 2}}}}{{1 + r \cdot a \cdot {{\textit{z}}^{ - 1}} + {r^2} \cdot {{\textit{z}}^{ - 2}}}}$$ (1) 式中:a为与陷波归一化频率相关的参数; r为小于且接近1的常数,可以影响陷波3 dB带宽。其中,r值越大,陷波张口越窄,即3 dB带宽越小,效果越理想。陷波幅值(陷波器的放大系数)随r值变化图如图1(b)所示。
实现陷波需要满足2个基本条件:一是,其传递函数零点须在单位圆上,使在陷波频率处陷阱深度无穷大;二是,其传递函数极点应与零点相匹配,使其他频率成分不受影响,其零极点分布如图1(c)所示,一对极点
${p_i}$ 和零点${{\textit{z}}_i}$ 的关系为$${{\textit{z}}_i} = r \cdot {p_i}$$ (2) 设在z平面单位时间走过的角度
$ \pm {\omega _{\rm{0}}}$ 处有一对共轭零点${{\textit{z}}_{1,2}} = {{\rm{e}}^{ \pm {\rm{j}}{\omega _0}}}$ ,根据式(2)可得${p_{1,2}} = r \cdot {{\rm{e}}^{ \pm {\rm{j}}{\omega _0}}}$ ,其传递函数可以表示为$$ \begin{split} & H\left( {\textit{z}} \right) = \frac{{\left( {{\textit{z}} - {{\rm{e}}^{{\rm{j}}{\omega _0}}}} \right) \cdot \left( {{\textit{z}} - {{\rm{e}}^{ - {\rm{j}}{\omega _0}}}} \right)}}{{\left( {{\textit{z}} - r \cdot {{\rm{e}}^{{\rm{j}}{\omega _0}}}} \right) \cdot \left( {{\textit{z}} - r \cdot {{\rm{e}}^{ - {\rm{j}}{\omega _0}}}} \right)}} =\\&\quad \frac{{1 - (2 \cdot \cos {\omega _0} \cdot {{\textit{z}}^{ - 1}}) + {{\textit{z}}^{ - 2}}}}{{1 - (2 \cdot r \cdot \cos {\omega _0} \cdot {{\textit{z}}^{ - 1}} )+ {r^2} \cdot {{\textit{z}}^{ - 2}}}} \end{split} $$ (3) 将式(1)和式(3)进行系数配比,得到
$$a = - 2\cos {\omega _0}$$ (4) -
根据式(4)陷波器传递函数相关参数a与角频率
${\omega _0}$ 之间的关系,设定a便可实现对特定频率的识别并提取。为避免上阶陷波器估计误差对下阶频率估计造成影响,用各阶陷波器的自身误差${\varepsilon _i}\left( n \right)$ 对其传递函数的参数${a_i}\left( n \right)$ (n为迭代次数)进行调整,以保证对激励频率精确估计。N阶级联自适应陷波器结构如图2所示,图中k 为第 k 个陷波器,vi(n)为各级陷波器的输入信号。由自适应算法写出迭代公式如下:
$${a_i}\left( {n + 1} \right) = a\left( n \right) - {\mu _i}{\nabla _i}\left( n \right)$$ (5) 式中,
${{\mu _i}}$ 为迭代步长;${\nabla _i}\left( n \right) $ 为均方误差梯度算符,即$${\nabla _i}\left( n \right) = \frac{{\partial \left[ {{\varepsilon ^2}_i\left( n \right)} \right]}}{{\partial {a_i}}} = 2{\varepsilon _i}\left( n \right)\frac{{\partial {\varepsilon _i}\left( n \right)}}{{\partial {a_i}}}$$ (6) 根据式(1)输入输出关系,可写为
$${\varepsilon _i}\left( n \right) = \frac{{{\rm{1 + }}{a_i}\left( n \right){{\textit{z}}^{ - 1}} + {{\textit{z}}^{ - 2}}}}{{1 + r{a_i}\left( n \right){{\textit{z}}^{ - 1}} + {r^2}{{\textit{z}}^{ - 2}}}} \cdot {v_i}\left( n \right)$$ (7) 将式(7)展开,可得
$$ \begin{split} & {\varepsilon _i}\left( n \right) + r{a_i}\left( n \right){\varepsilon _i}\left( {n - 1} \right) + {r^2}{\varepsilon _i}\left( {n - 2} \right) = \\&\quad {v_i}\left( n \right) + {a_i}\left( n \right){v_i}\left( {n - 1} \right) + {v_i}\left( {n - 2} \right) \end{split} $$ (8) 对式(8)关于
${a_i}\left( n \right)$ 求偏导,得到$$ \begin{split} & \frac{{\partial {\varepsilon _i}\left( n \right)}}{{\partial {a_i}}} + r{a_i}\frac{{\partial {\varepsilon _i}\left( {n - 1} \right)}}{{\partial {a_i}}} + r{\varepsilon _i}\left( {n - 1} \right){\rm{ + }} \\&\qquad {r^2}\frac{{\partial {\varepsilon _i}\left( {n - 2} \right)}}{{\partial {a_i}}} = {v_i}\left( {n - 1} \right) \end{split} $$ (9) 将式(5)~式(9)联立,得到线谱频率估计公式
$${a_i}\left( {n + 1} \right) = a\left( n \right) - {\mu _i}{v_i}\left( {n - 1} \right)\tilde \varepsilon \left( n \right)$$ (10) $$\tilde \varepsilon \left( n \right) \approx \frac{{\rm{1}}}{{{\rm{1 + }}r{a_i}\left( n \right) + {r^2}}}\left( {1 - \frac{{{\varepsilon _i}\left( {n - 1} \right)}}{{{v_i}\left( {n - 1} \right)}}} \right)$$ (11) -
若参考信号难以采集或者信噪比较大,将不利于前馈算法的开展和振动的有效控制。根据激励信号通过LTI(linear time-invariant system)系统(初级通道
$P\left( {\textit{z}} \right)$ 为线性时不变系统)其响应具有频率保持性的特点,因此,可基于反馈式控制结构从误差信号提取频率信息来实现自适应控制。本文所提算法框图如图3所示,该算法是将误差信号经级联自适应滤波器输出
${a_i}\left( n \right)$ 进入信号合成器。该算法中,相位信息对于振动控制的实现同样至关重要。在与参考信号相位值存在误差的情况下,算法控制效果较差,甚至会失效,故应对合成信号进行相位补偿。基于二阶自回归模型对参考信号${{{\mathit{\boldsymbol{x}}}}_{i1}}\left( n \right)$ 进行调节[14],即下式来实现对未知参考信号的相位估计。$${{{\mathit{\boldsymbol{x}}}}_{i1}}\left( n \right){\rm{ = }}{a_i}\left( n \right){{{\mathit{\boldsymbol{x}}}}_{i1}}\left( {n - {\rm{1}}} \right) - {{{\mathit{\boldsymbol{x}}}}_{i1}}\left( {n - 2} \right)$$ (12) 同时,对
${{{\mathit{\boldsymbol{x}}}}_{i1}}\left( n \right)$ 这一路参考信号进行Hilbert变换,得到另一路参考信号${{{\mathit{\boldsymbol{x}}}}_{i{\rm{2}}}}\left( n \right)$ $${{{\mathit{\boldsymbol{x}}}}_{i{\rm{2}}}}\left( n \right) = H\left[ {{{{\mathit{\boldsymbol{x}}}}_{i1}}\left( n \right)} \right] = \int_{{\rm{ - }}\infty }^{{\rm{ + }}\infty } {\frac{{{{{\mathit{\boldsymbol{x}}}}_{i1}}\left( n \right)}}{{{\text{π}} \left( {n - \tau } \right)}}{\rm{d}}\tau } $$ (13) 式中,Hilbert变换
$H\left[ {{{{\mathit{\boldsymbol{x}}}}_{i1}}\left( n \right)} \right]$ 是冲激响应为$h\left( t \right) = $ $ {1}/{{{\text{π}} t}}$ 的LTI系统,是实现90°相移的全通滤波器,适用于窄带信号;$\tau $ 为积分变量。两路合成参考信号分别进入控制滤波器进行幅值更新,输出信号表示为
$${y_i}\left( n \right) = {{{\mathit{\boldsymbol{w}}}}_{i1}}\left( n \right){{{\mathit{\boldsymbol{x}}}}_{i1}}\left( n \right) + {{{\mathit{\boldsymbol{w}}}}_{i2}}\left( n \right){{{\mathit{\boldsymbol{x}}}}_{i2}}\left( n \right)$$ (14) 式中,
${{{\mathit{\boldsymbol{w}}}}_{i1}}\left( n \right)$ 为第i个控制器权系数矢量。各控制单元总输出为
$$y\left( n \right){\rm{ = }}\sum\limits_{i = 1}^N {\left[ {{{{\mathit{\boldsymbol{w}}}}_{i1}}\left( n \right){{{\mathit{\boldsymbol{x}}}}_{i1}}\left( n \right) + {{{\mathit{\boldsymbol{w}}}}_{i2}}\left( n \right){{{\mathit{\boldsymbol{x}}}}_{i2}}\left( n \right)} \right]} $$ (15) 残余信号
$e\left( n \right)$ 为$$e\left( n \right) = d\left( n \right) - \sum\limits_{j = 0}^{M - 1} {{s_j} \cdot y\left( {n - j} \right)} $$ (16) 式中:
$d\left( n \right)$ 为期望信号;${s_j}$ 为真实次级通道的M阶系数。在实际中次级通道的传递特性是未知的,在控制前需对次级通道进行离线辨识。采用相同M阶数的FIR(finite impulse response)建模滤波器进行估计,得到的辨识模型如下:
$${\hat{{\mathit{\boldsymbol{s}}}}}\left( n \right)={{\left[ {{{\hat{s}}}_{0}}\left( n \right),{{{\hat{s}}}_{{1}}}\left( n \right),\cdots ,{{{\hat{s}}}_{M-1}}\left( n \right) \right]}^{\rm{T}}}$$ (17) 式中,
${\hat{{\mathit{\boldsymbol{s}}}}}\left( n \right)$ 为辨识次级通道特性。滤波参考信号${\hat{{\mathit{\boldsymbol{r}}}}_i}\left( n \right)$ 为:$$ \begin{split} & {{\hat{{\mathit{\boldsymbol{r}}}}}_{i1}}\left( n \right) = \sum\limits_{j = 0}^{M - 1} {{{\hat s}_j}{{{\mathit{\boldsymbol{x}}}}_{i1}}\left( {n - j} \right)} \\& {{\hat{{\mathit{\boldsymbol{r}}}}}_{i{\rm{2}}}}\left( n \right) = \sum\limits_{j = 0}^{M - 1} {{{\hat s}_j}{{{\mathit{\boldsymbol{x}}}}_{i2}}\left( {n - j} \right)} \end{split} $$ (18) 根据最小方差(least mean square, LMS)算法,得到控制滤波器权系数
${{{\mathit{\boldsymbol{w}}}}_i}\left( n \right)$ 迭代公式为:$$\begin{split} & {{{\mathit{\boldsymbol{w}}}}_{i1}}\left( {n + 1} \right) = {{{\mathit{\boldsymbol{w}}}}_{i1}}\left( n \right) - 2{\mu _{\rm{1}}}e\left( n \right){{\hat{{\mathit{\boldsymbol{r}}}}}_{i1}}\left( n \right) \\& {{{\mathit{\boldsymbol{w}}}}_{i{\rm{2}}}}\left( {n + 1} \right) = {{{\mathit{\boldsymbol{w}}}}_{i{\rm{2}}}}\left( n \right) - 2{\mu _{\rm{2}}}e\left( n \right){{\hat{{\mathit{\boldsymbol{r}}}}}_{i{\rm{2}}}}\left( n \right) \end{split} $$ (19) -
对频率估计窄带控制算法的性能开展了仿真验证。将单通道主动隔振试验测得的数据转化为600,200阶横向FIR滤波器系数,作为初级通道和次级通道脉冲响应的函数,图4所示为其幅值特性。
图 4 初级通道和次级通道传递函数幅值图
Figure 4. Amplitude diagram of transfer function of the primary and secondary channels
采用频率分别为30,37,60和110 Hz幅值均为1的四谐波信号,叠加信噪比20 dB白噪声作为参考信号,采样率为10 kHz,各LMS滤波器为300阶横向FIR滤波器。两种算法的迭代步长均调整到最佳值,级联自适应陷波器阶数N设为4。以残余信号加速度响应作为评价标准,开启主动控制前、后的信号落幅即为隔振效果。图5~图9所示为算法仿真结果相关图。
图 8 原始信号与合成信号相位补偿前\后时历曲线
Figure 8. Time history of original signal and synthetic signal before and after phase compensation
图 9 本文算法与FXLMS算法输出信号图
Figure 9. Output signal diagram of the algorithm in this paper and the FXLMS algorithm
图5中FXLMS算法迭代到6万次时才完成收敛,在时域上残余信号降低75%,在频域上对30,37,60和110 Hz振动线谱分别仅有7.17,11.95,0和8.81 dB的幅值衰减,控制效果较差。图6可见,反馈式线谱控制算法针对多频信号控制效果较好,残余信号迭代到2万次时就完成收敛,且幅值降低90%以上。频域上对30,37,60和110 Hz线谱幅值分别降低了23.45,30.52,22.87和38.08 dB,各振动线谱衰减情况如表1所示。
表 1 仿真中两种算法振动衰减
Table 1. Vibration attenuation of two algorithms in simulation
振动线谱/Hz FXLMS算法衰减值/dB 反馈算法衰减值/dB 控制前 控制后 控制前 控制后 30 −6.34 −13.51 −6.34 −29.79 37 −6.12 −18.07 −6.12 −36.64 60 −22.5 −22.5 −22.5 −45.37 110 −27.97 −36.78 −27.97 −66.05 从图7可以看出,算法迭代到5千次左右就基本完成了对频率信息的精确识别,且实现了对频率的实时跟踪。图8所示为原始信号与合成信号的时历曲线,由图可见:在未进行相位补偿时,合成信号与原始信号之间存在相位差,在补偿后,合成信号可以较好地与原始信号吻合。图9所示为本文所提算法和FXLMS算法控制器输出信号图。由图可见:FXLMS算法针对多频激励,其输入信号的自相关矩阵特征值分布较广,收敛效果较差;而反馈式多线谱控制算法将不同频率信号进入并行控制器处理,从仿真结果看控制效果较好。
-
为进一步验证算法可行性,本研究采用了NI PXI搭载RT系统进行主动隔振试验。试验对象为电磁激振器及附加配重,用以模拟中小型机械设备。设备主要有NI LabVIEW控制系统、Copley Xenus伺服驱动器、PCB加速度传感器、电磁激振器、以及主被动一体化隔振器、WYK-20040K直流稳压稳流电源等。主被动一体化隔振器为自主研制,在低频激励范围内输出力与输入电流之间呈线性关系。试验现场图及示意图如图10所示。试验台架分为上、下两层,其中上层装有电磁激振器,下层板固定在地面上。下层板4个角位置设有立柱导轨,上层板通过滚珠轴承与立柱形成滑动配合。一体化隔振器与上、下层板的安装由螺栓连接固定。
加速度传感器采集振动源激励信号和下层误差信号(其中,振动源信号为FXLMS算法的参考信号,下层误差信号为残余振动),经PXI 4498 AD采集板卡进入PXI 8840自适应控制器。输出控制信号经PXI 6733 DA输出板卡进入功率放大器,然后驱动主被动一体化隔振器产生次级振动来抑制初级振动。上位机通过TCP/IP通信线与PXIe-1042Q机箱相连,将算法程序配置到下位机来控制作动器输出力。PXI-1062Q控制激振器的运行及振动信号的采集,以便实时监测及后期处理。
图 10 主动隔振试验现场图及示意图
Figure 10. Test site diagram and schematic diagram of active vibration isolation experiment
由信号发生器产生30,37,60,110 Hz四谐波信号作为初始激励信号。控制滤波器长度均为300阶。算法迭代步长各调整到最佳值。采样频率设定为10 kHz,既满足采样定理,又保证了对信号的实时跟踪。针对AD板卡采集信号存在的50和100 Hz工频干扰,试验中设计了Butterworth带阻滤波器进行滤除。以下层误差信号传感器采集的加速度响应作为评价标准,开启主动控制前、后的信号落幅即为隔振效果。得到时频域控制效果如图11和图12所示。
从图11(a)得到的残余振动信号在开启控制后6 s实现收敛,速度较慢,幅值仅降低70%;图11(b)频域中对30,37,60和110 Hz振动线谱分别有11.28,5.03,5.33和29.92 dB的幅值衰减,其中对30和110 Hz线谱控制效果较好,而对37和60 Hz线谱控制效果较差。图12(a)得到残余振动信号在开启控制后2 s就实现收敛,速度较快,幅值降低了90%左右;图12(b)频域中对30,37,60和110 Hz振动线谱分别有19.91,25.93,20.28和31.39 dB的幅值衰减,控制效果较好。各振动线谱衰减情况如表2所示。通过开展多线谱主动隔振试验,验证了所提算法的有效性。
表 2 试验中两种算法振动衰减的比较
Table 2. Comparison of vibration attenuation by two algorithms in experiment
振动线谱/Hz FXLMS算法衰减值/dB 反馈算法衰减值/dB 控制前 控制后 控制前 控制后 30 −35.1 −46.38 −35.1 −55.01 37 −36.46 −41.49 −36.46 −62.39 60 −38.88 −44.21 −38.88 −59.16 110 −37.51 −67.43 −39.48 −70.87 -
本文针对多线谱振动和参考信号失配的问题,提出了一种反馈式多线谱控制算法。首先,误差信号通过级联自适应陷波器,根据自适应算法更新陷波器参数,以估计多个信号频率,然后进行相位补偿合成各参考信号,再将信号经Hilbert变换得到另一路参考信号,最后进入并行控制器完成线谱控制。
通过对该算法的仿真和试验,实现了对30,37,60和110 Hz激励信号的精确辨识,合成了可靠的参考信号进入算法,并对各线谱有20~40 dB的幅值衰减;
本文所提算法较好地解决了振动控制中参考信号失配和多线谱振动的问题,有效减弱和抑制了低频振动能量的传递。
Multi-line spectrum feedback control algorithm for active vibration isolation
-
摘要:
目的 针对振动控制中多频激励的传统自适应滤波算法控制效果不佳,以及工程中传感器不易安装和通道耦合等原因导致参考信号失配的问题,提出了一种反馈式多线谱控制算法。 方法 首先,将误差信号通过级联自适应陷波器并根据自适应算法更新陷波器参数来估计多个信号频率;然后,合成各参考信号,进行相位补偿,通过Hilbert变换得到另一路参考信号; 最后,进入并行控制器完成幅值更新,实现振动控制。 结果 仿真和试验结果表明,该算法能够精确估计频率信息,合成可靠参考信号,对30,37,60和110 Hz线谱均取得了20~40 dB能量衰减。 结论 该算法能较好解决振动控制中参考信号失配和多线谱振动的问题,有效减弱和抑制低频振动能量传递。 Abstract:Objectives Aiming at the poor control effect of the traditional adaptive filtering algorithm for multi-frequency excitation in vibration control, and the engineering problems of difficult sensor installation and channel coupling that cause reference signal mismatch, this paper proposes a feedback multi-line spectrum control algorithm. Methods First, the error signal is passed through the cascaded adaptive notch filter, and the notch filter parameters are updated according to the adaptive algorithm to estimate multiple signal frequencies. Next, each reference signal is synthesized and phase compensation performed, then another reference signal is obtained through Hilbert transform and the parallel controller is finally entered to complete the amplitude update and realize vibration control. Results Through simulation and experimental verification, the results show that the proposed algorithm can accurately estimate frequency information, synthesize reliable reference signals and achieve 20–40dB energy attenuation for 30 Hz, 37 Hz, 60 Hz and 110 Hz line spectrums Conclusions This algorithm provides a better solution to the problems of reference signal mismatch and multi-line spectrum vibration in vibration control, and effectively reduces and suppresses low-frequency vibration energy transmission. -
表 1 仿真中两种算法振动衰减
Table 1. Vibration attenuation of two algorithms in simulation
振动线谱/Hz FXLMS算法衰减值/dB 反馈算法衰减值/dB 控制前 控制后 控制前 控制后 30 −6.34 −13.51 −6.34 −29.79 37 −6.12 −18.07 −6.12 −36.64 60 −22.5 −22.5 −22.5 −45.37 110 −27.97 −36.78 −27.97 −66.05 表 2 试验中两种算法振动衰减的比较
Table 2. Comparison of vibration attenuation by two algorithms in experiment
振动线谱/Hz FXLMS算法衰减值/dB 反馈算法衰减值/dB 控制前 控制后 控制前 控制后 30 −35.1 −46.38 −35.1 −55.01 37 −36.46 −41.49 −36.46 −62.39 60 −38.88 −44.21 −38.88 −59.16 110 −37.51 −67.43 −39.48 −70.87 -
[1] 朱石坚, 楼京俊, 何其伟, 等. 振动理论与隔振技术[M]. 北京: 国防工业出版社, 2006. ZHU S J, LOU J J, HE Q W, et al. Vibration theory and vibration isolation[M]. Beijing: National Defense Industry Press, 2006 (in Chinese). [2] 陈绍青, 王永, 魏璀璨, 等. 一类多频线谱激励下的主动隔振控制方法[J]. 振动与冲击, 2012, 31(23): 128–131, 190. CHEN S Q, WANG Y, WEI C C, et al. An active vibration isolation control method under excitation of multi-line spectra[J]. Journal of Vibration and Shock, 2012, 31(23): 128–131, 190 (in Chinese). [3] LI Y, HE L, SHUAI C G, et al. Improved hybrid isolator with maglev actuator integrated in air spring for active-passive isolation of ship machinery vibration[J]. Journal of Sound and Vibration, 2017, 407: 226–239. doi: 10.1016/j.jsv.2017.07.007 [4] 牛宁, 刘钟, 孙玲玲. 考虑作动器输出约束的分散中间质量混合隔振系统建模分析[J]. 中国舰船研究, 2019, 14(6): 155–163. doi: 10.19693/j.issn.1673-3185.01505 NIU N, LIU Z, SUN L L. Modeling and analysis of decentralized intermediate mass hybrid vibration isolation system considering actuator output constraints[J]. Chinese Journal of Ship Research, 2019, 14(6): 155–163 (in Chinese). doi: 10.19693/j.issn.1673-3185.01505 [5] ELLIOT S J. Signal processing for active control[M]. London: Academic Press, 2001. [6] 赵洪亮, 徐健, 李晓东, 等. 一种改进的选频有源噪声控制算法[J]. 声学学报, 2004, 29(6): 499–506. ZHAO H L, XU J, LI X D, et al. A modified algorithm of the active noise control system based on frequency-selective filter[J]. Acta Acustica, 2004, 29(6): 499–506 (in Chinese). [7] 李彦, 何琳, 帅长庚. 多通道窄带Fx-Newton时频算法及动力机械有源隔振实验[J]. 声学学报, 2015, 40(3): 391–403. LI Y, HE L, SHUAI C G. MIMO Fx-Newton narrowband algorithm and experiment of active vibration isolation on diesel generator[J]. Acta Acustica, 2015, 40(3): 391–403 (in Chinese). [8] 高伟鹏, 贺国, 刘树勇, 等. 小波包变换自适应振动控制方法[J]. 控制与决策, 2021, 36(1): 83–89. GAO W P, HE G, LIU S Y, et al. Adaptive vibration control method based on wavelet packet transform[J]. Control and Decision, 2021, 36(1): 83–89 (in Chinese). [9] 张志谊, 王俊芳, 周建鹏, 等. 基于跟踪滤波的自适应振动控制[J]. 振动与冲击, 2009, 28(2): 64–67. ZHANG Z Y, WANG J F, ZHOU J P, et al. Adaptive vibration control with tracking filters[J]. Journal of Vibration and Shock, 2009, 28(2): 64–67 (in Chinese). [10] 翟晓军, 周波. 一种改进的插值FFT谐波分析算法[J]. 中国电机工程学报, 2016, 36(11): 2952–2958. ZHAI X J, ZHOU B. An improved interpolated FFT algorithm for harmonic analysis[J]. Proceedings of the CSEE, 2016, 36(11): 2952–2958 (in Chinese). [11] WANG H, SUN H L, SUN Y P, et al. A narrowband active noise control system with a frequency estimation algorithm based on parallel adaptive notch filter[J]. Signal Processing, 2019, 154: 108–119. doi: 10.1016/j.sigpro.2018.08.012 [12] DAI W H, QIAO C J, WANG Y K, et al. Adaptive cascaded notch filter for frequency estimation of multiple sinusoids[C]//Proceedings of the 2nd International Conference on Wireless Communication and Sensor Network. Changsha: World Scientific, 2016: 138-144. [13] SIRDI N K M, MONNEAU A, NAAMANE A. Adaptive notch filters for prediction of narrow band signals[C]//Proceedings of the 7th IEEE International Conference on Systems and Control (ICSC). Valencia, Spain: IEEE, 2018. [14] XIAO Y, MA L, KHORASANI K, et al. A new robust narrowband active noise control system in the presence of frequency mismatch[J]. IEEE Transactions on Audio, Speech, and Language Processing, 2006, 14(6): 2189–2200. doi: 10.1109/TASL.2006.872604 -