基于总体平均经验模态分解的主动噪声控制系统研究

罗磊　黄博妍　孙金玮　温良

基于总体平均经验模态分解的主动噪声控制系统研究

罗磊1黄博妍1孙金玮1温良1

为了提高宽窄带混合噪声的消噪效果，本文提出一种基于总体平均经验模态分解（Ensemble empirical mode decomposition，EEMD）的主动噪声控制（Active noise control，ANC）系统，利用实时EEMD算法逐段将混合噪声分解成若干个固有模态函数（Intrinsic mode functions，IMF）分量.因为这些IMF分量的频带各不相同，所以实现了混合噪声中宽带分量和窄带分量的有效分离，独立进行ANC处理后成功解决了处理混合噪声时带来的“火花”现象，而且避免了传统混合ANC（Hybrid ANC，HANC）系统中频率失调的影响.EEMD算法也是对混合噪声的平稳化处理过程，因此当混合噪声中出现非平稳变化时，本文提出的系统也能保持较好的系统稳定性.通过不同噪声环境下进行仿真分析，提出的ANC系统比HANC系统具有更好的系统稳定性和更小的稳态误差.

混合噪声，主动噪声控制，总体平均经验模态分解，固有模态函数，非平稳变化

引用格式罗磊，黄博妍，孙金玮，温良.基于总体平均经验模态分解的主动噪声控制系统研究.自动化学报，2016，42（9）：1432-1439

随着科技的巨大发展，更多的产业开始实现半机械化甚至全机械化，创造更多效益的同时却给人们带来了更多噪音困扰，这些噪音已经严重影响了人们的身心健康，特别是对老人和婴儿，在某种程度上，它已经成为了一种隐形的杀手.

大部分机械噪声都是由其内部的旋转机构引起的，这种噪声的频率主要集中在低频段，传统的被动噪声控制（Passive noise control，PNC）系统不能有效地对其进行抑制，这时主动噪声控制（Active noise control，ANC）系统应运而生，它对低频噪声作用效果显著［1］.ANC系统是一种基于声波相交干涉原理的电声学装置，通过引入一个与待处理噪声（初级噪声源）幅度相同、频率相同、相位相反的噪声（次级噪声源），使其达到降噪或者消噪的目的［2］.

随着ANC系统的诞生，大量的结构和算法也相应地被提出，其中基于滤波-X最小均方（Filtered-X least mean square，FXLMS）算法的结构在ANC研究领域中最受欢迎而且也最有效［3-4］.ANC系统按照噪声的频谱特性可以分为窄带ANC（Narrowband ANC，NANC）和宽带ANC（Broadband ANC，BANC）.这两种系统对其相应带宽的噪声消噪效果明显，但是应对宽窄带混合的噪声时，会出现“火花”现象，造成系统不稳定（这里是利用BANC系统处理宽窄带混合噪声，因为NANC系统不能处理宽带噪声）.生活中常见的噪声基本都具有宽窄带混合特性，例如风扇噪声、电机噪声以及切割机噪声等，为此，文献［5］中提出了一种混合ANC（Hybrid ANC，HANC）系统，理想情况下可以完整分离混合噪声中的宽带分量和窄带分量，然后分别利用BANC和NANC对其相应分量进行单独处理，这种系统在处理宽窄带混合噪声上取得了巨大的进步，同时也是处理宽窄带混合噪声最典型的结构.但是这种HANC系统主要存在两个问题：1）该系统中NANC子模块的输入频率是由非声学传感器（转速计、加速度计等）获取得到，因此很容易与真实的频率出现偏差，即频率失调（Frequency mismatch，FM），这种现象会造成混合系统的不稳定，使得噪声处理效果变差［6-7］.文献［8-10］给出了FM的几种解决方法，这些方法主要应对10%以内的FM，而且随着FM的增大，系统收敛速度会变慢，稳态误差也会随之增大，当FM大于10%时，这些方法不再适用；2）该系统应对非平稳的混合噪声时效果不佳.

根据HANC系统处理宽窄带混合噪声的思想：预先分离出混合噪声中的宽带分量和窄带分量，然后对各分量进行独立的ANC处理，本文提出了一种基于总体平均经验模态分解（Ensemble empirical mode decomposition，EEMD）的ANC系统.本文是第一次将EEMD算法应用到ANC领域，并且提出了实时EEMD算法，通过实时EEMD算法和ANC系统的有效结合，自适应地分离出混合噪声中的宽带分量和窄带分量，独立进行ANC处理后不仅可以消除平稳的混合噪声，而且对非平稳的混合噪声也有明显的抑制效果.与HANC系统相比，该系统没有FM的风险，而且混合噪声的处理效果更好.

EEMD是经验模态分解（Empirical mode decomposition，EMD）的改进形式.EMD是1998年由Huang等提出的一种时频数据分析方法［11］，可以认为是一种自适应的二进制带通滤波器组（由多个带通滤波器组成，这些带通滤波器的形状固定，但是相邻两个滤波器的频率范围相互重叠，而且一个是另外一个的一半或者二倍）［12］.EMD算法具有自适应多尺度分解特性，利用信号本身的时间尺度特征，通过一定的筛选条件将一个复杂信号中不同时间尺度的波动或趋势逐级分解为若干个固有模态函数（Intrinsic mode function，IMF）和一个余量［13-14］，每个分量的主要频带范围都不同，因此该方法可以实现混合噪声中宽窄带分量的分离，而且不需要估计窄带频率，因此不受FM的影响.EMD的提出主要是为了解决非线性、非平稳问题，可以实现对原始信号和数据的平稳化处理.EMD经常和Hilbert变化结合使用，通常称为希尔伯特黄变换（Hilbert-Huang transform，HHT）.传统的EMD有时会有模态混叠（即没有合适的时间特征尺度对不同的模态成分进行有效地分离，使得在一个固有模态函数中含有不只一个模态分量，从而不能清楚地反映信号的内在固有性质）.模态混叠现象一旦产生，会对信号的后续分解产生影响，使经过EMD分解后得到的结果不能完整地反映真实的物理过程［15-17］.受Flandrin和Gledhill在噪声辅助分析方面研究的启发，2009年，Huang等提出了总体平均经验模态分解（EEMD）［15］，该方法利用白噪声频谱的均匀分布，使不同尺度的信号自动地映射到合适的参考尺度上，可以得到合适的时间尺度对信号进行分解，成功地避免模态混叠的影响.EEMD算法的应用条件和EMD算法相同，但是需要的计算量更大［11］.

通过仿真表明，提出的系统在不引入FM的情况下也能很好分离出混合噪声中的宽带分量和窄带分量，从而解决混合噪声问题，同时由于EEMD对信号的平稳化处理（消除骑行波，平滑不均匀幅度，使其更加对称），当噪声中出现非平稳变化时，该系统也能保持较高的系统稳定性.

1　传统宽窄带混合噪声控制系统

由宽窄带分量组成的原始混合噪声信号可以表示为：

其中xn（n）是一个包含q种频率成分的窄带噪声信号，是其傅里叶系数.xw（n）表示宽带噪声信号，一般是均值为零、方差为的白噪声或有色噪声.

为了解决传统BANC在处理宽窄带混合噪声时的不稳定现象，Xiao等提出了一种混合主动噪声控制（HANC）系统［5］，其框图如图1所示，主要包含NANC子系统、BANC子系统以及正弦噪声消除器（Sinusoidal noise canceller，SNC）子系统三个部分，其中NANC子系统和BANC子系统分别针对噪声中的窄带分量和宽带分量，SNC子系统用来去除混合噪声中的窄带部分，使进入BANC中的噪声都为宽带噪声.P（z）为初级路径的传递函数

混合噪声经过SNC子系统后，理想情况下会变成一个纯净的宽带信号x（n），然后利用BANC子系统对其进行处理，产生次级噪声中的宽带分量yB（n）.“Sync signal”通过实时的信号传感器（如加速度计）获得，用来计算窄带信号的估计频率.次级噪声中的窄带分量yN（n）由NANC子系统生成.初级噪声p0（n）则由混合噪声经过初级通道P（z）直接产生，此时误差信号表示为

其中vp（n）为消噪点的背景噪声.

图1　HANC系统Fig.1HANC system

HANC系统虽然解决了混合噪声造成的系统不稳定现象，但是同时也引入了FM问题，给该系统带来了新的不稳定因素，FM的定义如下：

其中ωi为真实频率，为其估计值.

由于FM的影响，HANC系统中的SNC子系统不仅无法准确地去除混合噪声中的窄带分量，而且还会引入新的窄带分量，此时所有的窄带成分和噪声中的宽带成分最后都被送入到BANC子系统中同时进行处理，由于混合噪声中窄带分量的增多，使得整个系统变得更加不稳定，“火花”现象也更加严重.

2　本文提出的系统

解决宽窄带混合噪声的关键是分离出其中的宽带分量和窄带分量，并对它们单独进行处理，才能避免相互之间的影响，传统的HANC系统是一种精确的分离方式，即将窄带分量和宽带分量进行完整分离，然后对其单独处理.该方法虽然可以解决混合噪声问题，但是却会引入FM现象，严重影响ANC系统的处理效果.为此本文提出一种宽窄带分量模糊分离的方法，利用EEMD算法对宽窄带混合噪声进行分解，得到其所有IMF分量，即将混合噪声按频率从高到低分解到各个IMF分量中，此时宽带分量和窄带分量被分解到不同的IMF中.由于每个分量包含的是某一段频率的信号，而不是某个确定频率信号，所以称该分离方式为模糊分离.该方法不需要估计窄带分量的频率，因此也不受这种估计误差带来的FM影响.

本文中提出的基于总体平均经验模态分解的主动噪声控制（EEMDANC）系统如图2所示，主要包含两个部分：EEMD算法和若干个BANC子系统（具体数量由被处理噪声的复杂程度决定），利用EEMD算法将混合噪声分解成若干个分量，然后通过BANC子系统对每个分量进行独立消噪.

图2　EEMDANC系统Fig.2EEMDANC system

2.1实时EEMD算法处理

EEMD算法的基本原理为：给信号中加入高斯白噪声将信号和噪声组成一个“总体”，加入信号的白噪声将遍布整个时频空间，与被滤波器组分离的不同尺度分量相一致.当信号加在这些一致分布的白噪声背景上时，不同尺度的信号，自动地映射到合适的参考尺度上.而且由于零均值白噪声的特性，加入的次数越多，将这些多次分解的结果取平均后，白噪声最终将被互相抵消而消除.总体平均的结果即被当作真实信号.加噪声的次数越多，平均后的结果所含的剩余白噪声越小，结果就越接近真实值.EEMD算法是EMD算法的改进形式，同样假设任何信号都是由若干个IMF组成的，每一个IMF必须满足下面两个条件［13-14］：1）在整个信号上，极值点的个数和过零点的个数相差不大于1；2）在任意点处，上下包络的均值为0.

EEMD算法的实现步骤如下：

1）对待分解的宽窄带混合噪声xr（n）重复加入白噪声NE次，每次加入的白噪声都不同，幅度用νi表示，一般为0.1～0.4.通过在混合噪声中添加白噪声来平滑脉冲干扰，使脉冲干扰在EMD分解过程中融入到白噪声模式中，抑制脉冲干扰与其他固有模态之间的混叠.

2）对每次加入白噪声后的混合噪声进行EMD分解.

a）首先确定加入白噪声后的混合噪声xi（n）的极大值和极小值，然后通过三次样条插值函数拟合出极大值包络emax（n）和极小值包络emin（n）.

b）计算整个包络线的均值，用mi，1（n）表示，并得到xi（n）与mi，1（n）的差值hi，1（n），hi，1（n）为去掉低频的新信号.

c）hi，1（n）一般不满足IMF定义的两个条件，需要重复上述过程多次，直到其满足IMF的定义，此时得到第一个IMF分量ci，1（n）.但是实际中，上下包络的均值无法为零，Huang把标准偏差系数（Standard deviation，SD）作为停止准则，即当SD满足下式时，就认为包络的均值满足IMF的均值为零的条件［14］.

式中NT为待分解信号的总长度，ε为筛分门限，一般取值在0.2到0.3之间.

d）用加入白噪声后的混合噪声信号减去ci，1（n），得到一个去掉高频成分的余量ri，1（n）.判断ri，1（n）是否满足EMD的停止准则，如果满足则EMD结束，否则将ri，1（n）看成一个新的信号重复步骤a）～d）.EMD停止准则的判断标准通常有两种：一种是当余量为常数或单调函数时，从中不能再筛分出其他IMF时停止；另一种是当最后一个IMF分量或余量的幅度小于预设值时停止.最后得到每组加入白噪声后信号的EMD结果为

其中N为IMF的总个数，ri，N（n）是一个趋势项，表示信号的平均趋势或均值.

3）将上述每次分解得到的IMF分量及其余量对应求平均，最终得到的IMF分量和余量分别为

此时，就得到原始混合噪声xr（n）的EEMD结果

为保证信号分解过程中的实时性，这里采用实时的EEMD算法，即通过一个窗口W对输入信号从左到右进行实时分割，分割得到的部分利用EEMD对其分解.同时为了确保分解得到的每组IMF具有连续性，每次分割信号的窗口是连续移动的，整个过程的示意图如图3所示.

图3　实时EEMD示意图Fig.3The sketch map of real-time EEMD

2.2BANC子系统的FXLMS算法处理

利用实时EEMD对输入的混合噪声进行分段分解，每段的长度都为W，每段都可以被分解成N个IMF分量和一个余量，再把得到的IMF分量和余量分别看成新的噪声信号，利用FXLMS算法对其进行计算处理，得到其各自的次级噪声.总的次级噪声为

式（27）中的α是控制函数形状的常数，β是控制函数范围的常数，δj（n）表示每一个分量在每段信号中占的能量比，δj（n）e（n）则近似表示每个分量产生的误差.

2.3EEMDANC系统稳定性分析

EEMDANC系统是由多个BANC子系统组成，只有当每个子系统都稳定时，EEMDANC系统才能稳定.将式（23）代入式（22），并对其两端取期望

3　仿真

在仿真1和仿真2中，信号的采样频率都为2000Hz.初级路径P（z）和次级路径S（z）都是截止频率为0.4π（400Hz）的低通滤波器，次级路径训练的信号为白噪声，最终取其稳态时的1000个点的均值作为次级路径估计的系数.次级路径估计时的参数设置如表1所示.

表1　次级路径估计参数设置Table 1Parameters of secondary path estimation

3.1仿真1：平稳变化的混合噪声

在该仿真过程中，混合噪声中窄带分量的频率为ω=［0.03 0.06 0.09］π（因为采样率为2000Hz，该频率也可表示为（30Hz 60Hz 90Hz）），其中各频率的傅里叶系数分别为ar，i=［2 1 0.5］，br，i=［-1-0.5 0.1］（1≤i≤3），宽带分量为一均值为0、方差为1的白噪声.分别利用HANC系统以及提出的EEMDANC系统对该混合噪声进行处理.

根据实际情况，HANC在窄带频率估计时会出现FM现象，本次仿真中，FM分别为1%、5%、10%，对应的BANC子系统控制器步长都为0.0002，其他参数与文献［5］相同，在本文提出的EEMDANC系统中，实时EEMD算法每次分解的分量个数都为5（4个IMF分量和1个余量），分割信号的窗口长度为W=101，各分量控制器的固定步长分别为µc=［0.0004 0.0009 0.00009 0.00006 0.0003］，步长中其他参数分别为α=300，β=0.4，ξ=0.1，各分量主控制器的阶数都为51，即，在该仿真条件下，不同FM时 HANC系统的消噪结果分别如图4中的（a）、（b）、（c）所示.由于本文提出的EEMDANC系统没有频率估计，不受FM的影响，最终的消噪结果如图4（d）所示.该次仿真的点数都为30000，独立运行次数为40.

图4　平稳宽窄带混合噪声系统性能对比Fig.4Comparisons between systems＇performance for stationary mixture noise of wideband and narrowband

3.2仿真2：非平稳变化的混合噪声

分别利用HANC系统和提出的EEMDANC系统对非平稳的混合噪声进行消噪处理.仿真的参数设置为：总仿真点数为50000，独立运行次数为40次.在25000点时混合噪声中窄带分量的幅值和频率同时发生突变：幅度由ar，i，br，i变为-ar，i，-br，i，频率由［0.03 0.06 0.09］π（30Hz 60Hz 90Hz）变为［0.05 0.08 0.12］π（50Hz 80Hz 120Hz）.宽带分量保持不变.在本次仿真中，HANC的BANC子系统主控制器的步长为0.00025，提出的EEMDANC系统参数和HANC系统的其他参数设置与仿真1中一样.仿真结果如图5所示.

3.3仿真3：实际录制的混合噪声

在该仿真中，被处理噪声为实时录制的某地下室某房间排风扇产生的噪音，录制设备为NI公司的PXI-4472数据采集卡以及SKC公司的NP21麦克风.实际测量噪声的时域图和频域图分别如图6（a）和图6（b）所示，信号的采样频率为22kHz，测得窄带分量的频率成分主要包含11Hz（0.001π）和44Hz（0.004π）.EEMDANC系统控制器的固定步长为µc=［0.03 0.01 0.01 0.00009 0.001］，HANC系统中BANC子系统的控制器步长为0.0007，其他参数设置与仿真1中一样.消噪后的结果分别如图7和图8所示，其中图7为时域消噪结果，图8为频域消噪结果.

3.4仿真分析

1）当HANC系统引入FM后，该系统变得不再稳定，“火花”现象又会重新出现，而且变得更加严重，随着FM增大，“火花”的个数也逐渐增多.出现这种现象的原因是真实频率和估计频率之间存在误差，使得真正的窄带频率和估计的窄带频率都被送到了BANC子系统中进行处理，由于窄带频率个数的增多，导致“火花”现象更加严重.当FM减小时，估计频率和真实频率接近，产生的“火花”基本发生在同一处，所以“火花”个数减少，随着FM增大，真实频率和估计频率差距变大，各自产生的“火花”现象也逐渐分离，“火花”个数也随之增多.为了便于分析EEMDANC系统，这里给出平稳的混合噪声24950～25050点经过EEMD后时域图和频域图，分别如图9和图10所示.从图10中可以看出，宽窄带混合噪声被分解成5个不同频率段的分量：IMF1的带宽为300Hz以上；IMF2的带宽为100～600Hz；IMF3的带宽为50～300Hz；IMF4的带宽为20～150Hz；余量r4（n）的带宽为0～60Hz.这也充分说明EEMD类似于一个自适应二进制带通滤波器组.仿真中窄带信号的频率分别为30Hz，60Hz和90Hz，因此经过EEMD后窄带分量集中在IMF4和余量r4（n）中，宽带分量集中在IMF1～IMF3中，通过这种分解，同样有效分离了宽带分量和窄带分量，避免了宽窄带混合噪声同时处理带来的“火花”现象，而且该系统不需要估计窄带频率，因此也不存在FM带来的影响.仿真结果图4（d）也充分证明提出的EEMDANC系统在不需要估计窄带频率的条件下依然可以很好地处理宽窄带混合噪声.由于实际中FM现象很难避免，而且随着非声学传感器的老化，这种现象会更加严重，所以提出的EEMDANC系统比HANC系统具有更好的实用性.

图5　非平稳宽窄带混合噪声的系统性能对比Fig.5Comparisons between systems＇performance for non-stationary mixture noise of wideband and narrowband

图6　排风扇噪声Fig.6Exhaust fan noise

图7　排风扇噪声的消噪结果Fig.7Noise cancellation results of exhaust fan noise

图8　初级噪声和误差噪声的功率谱密度Fig.8PSD of primary noise and error noise

图9　平稳宽窄带混合噪声的EEMD时域图Fig.9Time domain EEMD view of stationary mixture noise of wideband and narrowband

图10　平稳宽窄带混合噪声的EEMD频域图Fig.10Frequency domain EEMD view of stationary mixture noise of wideband and narrowband

2）通过观察图5发现：当混合噪声发生非平稳变化时，HANC系统的收敛曲线发生迟滞现象，不能及时跟踪控制，需要经过重新调整才能再一次达到稳定状态，而提出的EEMDANC系统却可以正常工作，不受非平稳变化的影响.图11给出了非平稳变化混合噪声24950～25050点的EEMD图，原始混合噪声在25000点时出现突变，EEMDANC系统将这种突变分解到4个IMF分量和一个余量中. EEMD算法主要有两个优点：1）消除骑行波；2）平滑不均匀幅度，使其更加对称.这两个优点被称为EEMD的平稳化处理，本文就是利用这种平稳化处理过程来解决ANC系统中的非平稳现象.图12为混合噪声的Hilbert谱（上图）和经过EEMD分解后各分量的Hilbert谱（下图），从图中可以看出分解得到的各分量频率随着时间基本都是连续变化的，而且每个分量的能量变化都比较均匀，这说明EEMD分解后的各分量基本都趋于平稳，然后对这些平稳分量进行独立的ANC处理，就会减少甚至消除原始混合噪声非平稳变化带来的影响.

图11　非平稳宽窄带混合噪声的EEMD图Fig.11EEMD view of non-stationary mixture noise of wideband and narrowband

图12　混合噪声和EEMD后各分量的Hilbert谱Fig.12Hilbert spectrum of the mixture noise and the components after EEMD

3）通过观察图7和图8，提出的EEMDANC系统对排风扇噪声的消噪效果要明显优于HANC系统.这是由于实际排风扇噪声中窄带分量的频率复杂，而传统的HANC系统不能同时对所有的频率成分进行跟踪控制，只能对其主要频率进行作用，而且部分窄带分量的频率会随着风扇的转速不均匀实时变化，使得HANC系统不能对其频率做到精确估计.本文提出的EEMDANC系统是对信号频率的一种模糊分解，不需要对频率进行精确估计，因此也不受频率估计不完整和频率估计误差的影响，同时这种信号分解方式也简化了宽带分量，优化了系统对宽带噪声的消噪效果.

4　总结

本文主要研究参考噪声信号具有宽窄带混合特性的ANC系统，首先介绍了传统BANC系统以及Xiao等提出的混合ANC系统，探讨了这两种方法处理混合噪声的可能性及其存在的风险.同时提出一种基于总体平均经验模态分解（EEMD）的ANC系统，分析了其降低宽窄带混合噪声的机理.通过大量的仿真分析，提出的EEMDANC系统不仅能有效消除宽窄带混合噪声，而且也不会引入FM的影响，同时当噪声出现非平稳变化时，也能保证较高的系统稳态性.

References

1 Kuo S M，Morgan D R.Active noise control：a tutorial review.Proceedings of the IEEE，1999，87（6）：943-973

2 Wu L F，Qiu X J，Guo Y C.A simplified adaptive feedback active noise control system.Applied Acoustics，2014，81：40-46

3 Xiao Y G，Ikuta A，Ma L Y，Khorasani K.Stochastic analysis of the FXLMS-based narrowband active noise control system.IEEE Transactions on Audio，Speech，and Language Processing，2008，16（5）：1000-1014

4 Kajikawa K，Gan W S，Kuo S M.Recent applications and challenges on active noise control.In：Proceedings of the 8th International Symposium on Image and Signal Processing and Analysis.Trieste，Italy：IEEE，2013.661-666

5 Xiao Y G，Wang J.A new feedforward hybrid active noise control system.IEEE Signal Processing Letters，2011，18（10）：591-594

6 Kuo S M，Puvvala A B.Effects of frequency separation in periodic active noise control systems.IEEE Transactions on Audio，Speech，and Language Processing，2006，14（5）：1857-1866

7 Jeon H J，Chang T G，Kuo S M.Analysis of frequency mismatch in narrowband active noise control.IEEE Transactions on Audio，Speech，and Language Processing，2010，18（6）：1632-1642

8 Xiao Y G，Ma L Y，Ward R K.Fast RLS Fourier analyzers capable of accommodating frequency mismatch.Signal Processing，2007，87（9）：2197-2212

9 Xiao Y G，Ma L Y，Khorasani K，Ikuta A，Xu L.A filtered-X RLS based narrowband active noise control system in the presence of frequency mismatch.In：Proceedings of the 2005 IEEE International Symposium on Circuits and Systems. Kobe，Japan：IEEE，2005.260-263

10 Huang Bo-Yan，Chang Lin，Ma Ya-Ping，Sun Jin-Wei，Wei Guo.A new narrowband ANC system against nonstationary frequency mismatch.Acta Automatica Sinica，2015，41（1）：186-193（黄博妍，常林，马亚平，孙金玮，魏国.一种应对非平稳频率失调的窄带主动噪声控制系统.自动化学报，2015，41（1）：186-193）

11 Huang N E，Shen Z，Long S R，Wu M C，Shih H H，Zheng Q A，Yen N C，Tung C C，Liu H H.The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis.Proceedings of the Royal Society A：Mathematical，Physical and Engineering Sciences，1998，454（1971）：903-995

12 Flandrin P，Rilling G，Goncalves P.Empirical mode decomposition as a filter bank.IEEE Signal Processing Letters，2004，11（2）：112-114

13 Qu Cong-Shan，Lu Ting-Zhen，Tan Ying.A modified empirical mode decomposition method with applications to signal de-noising.Acta Automatica Sinica，2010，36（1）：67-73（曲从善，路廷镇，谭营.一种改进型经验模态分解及其在信号消噪中的应用.自动化学报，2010，36（1）：67-73）

14 Kang Chun-Yu，Zhang Xin-Hua.Adaptive noise reduction method based on EMD and ANC.Systems Engineering and Electronics，2008，30（5）：810-812（康春玉，章新华.一种基于EMD和ANC技术的自适应降噪方法.系统工程与电子技术，2008，30（5）：810-812）

15 Wu Z H，Huang N E.Ensemble empirical mode decomposition：a noise-assisted data analysis method.Advances in Adaptive Data Analysis，2009，1（1）：1-41

16 Wang Y H，Yeh C H，Vincent Young H-W，Hu K，Lo MT.On the computational complexity of the empirical mode decomposition algorithm.Physica A：Statistical Mechanics and its Applications，2014，400：159-162

17 Tang Jian，Chai Tian-You，Cong Qiu-Mei，Yuan Ming-Zhe，Zhao Li-Jie，Liu Zhuo，Yu Wen.Soft sensor approach for modeling mill load parameters based on EMD and selective ensemble learning algorithm.Acta Automatica Sinica，2014，40（9）：1853-1866（汤键，柴天佑，丛秋梅，苑明哲，赵立杰，刘卓，余文.基于EMD和选择性集成学习算法的磨机负荷参数软测量.自动化学报，2014，40（9）：1853-1866）

罗磊哈尔滨工业大学电气工程及自动化学院博士研究生.主要研究方向为信号与信息处理，主动噪声控制理论.E-mail：llei@hit.edu.cn（LUO LeiPh.D.candidate at the School of Electrical Engineering and Automation，Harbin Institute of Technology.His research interest covers signal and information processing，and active noise control theory.）

黄博妍哈尔滨工业大学电气工程及自动化学院讲师.主要研究方向为语音增强，主动噪声控制理论.E-mail：byhuang@hit.edu.cn（HUANG Bo-YanLecturer at the School of Electrical Engineering and Automation，Harbin Institute of Technology.Her research interest covers speech enhancement and active noise control theory.）

孙金玮哈尔滨工业大学电气工程及自动化学院教授.主要研究方向为生物医学传感器，主动噪声控制理论.本文通信作者.E-mail：jwsun@hit.edu.cn（SUN Jin-WeiProfessor at the School of Electrical Engineering and Automation，Harbin Institute of Technology.His research interest covers biomedical sensors，and active noise control theory.Corresponding author of this paper.）

温良哈尔滨工业大学电气工程及自动化学院博士研究生.主要研究方向为自适应信号处理，主动噪声控制理论.E-mail：lwen@hit.edu.cn（WEN LiangPh.D.candidate at the School of Electrical Engineering and Automation，Harbin Institute of Technology. His research interest covers adaptive signal processing and active noise control theory.）

A New ANC System Based on Ensemble Empirical Mode Decomposition

LUO Lei1HUANG Bo-Yan1SUN Jin-Wei1WEN Liang1

In order to obtain a better de-noising performance of mixture noise containing both wideband components and narrowband components，a new active noise control（ANC）system based on ensemble empirical mode decomposition（EEMD）is proposed in this paper.Real-time EEMD algorithm is used to decompose the mixture noise into several intrinsic mode functions（IMF）which have a different frequency range each other，so this decomposition can separate wideband components and narrowband components from the mixture noise adaptively.Each component controlled independently can not only process mixture noise without“firework”，but also avoid the frequency mismatch occurring in conventional hybrid ANC（HANC）system.The EEMD algorithm can smooth the mixture noise to make sure the proposed system has better stability when nonstationary phenomenon happens in mixture noise.Compared with HANC system，the proposed ANC system has better system stability and smaller steady-state error in processing different noise.

Mixture noise，active noise control（ANC），ensemble empirical mode decomposition（EEMD），intrinsic mode function，non-stationary change

Manuscript July 7，2015；accepted January 15，2016

10.16383/j.aas.2016.c150433

Luo Lei，Huang Bo-Yan，Sun Jin-Wei，Wen Liang. A new ANC system based on ensemble empirical mode decomposition.Acta Automatica Sinica，2016，42（9）：1432-1439

2015-07-07录用日期2016-01-15

国家自然科学基金（61171183，61471140），2012航天支撑基金（01320214），哈尔滨工业大学科研创新基金（HIT.NSRIF201614）资助

Supported by National Natural Science Foundation of China（61171183，61471140），2012 Aerospace Support Fund（01320214），Natural Scientific Research Innovation Foundation in Harbin Institute of Technology（HIT.NSRIF201614）

本文责任编委孙健

Recommended by Associate Editor SUN Jian

1.哈尔滨工业大学电气工程及自动化学院哈尔滨150001

1.School of Electrical Engineering and Automation，Harbin Insti-tute of Technology，Harbin 150001