毕业论文基于麦克风阵列的声源定位技术43655

《毕业论文基于麦克风阵列的声源定位技术43655》由会员分享，可在线阅读，更多相关《毕业论文基于麦克风阵列的声源定位技术43655（44页珍藏版）》请在装配图网上搜索。

1、2012届毕业设计说明书毕业设计说明书基于麦克风阵列的声源定位技术学生姓名： 学号： 学 院： 专 业： 指导教师： 2012年 6 月基于麦克风阵列的声源定位技术摘 要声源定位技术是利用麦克风拾取语音信号，并用数字信号处理技术对其进行分析和处理，继而确定和跟踪声源的空间位置。声源定位技术在视频会议、语音识别和说话人识别、目标定位和助听装置等领域有着重要的应用。传统的单个麦克风的拾音范围很有限，拾取信号的质量不高，继而提出了用麦克风阵列进行语音处理的方法，它可以以电子瞄准的方式对准声源而不需要人为的移动麦克风，弥补单个麦克风在噪声处理和声源定位等方面的不足，麦克风阵列还具有去噪、声源定位和跟踪

2、等功能，从而大大提高语音信号处理质量。本文主要对基于多麦克风阵列的声源定位技术领域中的基于时延的定位理论进行了研究，在此基础上研究了四元阵列、五元阵列以及多元阵列的定位算法，并且分别对其定位精度进行了分析，推导出了影响四元、五元阵列目标方位角、俯仰角及目标距离的定位精度的一些因素及相关定位方程，并通过matlab仿真软件对其定位精度进行了仿真；最后在四元、五元阵列的基础上，采用最小二乘法对多元阵列定位进行了计算；通过目标计算值和设定值对比，对多元阵列的定位精度进行了分析，并得出了多元阵列的目标定位的均方根误差。关键词：麦克风阵列，声源定位，时延，定位精度，均方根误差Based on Micro

3、phone Array for Sound Source Localization ResearchAbstractSound source positioning technology is to use the microphone to pick up voice signals, and digital signal processing technology used for their analysis and processing , Then identify and track the spatial location of sound source. Acoustic so

4、urce localization techniques have a variety of important uses in videoconferencing, speech recognition and speaker identification, targets direction finding, and biomedical devices for the hearing impaired. The pick up range of traditional single microphone is limited, the signal quality picked up i

5、s not high, then a voice processing methods with the microphone array has been proposed . It may be electronically aimed to provide a high-quality signal from desired source localization and does not require physical movement to alter these microphones direction of reception. Microphone array has th

6、e functions of de-noising, sound source localization and tracking functions, which greatly improved the quality of voice signal processing. The article discusses some issues of sound source localization based on microphone array, On the basis ，it studies a four element array,five element array and a

7、n multiple array positioning algorithm, then the positioning precision is analyzed. Derived some factors of the azimuth and elevation angle targets the target range of the estimation precision affected and positioning equation. And through MATLAB simulation software for its positioning accuracy of s

8、imulation. finally ，based on four yuan, five yuan of array, using the least square method ,the multiple array localization were calculated. Through the contrast of the target value and set value, multiple array positioning accuracy is analyzed, and the of diverse array target positioning.Keywords: M

9、icrophone Array, Sound Source Localization, Time Delay, Positioning precision , root mean square error目 录1 引言11.1 研究背景和意义11.2 声源定位技术的研究现状及发展 11.2.1 研究历史和现状11.2.2 发展趋势21.3 麦克风声源定位技术31.4 本文所要研究的内容42 声学理论基础知识52.1 空气的物理特性52.2 声波的物理特性52.3 声波在空气中的传播特性82.4 声源定位原理92.5 本章小结103 基于时延估计的声源定位算法及其精度分析113.1 时延估计算法

10、概述113.2 基于时延估计的声源定位的研究分析113.2.1 四元阵列定位算法123.2.2 五元阵列定位算法153.3 四元阵列的定位精度分析及其仿真153.3.1 方位角精度分析及仿真173.3.2 俯仰角精度分析及仿真193.3.3 距离估计精度分析及仿真213.4 五元阵列的定位精度分析及其仿真233.4.1 方位角精度分析及仿真233.4.2 俯仰角精度分析及仿真253.4.3 距离估计精度分析及仿真263.5 本章小结284 多元麦克风阵列声源定位分析294.1 多元麦克风阵列定位方程294.2 最小二乘法求声源位置304.3 定位精度分析314.4 结果及计算分析324.4 本

11、章小结335 总结与展望345.1 全文总结345.2 本文的不足之处及后续工作展望34参考文献36致 谢38第 II 页 共 页1 引言1.1 研究背景和意义在各种电子设备高度智能化的今天，语音增强与声源定位技术成为语音通信领域中两种不可缺少的技术。例如，在视频会议中，通过声源定位技术控制摄像头，使其自动转向感兴趣的说话人方向；对于高速行驶的车辆，为避免驾驶员用手去接听电话，车载免提电话应运而生1。然而，当车中坐有多个说话人时，语音识别系统就无法辨别实际命令的来源，此时就需要一种定位系统来提取驾驶员方位的语音，进而对其命令作出响应；助听器的出现为有听觉障碍的残疾人提供了帮助，基于阵列的语音增

12、强技术利用声源的位置信息进行空间滤波，可以进一步抑制除说话人以外的其他方向的噪声，使得助听器话音更加清晰。声源的定位在各个领域都有着广泛的应用，尤其是利用传声器阵列的方法，在电话会议、视频会议、可视电话等系统中控制摄像头和传声器阵列波束方向对准正在说话的人，发挥了强大的功能。早在20 世纪七八十年代，声源定位系统就开始被广泛地研究。用于声源定位的传声器阵列具有很强的空间选择性，不需要移动传声器就可以获得移动目标的声音信号，并且能够在一定范围内实现声源的自适应检测、定位及跟踪。这种声源定位技术主要可以分成3 类:第一类是基于最大输出功率的可控波束形成技术，当传声器阵列探测到信号时，对各路信号进行

13、加权求和形成波束，直到得到具有最大输出功率的波束为止；第二类是高分辨率谱估计技术，这类技术需要利用传声器阵列所获取的信号计算空间谱的相关矩阵，运算量大，定位不准确，在实际操作过程中很少采用；第三类是基于声达时间差(time difference of arrival , TDOA) 的定位技术，利用到达阵列上各传声器的声音信号间的时间差未定位声源，这类技术计算量小，比较适用于实时处理，在实际应用中占有很大的比重2。1.2 声源定位技术的研究现状及发展1.2.1 研究历史和现状近年来，传感器和探测技术、微电子技术、信号处理技术以及人工智能技术的飞速发展，均为声探测技术用于直升机等军事目标的定位、

14、跟踪和识别幵辟了新的应用前景，使声探测技术成为一种重要的军事侦察手段和防空作战中反电子干扰和反低空突防的一种有效途径。声探测技术就是利用直升机飞行及坦克行进时的噪声和振动信号，实现对目标的自动探测、跟踪和识别。美英等国研制的声智能雷弹系统，就是采用了这一思想来对付超低空飞行的直升机。采用声复合引信的反直升机、反坦克智能雷弹系统具有独立自主的作战能力。单枚雷弹的杀伤半径将超过100km，与普通雷弹系统相比，极大地提高了武器系统的有效性，如法国的MAZAC反坦克雷弹，单枚雷弹的效能相当于60枚普通地雷的综合效能3。因此，智能雷弹系统在构成雷场时所需的弹药数量最少，可以最大限度地减轻战时与和平时期的

15、后勤保障负担，具有显著的作战效益和经济效益。发展这种系统的关键在于，研制一种声或声复合引信，而声阵的布设方法及目标定位、跟踪方法的研究是声引信的核心。和平与发展是当今世界的两大主题，但世界并不太平，一些局部战争不断发生。历史遗留下来的边界和领土争端、南海诸岛的主权争议以及台湾问题等都对我国构成了不安定的因素。一些主要的邻国装备了相当数量性能优良的武装直升机和坦克，并在不断提高其武器装备的现代化水平，已对我国的国家安全构成了潜在的威胁。研究声源定位技术，加强我军反直升机、反坦克武器系统的装备水平，对我国的国防现代化建设具有重要意义。1.2.2 发展趋势声源定位技术是利用声学与电子装置接收目标声场

16、信息以确定目标声源位置的一种技术，能够用于探测武装直升机、坦克、火炮等军事目标的位置，是一种重要的军事侦察手段。其方法可分为被动定位系统和主动定位系统两大类。被动声定位系统属于辐射源无源定位，它具有作用距离远，隐蔽性好等优点。主动声定位系统容易暴露自己，被对方发现，影响系统自身的安全。进行利用声探测技术进行目标识别、跟踪和定位最早见于声呐，由于电磁波在海水中的传播衰减很大，作用距离受到限制，因此，声波被广泛用于水下目标的探测。地面声测产生于第一次世界大战。地面声测侦察在炮兵作战中曾立下汗马功劳，用于确定堑壕战中敌人武器的方位。据统计，在第二次世界大战和朝鲜战争中有75%的火炮侦察任务是利用声测

17、手段完成的。随着雷达、红外、激 光侦察技术的兴起，声测技术曾一度受到冷落。近年来，由于雷达面临着电子干扰、反辐射导弹、低空突防和隐身技术这四大威胁，越来越容易遭受攻击，因此，人们又开始重视被动式声探测系统，重新激起对声探测技术的兴趣4。反坦克等智能武器的研究幵发是必要的、紧迫的，这是对地面防空力量的重要补充。声测系统作为一种传统的侦察手段，近年来通过釆用新技术，提高了性能，满足了现代化战争的需要，其主要特点是：（1)不受通视条件的限制。可见光，激光和无线电侦察器材则需要通视目标，在侦察器和目标之间不能有遮蔽物，而声测系统可以侦察遮蔽物如山，树林等后面的声源。（2）隐蔽性强。声测系统不受电磁波干

18、扰也不会被无线电测向及定位，工作隐蔽性强。（3）不受能见度限制。其它侦察器材受环境气候影响较大。在恶劣气候条件下工作时，性能下降，甚至无法工作，声测系统可以在夜间，阴天，雾天和下雪天工作，具有全天候工作的特点。在几十年的发展过程中，声测系统的基本原理和设计变化不大，主要是通过釆用釆用各种现代化技术实现了自动化、集成电路化，提髙了性能。现有的声测 系统一般用热线，碳粒或电容传声器接收声音信号，用无线电技术或光纤技术传输信号,用电蚀纸带或墨水纸带记录测量结果，并把记录仪与计算机相连,用计算机处理测量结果。声测基线采用直线，弧线，正方形等多种形式。目前,声测系统的侦察距离可达30km，精度达到1%,

19、无线电传输方式的展幵时间仅需半小时，可同时测量50多个目标5。目前，这些国家研制的声探测系统都具有战场侦察、目标定位跟踪和敌我识别等综合作战能力，而且它们的研制已经进入了实质性的研制阶段，而我国在这方面的研究仍处于预研阶段，拟釆用的技术手段和途径仍处于可行性探索阶段，与其他国家相比还有很大的差距。1.3 麦克风声源定位技术麦克风阵列是获取高质量语音信号的重要途径。麦克风阵列语音采集系统能在多种环境下（混响、噪声、竞争声源、说话者运动）为后继的语音识别系统和编码系统等提供高质量的语音数据，而且还能定位声源位置和跟踪声源运动。准确的声源定位能引导麦克风阵列增强感兴趣的语音，同时抑制干扰语音和噪音。

20、声源定位技术是利用声学与电子装置接收目标声场信息以确定目标声源位置的一种技术。麦克风阵列的声源定位是指用麦克风阵列拾取声音信号，通过对多路声音信号进行分析与处理，在空间域中定曲一个或多个声源的平面或空间坐标，即得到声源的位置，以进一步控制摄像机和麦克风阵列波束对准正在说话的人。基于麦克风阵列的声源定位在视频会议、声控机器人、车载电话系统和语音导航系统、智能空间中的人机交互、声源的实时跟踪和监控等诸多领域具有广泛的应用前景。我们提出了有效的高分辨率多声源定位方法，适用于存在较强混响的环境。目前，基于麦克风阵列的声源定位方法大体上可分为三类：(1)基于子空间的定位技术；(2)基于可控波束形成的定位

21、技术；(3)基于到达时延(TDOA)的定位技术6。TDOA方法首先求出声音到达不同位置麦克风的时延，再利用这些时延求得声音到达不同位置麦克风的距离差，最后用搜索或几何知识确定声源位置。在现有的麦克风阵列声源定位方法中，基于到达时间差(TDOA)估计定位法计算量较小，硬件成本较低，定位精度较高，同时也易于实时实现，是目前声源定位法中常用的方法。1.4 本文所要研究的内容本文的研究内容是学习声学理论基础知识和被动声定位的原理，参与设计基于多麦克风的声源定位系统，在此基础上研究四元阵列、五元阵列以及多元阵的定位算法，分别对定位精度进行分析，采用matlab进行仿真，并指出影响定位精度的因素。具体研究

22、工作如下： 第一章概述了麦克风声源定位技术研究的背景和意义、历史发展以及国内外研究现状。第二章主要概述了声学理论的基础知识及声波的物理传播特性，并简述了声源定位原理，继而为后面研究声源定位技术打好基础。第三章论述了基于时延估计的声源定位算法，推导出了四元、五元十字阵的声源定位方程；并对四元和五元两种十字阵的定位精度进行了理论分析和仿真对比。 第四章在第三章的基础上，分析了多元麦克风阵列的声源定位及定位精度；并用实际数据对上述算法及麦克风阵列声源定位精度进行了计算求均方根误差。第五章对全文作了总结。2 声学理论基础知识2.1 空气的物理特性当声波在空气媒质中传播时，空气的物理特性，如密度、温度、

23、压强、比热 和粘滞性系数等，会直接影响到声波的物理特性。(1) 空气的密度空气密度可根据表达式（2.1）来计算 (2.1)式中 T绝对温度（K）； P大气压强（Pa）；(2) 声速对于理想气体而言，声波的传播速度可表示为 (2.2)式中 气体摩尔质量 比热比对空气来说，=1.402，=29千克/摩尔；R为气体常数，可表示为R=8.31焦耳/(开尔文摩尔)；因此，对于空气介质来说，声波在其中的传播速度可表示为： 或 (2.3)其中t表示摄氏温度（）。在常温下，空气中的声速大约为：340m/s。2.2 声波的物理特性当外力对介质的某一部分产生初始扰动时，这种扰动就将由一个质点传播到另一个质点，交替

24、形成密层和疏层；如此继续传播下去，就形成了波。当这种扰动表现为弹性变形时，就视其为以弹性波的形式传播信息。声音本质上是物体的机械振动对气体介质的一种扰动，这种扰动使气体介质的压强，以及密度、温度和速度等发生周期性的变化，忽强忽弱，疏密相间。这 种周期性的变化在各向同性的介质中，是借助介质的弹性向四面八方传播。由于声波的频率范围较广，从Hz到1012Hz，因此当那种周期性变化的频率在声频 范围2092201012内时，就将其称为是声音或声波：变化频率低于声频范围的弹性波被称为次声波；而变化频率髙于声频范围的弹性波则称为超声波。在声波探测技术中,习惯上将声波和超声波统称为声波。声波不仅可以在气体媒

25、质中传播，同时也可以在液体和固体中传播。当其在空气中传播时，只能发生压缩和膨胀，即空气质点的运动方向与波的传播方向一致，因此在空气中的传播的声波属于纵波。衡量声波的物理量主要有两个方面：一方面是声音强度的度量，即反应声的大小；另一方面是声频的高低度量。具体体现在：声压、声强、声功率以及声能密度等。(1) 声压声波对传播媒质作用时，使媒质质点受到挤压而产生压力变化，并发生周期性的压缩和膨胀，从而引起媒质中压强的变化。在空气中，由于声波扰动的影响,使空气压强发生起伏变化，而出现压强增量，该压强增量就称为声压，单位为帕 (Pa)。声压是用来表示声音强弱的物理量。通常釆用的声压有瞬时声压、峰值声压和有

26、效声压。瞬时声压是指媒质中某点瞬时压强和静压强的差值；某一时刻内最大的瞬时声压称为峰值声压；当瞬时声压对时间取方均根值时就称为有效声压。 通常所说的声压都是指有效声压。声压一般是时间和空间的函数，即： P=P(x,y,z,t) (2.4) P(t,x)= (2.5)式中 圆频率，= ； 波数，； 声压幅值；此声场是一个波阵面为平面，沿正x方向以速度传播的平面行波。(2) 声强和声功率声强是指在垂直于声波传播方向上，单位时间内通过单位面积的声能；其单位为：；而且声强越大，说明声音越强。在平面波自由声场中，声波仅来自声源方向，此时的声强可表示为： (2.6)式中 有效声压； 空气密度； 空气中声波

27、传播速度，即声速；在扩散声场中，声强可表示为： (2.7)通过上式，可知某点的声强和该点的声压平方成正比。其中：和的乘积也被称做是媒质的特性阻抗；对空气来说，20时的特性阻抗为：407瑞利（Pas/m）。声功率，是声源声输出的一种基本度量，也是声源本身的一种基本物理特性；它是指声源在单位时间内辐射出的总的声能量，单位为：W； 在自由声场中，声源中声功率与声强有如下关系： (2.8)式中 离开声源的距离； 声源辐射的声功率；(3) 声能密度声能密度是指在传播媒质中，媒质单位体积内所包含的声能。在平面波声场中，声能密度可表示为: (2.9)式中 有效声压； 空气度； 空气中声波传播速度，即声速；(

28、4) 声学量的级在声学中，直接使用声压、声强和声功率是极不方便的，为此，引入了 “级” 的概念，来表示声音的强弱。声学量的级就是指某个声学量与其同类基准值之比的对数。基准值是用来规定声学量级中的零分贝值。级的类别用名称表示有声压级、声强级和声功率级等7。通常所使用的级的单位是分贝，符号表示为dB。 声压级的数学表达式为： (2.10)声强级的数学表达式为： (2.11)声功率级的数学表达式为： (2.12)上述表达式中、分别表示基准声压、基准声强和基准声功率，其中和的值分别为： 和；2.3 声波在空气中的传播特性声波在空气中传播时，其形状可根据波阵面的形状来划分，主要有平面波和球面波。波阵面为

29、平面的波称为平面波，并且此波阵面具有与声传播方向的垂直平面相平行的特征；在不考虑介质吸收波能量的情况下，即在理想媒质中，声压不随着声源距离的变化而变化，即此时的声压为恒量。球面波则是波阵面为同心球面的声波；球面波在介质中传播时，球面波的声强与距声源距离的平方成反比，也就是说声压与传播距离成反比；在理想媒质中，声压与球面声波的半径成反比。声波在媒质中传播时，其振幅随传播距离的增大而减小的现象，就称为声波的衰减。声波在空气中传播时，声波的传播方向会由于大气温度梯度和传播速度梯度的存在而发生变化，同时大气的扰动会导致声波发生崎变，以及空气的粘滞性造成声能的吸收等，这些都造成了声波在空气中传播时发生衰

30、减，此衰减与空气的温度、湿度和声波的频率等有关8。除这些原因外，声波在传播媒质中的悬浮粒子上发生的散射，也是形成衰减的重要原因。对于声波衰减现象可用数学表达式表示为：当声波传播一段较小的距离办后，振幅的减小量与原来的振幅成正比，也与这段距离成正比，即 (2.13)当x=0时，因此有 (2.14)式中 媒质对声波的衰减系数，由两部分组成，即吸收衰减系数和散射衰减系数；其单位为：奈培.米-1在均匀介质中，对于被动声设备其接收点信号强度的衰减可表示为： (2.15)式中 、分别为接收点和发射点的信号强度；收发间距（或是带有接收换能器的声源装置与目标间的距离）；吸收系数； 由于介质往往是不均匀的，因此

31、对上式进行修正： (2.16)式中 空间衰减系数；异常值，即：在相同距离上，不均匀介质中信号强度与均匀介质中的信号强度之比；其不是距离的单调函数；2.4 声源定位原理声源定位技术是利用声学与电子装置接收目标声场信息以确定目标声源位置的一种技术。被动声目标定位就是指传声器工作在被动状态方式下，利用目标发射出的噪声信号对目标位置进行估测。如何布设性能优良、结构简单合理的传声器阵列是被动声目标定位和跟踪系统中的关键技术之一；通常传声器阵列布设的结构可以分为线型阵列、平面阵列和立体阵列等。对线型阵列来说，它只能对以阵列所在直线为界的半个平面进行目标定位，确定目标的二维参量，否则无唯一解；并且当目标位于

32、线型阵列的端射方向时，线型阵列将失去测距的能力；平面阵列不仅可以对整个平面进行目标定位，同时也可以对阵列所在平面为界的半个空间进行定位，确定目标的三维参量；立体阵列则可以对整个空间进行定位，但其定位算 法比较复杂。被动声目标定位算法按照测量依据来分主要有两种:一种是基于时延估计(信号到达不同传声器间的时间差)的目标参数估计测量法，而另一种则是基于瞬时频率估计的参数估计测量方法9；由于前者具有较高的测向、测距精度，并有较强的抗干扰性能，是目前被动声测系统中广泛釆用的方法。若按照传声器阵列的排布方式，有直线定位法、平面三角定位法、平面圆形定位法、平面正方形定位法和球面三角形定位法等。2.5 本章小

33、结本章主要概述了声学理论的基础知识及声波的物理传播特性，并简述了声源定位原理，列举了声源定位的两种方法及几种类型，继而为后面研究时延估计的声源定位技术打好基础。3 基于时延估计的声源定位算法及其精度分析3.1 时延估计算法概述所谓时延，就是指具有同一计时起点的两传声器所接收到的信号到达时刻差，即两信号间的时间延迟。由此可见，由N个传声器组成的阵列可以获得N-1个相互独立的时间延迟值；根据目标位置的自由度，可知要确定平面目标位置至少需要3个传声器，即2个时延值；而确定空间中的目标位置则至少需要4个传声器， 即3个时延值；但由于受到区域条件、阵列设计方法等的限制，传声器阵列中传声器元件的个数不宜过

34、多。利用声信号到达传声器间的时间差(时延值)进行目标定位的基本原理见表达式所示：式中 一代表目标信号的位置坐标；，代表阵列中传声器和传声器的位置坐标； 声波在媒质中的传播速度；表示声波传播到传声器和传声器之间的时间差，即时延值。3.2 基于时延估计的声源定位的研究分析声传感器阵列的设计在直升机被动声定位中具有十分重要的意义。采用的传声器阵列可分为线阵、面阵和立体阵。对于固定式阵列，线阵只能对以阵列所在直线为界的半个平面进行定位，否则解不唯一。立体阵可以对整个空间进行定位，但其算法要复杂的多。面阵可以在整个平面对目标进行定位，也可以对阵列所在平面为界的半个空间进行定位。由于定位系统布设于地面，目

35、标为低空或超低空飞行的武装直升机，因此采用面阵是可行的。由N个声传感器阵元组成的阵列，可以得到N-1个独立的时延，空中的直升机对于被动声定位系统来说可以看成点目标，有三个自由度。因此，要对目标进行定位，至少需要四个阵元组成的声传感器阵列。由于十字形阵具有分维特性，且阵列冗余度较小，因此，十字阵是较为合适的阵形10。本章推导了四元十字阵的目标定位方程以及五元十字阵的定位算法。3.2.1 四元阵列定位算法如图3.1所示，十字形声传感器阵列，由两个相互正交的线阵S1、S3和S2、S4组成。线阵的阵元间距为D，以两线阵的交点为坐标原点，建立如图3.1所示的直角坐标系和球坐标系。四阵元的直角坐标分别为：

36、S1(D/2,0,0), S3(-D/2,0,0), S2(0,D/2,0), S4(0,-D/2,0),设目标声源T的直角坐标为(x,y,z),球坐标为(r, ,)。即目标T到坐标原点的距离为r,方位角为，俯仰角为。图3.1 四元十字阵定位示意图当目标声源T离阵中心的距离比阵元间距大的多时，可以假设目标T为点声源，并以球面波形式进行传播。设声源到达阵元S1的传播时间为t1，到达阵元S2,S3,S4与相对于到达阵元S1的时间延迟(以下简称时延)分别为, , ,则声源传播到S2、S3、S4与传播到S1的声程差分别为， ，： (3.1)式中，C为空气速度，C=340m/s.设目标声源T与阵元S1的

37、距离为r1,则： (3.2)目标声源T以球面波进行传播，所以阵元S1、S2、S3、S4分别位于以T为球心以r1、r1十、r1+ ,r1+ 为半径的四个球面上,因而可联立列出方程组如下： (3.3)测量时延、，即可知程差、，再通过解方程组(3.3), 可得目标的位置坐标（x,y,z）。下面推导目标的球坐标(r, ,)与时延、和的关系。将3、4、5式分别与2式相减，得： (3.4)如图3.1所示，直角坐标系下目标声源T的位置坐标(x,y,z)与球坐标系下的位置坐标(r, ,)的关系式为： (3.5)式中， ,。由（3.5）式，可解得由直角坐标（x,y,z）表示的方位角，俯仰角以及目标声源T距阵中心

38、的距离r。 由于，=2,3,4.所以，上式可以简化为： (3.6) (3.7)2式减去1式，得： (3.8) 式中，。由于，釆用的是小尺寸的声学阵列，因此，上式可以简化为： (3.9) (3.10)由(3.8)式得： (3.11) 因此，通过估计时延、和，即可确定目标声源T的位置。3.2.2 五元阵列定位算法五阵元十字形阵列是在四元阵的中心增加一个声传感器构成的。建立如图3.2 所示的直角坐标系，五阵元所在位置的直角坐标分别为：,目标声源T的直角坐标为(x,y,z)。 目标T到坐标原点的距离为r，方位角为，俯仰角为。图3.2 五元十字阵定位示意图假设目标声源T以球面波形式进行传播，到达阵元S1

39、、S2、S3和S4相对于到达阵元S0的时延分别为、。同理，可以得到其目标定位方程为： (3.12) (3.13) (3.14) 3.3 四元阵列的定位精度分析及其仿真由十字阵的目标定位方程可知，十字阵被动声定位系统的定位精度跟空气声速C、阵元间距D和时延估计精度有关，其中时延估计精度是影响被动声定位精度的关键因素。下面分别就四阵元十字阵和五阵元十字阵的时延估计误差对系统的定向精度的影响进行分析。为了便于分析定向精度，首先讨论误差的合成公式。误差的合成是研究如何根据分项误差求总误差的问题。下面讨论一个普遍适用的公式误差传递公式11，并利用该公式推导出时延估计误差对定向精度的影响关系式。设变量y由

40、两个分项x1,x2合成，即y是x1,x2的函数： (3.15)设、分别是、的真值，如果函数在附近的各阶偏导数都存在，把上式在展幵成泰勒级数可以得到： (3.16)式中，表示展开式的余项。设，表示分项的误差，由于，则（3.16）式中的高阶小量可以舍去，总的合成误差可以表示为： (3.17)同理，当y由m个分项合成时，可以得到： (3.18)上式就是所需要的误差传递公式。 将式(3.18)的两边平方，得： 对上式进行n次求和，则： (3.19)如果x1、x2为相互独立的变量，则与互不相关，它们都是随机变量。因此，也是随机变化的，当时，(3.19)式的第二项趋近于零。不考虑第二项，将等式两端同除以n

41、，得到： (3.20)整理得到： (3.21)式(3.21)即是在m项分量相互独立时的方差合成公式。 3.3.1 方位角精度分析及仿真由于时延、的统计误差特征相同，不妨设（）的方差为,则根据(3.21)式，可得由时延估计误差引起的方位角误差为： (3.22)由3.2节推导的四元十字阵的目标定位方程12，可得： (3.23) (3.24)分别求偏导数，得： (3.25) (3.26) (3.27)将式(3.25)(3.27)联立，代入（3.22）式，得: (3.28)上式表明，对于四元十字阵，时延估计引起的方位角误差跟声速C、阵元间距D、俯仰角以及方位角有关。设声速为340m/s图3.3 不同阵

42、元间距下的方位角估计误差上图是对式(3.28)进行的仿真，由图可知，在时延估计误差一定时，方位角估计的精度和阵元间距以及目标的俯仰角、方位角有关。对于给定的时延估计误差，阵元间距越大，方位角估计的标准差越小，精度越高； 方位角越大，其估计的标准差越小，定位精度越高；俯仰角越大，定位精度也越高。图3.4 四元阵列的方位角定位精度由上图可知，方位角估计的精度还与时延估计误差有关，时延估计误差越大，方位角估计的标准差也越大，定位精度就越低。 3.3.2 俯仰角精度分析及仿真对于俯仰角，跟上述类似可得公式： (3.29)对公式(3.29)求偏导数得: (3.30)将方程组(3.30代入式(3.29)有

43、: (3.31)上式表明，对于四元十字阵，时延估计引起的俯仰角误差跟声速C、阵元间距D、方位角和俯仰角均有关。设声速C=340m/s图3.5 不同阵元间距下的俯仰角估计误差上图是对式(3.31)进行的仿真，由图可知，在时延估计误差一定时，俯仰角估计的精度和阵元间距以及目标的俯仰角、方位角有关。对于给定的时延估计误差，阵元间距越大，方位角估计的标准差越小，精度越高； 方位角越大，俯仰角估计的标准差越大，定位精度越低；同样俯仰角越大，定位精度也越低。图3.6 四元十字阵的俯仰角估计精度由上图可知，目标俯仰角的估计精度还受时延估计误差影响，时延估计误差越大，方位角估计的标准差也越大，定位精度就越低。

44、 3.3.3 距离估计精度分析及仿真由式(3.11 )可知，目标距离,结合和。由于对于对称阵来说，测距在三维空间和平面内的分析类似，考虑对称性可近似化简得到:目标距离估计的均方根误差为: (3.32)考虑目标距离估计的相对误差为: (3.33)式(3-32)表明:目标距离估计的均方根误差与阵元间距、目标距离、目标俯仰角、方位角以及时延估计精度有关。设声速C=340m/s图3.7 不同阵元间距下的目标距离估计误差由上图可知，目标距离估计的精度和阵元间距、目标距离以及目标俯仰角有关。在目标方位角及时延估计误差一定时，阵元间距越大，方位角估计的标准差越小，精度越高； 俯仰角越大，目标距离估计的标准差

45、越高，而定位精度越低。图3.8 不同目标距离时的时延估计误差由上图可知，目标距离估计误差还与时延估计误差有关，上图是在时延估计误差分别为1us,5us时，不同的目标距离下对目标距离估计误差的仿真。目标距离越大时，目标距离估计的标准差越大，定位精度越低；俯仰角越大时，定位精度也越低。图3.9 D=3m,r=100m，us时，五元十字阵的目标距离估计精度由上图可知，在时延估计误差、阵元间距及目标距离一定时，目标距离估计精度还与方位角有关，方位角越大时，目标距离定位精度稍有降低，但对其影响不大。3.4 五元阵列的定位精度分析及其仿真3.4.1 方位角精度分析及仿真设的方差为，则时延估计引起的方位角误

46、差为： (3.34)由(3.12)、(3.13)，可得： (3.35) (3.36)分别对上式求偏导数，得： (3.37) (3.38) (3.39) (3.40)将上式联立，代入（3.22）式，得： (3.41)由式(3.41)可看出，对于五阵元十字阵，时延估计引起的方位角误差跟声速 C阵元间距D和俯仰角有关，而跟目标声源的方位角无关。釆用四阵元十字阵时，阵元数目减少,降低了系统的复杂性，但仿真结果表明，目标方位角为0 度或180度附近时，方位角估计的误差较大。而五阵元十字阵系统的定向精度与目标的方位角无关。设时延估计的均方根误差为=5us，声速C=340m/s,由式(3.41)，可以得到五

47、元十字阵的目标方位角定位精度如图3.3所示。图3.10 不同阵元间距下的方位角估计误差由上图可知，五元阵中方位角估计的精度和阵元间距及目标的俯仰角有关，而与目标的方位角无关，克服了四元十字阵测向时受目标方位角影响的缺点。对于给定的时延估计精度，随阵元间距的增大，方位角估计精度提高；且目标俯仰角越大，方位角估计的均方根误差越小，精度也越高。 图3.11 五元阵的方位角估计精度由上图可知，方位角估计的精度还与时延估计误差有关，时延估计误差越大，方位角估计的标准差越大，定位精度越低；对于给定的时延估计精度，随着目标俯仰角的增大，精度提高。3.4.2 俯仰角精度分析及仿真同理，根据，可求得俯仰角对各时

48、延的偏导数为： (3.42)因此俯仰角方差同时延方差之间的关系式可表示为： (3.43)图3.12 不同阵元间距下的俯仰角估计误差由上图可知,目标俯仰角的估计精度与阵元间距及目标俯仰角有关，而与方位角无关。在时延估计误差一定时，增大阵元间距可以提高目标的俯仰角估计精度；与目标方位角估计精度情况相反，随着目标俯仰角的增大，俯仰角估计的标准差上升，而定位精度下降。 图3.13 五元阵的俯仰角估计精度由上图可知，目标俯仰角的估计精度还与时延估计误差有关，时延估计误差越大，俯仰角估计的标准差越大，定位精度越低。在时延估计误差一定时，随着目标俯仰角的增大，定位精度下降。 3.4.3 距离估计精度分析及仿

49、真 由距离r对各时延的偏导数可得： （） (3.44) 因此距离方差同时延方差之间的关系式可表示为： (3.45)图3.14 时延估计误差为5us时，不同阵元间距下的目标距离估计精度由上图可知,目标距离的定位精度与阵元间距、目标距离、目标的俯仰角有关，而与目标的方位角无关。阵元间距越大，目标距离定位精度越高；而在一定的阵元间距下和时延估计误差下，随着目标距离增大，定位精度有所下降；随着目标俯仰角的增大，定位精度也下降。 图3.15 不同目标距离时的时延估计误差由上图及式（3.45）可知，目标距离估计误差还与时延估计误差有关，上图是在目标距离估计误差分别为0.2,0.3时，不同的目标距离下对时延

50、估计误差的仿真。目标距离越大时，时延估计误差越小；目标距离估计误差越大，时延估计误差也越大；而俯仰角越大，时延估计误差越小。图3.16 D=3m,r=100m，us时，五元十字阵的目标距离估计精度由上图可知，在时延估计误差、阵元间距及目标距离一定时，目标距离估计精度与俯仰角有关，而与方位角无关。目标俯仰角的越大，定位精度越高。3.5 本章小结本章论述了基于时延估计的声源定位原理，研究了基于麦克风阵列的声源定位的算法，推导出了四元、五元十字阵的声源定位方程，并对四元和五元两种十字阵的定位精度进行了理论分析和对比，最后利用matlab仿真软件对其精度进行了仿真分析。4 多元麦克风阵列声源定位分析4

51、.1 多元麦克风阵列定位方程随着现代科技的发展，定位技术在航空、航天、交通、海 洋资源勘探等领域得到了广泛地应用。近年来，针对某些特殊的阵形(如平面三元阵、平面四元阵)的研究较多，而对任意阵列模型的研究较少13。由于实验设备几何外形尺寸的限制以及在野外传感器布设时受地形条件的影响(如基于智能浮漂阵列的水声定位系统中，浮漂单元易受海洋波浪、风等因素的影响而发生随机的移动)，很多情况下平面阵列不能满足实际应用的需要，阵列布设往往需要任意阵。针对试验中阵列传感器布设中存在的问题，本文提出了任意多元阵列定位模型14，利用最小二乘法的估计特性，解决迭代法的初始值问题，以提高定位精度。以声阵列定位系统为例

52、，通过试验验证，证明了算法的有效性。假设空间任意分布的N元传感器阵列()， 其空间相互位置已知并且在同一坐标系中，如图（4.1）所示。其中，为坐标原点（0,0,0），其它传感器的位置坐标为（=1,2n）,点为目标位置，信源在介质中的传播速度为, 表示目标位置P到各传感器(i=0,1,2,3n)的距离, 分别为目标到第i接收传感器与到第0接收传感器时间差，则有，因此可建立n-1个定位方程15。 图4.1 任意多元阵列定位原理图（=1，2,3n） (4.1) 式中: 整理上式(4.1)可得如下线形方程组： (4.2)式中： （=1,2n）令： 则上面方程式可简化为：求解上述方程组可得到目标位置。当

53、方程的个数大于未知数的个数时，等价于非线性最优化问题，可采用改进算法得到最优解。 (1)理论上当空间布设的传感器的个数为5时，可依据线性方程组(2)求解目标位置和目标到坐标原点的距离，而实际中为提高系统的定位精度和定位范围，传感器的个数要远超过5，即列出的方程的数目大于需求解的未知数的数目，因此采用最小二乘法使残差平方和最小，以提高定位精度。4.2 最小二乘法求声源位置最小二乘法求解是基于由多个传感器获得的到达时间所建立的式（4.3）所给出的固定方程组(4.2)得到声发射源位置坐标。 (4.3)对于线性组合的方程组(4.2)(式中n5)，利用最小二乘法求解16。假设x,y,z, 表示各测定值的最可信赖值，且以表示各测定值对应的残差,则有残差方程组： (4.4)根据最小二乘法的意义，使残差平方和最小,也即：有最小值.即： 由以上各式可得到如下正规方程组: