CN103916733A - 基于均方误差最小准则的声能量对比控制方法及系统 - Google Patents

Abstract

本发明涉及一种基于均方误差最小准则的声能量对比控制方法，包括：布放扬声器阵列，设定明区和暗区的控制点；在一系列频率点上分别计算扬声器阵列的输入复权系数向量；设定各个频率点上的目标向量，按照均方误差最小准则，分别调整各个频率点上的输入复权系数向量的幅度与相位；分别将每个扬声器通道在一系列频率点上的复权系数组成频率响应，对所述频率响应加窗并进行逆傅里叶变换，得到各个扬声器通道的时域冲激响应信号。

Description

基于均方误差最小准则的声能量对比控制方法及系统

技术领域

本发明涉及声学领域，特别涉及一种基于均方误差最小准则的声能量对比控制方法及系统。

背景技术

声能量对比控制方法被广泛应用于电视机、电脑和手机等个人音频系统领域。声能量对比控制方法能够在不打扰其它用户(暗区)的情况下将大部分声能量聚集在用户周围(明区)，从而为用户打造个性化聆听空间，提高生活品质。

现有技术中的声能量对比控制方法的具体实现步骤如下：

(1)确定扬声器阵列，明区和暗区控制点的位置；

(2)确定频率点f,采用B&K PULSE等音频测试仪器或者通过建模仿真获得扬声器阵列分别到明区和暗区的控制点在频率f处的传递函数矩阵G_B,f和G_D,f，其大小分别为M_B×L和M_D×L，其中M_B和M_D分别是明区和暗区的控制点个数，L为扬声器个数；

(3)计算出扬声器阵列在频率f处的大小为L×1的输入复权系数向量w_M,f，其表达式为：

$w_{M,f}=\arg \underset{w}{\max }\frac{w^H{G}_{B,f}^H{G}_{B,f}w}{w^H(G_{D,f}^H{G}_{D,f}+{\mathrm{\δ}}_{f}I)w},$

$=P_{\max }\left\{{(G_{D,f}^H{G}_{D,f}+{\mathrm{\δ}}_{f}I)}^{-1}{G}_{B,f}^H{G}_{B,f}\right\}$

其中P_max{·}是指取矩阵对应最大特征值的单位特征向量，I为单位矩阵，δ_f为鲁棒性参数。

(4)分别将每个扬声器通道上的各个频率点的复权系数组成频率响应，并进行逆傅里叶变换，得到各个通道的时域冲激响应信号。

从对现有技术的声能量对比控制方法的上述描述可以看出，在现有方法中，由于时域冲激响应信号是通过直接对频域信号进行逆傅里叶变换得到的，因此时域冲激响应信号的因果性无法得到保证。这使得现有的声能量对比控制方法在非选定频率点上的对比聚焦效果下降，并且在不同的频率点上有不同的响应，从而无法在明区得到较好的音质。在图1和图2中，给出了全消声室内8元阵列的常规声能量对比控制方法在明区中心处控制点上的冲激响应信号和的对比聚焦的结果，其系统采样频率f_s为8kHz，冲激响应信号长度为100。从图1中可以看出，冲激响应信号有较大程度的边缘拖尾现象，这会导致听到的声音有明显的混响效果，从而降低在明区的音质；从图2中可以看出，在选定的频率点上，现有技术中的声能量对比控制方法具有较好的对比聚焦效果，但在非选定频率上，对比聚焦效果下降很快，最差的点上可下降40dB，几乎不存在对比聚焦效果。现有技术中的声能量对比控制方法所存在的这些问题造成了在个人音频系统中并不适合处理宽带信号输入，降低了用户区域的音质，在暗区中泄露更多的声能量，无法为用户营造个性化的私人聆听空间。

针对现有技术中的声能量对比控制方法中存在的音质和宽带对比聚焦能力下降的问题，需要寻找更为简便有效的方法来改善其音质和宽带对比聚焦能力。

发明内容

本发明的目的在于克服现有技术中的声能量对比控制方法所存在的音质和宽带对比聚焦能力下降的问题，从而提供一种能够有效改善音质和宽带对比聚焦能力的方法。

为了实现上述目的，本发明提供了一种基于均方误差最小准则的声能量对比控制方法，包括：

步骤1)、布放扬声器阵列，设定明区和暗区的控制点；其中，所述明区为具有高声能量的区域，所述暗区为具有低声能量的区域；

步骤2)、在一系列频率点上分别计算扬声器阵列的输入复权系数向量；

步骤3)、在所述一系列频率点上设定各个频率点的目标向量，按照均方误差最小准则，分别调整步骤2)所得到的各个频率点上的输入复权系数向量的幅度与相位；

步骤4)、分别将每个扬声器通道在步骤3)得到的所述一系列频率点上的复权系数组成频率响应，对所述频率响应加窗并进行逆傅里叶变换，得到各个扬声器通道的时域冲激响应信号。

上述技术方案中，所述的步骤3)包括：

步骤3-1)、从一系列频率中选择频率点f,设定其目标向量d_f大小为L×1,L为扬声器个数；

步骤3-2)、在步骤2)得到的输入复权系数向量w_M,f的基础上，求取在频率点f处基于均方误差最小准则的输入最优复权系数向量w_O,f；

步骤3-3)、从选定的一系列频率点中选择其它未经处理的频率点，重复步骤3-1)-步骤3-2)，直至计算出所有频率点的输入复权系数向量w_O,f。

上述技术方案中，在所述的步骤3-2)中，求取在频率点f处基于均方误差最小准则的输入最优复权系数向量w_O,f包括：

假定在频率点f处基于均方误差最小准则的输入最优复权系数向量w_O,f的表达式为：

$w_{O,f}={\mathrm{\β}}_f{e}^{{\mathrm{j\θ}}_f}{w}_{M,f}$

其中β_f为任意大于0的缩放系数，θ_f为任意旋转角度；

使得所述β_f和θ_f满足均方误差最小准则，即

$\underset{{\mathrm{\β}}_f,{\mathrm{\θ}}_f}{\min }{\left|\right|G_{B,f}{w}_{O,f}-d_f\left|\right|}_2^2={\mathrm{\β}}_f^2{w}_{M,f}^H{G}_{B,f}^H{G}_{B,f}{w}_{M,f}-2{\mathrm{\β}}_f\mathrm{Re}\left\{e^{-j{\mathrm{\θ}}_f}{w}_{M,f}^H{G}_{B,f}^H{d}_f\right\}+d_f^H{d}_f$

其中Re{·}取复数的实数部分，G_B,f为扬声阵列到明区控制点之间的传递函数矩阵，

对β_f和θ_f分别进行求偏导，并置0得：

${\mathrm{\θ}}_f=\∠\left(w_{M,f}^H{G}_{B,f}^H{d}_f\right)$

${\mathrm{\β}}_f=\frac{\mathrm{abs}\left(w_{M,f}^H{G}_{B,f}^H{d}_f\right)}{w_{M,f}^H\left(G_{B,f}^H{G}_{B,f}\right)w_{M,f}}$

其中∠(·)是取复数的相位部分，abs(·)是取复数的幅度部分。

所要求取的最优复权系数向量w_O,f的最终表达式为：

$w_{O,f}=\frac{\mathrm{abs}\left(w_{M,f}^H{G}_{B,f}^H{d}_f\right)}{w_{M,f}^H\left(G_{B,f}^H{G}_{B,f}\right)w_{M,f}}{e}^{j\∠\left(w_{M,f}^H{G}_{B,f}^H{d}_f\right)}{w}_{M,f}.$

上述技术方案中，在步骤1)中，所布放的扬声器阵列为线性阵列或圆形阵列，或随机阵列。

上述技术方案中，在步骤1)中，所述明区或暗区的区域形状是方形或圆形或线型。

上述技术方案中，在步骤4)中，对频率响应加窗时所要用到的窗函数为矩形窗或汉宁窗或者海明窗。

本发明还提供了一种基于均方误差最小准则的声能量对比控制系统，该系统包括扬声器阵列设定模块、输入复权系数向量计算模块、输入复权系数向量调整模块、时域冲激响应信号生成模块；其中，

所述的扬声器阵列设定模块用于布放扬声器阵列，并设定明区和暗区的控制点；其中，所述明区为具有高声能量的区域，所述暗区为具有低声能量的区域；

所述的输入复权系数向量计算模块用于在一系列频率点上分别计算扬声器阵列的输入复权系数向量；

所述的输入复权系数向量调整模块用于在所述一系列频率点上设定各个频率点的目标向量，按照均方误差最小准则，分别调整所述输入复权系数向量计算模块所得到的各个频率点上的输入复权系数向量的幅度与相位；

所述时域冲激响应信号生成模块分别将每个扬声器通道在一系列频率点上的复权系数组成频率响应，对所述频率响应加窗并进行逆傅里叶变换，得到各个通道的时域冲激响应信号。

上述技术方案中，所述的输入复权系数向量调整模块还包括目标向量设定单元、输入最优复权系数向量求取单元；其中，

所述的目标向量设定单元从一系列频率中选择频率点f,设定其目标向量d_f大小为L×1；L为扬声器个数；

所述的输入最优复权系数向量求取单元在所述的输入复权系数向量计算模块得到的输入复权系数向量w_M,f的基础上，求取在频率点f处基于均方误差最小准则的输入最优复权系数向量w_O，f。

本发明的优点在于：

A、本发明通过预设目标向量，采用均方最小误差准则，使得明区的频率响应逼近于理想的Dirac脉冲的频率响应，从而改善了明区的音质。

B、只通过简单的对频率加窗操作，可以保证时域冲激响应的因果性。

C、与常规声能量对比控制方法相比较，在相同配置情况下，本发明的宽带信号对比聚焦能力会远大于常规声能量对比控制方法的宽带信号对比聚焦能力。

D、相比于常规声能量对比控制方法，只增加了少量的运算量，而且计算简单，便于实时实现。

E、本发明不仅能在窄带信号上取得与常规声能量对比控制方法一样的对比聚焦能力，而且能在宽带信号上取得比常规声能量对比控制方法更好的对比聚焦能力。

F、本发明的声能量对比控制方法能够广泛应用于电视机，电脑和手机等个人音频领域，在宽带信号输入时，能够提高明区的音质，并取得较好的聚焦能力，从而营造出较好的私人聆听空间。在常规声能量对比控制过程中，由于并未考虑明区的频率响应，同时，由于只是对离散频率点进行控制，并且转换为时域冲激响应信号时可能存在非因果性，这使得常规声能量对比控制方法对宽带信号处理时无法取得较好的音质和对比聚焦效果。

附图说明

图1是现有技术中的声能量对比控制方法在明区中心控制点的冲激响应信号的示意图，其系统采样频率f_s为8kHz，冲激响应信号长度为100；

图2是现有技术中的声能量对比控制方法在选定频率点上和非选定频率点上的对比聚焦性能，其系统采样频率f_s为8kHz；

图3是本发明方法的流程图；

图4是本发明的实施例中实验配置示意图；

图5是本发明的实施例中扬声器第一通道到明区中心控制点的频率响应示意图；

图6是本发明的实施例中采用不同方法在明区中心控制点处的冲激响应，其中，(a)表示本发明方法，(b)表示现有技术中的声能量对比控制方法；

图7是本发明的实施例中分别采用本方法和现有技术中的声能量对比控制方法在宽带信号上的对比聚焦性能的示意图。

具体实施方式

现结合附图对本发明作进一步的描述。

本发明的基本构思是利用设计均方误差最小准则，预设目标向量，使得明区控制点的频率响应逼近于理想的Dirac-delta脉冲函数的频率响应，从而使本发明的声能量对比控制方法在明区能够改善音质，同时提高了在非选定频率点的对比聚焦能力。本发明消除了现有技术中的声能量控制方法存在的明区冲激响应波形拖尾和在非选定频率点上聚焦能力剧烈下降问题。

如图3所示，本发明的方法包括以下步骤：

步骤1)、布放扬声器阵列，设定明区和暗区的控制点；其中，所述明区为具有高声能量的区域，所述暗区为具有低声能量的区域；

步骤2)、在一系列频率点上分别计算扬声器阵列的输入复权系数向量；

步骤3)、在所述一系列频率点上设定各个频率点的目标向量，按照均方误差最小准则，分别调整各个频率点上的输入复权系数向量的幅度与相位；

步骤4)、分别将每个扬声器通道在一系列频率点上的复权系数组成频率响应，对频率响应加窗，并进行逆傅里叶变换，得到各个通道的时域冲激响应信号。

下面对本发明方法中的各个步骤做进一步说明。

在步骤1)中，所布放的扬声器阵列为线性阵列或圆形阵列，也可以为随机阵列。所述明区或暗区的区域形状可以是方形或圆形，也可以是线型。

在步骤2)中，在一系列频率点上计算扬声器阵列的输入复权系数向量可以采用现有技术中的方法实现，该步骤具体包括：

步骤2-1)、从一系列频率中选择频率点f,采用诸如B&K PULSE的音频测试仪器或者通过建模仿真获得扬声器阵列分别得到明区和暗区的控制点在频率f处的传递函数矩阵G_B,_f和G_D,f，其大小分别为M_B×L和M_D×L，其中M_B和M_D分别是明区和暗区的控制点个数，L为扬声器个数；

步骤2-2)、确定鲁棒性参数δ_f数值，鲁棒性参数δ_f选取由系统测量传递函数矩阵G_B,f与真实传递函数的误差ΔG_B,f决定，假定矩阵ΔG_B,f的F范数最大值不大于0.01，即||ΔG_B,f||_F≤0.01，则δ_f的数值选为0.01；

步骤2-3)、根据现有技术中的声能量对比控制方法，计算出扬声器阵列在频率f处的大小为L×1的输入复权系数向量w_M,f，其表达式为：

$w_{M,f}=\arg \underset{w}{\max }\frac{w^H{G}_{B,f}^H{G}_{B,f}w}{w^H(G_{D,f}^H{G}_{D,f}+{\mathrm{\δ}}_{f}I)w}---\left(1\right)$

$=P_{\max }\left\{{(G_{D,f}^H{G}_{D,f}+{\mathrm{\δ}}_{f}I)}^{-1}{G}_{B,f}^H{G}_{B,f}\right\}$

其中P_max{·}是指取矩阵对应最大特征值的单位特征向量，I为单位矩阵。

步骤2-4)、从选定的一系列频点中选择其它未经处理的频率点，重复上述步骤2-1)-步骤2-4)，直至计算出所有频率点的输入复权系数向量w_M,f。

所述的步骤3)具体包括以下步骤：

步骤3-1)、从一系列频率中选择频率点f,设定其目标向量d_f大小为L×1。

步骤3-2)、在之前步骤得到的输入复权系数向量w_M,f的基础上，求取在频率点f处基于均方误差最小准则的输入最优复权系数向量w_O,f。

假定在频率点f处基于均方误差最小准则的输入最优复权系数向量w_O,f的表达式为：

$w_{O,f}={\mathrm{\β}}_f{e}^{{\mathrm{j\θ}}_f}{w}_{M,f}---\left(2\right)$

其中β_f为任意大于0的缩放系数，θ_f为任意旋转角度；

所述β_f和θ_f应满足均方误差最小准则，即

$\underset{{\mathrm{\β}}_f,{\mathrm{\θ}}_f}{\min }{\left|\right|G_{B,f}{w}_{O,f}-d_f\left|\right|}_2^2={\mathrm{\β}}_f^2{w}_{M,f}^H{G}_{B,f}^H{G}_{B,f}{w}_{M,f}-2{\mathrm{\β}}_f\mathrm{Re}\left\{e^{-j{\mathrm{\θ}}_f}{w}_{M,f}^H{G}_{B,f}^H{d}_f\right\}+d_f^H{d}_f---\left(3\right)$

其中Re{·}取复数的实数部分

在公式(3)中，对β_f和θ_f分别进行求偏导，并置0可得：

${\mathrm{\θ}}_f=\∠\left(w_{M,f}^H{G}_{B,f}^H{d}_f\right)$

${\mathrm{\β}}_f=\frac{\mathrm{abs}\left(w_{M,f}^H{G}_{B,f}^H{d}_f\right)}{w_{M,f}^H\left(G_{B,f}^H{G}_{B,f}\right)w_{M,f}}---\left(4\right)$

其中∠(·)是取复数的相位部分，abs(·)是取复数的幅度部分。

所要求取的最优复权系数向量w_O,f的最终表达式为：

$w_{O,f}=\frac{\mathrm{abs}\left(w_{M,f}^H{G}_{B,f}^H{d}_f\right)}{w_{M,f}^H\left(G_{B,f}^H{G}_{B,f}\right)w_{M,f}}{e}^{j\∠\left(w_{M,f}^H{G}_{B,f}^H{d}_f\right)}{w}_{M,f}---\left(5\right)$

步骤3-3)、从选定的一系列频点中选择其它频率点，重复步骤3-1)-步骤3-2)，直至计算出所有频率点的输入复权系数向量w_O,f。

在步骤4)中，所述的窗函数可以为矩形窗，也可以为汉宁窗或者海明窗。

为了更好的理解本发明，下面结合附图和具体实施例对本发明作进一步详细描述。

在一个实施例中，如图4所示，在中国科学院声学研究所全消声室内放置一线性扬声器阵列，该扬声器阵列由8个动圈式惠威扬声器组成，相邻扬声器之间的距离为12cm，扬声器单元最大功率为15W，额定阻抗8Ω；明区和暗区处在扬声器阵列中垂线上左右各45度方向，距离扬声器阵列间距离都为1m，并与扬声器阵列处于同一水平面上；明区和暗区的控制点都由5元传声器阵列定义，间隔为8cm，传声器为丹麦B&K公司的4189传声器；实验选用的硬件装置为功放、PULSE测量仪和PC机。下面以这一实验场景为例，对本发明方法的实施过程做详细说明。

本实施例的具体实施过程包括以下几步：

(1)将系统采样率f_s设为8kHz，时域冲激响应信号的长度设为100，选择的频率点起止范围为80Hz到3920Hz，频率间隔为80Hz。首先测量扬声器阵列在各个频率点分别到明区和暗区控制点的传递函数矩阵，这可以通过采用PULSE仪器直接测量扬声器阵列到各个控制点的频响应函数，然后将频率响应导入PC平台上的Matlab软件进行处理，从而得到所有频率的传递函数矩阵G_B,f和G_D,f。图6为扬声器第一通道到明区中心控制点的频率响应示意图，从图中可以看出，在考虑的频带频响应比较平坦。因为在本实施例中，系统采样率f_s设为8kHz，只考虑人声频段，因此系统里接了一个固定的低通滤波器，截止频率设在了3.5kHz。

(2)其次，选定所有频率点的鲁棒性参数δ_f都为0.01，按常规的声能量对比控制方法计算出所有选择频率点的输入复权系数向量w_M,f。

(3)在频率点f处设定目标向量为输入最优复权系数向量w_O,f的表达式为并根据均方误差最小准则计算所有频率的β_f和θ_f，并换算成w_O,f。

(4)假定r_f，i是最优复权系数向量w_O,f的第i个系数，1≤i≤L，窗函数直接选择为矩形窗，则第i通道的频率响应为：

$\mathrm{Fre}=[0,r_{80,i},r_{160,i},\·\·\·,r_{3920,i},0,r_{3920,i}^{*},\·\·\·,r_{160,i}^{*},r_{80,i}^{*}]$

再对Fre做逆傅里叶变换得到第i通道的时域冲激响应信号，并依次求出所有通道的时域冲激响应信号。

图6分别给出了本发明方法(图6(a))和现有技术中的声能量对比控制方法(图6(b))在明区中心控制点处的冲激响应。比较两幅图可以看出，现有技术中的声能量对比控制方法则振荡幅度大，拖尾现象严重，而采用本发明方法的冲激响应拖尾衰减速度很快，更逼近于理想的Dirac-delta脉冲，因此能极大改善音质。

改善音质可以通过PESQ(Perceptual evaluation of speech quality)分数来评估，PESQ的分值为0到4.5分之间，0代表最差，4.5代表最好，处理后的语音与原语音只有纯延迟。这里采用的两个语音信号采样率为8kHz，分别命名为e1和e2，e1和e2的内容分别是“the birch canoe slid on the smooth planks”和“glue the sheet to thedark blue background”，由男声朗读。表1给出了明区中心控制点处产生的语音音质PESQ分数，从表中可以看出，现有技术中的声能量对比控制方法对语音音质有损害，而相比于现有技术中的声能量对比控制方法，本方法极大改善了语音音质，几乎可以将语音音质恢复成原始语音的音质。

表1明区中心控制点处产生的语音音质评价，采用PESQ分数

	本发明的方法	现有技术中的方法
			e1	4.412	3.907
e2	4.398	3.977

图7给出了本发明方法在宽带信号上的对比聚焦性能，从图中可以看出，在选定的频率点上，本发明方法与现有技术中的声能量对比控制方法的对比聚焦性能一样，而在其它频率点上，本发明方法的对比聚焦性能要远远好于现有技术中的声能量对比控制方法。这表明相比于现有技术中的声能量对比控制方法，本发明方法更适合处理宽带信号输入。

本实施例中虽然限定了采样频率为8kHz，并选定明区和暗区为线型区域，但这仅仅是对本发明所提供方法的一个举例说明，并不限定本发明所提供方法仅适用人说话声频率范围内，或者明区，暗区只能选择线型。事实上，本发明所提供方法能够拓展到整个可听声频率段的宽带信号并实现声能量对比聚焦。

除了上述方法外，本发明还提供了一种基于均方误差最小准则的声能量对比控制系统，该系统包括扬声器阵列设定模块、输入复权系数向量计算模块、输入复权系数向量调整模块、时域冲激响应信号生成模块；其中，

所述的扬声器阵列设定模块用于布放扬声器阵列，并设定明区和暗区的控制点；其中，所述明区为具有高声能量的区域，所述暗区为具有低声能量的区域；

所述的输入复权系数向量计算模块用于在一系列频率点上分别计算扬声器阵列的输入复权系数向量；

所述的输入复权系数向量调整模块用于设定各个频率点上的目标向量，按照均方误差最小准则，分别调整所述输入复权系数向量计算模块所得到的各个频率点上的输入复权系数向量的幅度与相位；

所述时域冲激响应信号生成模块分别将每个扬声器通道在一系列频率点上的复权系数组成频率响应，对所述频率响应加窗并进行逆傅里叶变换，得到各个通道的时域冲激响应信号。

其中，所述的输入复权系数向量调整模块还包括目标向量设定单元、输入最优复权系数向量求取单元；其中，

所述的目标向量设定单元从一系列频率中选择频率点f,设定其目标向量d_f大小为L×1；L为扬声器个数；

所述的输入最优复权系数向量求取单元在所述的输入复权系数向量计算模块得到的输入复权系数向量w_M,f的基础上，求取在频率点f处基于均方误差最小准则的输入最优复权系数向量w_O，f。

最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (8)

1.一种基于均方误差最小准则的声能量对比控制方法，包括：

步骤1)、布放扬声器阵列，设定明区和暗区的控制点；其中，所述明区为具有高声能量的区域，所述暗区为具有低声能量的区域；

步骤2)、在一系列频率点上分别计算扬声器阵列的输入复权系数向量；

步骤3)、在所述一系列频率点上设定各个频率点的目标向量，按照均方误差最小准则，分别调整步骤2)所得到的各个频率点的输入复权系数向量的幅度与相位；

步骤4)、分别将每个扬声器通道在步骤3)得到的所述一系列频率点上的复权系数组成频率响应，对所述频率响应加窗并进行逆傅里叶变换，得到各个扬声器通道的时域冲激响应信号。

2.根据权利要求1所述的基于均方误差最小准则的声能量对比控制方法，其特征在于，所述的步骤3)包括：

步骤3-1)、从一系列频率点中选择频率点f,设定其目标向量d_f大小为L×1,L为扬声器个数；

步骤3-2)、在步骤2)得到的输入复权系数向量w_M,f的基础上，求取在频率点f处基于均方误差最小准则的输入最优复权系数向量w_O,f；

步骤3-3)、从选定的一系列频率点中选择其它未经处理的频率点，重复步骤3-1)-步骤3-2)，直至计算出所有频率点的输入复权系数向量w_O，f。

3.根据权利要求2所述的基于均方误差最小准则的声能量对比控制方法，其特征在于，在所述的步骤3-2)中，求取在频率点f处基于均方误差最小准则的输入最优复权系数向量w_O，f包括：

假定在频率点f处基于均方误差最小准则的输入最优复权系数向量w_O,f的表达式为：

$w_{O,f}={\mathrm{\β}}_f{e}^{{\mathrm{j\θ}}_f}{w}_{M,f}$

其中β_f为任意大于0的缩放系数，θ_f为任意旋转角度；

使得所述β_f和θ_f满足均方误差最小准则，即

$\underset{{\mathrm{\β}}_f,{\mathrm{\θ}}_f}{\min }{\left|\right|G_{B,f}{w}_{O,f}-d_f\left|\right|}_2^2={\mathrm{\β}}_f^2{w}_{M,f}^H{G}_{B,f}^H{G}_{B,f}{w}_{M,f}-2{\mathrm{\β}}_f\mathrm{Re}\left\{e^{-j{\mathrm{\θ}}_f}{w}_{M,f}^H{G}_{B,f}^H{d}_f\right\}+d_f^H{d}_f$

其中Re{·}取复数的实数部分，G_B，f为扬声阵列到明区控制点之间的传递函数矩阵，

对β_f和θ_f分别进行求偏导，并置0得：

${\mathrm{\θ}}_f=\∠\left(w_{M,f}^H{G}_{B,f}^H{d}_f\right)$

${\mathrm{\β}}_f=\frac{\mathrm{abs}\left(w_{M,f}^H{G}_{B,f}^H{d}_f\right)}{w_{M,f}^H\left(G_{B,f}^H{G}_{B,f}\right)w_{M,f}}$

其中∠(·)是取复数的相位部分，abs(·)是取复数的幅度部分。

所要求取的最优复权系数向量w_O,f的最终表达式为：

$w_{O,f}=\frac{\mathrm{abs}\left(w_{M,f}^H{G}_{B,f}^H{d}_f\right)}{w_{M,f}^H\left(G_{B,f}^H{G}_{B,f}\right)w_{M,f}}{e}^{j\∠\left(w_{M,f}^H{G}_{B,f}^H{d}_f\right)}{w}_{M,f}.$

4.根据权利要求1所述的基于均方误差最小准则的声能量对比控制方法，其特征在于，在步骤1)中，所布放的扬声器阵列为线性阵列或圆形阵列，或随机阵列。

5.根据权利要求1所述的基于均方误差最小准则的声能量对比控制方法，其特征在于，在步骤1)中，所述明区或暗区的区域形状是方形或圆形或线型。

6.根据权利要求1所述的基于均方误差最小准则的声能量对比控制方法，其特征在于，在步骤4)中，对频率响应加窗时所要用到的窗函数为矩形窗或汉宁窗或者海明窗。

7.一种基于均方误差最小准则的声能量对比控制系统，其特征在于，该系统包括扬声器阵列设定模块、输入复权系数向量计算模块、输入复权系数向量调整模块、时域冲激响应信号生成模块；其中，

所述的扬声器阵列设定模块用于布放扬声器阵列，并设定明区和暗区的控制点；其中，所述明区为具有高声能量的区域，所述暗区为具有低声能量的区域；

所述的输入复权系数向量计算模块用于在一系列频率点上分别计算扬声器阵列的输入复权系数向量；

所述的输入复权系数向量调整模块用于在所述一系列频率点上设定各个频率点的目标向量，按照均方误差最小准则，分别调整所述输入复权系数向量计算模块所得到的各个频率点上的输入复权系数向量的幅度与相位；

所述时域冲激响应信号生成模块分别将每个扬声器通道在一系列频率点上的复权系数组成频率响应，对所述频率响应加窗并进行逆傅里叶变换，得到各个通道的时域冲激响应信号。

8.根据权利要求7所述的基于均方误差最小准则的声能量对比控制系统，其特征在于，所述的输入复权系数向量调整模块还包括目标向量设定单元、输入最优复权系数向量求取单元；其中，

所述的目标向量设定单元从一系列频率点中选择频率点f,设定其目标向量d_f大小为L×1；

所述的输入最优复权系数向量求取单元在所述的输入复权系数向量计算模块得到的输入复权系数向量w_M,f的基础上，求取在频率点f处基于均方误差最小准则的输入最优复权系数向量w_O，f。