专利名称：一种用于消减噪声的多通道有源噪声控制方法
技术领域：
本发明涉及一种消减噪声的方法，尤其是涉及一种用于消减噪声的多通道有源噪声控制方法。
背景技术：
现有的有源噪声控制方法是利用声波的叠加原理，通过次级声源产生一组与初级噪声幅度相同相位相反的声波，二者在误差传感器相加，从而达到消减噪声的效果。几十年来，人们提出许多噪声控制算法，其中最著名的是Widrow提出的FXLMS算法。
鉴于FXLMS算法收敛速度慢，对宽带噪声降噪效果低的情况，人们根据噪声频谱宽带特性，提出了基于子带自适应滤波器以及频域FXLMS类的算法，但它们在初级参考噪声路径上引入了较多的时延，因而在宽带噪声控制上效果并不理想。为克服子带滤波器对初级噪声时延的影响，Morgan(1995)提出了一种无延时子带滤波器结构，该结构利用FFT，自适应滤波器系数在各子带内分别调整，通过特定组合，再转换到时域滤波器，然后产生控制信号。因其消除了对初级参考信号的时延，所以提高了算法的收敛速度。
然而，从本质上看，有源噪声控制是建立在声波的叠加原理之上的。上述单次级源单误差传感器的单通道算法，无论其性能如何，只能在误差传感器附近有效。观测点与误差传感器距离远到一定程度，噪声控制效果明显变差。这样，对需要一定空间控制范围的应用，单通道结构不再适用。因此，人们又提出了以多通道FXLMS算法为代表的多通道噪声控制算法。但研究发现，多通道FXLMS算法的工作过程不像单通道算法那样简单，各通道之间存在强烈的相互耦合效应，使得多通道FXLMS算法收敛速度大为降低，控制效果也不尽人意。尽管已经提出了不少解耦算法，但它们大多结构复杂，运算量过大，不利于实际应用。
不考虑附属电路，典型的前馈结构多通道ANC系统的数学模型如图1所示. 图中，共有I个参考传感器J个次级声源和K个误差传感器，即典型的I×J×K结构。其中，P(z)是从参考传感器到误差传感器之间的传输函数，称之为初级通道模型.S(Z)是从滤波器W(Z)输出到误差传感器之间的传输函数。W(Z)是控制器的传输函数，一般由自适应滤波器实现。
是对S(z)的估计，需事先确定。
定义 xi(n)来自参考传感器i的第i个参考信号 yj(n)第j个次级源的输出信号 dk(n)第k个误差传感器的期望信号 ek(n)第k个误差传感器的误差信号 Wij(n)第i个参考信号到第j个次级源输出的自适应滤波器系数向量 Sjk(n)第j个次级源到第k个误差传感器之间的传递函数
第j个次级源到第k个误差传感器之间的传递函数的估计，需要预先辨识 vijk(n)xi(n)通过次级通道模型
后的滤波参考信号 yjk(n)yj(n)经次级通道Sjk(n)到达第k个误差传感器后的信号设自适应滤波器阶数为L，次级通道模型阶数为M，则 Xi(n)＝[xi(n)，xi(n-1)，xi(n-2)...，xi(n-L+1)]T Yj(n)＝[yj(n)，yj(n-1)，yj(n-2)...，yj((n-M+1)]T Vijk(n)＝[vijk(n)，vijk(n-1)，vijk(n-2)，...vijk(n-L+1)]T 根据以上定义，可得多通道FLMS算法的步骤如下 根据上述步骤，可以得到传统多通道FLMS算法的计算量为L(IJ+IJK)+M(IJK)+K次乘法运算。一般而言L＞＞1，M＞＞1，所以整个系统的运算量可近似为IJK(2L+M)，即系统运算量与参考传感器个数I、次级源个数J和误差传感器个数K的乘积近似成正比关系。随着I，J，K的增加，运算量将指数上升，所以系统规模很难做大。如果I＝J＝K＝1，则可以得到单通道FLMS算法的运算量为2L+M+1。通过比较，可以发现，多通道FLMS算法运算量近似为单通道FLMS的IJK倍。
发明内容
本发明所要解决的技术问题是提供一种用于消减噪声的多通道有源噪声控制方法，在基本不降低系统性能的前提下，使得系统的运算量成倍降低，更加适宜实际应用。
本发明解决上述技术问题所采用的技术方案为一种用于消减噪声的多通道有源噪声控制方法，包括设置1个用于拾取初级噪声源信号的参考传感器、至少2个用以产生与初级噪声源信号相反的抵消噪声信号的次级源和与所述的次级源数量相同的用以拾取初级噪声源信号与抵消噪声信号相互抵消后的残留信号的误差传感器，将所述的次级源与所述的误差传感器成对布放形成一路信号传输通道，然后按以下步骤进行 1)通过控制电路产生辨识激励信号z(n)，将辨识激励信号z(n)依次传送到各路信号传输通道，使次级源发出激励声波，由对应的误差传感器拾取激励声波，并送回所述的控制电路，用普通的LMS算法对所述的次级源与所述的误差传感器之间的传递函数Sjj(n)进行辨识，依次得到各路信号传输通道的传递函数估计值
其中j表示信号传输通道的路序号，n是时间，为大于零的正整数； 2)控制电路读取参考传感器的信号得到初级噪声源信号x(n)，根据公式得到第j路信号传输通道的合成参考信号x′j(n)，其中J为总的信号传输通道路数，y′l(n)＝αlyl(n)，yl(n)是第l路信号传输通道自适应滤波器输出信号，l为大于等于1小于等于J且不等于j的正整数，0＜αl＜0.5为第l路输出yl(n)的耦合系数； 3)根据公式得到第j路信号传输通道的自适应滤波器输出信号 yj(n)，式中，WTj(n)是第j路自适应滤波器系数向量Wj(n)的转置，而Wj(n)表示为公式Wj(n)＝[wj(n)，wj(n-1)，wj(n-2)...，wj(n-L+1)]T；X′j(n)是第j路合成参考信号向量，表示为公式X′j(n)＝[x′j(n)，x′j(n-1)，x′j(n-2)...，x′j(n-L+1)]T；式中，上标T是转置 符号，L是自适应滤波器的阶数； 4)根据公式得到j路信号传输通道的滤波参考信号vjj(n)； 5)根据公式Wj(n+1)＝Wj(n)-μjVjj(n)ej(n)，对j路信号传输通道的自适应滤波器系数向量进行调整，μj是自适应滤波器收敛因子，表达为公式其中Vjj(n)＝[vjj(n)，vjj(n-1)，vjj(n-2)，...vjj(n-L+1)]T是第j路信号传输通道的滤波参考信号向量，δj的取值范围为0＜δj＜＜1，用以防止分母为零，
取的正常数；上标T是转置符号；ej(n)是第j路信号传输通道中的误差传感器得到的误差信号；重复步骤2)～5)，ej(n)逐渐减小，从而实现在第j个误差传感器处降低噪声的目的。
步骤1)中获得传递函数估计值
的具体过程为控制电路产生辨识激励信号z(n)，该辨识激励信号z(n)一路经次级源发出激励声波，该激励声波经所述的次级源与所述的误差传感器之间的传递函数Sjj(n)后，在误差传感器处得到期望信号ρ(n)；另一路经LMS算法所用的自适应滤波器U(n)，得到自适应滤波器输出信号uo(n)，自适应滤波器系数向量为U(n)＝[u(n)，u(n-1)，u(n-2)...，u(n-M+1)]T，其中M为自适应滤波器的阶数，根据公式 uo(n)＝UT(n)Z(n) e(n)＝ρ(n)-uo(n) U(n+1)＝U(n)+ηZ(n)e(n) 对自适应滤波器的系数向量U(n)进行调整，循环往复，直到U(n)收敛，此时U(n)就是所述的次级源与所述的误差传感器之间传递函数的估计
其中，Z(n)＝[z(n)，z(n-1)，z(n-2)...，z(n-M+1)]T是辨识激励信号向量，β的取值范围为0＜β＜＜1，用以防止分母为零；
取的正常数；e(n)是期望信号ρ(n)与自适应滤波器输出信号uo(n)之间的误差信号。
所述的次级源是扬声器。
与现有技术相比，本发明的优点在于通过合理布放次级声源和误差传感器，在基本不降低系统性能的前提下，将传统的多通道结构简化为多路单通道结构，从而使得系统的运算量成倍降低，更加适宜实际应用；与传统结构不同的是，第j个自适应滤波器的输入信号，除了直接来自参考传感器输入x(n)外，还外加其它次级源的输出信号，这样，对第j个自适应滤波器而言，其它次级源的输出信号yi(n)，1≤l≤J，l≠j等效为J-1个初级参考信号；对第j个误差传感器而言，来自其它次级源的信号yl(n)，1≤l≤J，l≠j，经次级通道传递函数Slj(n)，l≠j到达误差传感器j后的信号yl(n)*Slj(n)不作为干扰信号看待，而是作为期望信号看待，从而形成了具有“J”个参考信号，一个次级源和一个误差传感器的ANC结构。因这“J”个参考信号并不是独立的，它们均来自于初级参考信号x(n)，所以将它们简单的相加，得到一个合成参考信号x′(n)。经过这样安排，第j个自适应滤波器的更新步骤与单通道FLMS算法结构非常相近，不同之处只在于输入信号由x(n)变为x′j(n)，其它通道类似。根据前面内容，单通道FLMS算法的运算量为2L+M+1次乘法运算，因共有J个通道，所以新结构的总运算量为J(2L+M+1)次乘法运算，相比传统的多通道FLMS算法JK(2L+M)次乘法，新结构运算量降低了K-1倍。
图1为现有技术典型的前馈结构多通道ANC系统的数学模型； 图2为本发明装置的结构示意图； 图3为本发明新型的多通道ANC系统的数学模型； 图4为本发明实施例的装置的结构示意图。
具体实施例方式 以下结合附图实施例对本发明作进一步详细描述。
实施例一种用于消减噪声的多通道有源噪声控制方法，包括设置1个用于拾取初级噪声源信号的参考传感器R、2个用以产生与初级噪声源信号相反的抵消噪声信号的扬声器S1和S2和2个用以拾取初级噪声源信号与抵消噪声信号相互抵消后的残留信号的误差传感器M1和M2，将扬声器S1与误差传感器M1成对布放形成第一路信号传输通道，将扬声器S2与误差传感器M2成对布放形成第二路信号传输通道，用合众达TMS320c6727高级音频开发板组成一个控制电路1，然后按以下步骤进行 1)通过控制电路1产生辨识激励信号z(n)，将辨识激励信号z(n)经扬声器S1发出激励声波，由误差传感器M1拾取该激励声波，并送控制电路1，用普通的LMS算法对扬声器S1与误差传感器M1之间的传递函数S11(n)进行辨识，得到扬声器S1与误差传感器M1的传递函数估计值
即第一路信号传输通道传递函数估计值；然后，控制电路1再次产生辨识激励信号，经次级源S2发出激励声波，由误差传感器M2拾取该激励声波，并送控制电路1，用普通的LMS算法对扬声器S2与误差传感器M2之间的传递函数S22(n)进行辨识，得到扬声器S2与误差传感器M2的传递函数估计值
即第二路信号传输通道传递函数估计值；n是时间，为大于零的正整数； 2)控制电路1读取参考传感器R的信号得到初级噪声源信号x(n)，根据公式x＇1(n)＝x(n)+α2y2(n)得到第一路信号传输通道的合成参考信号x′1(n)，y2(n)是第二路信号传输通道自适应滤波器输出信号，α2为第二路输出信号传输通道自适应滤波器输出信号y2(n)的耦合系数，在这里取0.1；根据公式x′2(n)＝x(n)+α1y1(n)得到第二路信号传输通道的合成参考信号x′2(n)，y1(n)是第一路信号传输通道自适应滤波器输出信号，α1为第一路输出信号传输通道自适应滤波器输出信号y1(n)的耦合系数，在这里取0.1； 3)根据公式得到第一路信号传输通道的自适应滤波器输出信号y1(n)，式中，W1T(n)是第一路自适应滤波器系数向量W1(n)的转置，而W1(n)表示为公式W1(n)＝[w1(n)，w1(n-1)，w1(n-2)...，w1(n-L+1)]T；X′1(n)是第一路合成参考信号向量，表示为公式X′1(n)＝[x′1(n)，x′1(n-1)，x′1(n-2)...，x′1(n-L+1)]T；根据公式得到第二路信号传输通道的自适应滤波器输出信号y2(n)，式中，WT2(n)是第二路自适应滤波器系数向量W2(n)的转置，而W2(n)表示为公式W2(n)＝[w2(n)，w2(n-1)，w2(n-2)...，w2(n-L+1)]T；X′2(n)是第二路合成参考信号向量，表示为公式X′2(n)＝[x′2(n)，x′2(n-1)，x′2(n-2)...，x′2(n-L+1)]T；上述公式中，上标T是转置符号，L是自适应滤波器的阶数； 4)根据公式得到第一路信号传输通道的滤波参考信号v11(n)；其中
是第一路信号传输通道传递函数估计值
的转置，X′1(n)＝[x1(n)，x′1(n-1)，x′1(n-2)...，x′1(n-L+1)]T是第一路信号传输通道滤波参考信号向量；根据公式得到第二路信号传输通道的滤波参考信号v22(n)；其中
是第二路信号传输通道传递函数估计值
的转置，X′2(n)＝[x′2(n)，x′2(n-1)，x′2(n-2)x′1，x′2(n-L+1)]T是第二路信号传输通道滤波参考信号向量； 5)根据公式W1(n+1)＝W1(n)-μ1V11(n)e1(n)，对第一路信号传输通道的自适应滤波器系数向量进行调整，μ1是自适应滤波器收敛因子，表达为公式其中V11(n)＝[v11(n)，v11(n-1)，v11(n-2)，...v11(n-L+1)]T是第一路信号传输通道的滤波参考信号向量，δ1的取值为0.000001，
取1；上标T是转置符号；e1(n)是第一路信号传输通道中的误差传感器M1得到的误差信号；重复步骤2)～5)，e1(n)逐渐减小，从而实现在第一路信号传输通道中的误差传感器M1处降低噪声的目的；根据公式W2(n+1)＝W2(n)-μ2V22(n)e2(n)，对第二路信号传输通道的自适应滤波器系数向量进行调整，μ2是自适应滤波器收敛因子，表达为公式其中V22(n)＝[v22(n)，v22(n-1)，v22(n-2)，...v22(n-L+1)]T是第二路信号传输通道的滤波参考信号向量，δ2的取值为0.000001，
取1；上标T是转置符号；e2(n)是第二路信号传输通道中的误差传感器M2得到的误差信号；重复步骤2)～5)，e2(n)逐渐减小，从而实现在第二路信号传输通道中的误差传感器M2处降低噪声的目的。
步骤1)中获得传递函数估计值
的具体过程为控制电路1产生辨识激励信号z(n)，该辨识激励信号z(n)经扬声器S1发出激励声波，由误差传感器M1拾取该激励声波，该激励信号一路经扬声器S1和误差传感器M1之间的传递函数S11(n)后，在误差传感器处得到期望信号ρ(n)；另一路经LMS算法所用的自适应滤波器U(n)，得到自适应滤波器输出信号uo(n)，自适应滤波器系数向量为 U(n)＝[u(n)，u(n-1)，u(n-2)...，u(n-M+1)]T，其中M为自适应滤波器的阶数，根据公式 uo(n)＝UT(n)Z(n) e(n)＝ρ(n)-uo(n) U(n+1)＝U(n)+ηZ(n)e(n) 对自适应滤波器的系数向量U(n)进行调整，循环往复，直到U(n)收敛，此时U(n)就是扬声器S1和误差传感器M1之间传递函数的估计
其中，Z(n)＝[z(n)，z(n-1)，z(n-2)...，z(n-M+1)]T是辨识激励信号向量，β的取值为0.000001，
取1；e(n)是期望信号ρ(n)与自适应滤波器输出信号uo(n)之间的误差信号。获得传递函数估计值
的具体过程与获得
的过程相同。
在上述实施例的基础上，设置一个初级扬声器2，扬声器2、扬声器S1和扬声器S2均为惠威1010MKII多媒体音箱；参考传感器R、误差传感器M1和误差传感器M2型号均为宁波Altron公司PM02；初级扬声器1距扬声器S1和扬声器S2的距离均为1.5m，扬声器S1距误差传感器M1的距离和扬声器S2距误差传感器M2的距离均为0.4m，扬声器S1和扬声器S2之间距离0.6m。数据采集和算法的运行在合众达TMS320c6727高级音频开发板上进行，系统的采样频率设为2Khz。
采用本发明提出的方法，对初级扬声器2发出152Hz的单频信号进行实验，各滤波器阶数均取128，耦合系数α1、α2均取0.1，δ1的取值为0.000001，
取1，δ2的取值为0.000001，
取1。分别对比算法收敛前后，误差传感器M1，M2处的噪声声级，发现算法收敛后，两处的噪声声级均下降分别为18、16dB左右。与在同样的实验条件下，采用传统的多通道FLMS算法得到的结果相比基本一致。说明本发明提出的方法在大幅降低计算量的情况下，可保持系统性能基本不变，证明了本发明的有效性。
权利要求
1.一种用于消减噪声的多通道有源噪声控制方法，其特征在于包括设置1个用于拾取初级噪声源信号的参考传感器、至少2个用以产生与初级噪声源信号相反的抵消噪声信号的次级源和与所述的次级源数量相同的用以拾取初级噪声源信号与抵消噪声信号相互抵消后的残留信号的误差传感器，将所述的次级源与所述的误差传感器成对布放形成一路信号传输通道，然后按以下步骤进行
1)通过控制电路产生辨识激励信号z(n)，将辨识激励信号z(n)依次传送到各路信号传输通道，使次级源发出激励声波，由对应的误差传感器拾取激励声波，并送回所述的控制电路，用普通的LMS算法对所述的次级源与所述的误差传感器之间的传递函数Sjj(n)进行辨识，依次得到各路信号传输通道的传递函数估计值
其中j表示信号传输通道的路序号，n是时间，为大于零的正整数；
2)控制电路读取参考传感器的信号得到初级噪声源信号x(n)，根据公式
得到第j路信号传输通道的合成参考信号x′j(n)，其中J为总的信号传输通道路数，y′l(n)＝αlyl(n)，yl(n)是第l路信号传输通道自适应滤波器输出信号，l为大于等于1小于等于J且不等于j的正整数，0＜αl＜0.5为第l路输出yl(n)的耦合系数；
3)根据公式
得到第j路信号传输通道的自适应滤波器输出信号yj(n)，式中，WjT(n)是第j路自适应滤波器系数向量Wj(n)的转置，而Wj(n)表示为公式Wj(n)＝[wj(n)，wj(n-1)，wj(n-2)...，wj(n-L+1)]T；X′j(n)是第j路合成参考信号向量，表示为公式X′j(n)＝[x′j(n)，x′j(n-1)，x′j(n-2)...，x′j(n-L+1)]T；式中，上标T是转置符号，L是自适应滤波器的阶数；
4)根据公式
得到j路信号传输通道的滤波参考信号vjj(n)，其中
是信号传输通道的传递函数估计值
X′j(n)是；
5)根据公式Wj(n+1)＝Wj(n)-μjVjj(n)ej(n)，对j路信号传输通道的自适应滤波器系数向量Wj(n)＝[wj(n)，wj(n-1)，wj(n-2)...，wj(n-L+1)]T进行调整，μj是自适应滤波器收敛因子，表达为公式
其中Vjj(n)＝[vjj(n)，vjj(n-1)，vjj(n-2)，...vjj(n-L+1)]T是第j路信号传输通道的滤波参考信号向量，δj的取值范围为0＜δj＜＜1，用以防止分母为零，
取
的正常数；上标T是转置符号；ej(n)是第j路信号传输通道中的误差传感器得到的误差信号；重复步骤2)～5)，ej(n)逐渐减小，从而实现在第j个误差传感器处降低噪声的目的。
2.如权利要求1所述的一种用于消减噪声的多通道有源噪声控制方法，其特征在于步骤1)中获得传递函数估计值
的具体过程为控制电路产生辨识激励信号z(n)，该辨识激励信号z(n)一路经次级源发出激励声波，该激励声波经所述的次级源与所述的误差传感器之间的传递函数Sjj(n)后，在误差传感器处得到期望信号ρ(n)；另一路经LMS算法所用的自适应滤波器U(n)，得到自适应滤波器输出信号uo(n)，自适应滤波器系数向量为U(n)＝[u(n)，u(n-1)，u(n-2)...，u(n-M+1)]T，其中M为自适应滤波器的阶数，根据公式
uo(n)＝UT(n)Z(n)
e(n)＝ρ(n)-uo(n)
U(n+1)＝U(n)+ηZ(n)e(n)
对自适应滤波器的系数向量U(n)进行调整，循环往复，直到U(n)收敛，此时U(n)就是所述的次级源与所述的误差传感器之间传递函数的估计
其中
，Z(n)＝[z(n)，z(n-1)，z(n-2)...，z(n-M+1)]T是辨识激励信号向量，β的取值范围为0＜β＜＜1，用以防止分母为零；
取
的正常数；e(n)是期望信号ρ(n)与自适应滤波器输出信号uo(n)之间的误差信号。
3.如权利要求1或2所述的一种用于消减噪声的多通道有源噪声控制方法，其特征在于所述的次级源是扬声器。
全文摘要
本发明公开了一种用于消减噪声的多通道有源噪声控制方法，通过设置1个用于拾取初级噪声源信号的参考传感器、至少2个用以产生与初级噪声源信号相反的抵消噪声信号的次级源和与次级源数量相同的用以拾取残留信号的误差传感器，并将次级源与误差传感器成对布放形成一路信号传输通道，与传统结构不同的是，第j个自适应滤波器的输入信号，除了直接来自参考传感器输入x(n)外，还外加其它次级源的输出信号，优点在于通过合理布放次级声源和误差传感器，在基本不降低系统性能的前提下，将传统的多通道结构简化为多路单通道结构，从而使得系统的运算量成倍降低，更加适宜实际应用；相比传统的多通道FLMS算法JK(2L+M)次乘法，新结构运算量降低了K-1倍。
文档编号G10K11/178GK101819766SQ20101004006
公开日2010年9月1日 申请日期2010年1月15日 优先权日2010年1月15日
发明者张瑞华, 谢智波 申请人:浙江万里学院