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

《基于四元十字阵的声源定位终期报告.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、.信号课外设计制作总结报告题目：信号-微弱声源定位系统组号：B28组长：14120763宛加明张海燕20%成员：14121596曾清清陈俊丽20%成员：14121572叶善云陈俊丽20%成员：14120946曾庆权张海燕20%成员：13121253马文辉孙晓岚20%联系方式：二零一六年五月二十日Word文档.一、实验要求设计要求一）设计要求  设计微弱声源定位系统，要求：（1） 声源功率尽量小；（2）写出数学模型，说明详细算法； (3) 能够现场演示。评分依据（1） 功能实现；（2） 性能指标；（3） 工程规范（软件规范、可靠性等）；（4） 理论水平；（5） 团队分工合作情况。在确定了项目小

2、组之后，我们组的组员一起讨论了该如何实现这个微弱声源定位系统并进行了分工，经过我们的讨论之后，我们大致把项目设计分为如下几个过程：1搞懂设计方案的原理，利用所学知识以及上网查找资料，初步确定设计的方向。2分析这个项目的功能要求，一起讨论并编写代码。3利用代码使用MATLAB进行实验，对一些地方进行修改加工。4进行最后的调整，然后纂写报告。二、实验原理1．概述本实验采用的环境为Matlab，核心算法为四元十字算法。声源定位基本原理是将传声器在空间布置成一定几何形状的阵列，以接收目标噪声的声场信息，通过检测或计算各传声器所测得的信号的时延来确定目标的位置。对于无方向性传声器，各个传声器接收到同一

3、声源的信号先后，反映了声源所在传声器连线的方向和距离。Word文档.声学阵列的设计在被动声定位中具有十分重要的意义。采用的传声器阵列可分为线阵、面阵和立体阵。对于固定式阵列，线阵由于其轴对称性，故在定位时会造成空间模糊。立体阵可以对整个空间进行定位，但其算法要复杂的多。面阵可以在整个平面对目标进行定位，也可以对阵列所在平面为界的半个空间进行定位。一般要对目标进行定位，至少需要四个阵元组成的声传感器阵列。由于十字形阵列具有分维持性(指二维参量可分开估计)，且阵列冗余度也较小，因此，平面十字阵是较为合适的阵形。2.四元十字算法原理平面四元十字阵进行分析。对十字形声传感器阵列，可建立如图3.1.1

4、所示的直角坐标系。四个阵元的坐标分别为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的关系式可推导如下:(设

5、C为声速)球坐标和直角坐标之间的互化。如式：求解方程式，消去x，y，z可得出目标距离，方位角和俯仰角，得出结果。Word文档.3.时延估计常用的时延估计方法有广义互相关法(Gcc)、互功率谱相位法(csP)、自适应滤波法和参量模型法等多种方法,广义相互函数法、互功率谱相位法、自适应滤波法和参量模型法等几种方法中，互功率谱相位法由于目标信号频率较低，实时性要求高，最小二乘拟合次数不能太高，故估计精度一般不高。自适应滤波法运算速度快，但点声源干扰对其估计精度影响较大。广义互相关法具有较好的综合性能，它是目前人们关注较多的时延估计方法。广义互相关法(GCC)互相关，顾名思义是通过求两路信号的互相关

6、来估计时延D。x1(t),x2(t)互相关函数为:由式(4.1.3)可知，x1(t),x2(t)互相关函数取最大值时，Rss(τ-D)也取最大值，又因为Rss(τ-D)<=Rss(0)，所有取得最大值的时候的τ即为时延D.由于观察时间受限，所以求得的互相关只是一个估计值，记为Rx1x2(τ)在各态历经的过程中，定义为:4.定位算法Word文档.Word文档.3.伪随机白噪声产生仿真过程中要人为的加入噪声干扰，以模拟目标的运动噪声和观测噪声。在仿真中加入的噪声一般为高斯随机白噪声。其产生方法有三种:第一种利用专门的随机数表:第二种利用物理装置产生随机噪声;第三种利用计算机通过数学方法获得，即运

7、用专门的程序产生伪随机噪声。前面两种由于其固有的缺陷而降低了其使用价值，后一种由于其灵活方便，在计算机仿真中被广泛应用[(s]。严格说来计算机生成的随机数显然不是真正的随机数，因为它生成的方式是完全确定的，所以常把这种随机数称为伪随机数。一般伪随机数个数远小于它的周期长度，可以当作随机数来使用。在matlab中可以直接用命令来产生伪随机白噪声。Word文档.三、代码展示clearfs=10000;%频率t=0