基于四元十字阵的声源定位终期报告选编.docxVIP

 PAGE \* MERGEFORMAT 15 信号课外设计制作 总结报告 题 目： 信号-微弱声源定位系统 组 号： B28 组 长： 加明 张海燕 20% 成 员： 曾清清 陈俊丽 20% 成 员： 叶善云 陈俊丽 20% 成 员： 曾庆权 张海燕 20% 成 员： 马文辉 孙晓岚 20% 联系方式： 二零一六年 五月二十日 一、 实验要求 设计要求 一）设计要求?? 设计微弱声源定位系统，要求： （1）?声源功率尽量小； （2） 写出数学模型，说明详细算法； ? (3)? 能够现场演示。 评分依据（1）?功能实现；（2）?性能指标；（3）?工程规范（软件规范、可靠性等）；（4）?理论水平；（5）?团队分工合作情况。 在确定了项目小组之后，我们组的组员一起讨论了该如何实现这个微弱声源定位系统并进行了分工，经过我们的讨论之后，我们大致把项目设计分为如下几个过程： 1搞懂设计方案的原理，利用所学知识以及上网查找资料， 初步确定设计的方向。 2 分析这个项目的功能要求，一起讨论并编写代码。 3利用代码使用MATLAB进行实验，对一些地方进行修改加工。 4进行最后的调整，然后纂写报告。 二、实验原理 概述 本实验采用的环境为Matlab，核心算法为四元十字算法。 声源定位基本原理是将传声器在空间布置成一定几何形状的阵列，以接收目标噪声的声场信息，通过检测或计算各传声器所测得的信号的时延来确定目标的位置。对于无方向性传声器，各个传声器接收到同一声源的信号先后，反映了声源所在传声器连线的方向和距离。 声学阵列的设计在被动声定位中具有十分重要的意义。采用的传声器阵列可分为线阵、面阵和立体阵。对于固定式阵列，线阵由于其轴对称性，故在定位时会造成空间模糊。立体阵可以对整个空间进行定位，但其算法要复杂的多。面阵可以在整个平面对目标进行定位，也可以对阵列所在平面为界的半个空间进行定位。一般要对目标进???定位，至少需要四个阵元组成的声传感器阵列。由于十字形阵列具有分维持性(指二维参量可分开估计)，且阵列冗余度也较小，因此，平面十字阵是较为合适的阵形。 2.四元十字算法原理 平面四元十字阵进行分析。 对十字形声传感器阵列，可建立如图3.1.1所示的直角坐标系。四个阵元的坐标分别为S1 (D /2,0,0) , S2(0, D / 2,0) , S3 (-D / 2,0,0) , S4 (0,-D / 2,0)，设目标声源的直角坐标为((x, y, z)、球坐标为((r,θ ,ψ )，即目标到坐标原点的距离为r，方位角为θ，俯仰角为 ψ，D为十字阵对角线阵元之间的距离 。 假设目标为点声源，目标产生的声源以球面波形式传播，声源到达阵元S1的传播时间为t1，相对于S1声源到达阵元S2 , S3 ， S4的时延分别为 τ12， τ13 ，τ14。目标到坐标原点的距离为r，俯仰角为θ(0=θ=90)，方位角为ψ(0=ψ=90)与时延τ 12， τ13 ， τ 14及阵元间距D的关系式可推导如下:(设C为声速) 球坐标和直角坐标之间的互化 。如式： 求解方程式，消去x，y，z可得出目标距离，方位角和俯仰角，得出结果。 时延估计 常用的时延估计方法有广义互相关法(Gcc )、互功率谱相位法(csP )、自适应滤波法和参量模型法等多种方法,广义相互函数法、互功率谱相位法、自适应滤波法和参量模型法等几种方法中，互功率谱相位法由于目标信号频率较低，实时性要求高，最小二乘拟合次数不能太高，故估计精度一般不高。自适应滤波法运算速度快，但点声源干扰对其估计精度影响较大。广义互相关法具有较好的综合性能，它是目前人们关注较多的时延估计方法。 广义互相关法(GCC) 互相关，顾名思义是通过求两路信号的互相关来估计时延D 。 x1 (t) , x2 (t)互相关 函数为: 由式(4. 1. 3)可知，x1(t) , x2(t)互相关函数取最大值时，Rss(τ - D)也取最大值，又因为Rss(τ-D) =Rss(0)，所有取得最大值的时候的τ即为时延D. 由于观察时间受限，所以求得的互相关只是一个估计