基于混合引力搜索算法的立体传声器阵列优化设计方法
【专利说明】基于混合引力搜索算法的立体传声器阵列优化设计方法 【技术领域】
[0001] 本发明属于传声器技术领域，特别涉及一种立体传声器阵列优化设计方法。 【【背景技术】】
[0002] 实际情况下，变电站的噪声源较复杂，平面阵的定位结果无法分辨出前后波达方 向，而体积阵能够实现对全空间的扫描，抑制后向来波。体积阵性能的优劣很大程度上受阵 型结构，阵元数量以及阵列孔径的影响。在相同的阵元数量和阵列孔径的前提下，不同形式 阵列结构的波束方向图的主瓣宽度和最大旁瓣级不一样，产生的波束响应也不同，对噪声 源识别定位的性能影响很大。主瓣宽度和旁瓣级，是相互矛盾的性能指标，无法得到同时具 有最窄主瓣宽度和最低旁瓣级的最优阵形，因此需要综合考虑，取得二者的折中。如何结合 声源特征，优化设计体积阵列结构是多传感器阵列声成像定位的一个难点，具有重要的意 义。
[0003] 就智能优化算法来说，目前应用比较广泛的有遗传算法、蚁群算法、神经网络、粒 子群算法。这些算法大多启蒙于自然现象或自然规律，各具优缺点，适用面也不同。遗传算 法对参数的依赖性较强，和粒子群算法一样容易早熟，陷入局部最优；蚁群算法和模拟退火 虽然不受初始条件的约束，但计算量大，优化过程较长。 【
【发明内容】
】
[0004] 本发明的目的在于提出一种基于混合引力搜索算法的立体传声器阵列优化方法， 以克服现有技术所存在的问题；本发明方法优化设计的立体传声器阵列，能弥补传统平面 阵定位结果无法分辨出前后波达方向的不足，抑制后向来波，实现对复杂声环境下目标噪 声源的被动定位。
[0005] 为了实现上述目的，本发明采用如下技术方案：
[0006] 基于混合引力搜索算法的立体传声器阵列优化设计方法，包括以下步骤：
[0007] 1)、构建立体传声器阵列优化参数模型：
[0008] 待优化多臂立体传声器阵列有N个阵臂，类似伞骨架呈角度识张开，每个臂上有M 个麦克风，共计N*M个麦克风；阵臂之间的夹角α = 360° /N，固定每个阵臂第一个阵元到 阵中心的距离I1= 0. 2m，最后一个阵元到阵中心的距离~ =1/sin% 只考虑单个臂上 阵元的距离L = [12, ...，1M J的变化，其余阵臂依次复制；与此同时，为防止阵元重叠，限 制每条臂上相邻两个阵元的最小间距为lc，即
[0009] lm+1-lm> I 0，m = 1，2,
[0010] 2)、设定引力系统粒子形式及其优化目标：
[0011] 假设一个包含有N。个粒子的引力系统，系统维度为D，这个引力系统的所有粒 子用向量表示为I =[心A经过步骤1)处理之后，引力系统粒子的形式设置为 JC - [/0 ? · * M s ?
[0012] 设频率为f。的声信号所对应的常规波束旁瓣级指标为SLL。，要求旁瓣级小于该指 标，若阵列结构为P，则适应度函数Fit (P，f。）可以表示为
[0014] 当旁瓣级满足指标要求时，则以主瓣宽度为适应度函数返回值；若旁瓣级不满足 指标要求，则返回给适应度函数的值为无穷大；因此在一定旁瓣级的要求下，优化目标函数 为
[0015] min Fit (P, f〇), s. t. SLL (P, f〇) ^ SLL0
[0016] 3)、具体的迭代优化过程：
[0017] 将粒子的位置用向量表示为
；D为引 力系统的维度；根据牛顿万有引力定律：任意两个质点在连心线方向上有力的相互吸引， 设第i个粒子受到第j个粒子的引力
则第d维方 向上的引力分量
[0019] 其中，M1UhMj (t)是第i个粒子和第j个粒子被动引力质量；;的、考(〇是第j 个粒子和第i个粒子在d维方向上的坐标，I I X1 (t)，X] (t) I 12是这两个粒子的欧式距离；ε 是一个很小的常数，一般设为〇 ;G(t)是随迭代时间递减的万有引力常系数，且
[0021] G。是引力常系数的初始值，k是一个递减的系数，令G。= l，k = 20 ;t为当前迭代 次数，T为设定的最大迭代次数；
[0022] 给任意两个粒子在d维方向上的引力分量加一个在（0, 1)范围内变化的随机权值 randO，则第i个粒子在d维方向上受到的其他粒子的引力随机加权总和
[0024] 这一迭代时刻，粒子在d维方向的引力加速度
[0026] 其中，M11⑴是第i个粒子的惯性质量；粒子被动引力质量和惯性质量的计算依赖 于适应度函数Fit返回的函数值；假设被动引力质量、惯性质量和个体质量三者相等，记为 M1, M1根据以下几个公式进行计算：
[0029] 其中 Fitbest= min(Fit Jt))，i e [1，N]，Fitworst= max(Fit Jt))，i e [1，N]; Fit1U)是第i个粒子在t时刻的适应度函数值。
[0030] 混合引力搜索算法的粒子速度和位置按以下两个公式进行更新：
[0033] 权重因子c3、(：4取0. 5和1. 5 ;惯性权值w取（0, 1)内的随机数；以.a(i)是指当前 迭代时刻粒子在d维方向上的最优解；当循环达到最大次数T时，终止迭代。
[0036] 阵列投影直径为2m，阵列中心高度为I. 5m到2. 5m可调，阵列中心配置高分辨率摄 像头1个，阵列工作时张开并锁定，不工作时能够折叠存放。
[0037] 步骤2)中粒子个数为25个，系统维度D = 8。
[0038] 基于混合引力搜索算法的立体传声器阵列优化设计方法，具体包括以下步骤：
[0039] 1)、目标声源的频率为1000Hz，要求优化后阵列的旁瓣级不高于-8dB，优化目标 为满足旁瓣级的情况下获得最窄的主瓣宽度；阵列是投影直径为2米的7臂阵，固定每个阵 臂第一个阵元到阵中心的距离I 1= 〇.2m，每个臂阵元数为9,最后一个阵元到阵中心的距 离/9 二 l/(sinp:h
[0040] 2)、设定引力系统粒子形式及其优化目标：
[0041] 将体积阵的参数设置为引力系统粒子，粒子个数为25个，系统维度D = 8 ;这个引 力系统的所有粒子可以用向量表示为X = [X1, X2，…，x25];设频率为1000Hz的声信号所对 应的常规波束旁瓣级指标为_8dB，要求旁瓣级小于该指标；当旁瓣级满足指标要求时，则 以主瓣宽度为适应度函数返回值；若旁瓣级不满足指标要求，则返回给适应度函数的值为 无穷大；因此在一定旁瓣级的要求下，优化目标函数为
[0042] min Fit (P，1000)，s. t. SLL (P，1000)彡-8
[0043] 3-1)、初始化引力系统粒子；
[0044] 3-2)、通过初始粒子计算阵列的波束指标，进而得到适应度函数值；
[0045] 3-3)、根据适应度函数值计算粒子的引力质量和惯性质量；计算粒子引力，引力加 速度，按照粒子的速度和位置进行更新；
[0046] 3-4)、更新后的粒子按照步骤3-2)得到适应度函数值，挑选适应度函数值最小 的引力子为当前最优参数，若满足优化目标函数，迭代停止；若不满足，则重复步骤3-3)、 3-4);若达到最大迭代次数2000,目标函数仍不满足，则终止迭代，更换初始化粒子，重复 步骤 3-2)、3-3)、3-4)。
[0047] 相对于现有技术，本发明具有以下有益效果：
[0048] 本发明结合了粒子群算法的混合引力搜索算法既改善了传统引力搜索算法在求 解复杂多峰问题时容易陷入局部最优的现象，又比原始的粒子群算法具有更快的迭代收敛 速度，大大提高了阵列优化的有效性；所设计的立体传声器阵列成功抑制了后向噪声的干 扰，实现了对复杂声环境下目标噪声源精确定位的要求。 【【附图说明】】
[0049] 图1 :待优化的7臂体积阵；
[0050] 图2 :优化后的7臂体积阵；
[0051] 图3 :优化前后的定位结果比较；图3(a)为优化前随机产生的阵列参数得到的定 位结果示意图；图3(b)为优化后的定位结果示意图。
[0052] 图4为本发明方法的流程图。 【【具体实施方式】】
[0053] -、伞形阵模型
[0054] 假设一个近场常规波束模型，将M元伞形阵垂直放置在空间中，建立三维直角坐 标系，各个阵元的三维空间坐标可以表示为，且P 1 (Xl，yi，Z1)，i = 1，2,…M，且
[0056] 其中，I1是第i个阵元到基阵原点的距离，朽是第i个阵元所在阵臂与z轴的夹 角，S 1是该阵臂在X〇y平面的投影与X轴的夹角。阵列结构可以用矩阵表示为
[00