《基于PSOMVDR的声场重构与鬼影抑制方法.pdf》由会员分享，可在线阅读，更多相关《基于PSOMVDR的声场重构与鬼影抑制方法.pdf（16页完整版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 104008287 A (43)申请公布日 2014.08.27 C N 1 0 4 0 0 8 2 8 7 A (21)申请号 201410222150.7 (22)申请日 2014.05.23 G06F 19/00(2011.01) G01H 17/00(2006.01) (71)申请人北京科技大学 地址 100083 北京市海淀区学院路30号 (72)发明人魏龙 黎敏 阳建宏 杨德斌 付强 秦胜 杨海波 孙冬柏 (74)专利代理机构北京金智普华知识产权代理 有限公司 11401 代理人皋吉甫 (54) 发明名称 基于PSO-MVDR的声场重构与鬼影抑制方法 (5。
2、7) 摘要 发明一种基于PSO-MVDR的声场重构与鬼影 抑制方法。该方法首先利用PSO算法实现了不等 间距阵列的阵元优化布置，抑制栅瓣值的同时，并 尽可能减小原旁瓣值的增加量；在优化不等间距 阵列的基础上，利用MVDR算法进行声场重构，抑 制旁瓣对声场重构的影响，实现鬼影的抑制，提高 声源定位的准确性。该方法的优点在于：结合了 不等间距阵列抑制栅瓣和MVDR算法较强的空域 滤波特性的优势，能够利用阵元数量较少的传声 器阵列，实现声场重构中的鬼影抑制。 (51)Int.Cl. 权利要求书3页 说明书7页 附图5页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书3页 说明。
3、书7页 附图5页 (10)申请公布号 CN 104008287 A CN 104008287 A 1/3页 2 1.一种基于PSO-MVDR的声场重构与鬼影抑制方法，其特征在于：该方法是在对等间 距十字传声器阵列进行优化的基础上，采用MVDR波束形成方法进行声场重构，得到声压云 图；具体包括以下步骤： 步骤1，将含有一定阵元的等间距十字传声器阵列PSO粒子群算法，寻找最优的各阵元 移动距离的组合，获得没有栅瓣，旁瓣影响最小的不等间距十字阵；在通过PSO算法得到优 化的不等间距十字阵之后，该阵列的栅瓣已尽可能多地被抑制，并已退化成为旁瓣； 步骤2，为了进一步削弱旁瓣对声场重构的影响，在获得优化不。
4、等间距十字阵的基础 上，使用借助MVDR波束形成算法来进行声场的重构，利用MVDR算法较强的空域滤波特性， 最大限度降低旁瓣对声场重构的影响，提高声场重构的准确性。 2.根据权利要求1所述的一种基于PSO-MVDR的声场重构与鬼影抑制方法，其特征在 于：所述步骤1中的含有一定阵元的等间距十字传声器阵列PSO粒子群算法，包括：首先， 确立阵列优化的最终目标，即在保持十字阵列形状不变的前提下，寻找最优的各阵元移动 距离的组合，使得在不等间距阵列情况下的声场重构的最大旁瓣级SL最小；然后，初始化 粒子群，进行多次迭代运算，每一次迭代运算均需要更新粒子群算法的全局最优解，并记录 全局最优解的最大旁瓣值。
5、。当最大旁瓣值SL趋于收敛时，算法迭代结束，得到一组最优的 阵元位置坐标解。 3.根据权利要求1或2所述的一种基于PSO-MVDR的声场重构与鬼影抑制方法，其特征 在于：所述步骤1中的含有一定阵元的等间距十字传声器阵列PSO粒子群算法具体包括以 下步骤： 首先，用声场重构结果中最大旁瓣级SL的大小作为优化目标，以此来衡量阵列的优化 效果；在传声器阵元数量有限的情况下，为确保重构后的声场云图具有足够的空间分辨率， 选取十字传声器阵列来进行位置的优化；为保持阵列的十字形状，十字交叉点处的阵元位 置在优化后保持不变；由上述分析可知，阵列优化的目标为：寻找最优的各阵元移动距离 的组合，使得在不等间距阵。
6、列情况下的声场重构的最大旁瓣级SL最小； 然后，进行算法的初始化和迭代运算： 设Z i z i1 ,z i2 ,z ij ,z i13 为某一种可能的阵元位置坐标解，称Z i 为一个粒子，其 中Z ij 表示第i个粒子中第j个阵元的坐标位置；由m个粒子组成一个组群ZZ 1 ,Z 2 , Z i ,Z m ，将粒子群初始化； 在迭代运算中，每个粒子的速度，记做V i v i1 ,v i2 ,v id v i13 ，其中v id 表示在第i 个粒子中第d个阵元坐标的速度；在每一次迭代计算中，每个粒子可以根据式(1)来更新速 度： v id (t+1)wv id (t)+ 1 rand()(p id。
7、 -z id (t)+ 2 rand()(p gd -z id (t) (1) 式(1)中，v id (t)表示第i个粒子第d个阵元坐标在第t次迭代中的速度，z id (t)则表 示经过第t次迭代后，第i个粒子中第d个阵元的坐标解，w为惯性权重， 1 、 2 为加速常 数，rand()为0到1之间的随机数；p id 和p gd 分别表示自搜索开始到经过t次迭代后每个粒 子和整个粒子群所搜索到的最优坐标解；在更新了粒子速度的基础上，则可得到第(t+1) 次迭代后第i个粒子中第d个阵元的坐标解z id (t+1)： z id (t+1)z id (t)+v id (t+1) (2) 粒子更新后，计。
8、算各粒子的最大旁瓣级，同已有的各粒子和整个粒子群的最优坐标 权 利 要 求 书CN 104008287 A 2/3页 3 解的最大旁瓣级进行比较，以更新各粒子和整个粒子群的最优坐标解，继续下一次迭代运 算； 最后，当最大旁瓣级SL达到最小，并趋于收敛时，整个算法迭代结束，即可得到一组最 优的阵元位置坐标解； 综上所述，为了减小声场重构中的栅瓣值，同时也减小原有旁瓣值的增加量，利用PSO 算法，在阵元个数有限的情况下，对不等间距十字阵进行了位置的优化，以最大旁瓣级SL 最小作为优化目标，最终获得了阵列的最优布置方案。 4.根据权利要求1所述的一种基于PSO-MVDR的声场重构与鬼影抑制方法，其特。
9、征在 于：所述步骤2中所述的针对不等间距十字阵列的MVDR波束形成方法，其原理如下所述： 波束形成算法中，阵列的声压响应p为： 其中，w为加权向量x为各通道信号x 1 ，x 2 ，x M T ,其中，m1、2M； 在MVDR算法中，加权向量w的形式为： 其中，为声源恰好从重构方向入射时的阵列流形向量， n 为干扰信号的协方 差矩阵。而在实际应用中，由于接收数据中的干扰信号的协方差矩阵 n 是无法估计的，往 往用接收数据的互功率谱矩阵R x 来代替： R x Ex()x H ( (5) 于是，式(4)便转化为式(6) 若使声源的入射方向角为一定值 0 ，并将式(6)代入式(1)中，则得到基于MV。
10、DR 算法的声场重构响应公式： 为便于求解阵列的声压响应在此通过MVDR波束输出功率进行计算。MVDR波 束输出功率谱为： 结合式(4)、式(7)、式(9)可得 由式(9)得 权 利 要 求 书CN 104008287 A 3/3页 4 在实际计算过程中，将某一声场平面离散成多个坐标点，将该声场平面称为重建面，即 可根据式(10)计算出该重建区域中每一坐标点处的重构声压，进而得到重构的声压云图。 权 利 要 求 书CN 104008287 A 1/7页 5 基于 PSO-MVDR 的声场重构与鬼影抑制方法 技术领域 0001 本发明涉及一种声场重构和鬼影抑制方法，该方法在有限空间、有限传声器个。
11、数 的情况下，能够有效抑制声场重构过程中栅瓣和旁瓣带来的干扰，提高声场重构的精度。 背景技术 0002 波束形成(Beamforming)是阵列信号处理的一个重要任务，其主要功能包括：形 成基阵接收系统的方向性；进行空域滤波，抑制空间干扰与环境噪声，提高信噪比；估计信 号到达方向，为信号源定位创造条件等。通过波束形成处理，可实现对声场声源的重构与定 位。 0003 声源识别的过程就是将重构声场空间离散成一系列的重构点，利用波束形成算法 计算出重构点的重构声压值，绘制出声压云图，由图中热点的位置来判断声源位置。传声器 阵列的声源识别效果与其结构形状，阵元数目及处理算法等因素决定。一个合适的传声器。
12、 阵列能够抑制栅瓣和旁瓣，提高声源的辨识精度，提高声场的重构质量。在基于传声器阵列 的声场重构研究中，前人搭建网格阵、十字阵、圆形阵等一些等间距的阵列来进行声场实测 与分析。但等间距阵列有其固有缺陷，即：为了保证较小的主瓣宽度，提高声场重构分辨率， 需要维持一个较大的阵列尺寸，使得相邻阵元的间距增加，并导致栅瓣的出现，极大地削弱 了阵列的声场重构能力。针对等间距阵列重构能力的缺陷，有学者提出了利用不等间距阵 列进行声场重构的方法。该方法通常是在等间距阵列的基础上，利用寻优算法，对阵元坐标 分布进行优化设计，获得优化不等间距阵。不等间距阵能够抑制栅瓣 8 ，在一定程度上可提 高声场重构的质量。但。
13、是，阵列中较大的旁瓣幅值依然会给真实声源的正确辨识带来干扰。 针对上述问题，若能够通过阵列优化，实现对栅瓣的抑制，同时在此基础上运用有效地算法 抑制较大旁瓣的产生，获得准确声场重构云图，实现鬼影抑制。本发明正是为了实现这一目 的而产生的。 0004 2012-09-19公开的专利：CN102680071A，题目：采用振速测量和局部近场声全息 法的噪声源识别方法，发明人：杨德森。该发明提供的是一种采用振速测量和局部近场声 全息法的噪声源识别方法。测量位于近场的有限孔径测量面上的法向质点振速，对其进行 补零扩展，最后计算得到声源面上的声压和法向质点振速。该发明的特点是采用法向质点 振速作为输入量进。
14、行声场重建，与传统迭代局部近场声全息法相比，该方法计算简单，效率 较高。 0005 2013-01-30公开的专利：CN102901950A，题目：平面阵列识别声源三维坐标的方 法，发明人：卢奂采；李春晓；金江明；丁浩；梅东挺。该发明提供的是一种采用平面阵列 识别定位声源三维坐标的方法。基于球面波声场模型,采用波束形成方法,以三维空间网 格化后的网格交点作为重构点,进行不同距离的平面声场重构,比较沿距离方向(即z方 向)上不同重构面的波束形成器的输出值,最大输出值对应的z轴坐标,即确定为声源的 距离。最后根据声源距离对应的重构面上波束形成器输出的最大值点,确定声源在x和y 轴坐标。该发明的优点。
15、在于可用二维的平面阵列,识别声源的三维坐标,可直接应用于现 说 明 书CN 104008287 A 2/7页 6 有的仅能识别声源两维坐标的平面阵列声全息设备上,实验已验证了该发明方案的有效 性和可靠性。 0006 本发明与现有发明专利的不同在于：将PSO粒子群算法与MVDR波束形成算法相结 合，能够抑制声场重构中的鬼影，提高声场重构准确性。具体来说，该方法分为两步：首先， 利用PSO算法对一个十字等间距传声器阵列的布局进行优化，获得阵元位置布局最优的不 等间距阵列；在此基础上，对各阵元采集到的信号进行MVDR波束形成运算，获得声场重构 云图。本发明结合了不等间距阵列抑制栅瓣和MVDR波束形成。
16、算法具有较强的空域滤波特 性的优势，能够利用阵元数量较少的传声器阵列，实现声场重构中的鬼影抑制。 发明内容 0007 本发明的目的是提出一种能够有效抑制声场重构中鬼影现象的方法。该方法在有 限空间、有限传声器个数的情况下，能够有效避免声场重构云图中栅瓣和旁瓣的存在，提高 声场重构的声源辨识精度，实现鬼影的抑制。 0008 本发明的技术方案是：一种基于PSO-MVDR的声场重构与鬼影抑制方法，该方法是 在对等间距十字传声器阵列进行优化的基础上，采用MVDR波束形成方法进行声场重构，得 到声压云图。具体包括以下步骤： 0009 步骤1，将含有一定阵元的等间距十字传声器阵列PSO粒子群算法，寻找最优。
17、的各 阵元移动距离的组合，获得没有栅瓣，旁瓣影响最小的不等间距十字阵。在通过PSO算法得 到优化的不等间距十字阵之后，该阵列的栅瓣已尽可能多地被抑制，并已退化成为旁瓣； 0010 步骤2，为了进一步削弱旁瓣对声场重构的影响，在获得优化不等间距十字阵的基 础上，使用借助MVDR波束形成算法来进行声场的重构，利用MVDR算法较强的空域滤波特 性，最大限度降低旁瓣对声场重构的影响，提高声场重构的准确性。 0011 优选的，所述步骤1中的含有一定阵元的等间距十字传声器阵列PSO粒子群算法， 包括：首先，确立阵列优化的最终目标，即在保持十字阵列形状不变的前提下，寻找最优的 各阵元移动距离的组合，使得在不。
18、等间距阵列情况下的声场重构的最大旁瓣级SL最小；然 后，初始化粒子群，进行多次迭代运算，每一次迭代运算均需要更新粒子群算法的全局最优 解，并记录全局最优解的最大旁瓣值。当最大旁瓣值SL趋于收敛时，算法迭代结束，得到一 组最优的阵元位置坐标解。 0012 进一步，如所述步骤1中的基于等间距十字阵列的粒子群优化算法，具体包括以 下步骤： 0013 首先，用声场重构结果中最大旁瓣级SL的大小作为优化目标，以此来衡量阵列的 优化效果。在传声器阵元数量有限的情况下，为确保重构后的声场云图具有足够的空间分 辨率，选取十字传声器阵列来进行位置的优化。为保持阵列的十字形状，十字交叉点处的阵 元位置在优化后保持。
19、不变。由上述分析可知，阵列优化的目标为：寻找最优的各阵元移动距 离的组合，使得在不等间距阵列情况下的声场重构的最大旁瓣级SL最小。 0014 然后，进行算法的初始化和迭代运算。 0015 设Z i z i1 ,z i2 ,z ij ,z i13 为某一种可能的阵元位置坐标解，称Z i 为一个 粒子，其中Z ij 表示第i个粒子中第j个阵元的坐标位置。由m个粒子组成一个组群Z Z 1 ,Z 2 ,Z i ,Z m ，将粒子群初始化。 说 明 书CN 104008287 A 3/7页 7 0016 在迭代运算中，每个粒子的速度，记做V i v i1 ,v i2 ,v id v i13 ，其中v i。
20、d 表示在 第i个粒子中第d个阵元坐标的速度。在每一次迭代计算中，每个粒子可以根据式(1)来 更新速度： 0017 v id (t+1)wv id (t)+ 1 rand()(p id -z id (t)+ 2 rand()(p gd -z id (t) (1) 0018 式(1)中，v id (t)表示第i个粒子第d个阵元坐标在第t次迭代中的速度，z id (t) 则表示经过第t次迭代后，第i个粒子中第d个阵元的坐标解，w为惯性权重， 1 、 2 为加 速常数，rand()为0到1之间的随机数。p id 和p gd 分别表示自搜索开始到经过t次迭代后 每个粒子和整个粒子群所搜索到的最优坐标解。
21、。在更新了粒子速度的基础上，则可得到第 (t+1)次迭代后第i个粒子中第d个阵元的坐标解z id (t+1)： 0019 z id (t+1)z id (t)+v id (t+1) (2) 0020 粒子更新后，计算各粒子的最大旁瓣级，同已有的各粒子和整个粒子群的最优坐 标解的最大旁瓣级进行比较，以更新各粒子和整个粒子群的最优坐标解，继续下一次迭代 运算。 0021 最后，当最大旁瓣级SL达到最小，并趋于收敛时，整个算法迭代结束，即可得到一 组最优的阵元位置坐标解。 0022 综上所述，为了减小声场重构中的栅瓣值，同时也减小原有旁瓣值的增加量，利用 PSO算法，在阵元个数有限的情况下，对不等间。
22、距十字阵进行了位置的优化，以最大旁瓣级 SL最小作为优化目标，最终获得了阵列的最优布置方案，从而为鬼影的抑制奠定基础。 0023 进一步，如步骤2中所述的针对不等间距十字阵列的MVDR波束形成方法，其原理 如下所述： 0024 波束形成算法中，阵列的声压响应p为： 0025 0026 其中，w为加权向量x为各通道信号x 1 ，x 2 ， x M T ,其中，m1、2M。 0027 在MVDR算法中，加权向量w的形式为： 0028 (4)其中，为声源恰好从重构方向入射时的阵列 流形向量， n 为干扰信号的协方差矩阵。而在实际应用中，由于接收数据中的干扰信号的 协方差矩阵 n 是无法估计的，往往用。
23、接收数据的互功率谱矩阵R x 来代替： 0029 R x Ex()x H () (5) 0030 于是，式(4)便转化为式(6) 0031 0032 若使声源的入射方向角为一定值 0 ，并将式(6)代入式(1)中，则得到基于 说 明 书CN 104008287 A 4/7页 8 MVDR算法的声场重构响应公式： 0033 0034 为便于求解阵列的声压响应在此通过MVDR波束输出功率进行计算。 MVDR波束输出功率谱为： 0035 0036 结合式(4)、式(7)、式(9)可得 0037 0038 由式(9)得 0039 0040 在实际计算过程中，将某一声场平面离散成多个坐标点，将该声场平面。
24、称为重建 面，即可根据式(10)计算出该重建区域中每一坐标点处的重构声压，进而得到重构的声压 云图。 0041 上述发明的优点在于：(1)通过PSO算法获得的优化不等间距阵列，不但能够抑制 栅瓣，也能够抑制旁瓣值的增加，有利于获得较为准确的声源识别结果；(2)与常规波束形 成算法相比，基于MVDR的波束形成算法具有较强的空域滤波特性，可进一步提高声源的辨 识精度，保证声场的重构质量。 附图说明 0042 图1是PSO-MVDR算法流程图； 0043 图2是原等间距十字传声器阵列位置分布图； 0044 图3是实验参数设置图； 0045 图4是优化后不等间距阵阵元坐标图； 0046 图5是优化阵元。
25、位置分布图； 0047 图6是全消音室试验环境示意图； 0048 图7-10分别为2000hz，4000hz，6000hz，8000hz条件下，未采用PSO-MVDR算法和 采用PSO-MVDR算法的声压重构云图。 具体实施方式 0049 为了验证本发明的可行性和正确性，在全消音室中进行实验验证。图1是 PSO-MVDR算法流程图；按照图1的流程图实施如下: 说 明 书CN 104008287 A 5/7页 9 0050 首先，在全消音室内进行声场测量与重构装置的搭建，图6为实验装置示意图，等 间距十字传声器阵列所在平面与声源重构面相距1m，正方形重建面边长为1.4m，所选用的 声源为中高频体。
26、积声源，发声的频率范围为200Hz至10000Hz，可近似为点声源。实验中所 用的声音频率为2000Hz、4000Hz、6000Hz和8000Hz的纯音。传声器为自由场传声器，全指 向性，并具有较高的灵敏度。所布置的传声器阵列分别为13阵元等间距十字阵和优化后的 不等间距十字阵。采集设备是基于PXI-e总线的数据采集系统，能实现多通道高速同步数 据采集，采样频率为44100Hz。重构实验参数如图3所示。 0051 以下详述利用本发明对声场进行重构和鬼影抑制的实现过程。 0052 首现，利用PSO粒子群算法对图2中13个阵元的等间距十字阵列进行优化。初始 阵元间距为d0.1m，1至13号各阵元初。
27、始位置坐标分别为：(0.1,0)，(0,0.1)，(-0.1,0)， (0,-0.1)，(0.2,0)，(0,0.2)，(-0.2,0)，(0,-0.2)，(0.3,0)，(0,0.3)，(-0.3,0)，(0,0)。 0053 为保持阵列的十字形状，十字交叉点处的第13号阵元位置在优化后保持不变。经 过优化处理后，设1至12号阵元坐标依次变为：(0.1+ d 1 ,0)，(0,0.1+ d 2 )，(-0.1+ d 3 ,0)， (0,-0.1+ d 4 )，(0.2+ d 5 ,0)，(0,0.2+ d 6 )，(-0.2+ d 7 ,0)，(0,-0.2+ d 8 )，(0.3+ d 9。
28、 ,0)， (0,0.3+ d 10 )，(-0.3+ d 11 ,0)，(0,-0.3+ d 12 )。结合PSO算法原理，以最大旁瓣值SL作为适 应度函数，寻找最优的各阵元移动距离的组合，使得在不等间距阵列情况下的声场重构的 最大旁瓣级SL最小，即使得式(1)成立： 0054 minSL( d 1 , d 2 d 12 ) (1) 0055 式(1)中，SL为最大旁瓣级， d 1 , d 2 d 12 为各阵元移动的距离， d 1 , d 3 , d 5 d 11 为1、3、511号阵元的水平位移， d 2 , d 4 , d 6 d 12 为2、4、612号阵元的竖直位 移。设Z i z。
29、 i1 ,z i2 ,z ij ,z i13 为某一种可能的阵元位置坐标解，称Z i 为一个粒子，其 中z ij 表示第i个粒子中第j个阵元的坐标位置。由20个粒子组成一个组群ZZ 1 ,Z 2 , Z i ,Z m 。每个粒子的速度，记做V i v i1 ,v i2 ,v id v i13 ，其中v id 表示在第i个粒子中 第d个阵元坐标的速度。在每一次迭代计算中，每个粒子可以根据式(2)来更新速度： 0056 v id (t+1)wv id (t)+1rand()(p id -z id (t)+ 2 rand()(p gd -z id (t) (2) 0057 式(2)中，v id (t。
30、)表示第i个粒子第d个阵元坐标在第t次迭代中的速度，z id (t) 则表示经过第t次迭代后，第i个粒子中第d个阵元的坐标解，w为惯性权重， 1 、 2 为加 速常数，rand()为0到1之间的随机数。p id 和p gd 分别表示自搜索开始到经过t次迭代后 每个粒子和整个粒子群所搜索到的最优坐标解。在更新了粒子速度的基础上，则可得到第 (t+1)次迭代后第i个粒子中第d个阵元的坐标解z id (t+1)： 0058 z id (t+1)z id (t)+v id (t+1) (3) 0059 当最大旁瓣值SL达到最小，并趋于收敛时，整个算法迭代结束，可以得到一组图4 所示的最优阵元位置坐标解。
31、，最终各阵元优化分布结果如图5所示。 0060 然后，为了进一步削弱旁瓣对声场重构的影响，更好地抑制鬼影，提高声场重构的 准确性，需借助MVDR算法来进行声场的重构。 0061 波束形成算法中，阵列的声压响应p为： 0062 说 明 书CN 104008287 A 6/7页 10 0063 其中，w为加权向量x为各通道信号x 1 ，x 2 ， x M T ,其中，m1、2M。 0064 在MVDR算法中，加权向量w的形式为： 0065 (5)其中，为声源恰好从重构方向入射时的阵列 流形向量， n 为干扰信号的协方差矩阵。而在实际应用中，由于接收数据中的干扰信号的 协方差矩阵 n 是无法估计的，。
32、往往用接收数据的互功率谱矩阵R x 来代替： 0066 R x Ex()x H () (6) 0067 于是，式(5)便转化为式(7) 0068 0069 MVDR波束输出功率谱： 0070 0071 结合式(4)、式(7)、式(9)可得 0072 0073 由式(9)得 0074 0075 在实际计算过程中，将某一声场平面离散成多个坐标点，将该声场平面称为重建 面，根据式(10)即可计算出重建面区域中每一坐标点处的重构声压，进而得到重构的声压 云图。 0076 图7-图10为全消音室声环境下重构获得的重构声压等高线云图。图中数字表示 云图相应热点位置处的声压值。图7和图8分别为2000Hz和。
33、4000Hz下通过原始的等间距 阵常规波束形成方法和PSO-MVDR方法获得的重构声压等高线云图。由于频率较低时不会 出现栅瓣，因此均未出现因栅瓣而导致的鬼影，而在图7(b)和8(b)中，由于PSO-MVDR方法 可以削弱旁瓣的影响，重构声源热点更加突出，位置更加精确。 0077 图9和图10分别为6000Hz和8000Hz下通过原始的等间距阵常规波束形成方法 和PSO-MVDR方法获得的重构声压等高线云图。在图9(a)和图10(a)的等间距阵列重构云 图中，重建面四边中点和拐角处均由于栅瓣的存在而产生鬼影；而在如图9(b)和图10(b) 的优化后不等间距阵列的重构云图中，上述鬼影区域幅值均因。
34、栅瓣的抑制而显著减小，同 说 明 书CN 104008287 A 10 7/7页 11 时旁瓣的干扰也得到有效抑制，声源定位也更加准确清晰。这是由于PSO-MVDR不仅对栅瓣 有抑制作用，使其退化为旁瓣；同时该算法具有较强的空域滤波特性，在抑制栅瓣的同时使 阵列在不同方位的输出值免受旁瓣干扰，能够尽可能无失真地输出，因此可以获得较好的 鬼影抑制效果。 0078 综上所述，在全消音室的声场重构中，利用本文所提出PSO-MVDR的方法，从低频 到高频，均能够有效抑制鬼影，提高声场重构的准确性。 说 明 书CN 104008287 A 11 1/5页 12 图1 图2 说 明 书 附 图CN 104008287 A 12 2/5页 13 图3 图4 图5 说 明 书 附 图CN 104008287 A 13 3/5页 14 图6 图7 说 明 书 附 图CN 104008287 A 14 4/5页 15 图8 图9 说 明 书 附 图CN 104008287 A 15 5/5页 16 图10 说 明 书 附 图CN 104008287 A 16 。