【摘要】 本发明涉及一种基于随机传声器阵列的运动物体声场测量方法，属于噪声分析和控制技术领域。通过基于阵列单元的随机传声器阵列生成、根据基线向量统计筛选、仿真计算评价三个步骤设计二维随机传声器阵列布局，按照阵列布局搭建随机传声器阵列，对运动物体声场进行测量。本发明方法减少了传声器布局的任意性，降低布局设计中的空间；应用基线统计指标预先滤除不满足要求的阵列布局，减少了性能验证的计算量，提高了效率；通过仿真计算对不同频率下的随机传声器阵列性能进行验证，保证所得到的阵列满足测量使用要求。本发明方法使运动物体声场测量系统能够在较短时间内确定随机传声器阵列布局，实现对运动声源更准确的定位和分析。 【专利类型】发明授权 【申请人】清华大学 【申请人类型】学校 【申请人地址】100084 北京市海淀区清华园 【申请人地区】中国 【申请人城市】北京市 【申请人区县】海淀区 【申请号】CN200810238887.2 【申请日】2008-12-04 【申请年份】2008 【公开公告号】CN101413824B 【公开公告日】2010-04-07 【公开公告年份】2010 【授权公告号】CN101413824B 【授权公告日】2010-04-07 【授权公告年份】2010.0 【IPC分类号】G01H17/00 【发明人】杨殿阁; 连小珉; 郑四发; 罗禹贡; 李克强; 王建强; 李兵; 邵林 【主权项内容】1.一种基于随机传声器阵列的运动物体声场测量方法，其特征在于该方法包括以下步骤： (1)设定测量用传声器之间的最短间隔距离为d， 其中c为声速，fmax为设定的最大测量分析频率； (2)设上述测量用传声器所在的阵列平面A的长度为L1、高度为H1，待测运动物体所在平面R的长度为L2、高度为H2，阵列平面A与平面R之间的垂直距离为D，阵列平面A的长度L1和高度H1是上述最短间隔距离d的整数倍； (3)将上述阵列平面A划分为n个大小形状相同的阵列单元，每个阵列单元的长度和宽度分别为k1×d、k2×d，其中k1、k2为正整数，d为上述最短间隔距离，将每个阵列单元划分为k1×k2个方格，每个方格的中心点为测量用传声器的预设置点，共有k1×k2预设置点； (4)在上述k1×k2预设置点中，随机选定一个预设置点作为传声器放置位置，得到一个由n个传声器构成的随机传声器阵列，重复该过程N1次，生成N1个随机传声器阵列； (5)计算上述N1个随机传声器阵列中每个随机传声器阵列的所有基线向量，将同一随机传声器阵列的所有基线向量放到同一点距坐标系下，得到基线向量的点距坐标图，所述的基线向量是从一个传声器位置指向另外一个传声器位置的有向线段，设定点距坐标图中坐标轴的单位为最短间隔距离d，形成U个由直线x＝±J和y＝±J所构成的正方形坐标环，其中J＝1，2，…U-1，U，U＝max(H1/d，L1/d)，计算基线向量终点在第1坐标环上的所有基线向量数M1，基线向量终点在第1坐标环上且互相不重合的所有基线向量总数K1，基线向量终点在第2坐标环上且互相不重合的所有基线向量总数K2，以及基线向量终点在第U坐标环上且互相不重合的所有基线向量总数KU； (6)根据上述基线向量总数，对上述N1个随机传声器阵列进行筛选，其过程为： (a)若上述阵列单元数n≤24，且上述基线向量数M1≤4，则进行步骤(b)，若上述阵列单元数n＞24，且上述基线向量总数K1≤4，则进行步骤(b)； (b)若上述阵列单元数n≥14，且上述基线向量数K2≥14，则进行步骤(c)，若上述阵列单元数10≤n≤13，且上述基线向量数K2≥12，则进行步骤(c)，若上述阵列单元数8≤n≤9，且上述基线向量数K2≥10，则进行步骤(c)； (c)若上述基线向量数KU≥4，则选择该随机传声器阵列； (7)设上述选择的随机传声器阵列数有N2个，分别对N2个随机传声器阵列进行仿真计算，得到随机传声器阵列中各个传声器的仿真声压信号 其中第i个传声器的仿真声压信号为： 其中， 为仿真被测声源特征函数 为仿真实测声源强度，f为仿真被测声源的频率， 为在t时刻Q点与第i个传声器间的距离，θi(t)为在t时刻Q点到第i个传声器连线与Q点运动方向间的夹角，M为马赫数； 根据上述传声器的仿真声压信号，计算出待测运动物体所在平面R上任意一点 处在t1至t2时间段内的仿真声场特征函数： 其中， (t，ε，η)为t时刻面上任意点 (ε，η)与第i个传声器之间的物理距离，遍历待测运动物体所在平面R，得到仿真声场特征函数分布图，设分布图中主瓣峰值为hp，最大旁瓣峰值为hv，则旁瓣抑制比γf为： 若γf≥n且上述阵列单元数8≤n＜16，或者γf≥16且n≥16，则该随机传声器阵列满足频率f下的筛选条件，进行下一频率的筛选； 在下一频率重复上述仿真步骤，若在对每一频率都满足上述筛选条件，则得到测量用的随机传声器阵列，否则对上述选择的N2个随机传声器阵列中的下一个随机传声器阵列重复上述筛选过程，直至得到一个满足上述筛选条件的随机传声器阵列； (8)根据上述筛选得到的随机传声器阵列组成测量系统，采集运动物体发出的声压信号，计算得到待测运动物体所在平面R上任意一点s(ε，η)处在测量时间t1至t2内的声场特征函数： 其中， 其中，pi(t)为t时刻第i个传声器接收到的信号声压，ri(t，ε，η)为t时刻，任意点s(ε，η)与第i个传声器之间的距离。 【当前权利人】清华大学 【当前专利权人地址】北京市海淀区清华园 【专利权人类型】公立 【统一社会信用代码】12100000400000624D 【引证次数】4.0 【自引次数】3.0 【他引次数】1.0 【家族引证次数】4.0 【家族被引证次数】16
