一种麦克风阵列频响校准方法技术

本发明专利技术公开了一种麦克风阵列频响校准方法，包括以下步骤：S1：电音频信号w(n)经扬声器和声传播信道到达麦克风阵列，麦克风阵列输入声信号为x(n)；S2：声信号x(n)经过k个不同的麦克风及前置放大器，分别得到不同的电信号x1(n)～xk(n)；S3：电信号x1(n)～xk(n)分别当作各路校准滤波器的输入信号用于调整各校准滤波器的滤波系数；S4：根据上述调整方式，计算各校准滤波器的滤波系数，以完成麦克风阵列的频响校准。该方法使同一阵列中不同麦克风的频响接近一致，从而提高了麦克风阵列的信号处理能力。

【技术实现步骤摘要】
一种麦克风阵列频响校准方法
本专利技术涉及信号处理
，尤其涉及一种麦克风阵列频响校准方法。
技术介绍
声音是人类广泛使用的信息载体，音频信号是信号处理技术的一个重要研究对象。麦克风作为声音传感器，在声音信号处理中发挥着重要的作用。由于单个麦克风具有声音拾取范围有限及噪声抑制能力弱，无法满足日益增长的音频信号质量要求，人们提出了麦克风阵列技术。由于其在时域、频域外还引进了空间域的信息，故对声音信号的处理能力得到增强，目前已成为众多高质量音频处理应用中的首选目标。作为声音信息采集工具，麦克风阵列对阵列音频处理领域中有着重要的作用。麦克风拾取声音过程一般由以下两个部分组成：(1)声信号转换成电信号；(2)电信号再经过前置放大器得到足够大的电信号。由于麦克风声-电参数的离散性或放大器中元器件参数的离散性等原因，不同麦克风之间不论是声/电信号转换模块还是前置放大器模块均会存在一些差异。这些差异在廉价麦克风或不同品牌的麦克风之间更为明显。这种差异导致了同一阵列中各麦克风之间的频率响应互不相同。具体表现为：同一位置不同的麦克风拾取同一段声信号时，所得到的电信号之间均存在着一定的差异。麦克风之间的频响差异会直接影响后续的音频信号处理结果，如声场感知、语音增强等算法的性能。然而，由于硬件电路等原因，同一阵列中的各个麦克风频率响之间往往存在一些差异。麦克风阵列的一些应用如声场感知、语音增强等都将受到上述差异的影响，导致处理性能下降。因此，校准同一阵列中各个麦克风的频响，使它们接近一致就成为麦克风阵列中的一个现实技术需求。现有技术中关于麦克风阵列的增益校准技术中由于只校准麦克风阵列中各麦克风之间的增益，并不能保证各麦克风在频响上相互一致，导致校准后各麦克风的幅频响应和相频响应仍然存在一定差异
技术实现思路
根据现有技术存在的问题，本专利技术公开了一种麦克风阵列频响校准方法，包括以下步骤：S1：电音频信号w(n)经扬声器和声传播信道到达麦克风阵列，麦克风阵列输入声信号为x(n)；S2：声信号x(n)经过k个不同的麦克风及前置放大器，分别得到不同的电信号x1(n)～xk(n)；S3：电信号x1(n)～xk(n)分别当作各路校准滤波器的输入信号用于调整各校准滤波器的滤波系数；所述准滤波器的滤波系数的调整方式为：通过调整各路滤波器的滤波系数使输出信号y1(n)～yk(n)均向目标信号d(n)逼近，逼近原则是使y1(n)～yk(n)与d(n)之间的均方误差分别最小；S4：根据上述调整方式，计算各校准滤波器的滤波系数，以完成麦克风阵列的频响校准。进一步的，当各路校准滤波器输出信号y1(n)～yk(n)向目标信号d(n)逼近实现阵列频响校准时采用如下方式：设一帧数字信号的长度为N，则第i路校准滤波器的输入信号xi(n)与目标信号d(n)写成以下向量形式xi＝[xi(0),xi(1),...,xi(N-1)]T(2)d＝[d(0),d(1),...,d(N-1)]T(3)其中，[]T表示向量或矩阵的转置，列出以下代价函数其中，Σ表示级数求和，使代价函数J最小，求出目标信号d(n)为d＝(x0+x1+...+xk)/k(5)即，目标信号的向量形式d等于对各滤波器的输入向量之和求平均。当目标信号d(n)确定后，以第i路滤波器输出信号yi(n)与目标信号d(n)之间的均方误差最小为准则，计算出第i路滤波器系数hi(n)采用如下方式：yi＝[yi(0),yi(1),...,yi(N-1)]T(6)则yi(n)与d(n)之间误差信号ei(n)写成以下向量形式ei＝[ei(0),ei(1),...,ei(N-1)]T＝yi-d(7)列出另一组代价函数Fi＝eiTei(8)其中，i＝1,2,...,k，k为待校准麦克风的个数，使Fi最小，计算出第i路校准滤波器的滤波系数hi(n)，计算规则为：hi＝[hi(0),hi(1),...,hi(Mi-1)]T＝(XiTXi)-1DTxi_m(9)其中，xi_m是由第i路校准滤波器输入向量的前N-Mi+1个元素组成的向量，即xi_m＝[xi(0),xi(1),...,xi(N-Mi)]T(10)而式(9)中的Xi与D分别为当i从1取到k时，用上述方法计算出每路校准滤波器的滤波系数，完成麦克风阵列频响校准。当分别对多帧输入数据时进行校准，对得到的多组滤波器系数进行融合、得到滤波器系数采用如下方式：设第i路麦克风经过Q帧数据进行校准，根据式(9)得到Q组滤波器系数向量分别为hi_1、hi_2、...hi_Q，其中，hi_Q表示第i路麦克风在第Q帧输入数据得到的滤波器系数向量；此外，Q个滤波器系数向量的平均向量定义为Hi＝(hi_1+hi_2+...+hi_Q)/Q(13)第i路第q帧滤波器系数向量的方差向量定义为δi_q＝(hi_q-Hi)·(hi_q-Hi)(14)其中，符号“·”表示向量点乘，q＝1,2,...,Q，第i路第q帧滤波器系数向量的方差倒数向量θi_q被定义为θi_q＝[1/δi_q(0),1/δi_q(1),...,1/δi_q(Mi-1)]T(15)即，向量θi_q中的每个元素分别是向量δi_q所对应元素的倒数，关于第i路校准滤波器，其滤波器系数向量的方差倒数向量和被定义为θi＝θi_1+θi_2+。。。+θi_Q(16)根据θi和θi_q(q＝1,2,...,Q)得出Q个滤波器系数权重向量Ci_1，Ci_2，...，Ci_q，...，Ci_Q，其中，Ci_q＝[θi_q(0)/θi(0),θi_q(1)/θi(1),...,θi_q(Mi-1)/θi(Mi-1)]T(17)根据hi_q和Ci_q得出第i路校准滤波器最终的滤波系数hi＝hi_1·Ci_1+hi_2·Ci_2+...+hi_Q·Ci_Q(18)其中，i＝1,2,...,k，k为麦克风数量。由于采用了上述技术方案，本专利技术提供的一种麦克风阵列频响校准方法,该方法使同一阵列中不同麦克风的频响接近一致，从而提高了麦克风阵列的信号处理能力。因此，一些麦克风阵列的应用，如声源定位或语音增强的性能都采用本专利技术公开的麦克风阵列频响校准方法来对信号进行处理。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术麦克风阵列频响校准方法的流程图；图2为本专利技术中麦克风与扬声器、麦克风与麦克风之间的距离示意图；图3为本专利技术中校准滤波器的工作原理图；图4为本专利技术中麦克风阵列频响校准方法的效果示意图；图5为本专利技术中麦克风阵列频响校准方法的效果示意图；图6为本专利技术中麦克风阵列频响校准方法的效果示意图；图7为本专利技术中麦克风阵列频响校准方法的效果示意图。具体实施方式为使本专利技术的技术方案和优点更加清楚，下面结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚完整的描述：如图1所示的一种麦克风阵列频响校准方法，具体包括以下步骤：S1：电音频信号w(n)经扬声器和声传播信道到达麦克风阵列，麦克风阵列输入声信号为x(n)；S2：声信号本文档来自技高网...

【技术保护点】
一种麦克风阵列频响校准方法，其特征在于：包括以下步骤：S1：电音频信号w(n)经扬声器和声传播信道到达麦克风阵列，麦克风阵列输入声信号为x(n)；S2：声信号x(n)经过k个不同的麦克风及前置放大器，分别得到不同的电信号x1(n)～xk(n)；S3：电信号x1(n)～xk(n)分别当作各路校准滤波器的输入信号用于调整各校准滤波器的滤波系数；所述准滤波器的滤波系数的调整方式为：通过调整各路滤波器的滤波系数使输出信号y1(n)～yk(n)均向目标信号d(n)逼近，逼近原则是使y1(n)～yk(n)与d(n)之间的均方误差分别最小；S4：根据上述调整方式，计算各校准滤波器的滤波系数，以完成麦克风阵列的频响校准。

【技术特征摘要】
1.一种麦克风阵列频响校准方法，其特征在于：包括以下步骤：S1：电音频信号w(n)经扬声器和声传播信道到达麦克风阵列，麦克风阵列输入声信号为x(n)；S2：声信号x(n)经过k个不同的麦克风及前置放大器，分别得到不同的电信号x1(n)～xk(n)；S3：电信号x1(n)～xk(n)分别当作各路校准滤波器的输入信号用于调整各校准滤波器的滤波系数；所述准滤波器的滤波系数的调整方式为：通过调整各路滤波器的滤波系数使输出信号y1(n)～yk(n)均向目标信号d(n)逼近，逼近原则是使y1(n)～yk(n)与d(n)之间的均方误差分别最小；S4：根据上述调整方式，计算各校准滤波器的滤波系数，以完成麦克风阵列的频响校准。2.根据权利要求1所述的一种麦克风阵列频响校准方法，其特征还在于：当各路校准滤波器输出信号y1(n)～yk(n)向目标信号d(n)逼近实现阵列频响校准时采用如下方式：设一帧数字信号的长度为N，则第i路校准滤波器的输入信号xi(n)与目标信号d(n)写成以下向量形式xi＝[xi(0),xi(1),...,xi(N-1)]T(2)d＝[d(0),d(1),...,d(N-1)]T(3)其中，[]T表示向量或矩阵的转置，列出以下代价函数其中，Σ表示级数求和，使代价函数J最小，求出目标信号d(n)为d＝(x0+x1+...+xk)/k(5)即，目标信号的向量形式d等于对各滤波器的输入向量之和求平均。3.根据权利要求1所述的一种麦克风阵列频响校准方法，其特征还在于：当目标信号d(n)确定后，以第i路滤波器输出信号yi(n)与目标信号d(n)之间的均方误差最小为准则，计算出第i路滤波器系数hi(n)采用如下方式：yi＝[yi(0),yi(1),...,yi(N-1)]T(6)则yi(n)与d(n)之间误差信号ei(n)写成以下向量形式ei＝[ei(0),ei(1),...,ei(N-1)]T＝yi-d(7)列出另一组代价函数Fi＝eiTei(8)其中，i＝1,2,...,k，k为待校准麦克风的个数，使Fi最小，计算出第i路校准滤波...

【专利技术属性】
技术研发人员：陈喆，殷福亮，呼德，
申请(专利权)人：大连理工大学，
类型：发明
国别省市：辽宁,21