平面阵列识别声源三维坐标的方法技术

本发明专利技术提供的是一种采用平面阵列识别定位声源三维坐标的方法。基于球面波声场模型，采用波束形成方法，以三维空间网格化后的网格交点作为重构点，进行不同距离的平面声场重构，比较沿距离方向（即z方向）上不同重构面的波束形成器的输出值，最大输出值对应的z轴坐标，即确定为声源的距离。最后根据声源距离对应的重构面上波束形成器输出的最大值点，确定声源在x和y轴坐标。本发明专利技术的优点在于可用二维的平面阵列，识别声源的三维坐标，可直接应用于现有的仅能识别声源两维坐标的平面阵列声全息设备上，实验已验证了该发明专利技术方案的有效性和可靠性。

【技术实现步骤摘要】
本专利技术属于结构声学、气动声学、水声学、信号处理、图像处理、噪声一振动控制等交叉领域，具体的说是一种采用平面阵列识别声源三维坐标的波束形成方法。
技术介绍
准确识别定位飞机、高铁、船舶及汽车等复杂系统的噪声源是实现对其有效降噪的前提，近年来，利用传感器阵列对噪声源精确地识别定位是一个重要的技术。目前常用的声源识别定位方法有近场声全息方法和波束形成的方法，近场声全息方法对于中低频率的声源，特别是结构噪声源具有较高的空间分辨率。然而，当声源为高频，且不能在声源近距离进行全息测量时，波束形成方法将成为一种有效的声源识别定位方法。本专利技术应用于声源的识别定位，现有波束形成方法和技术的局限性有(1)需要人工提供传感器阵列与声场搜索面之间的距离，才能确定声源在X和y方向的位置，由于提供的距离是人为估计的结果，无法准确定位事先未知声源在三维空间中的真实位置；(2)无法自动给出复杂对象及复杂结构表面的声源的分布状况。
技术实现思路
本专利技术的目的就是在现有的2维波束形成的基础上，智能地确定声源在第三维z·方向上的位置。本专利技术提出的识别声源三维坐标的波束形成方法，可应用于任何传声器分布的、任何外轮廓形状的平面阵列，由阵列面构成测量面，在测量面节点上放置传声器接收信号。将整个三维声场空间网格化，从某一初始的设定距离开始，由平行于测量面的、由X和y方向构成的声场重构面，根据各个传声器相对于重构点的空间位置对接收到的信号进行时延补偿后，加权求和输出，依次进行不同距离的声场重构，最终获得整个三维空间的波束形成器的输出，输出响应的最大值所对应的空间位置即声源坐标。下面对本专利技术方法做进一步说明。设定以M个传声器组成二维平面阵列，各传声器之间的间距满足λ/2准则，λ为声波的波长。以平面阵列轴心为坐标原点，轴心附近的某一传声器为参考传声器，其坐标为(x0, I0, Z0)，按如下步骤完成声源三维坐标的识别(I)设定初始声场重构面距平面阵列的距离Z=Ztl，截至距离ζΝ以及ζ方向的搜索步长Δζ=^^，其中N为整数；(2)将X和y坐标确定的声场重构面在感兴趣的矩形区间内进行声场网格化，假设共产生了 PXQ个节点；(3)将网格的节点即重构点作为假定声源的空间位置(X，y, z)，确定声波传播到不同传声器时产生的时延。重构点距参考传声器的距离为Ir I，其表达式为|r| = ^/(x-.r0)2+ (v — V0 )2+(z-z0)2第m个传感器的坐标为0^7111，2 1)，重构点距第111个传声器的距离为I rm |，其表达式为_4]KJ=V(x- xJ2+(y -凡)2+—)第m个传声器相对于参考传声器的时延表示为τ (rm)，其表达式如下r(r、=」(x-xm)2+(y-ym) + ~zm)2 ^λ](χ~χο)~+(y~}o) +(z~z0fq) mC其中c为介质中的声速；(4)波束形成的权向量为·c(x,>，,)= [l ... e—抑冲... e-^Mf(3)其中ω为声源的角频率，将整个搜索平面内所有节点的波束形成的权向量组成3维矩阵 C(W) ··· C(W) ... C(Xl5^g3Z) · * ■ · ■ * * · ·C= C(xp,}\,) ··· C(xp，},9,二)…c(xp,yg,z)(4)_c(w)…c{xP,yq,z)…c(xP,yQ,=)_其中c(xp，yq, z)为节点空间坐标为(V yq, z)所对应的波束形成的权向量。(5)将平面阵列上的传声器接收到的信号，按各个通道进行快速傅立叶变换，选取感兴趣的频点，以公式(4)确定的波束形成的权向量进行加权求和输出，获得声场重构面的波束形成输出，表示为一个PXQ维的矩阵Bz(6)设定下一个搜索面距平面阵列的距离ζ=ζ+Λζ，回到步骤(3)，获得相应的声场重构面上波束形成的输出Βζ+Λζ ；直到ζ = ζΝ，循环结束，共获得了 N个搜索面上的PXQ维的波束形成输出响应矩阵；(7)找出N个二维波束形成输出响应矩阵的最大值，并记录相应的三维坐标；(8)再从N个最大值中找出整个三维空间波束形成输出的最大值，对应的坐标(xs, ys, zs)即声源三维坐标。附图说明图I为60通道轮辐阵列图示。图2为本专利技术形成的四维功率谱图示。图3为本专利技术在X方向的定位效果图。图4为本专利技术在Y方向的定位效果图。图5为本专利技术在Z方向的定位效果图。具体实施方案下面通过具体实施例子对本专利技术作进一步地描述。本实施例子中，选择轮辐阵列作为测量阵，以图I所示的不规律方式布置60个传声器，传声器之间的间距不等，工作频段为O. 5 8kHz。设声源为点声源，位于(-0.5m，O. 5m，O. 8m)处，暂不考虑背景噪声，声源频率为4kHz，声源强度为IPa。感兴趣的搜索区域为长方体，其中心位于坐标原点0，Χ、Υ轴方向各从-2m搜索到2m, Z轴方向从Om搜索到3m,搜索步长均为O. 02m。平面阵列位于x、y坐标平面内，阵列面的中心位于坐标原点。根据前面设置的仿真参数，通过Matlab软件仿真可以得到所有B (p，q, η)的值在网格化的三维空间散点图，但是为了更清楚的显示识别定位结果，首先找到B(p，q，η)的最大值所在的空间位置，再选取与此空间位置共面同时又分别垂直于X、I、ζ轴的所有散点，最后经过筛选后的散点图如图2所示。由图2可知，本专利技术可以同时定位声源的X、Y、Z三维坐标，而且在感兴趣的区域内，旁瓣较低。图3 5是图2在各方向的剖面图。 图3显示当Y和Z轴固定在声源位置，X方向的波束模式图，由图可知，X轴方向的3dB主波束宽度为O. lm，第一旁瓣有-13dB。图4显示当X和Z轴固定在声源位置，Y方向的波束模式图，由图可知，Y轴方向的3dB主波束宽度为O. lm，第一旁瓣有-15dB。图5显示当X和Y轴固定在声源位置，Z方向的波束模式图，由图可知，Z轴方向的3dB主波束宽度为O. 2m，第一旁瓣有-18dB。由图3 5可知，本专利技术在X、Y和Z方向都有较高的分辨率和较好的抗干扰能力。本文档来自技高网...

【技术保护点】
一种采用平面阵列识别声源三维坐标的方法，其特征在于：按如下步骤进行：（1）设定初始声场重构面距平面阵列的距离z=z0，截至距离zN以及z方向的搜索步长其中N为整数；（2）将x和y坐标确定的声场重构面在感兴趣的矩形区间内进行声场网格化，假设共产生了P×Q个节点；（3）将网格的节点即重构点作为假定声源的空间位置(x,y,z)，确定声波传播到不同传声器时产生的时延。重构点距参考传声器的距离为|r|，其表达式为：|r|=(x-x0)2+(y-y0)2+(z-z0)2第m个传感器的坐标为(xm,ym,zm)，重构点距第m个传声器的距离为|rm|，其表达式为：|rm|=(x-xm)2+(y-ym)2+(z-zm)2第m个传声器相对于参考传声器的时延表示为τ(rm)，其表达式如下：τ(rm)=(x-xm)2+(y-ym)2+(z-zm)2-(x-x0)2+(y-y0)2+(z-z0)2c---(2)其中c为介质中的声速；（4）波束形成的权向量为c(x,y,z)=1...e-jωτ(rm)...e-jωτ(rm)T---(3)其中ω为声源的角频率，将整个搜索平面内所有节点的波束形成的权向量组成3维矩阵：C=c(x1,y1,z)...c(x1,yq,z)...c(x1,yQ,z)...............c(xp,y1,z)...c(xp,yq,z)...c(xp,yQ,z)...............c(xP,y1,z)...c(xP,yq,z)...c(xP,yQ,z)---(4)其中c(xp,yq,z)为节点空间坐标为(xp,yq,z)所对应的波束形成的权向量。（5）将平面阵列上的传声器接收到的信号，按各个通道进行快速傅立叶变换，选取感兴趣的频点，以公式（4）确定的波束形成的权向量进行加权求和输出，获得声场重构面的波束形成输出，表示为一个P×Q维的矩阵Bz（6）设定下一个搜索面距平面阵列的距离z=z+Δz，回到步骤（3），获得相应的声场重构面上波束形成的输出Bz+Δz；直到z=zN，循环结束，共获得了N个搜索面上的P×Q维的波束形成输出响应矩阵；（7）找出N个二维波束形成输出响应矩阵的最大值，并记录相应的三维坐标；（8）再从N个最大值中找出整个三维空间波束形成输出的最大值，对应的坐标(xs,ys,zs)即声源三维坐标。FDA00002165934300011.jpg...

【技术特征摘要】
1.一种采用平面阵列识别声源三维坐标的方法，其特征在于按如下步骤进行 (O...

【专利技术属性】
技术研发人员：卢奂采，李春晓，金江明，丁浩，梅东挺，
申请(专利权)人：浙江工业大学，
类型：发明
国别省市：