简化窄带无次级路径建模有源控制方法与流程

一.
技术领域：
本发明涉及无次级路径建模的有源控制方法，针对高采样率系统窄带信号提出了一种简化窄带无次级路径建模的有源控制方法。二.
背景技术：
：大多有源控制方法需要获取次级声源和误差传声器之间传递函数信息，即次级路径建模。前馈有源控制算法常采用滤波最小均方(fxlms)算法，当真实次级路径传递函数和建模次级路径传递函数之间误差较大时，fxlms算法可能不稳定。zhou等为避免对次级路径建模，提出一种无次级路径建模算法(d.zhouandv.debrunner，“anewactivenoisecontrolalgorithmthatrequiresnosecondarypathidentificationbasedonsprproperty，”ieeetrans.signalprocessing，vol.55，no.5，pp.1719-1729，may2007)。这种算法无需对次级路径建模，在更新自适应滤波器系数的过程中，通过判断残余噪声能量降低与否，在0°和180°两个方向选择更新过程中步长前的符号，以降低算法残差。但当信号频率对应的次级路径传递函数的相位与自适应滤波器系数更新方向间的相位差接近±90°时，算法收敛较慢。wu等对之做了改进，将算法改在频域实施，且引入±90°两个更新方向(m.wu，g.chen，andx.qiu，“animprovedactivenoisecontrolalgorithmwithoutsecondarypathidentificationbasedonthefrequency-domainsubbandarchitecture，”ieeetrans.onspeechandaudioprocessing，vol.16，pp.1409-1419，nov2008)，提高了这种情况下zhou等算法的收敛速度。当参考信号为宽带信号时，zhou等和wu等的算法均指出可采用子带结构将参考信号、误差信号以及自适应滤波器均做子带分解，并通过子带参考信号和子带误差信号调整对应子带自适应滤波器权系数，之后通过权变换的方法将各子带自适应滤波器系数堆砌为全频带自适应滤波器系数。debrunner等(v.debrunner，l.debrunner，andl.wang，“subbandadaptivefileringwithdelaycompensationforactivenoisecontrol，”ieeetrans.signalprocessing，vol.52，pp.2932-2937，2004)将该子带结构应用到有源控制系统，通过子带滤波参考信号和子带误差信号更新全带自适应滤波器系数，分析子带误差信号组成，从而补偿由次级路径和分析滤波器组对误差信号造成的时延。gao等(m.gao，j.luandx.qiu，“asimplifiedsubbandancalgorithmwithoutsecondarypathmodeling，”ieee/acmtrans.onaudio，speech，andlanguageprocessing，vol.24，pp：1164-1174，2017)定义代价函数为各子带误差信号模平方之和，在时域直接更新全带自适应滤波器系数。将参考信号通过直达声时延产生时延参考信号，同时，将参考信号经由希尔伯特滤波器产生正交时延参考信号，在时域实现子带内0°，180°或±90°四个方向搜索最佳更新方向。由于直达声时延补偿部分次级路径脉冲响应，只需通过搜索两类子带参考信号最佳更新方向近似残余次级路径脉冲响应相位信息，故减少子带个数，缩短搜索过程占用时间。然而，该算法需要采用多项fft方法获取子带信号信息。中国专利zl201410114263.5公开一种简化无次级路径建模的有源控制方法，该方法在时域实施，通过希尔伯特滤波器产生正交时延参考信号和时延参考信号一起在0°，180°两个方向搜索最佳更新方向，时域实现±90°两个搜索方向，且以估计直达声时延补偿次级路径信息，减少子带个数，结构简单。然而该方法需要通过多相fft方法获取子带信号信息，当在高采样率有源控制系统中，参考信号为窄带噪声时，参考信号频率范围内有效频率柜较少，该方法提取子带信号会较为复杂，且对如何直接提取参考信号和误差信号有效频率柜信息，并用一个频率柜相位和幅度信息近似周围频率点相应信息并未提及。据此，本发明的方法是在gao等算法上的一种改进，采用离散傅里叶变换(dft)直接提取参考信号和误差信号频率范围内有效频率柜信息，并均匀划分为若干子频率柜组，依次搜索子频率柜组更新方向，然后直接更新全带自适应滤波器系数。由于用一个频率点的幅度和相位信息近似周围频率点相应信息，不同于多相fft获取子带信号方式，从而简化子带结构，设计简单，且降低计算量。三.技术实现要素：1.发明目的：提出一种简化窄带无次级路径建模有源控制方法。该方法无需次级路径建模，采用离散傅里叶变换(dft)直接提取时延参考信号和误差信号有效频率柜信息，并用一个频率点的幅度和相位信息近似周围频率点相应信息，根据各频率柜组参考信号、误差信号以及相应更新方向一起直接更新全带自适应滤波器系数，系统结构简单，易于实现。2.技术方案本发明的目的通过以下技术方案实现一种简化窄带无次级路径建模有源控制方法，其步骤为：(1)计算实验中所用电路板电时延τ1＝q/fs，其中fs为采样频率，q为参考信号从控制器参考信号输入端到控制信号输出端延迟的采样点数。(2)计算实验中声时延τ2，计算公式为τ2＝l/c，其中，l为次级声源与误差传感器之间的距离，c为声速，一般环境中取340m/s。(3)基于步骤(1)电时延τ1和步骤(2)声时延τ2，计算总的时延τ＝τ1 τ2，总时延对应直达声时延点数为δ，计算公式为δ＝τfs，其中fs为采样频率。(4)产生时延修正参考信号时延修正参考信号r(n)为参考信号x(n)经由纯时延滤波器z-δ滤波得到，计算公式为r(n)＝x(n-δ)，其中δ为估计次级路径直达声时延采样点数。(5)提取频域时延修正参考信号和误差信号有效频率柜信息频域时延参考信号r(k)第i个频率柜r(i)可通过dft展开，采用同样的方法可得到第i个误差信号频率柜e(i)，计算公式分别为：式中，j为虚数单位。(6)子频率柜组更新方向的选择步骤1.不更新自适应滤波器系数，分别计算第t个频率柜组时延参考信号噪声能量χ1、误差信号能量ξ1以及最大误差幅度emax，计算公式：步骤2.分别从μ，-μ和±μi四个方向中搜索每个频率柜组最佳更新方向，重复下面步骤，寻找合适更新方向。a.采用式(1)更新自适应滤波器系数，计算p个采样点的频率柜组参考信号能量χ2和误差信号能量ξ2，如果|em(k)|＞(1 δ)emax，停止更新。b.如果|em(k)|＜(1 δ)emax，采用式(1)更新自适应滤波器系数。式中，μ(t)为步长，且为1，-1和±j，分别对应0°，180°和±90°四种不同更新方向。步骤3.如果满足|em(k)|＞(1 δ)emax或ξ2/χ2＞δ1ξ1/χ1、其中δ＜1，δ1＜1，则跳转至步骤5；步骤4.若μm＝μ，-μ和±μi，都满足|em(k)|＞(1 δ)emax或ξ2/χ2＞δ1ξ1/χ1，则μm＝μmα，α＜1，并返回步骤(1)；否则，μm＝μm/α，并返回步骤(1)；步骤5.重复执行步骤1至步骤4，直至m＝m-1；(7)所有频率柜组搜索最佳更新方向后，一起调整全带自适应滤波器系数，直至误差降低达到稳态，更新全带自适应滤波器w(z)系数采用如下迭代公式：式中，t为频率柜组个数，μ(t)为步长，且为1，-1和±j，分别对应0°，180°和±90°四种不同更新方向。每个子频率柜组初始更新方向为0°方向。3.有益效果本发明与现有技术相比，其显著优点：(1)通过dft方法提取出参考信号、误差信号和相应次级路径有效频率柜信息。(2)通过系统总时延，粗略估计次级路径直达声时延，并将直达声时延作为次级路径传递函数的一部分先进行补偿，然后将参考信号和误差信号的有效频率柜均匀划分为若干频率柜组，以减少无次级路径有源控制方法应用中子频率柜组个数。(3)无需次级路径建模，每个子频率柜组在0°、180°和±90°四个方向切换，搜索相应子频率柜组最佳更新方向，调整全带自适应滤波器系数。(4)针对高采样率频域窄带有源噪声控制系统，提取频率范围内参考信号和误差信号有效频率柜信息，避免全频带频域运算，系统结构简单，易于实现，减少计算量。四.具体实施方式下面以两端端口均封闭的管道声场中实现简化窄带无次级路径建模有源控制方法为例，对该方法的具体实施方式做详细说明。1.系统构成实验装置是端口直径为17cm的方形直管道，管道两端端口封闭，误差传声器位于次级声源的下游管道中，初级声源置于次级声源的上游管道中，测得初级声源与误差传声器之间的距离为136cm，次级声源与误差传声器之间的距离l为40cm。假设系统采样频率为16000hz。图1为简化窄带无次级路径建模有源控制方法框图，x(n)为210hz～250hz的带限信号，误差传感器用于获取误差信号e(n)。2.实施过程(1)将控制器的输入端口和输出端口短接，以白噪声作为控制器的输入信号，采用lms算法，测得控制器从输入端到输出端延迟采样点数。如ad21161评估板，降采样后，采样频率为2000hz，测得控制器从输入端到输出端共延迟13个采样点，可得控制器电时延τ1为6.5ms。(2)计算实验中声时延τ2，计算公式为τ2＝l/c，其中l为次级声源与误差传感器之间的距离，c为声速，一般环境中取340m/s。例如次级声源与误差传感器之间的距离l为40cm，计算声时延τ2近似为1.2ms。(3)基于步骤(1)电路板电时延τ1和步骤(2)声时延τ2，计算总的时延τ＝τ1 τ2，直达声时延δ＝τfs，fs为采样频率。例如控制器电时延τ1为6.5ms，τ2近似为1.2ms，计算总时延τ为7.7ms，fs为2000hz，对应的直达声时延δ为15。(4)产生时延修正参考信号时延修正参考信号r(n)为参考信号x(n)经由纯时延滤波器z-δ滤波得到，计算公式为r(n)＝x(n-δ)，其中δ为估计次级路径直达声时延采样点数。(5)提取参考信号和误差信号频率柜信息频域时延参考信号r(k)第i个频率柜r(i)可通过dft展开，采用同样的方法可得到第i个误差信号频率柜e(i)，其第i个频率柜对应的频率为fi，计算公式分别为：fi＝(i-1)×fs/2n式中，j为虚数单位。例如系统采样频率fs为16000hz，自适应滤波器长度n为512，参考信号频率范围为[210250]hz。当f＝210hz时，i＝14；当f＝250hz时，i＝17。有效频率柜范围为14至17，有效频率柜个数j＝4。(6)子频率柜组更新方向的选择p取值为2000，初始步长μ0为0.002，α取值为0.5，子频率柜组个数m为2，具体实施如下：步骤1.不更新自适应滤波器系数，获得最大误差信号幅度，迭代计算p个点的子频率柜误差信号功率和ξ1、子频率柜参考信号功率和χ1，计算公式为emax＝max(emax，|em(k)|)ξ1＝ξ1 |em(k)|2χ1＝χ1 |rm(k)|2步骤2.分别在μm＝μ，-μ和±μi两种更新方向情况下，重复下面步骤，寻找最佳更新方向。(a)迭代计算p个点的子频率柜误差信号功率和ξ2、子频率柜参考信号功率和χ2。(b)如果|em(k)|＜(1 δ)emax，采用下式更新自适应滤波器系数；否则停止更新。式中，r(i)和e(i)分别为时延参考信号r(k)和误差信号e(k)的第i个频率柜，μ(t)为步长，分别对应0°，180°和±90°四种不同更新方向。步骤3.如果满足|em(k)|＞(1 δ)emax或ξ2/χ2＞δ1ξ1/χ1、其中δ＜1，δ1＜1，则跳转至步骤5；步骤4.若μm＝μ0和-μ0，都满足|em(k)|＞(1 δ)emax或ξ2/χ2＞δ1ξ1/χ1，则μm＝μmα，α＜1，并返回步骤(1)；否则，μm＝μm/α，并返回步骤(1)；步骤5.比较步骤2中四种情况下所得子带误差信号功率ξ2，ξ2最小值对应合适的子带参考信号类型以及更新方向；步骤6.重复执行步骤1至步骤5，直至m＝m-1。(7)所有子频率柜组一起调整全带自适应滤波器系数，直至误差降低达到稳态，更新全带自适应滤波器w(z)系数采用如下迭代公式：式中，t为频率柜组个数，μ(t)为步长，且为1，-1和±j，分别对应0°，180°和±90°四种不同更新方向。3.实验结果本发明方法(简称ntanc方法)与fdfxlms方法的时域收敛曲线，如图2所示，达到稳态时，所提ntanc方法与fdfxlms方法最小均方误差(mse)一致，然而ntanc算法需要依次搜索每个子频率柜组的最佳更新方向，然后进入控制过程，因而占用一定时间，故与fdfxlms算法相比较，收敛速度稍慢。对两种方法降噪量进行比较，降噪前后频谱如图3所示，达到稳态时，本发明方法与fdfxlms方法降噪量相当。4.计算量分析窄带噪声有效频率柜范围为第r频率柜至第s频率柜，r＜s，且s-r 1＝j。当0≤i≤r-1和s 1≤i≤2n-1时，r(i)和e(i)均为0，又μ(k)r*(i)与e(i)乘积为复数，因此计算一个自适应滤波器系数需要复数乘为4j，共有n个自适应滤波器系数，故自适应滤波器系数更新式需要计算量为4jn复数乘，即8jn实数乘。比较fdfxlms算法和ntanc算法计算量如表1。表1比较fdfxlms算法和ntanc算法计算量步骤fdfxlmsntanc算法次级路径建模n0频域滤波参考信号2n×4 2n×log2(2n)0全频带频域信号4n×log2(2n)0参考信号有效频率柜02jn误差信号有效频率柜02jn控制器系数更新2nlog2(2n) 2n×48jn总计17n 8nlog2(2n)12jn假设自适应滤波器长度n＝512，则fdfxlms计算量为45568实数乘，ntanc算法计算量为6144j实数乘，j为窄带信号有效频率柜个数。当前设置下，当有效频率柜个数j≈8.1，ntanc算法与fdfxlms算法计算量相等；当有效频率柜个数j≤8，ntanc算法计算量比fdfxlms算法小；当频率柜个数j＞8时，ntanc算法计算量比fdfxlms算法计算量大。五.附图说明图1算法结构框图图2mse曲线比较图3仿真降噪前后频谱。当前第1页12