一种用于二维波达方向估计的声速测量方法及装置制造方法及图纸

技术编号:14836419 阅读:130 留言:0更新日期:2017-03-17 03:58
本发明专利技术公开了一种用于二维波达方向估计的声速测量方法及装置。测量方法包括如下步骤:确定水平均匀线阵和竖直均匀线阵的接收信号模型;计算水平均匀线阵和竖直均匀线阵的旋转算子;建立声波在正交均匀线阵上两个方向角之间的关系;根据两方向角之间的关系求出探测路径上声速的大小。装置包括数据采集处理与控制模块、发射模块、接收模块、输出模块和电源模块;数据采集处理与控制模块由一个处理器、一个D/A转换器和一个A/D转换器组成,并且与发射模块、接收模块、输出模块和电源模块相连;本发明专利技术的一种用于二维DOA估计的声速测量方法,实现对二维DOA估计中的探测路径声速大小进行精确的测量,本发明专利技术装置可行性强,安装简单。

【技术实现步骤摘要】

本专利技术涉及声速测量的
,尤其是指一种用于二维波达方向估计的声速测量方法及装置
技术介绍
空间信号的波达方向(简称:DOA)估计是阵列信号处理领域的一项重要技术,它利用阵列对空间信号进行接收和处理,提取有用的信号特征或信息,以此来估计信号的入射角,并常常用于水下目标定位中。二维DOA估计可以对目标的方向角和距离进行估计,但最终的估计结果受探测路径上的声速大小影响很大,如果探测路径声速测量得不精准,就会对二维DOA估计的结果产生影响。而在进行水下目标定位时,二维DOA估计通常是利用一片海水区域内的声速来代替每条探测路径上的声速,但是海水环境复杂,且受水流因素影响较大,每一条探测路径上的声速其实是不一样的,即二维DOA估计中对探测路径声速的确定是不准确的,这将导致二维DOA估计的估计误差较大。因此精准的水下探测路径声速测量对二维DOA估计至关重要。目前,海水声速测量方法主要分为两种:一种是直接法,另外一种是间接法。直接法是用水声设备直接在海洋现场测量海水声速,它可分为时差法、共振声谱法、驻波干涉法与相位比较法等方法。而间接法则是根据海水中的深度、温度和盐度三个主要影响海水声速大小的因素,经过成百上千次的分析和实践,提出一种经验公式,利用此经验公式来测量海水声速。根据经验公式的不同,间接法可分为DelGrosso声速算法、Wilson声速算法和Chen-Millero-Li声速算法三种声速算法。直接法和间接法虽然原理简单且操作较方便,但它们都需要在目标位置已知的情况下才能对某一条探测路径上的声速进行测量,而在二维DOA估计中,目标位置都未知的且是用来估计的;同时,间接法测量的是一个声速分布较均匀的海水环境中的声速,它测量的并不是某一条路径上的声速,即点对点声速,因此直接法和间接法对二维DOA估计的水下探测路径声速大小的测量并不具有实用性。为了克服在利用直接法和间接法测量二维DOA估计的声速过程中存在的局限性,本专利使用两个相互正交的均匀阵列作为接收阵列,利用从目标反射回来的声波在两个阵列中形成的方向角间的关系,得到二维DOA估计中探测路径声速与两个旋转算子之间的关系式,然后再计算出声速的大小,实现了在二维DOA估计中对水下探测路径声速的精准测量。因此,本专利方法对二维DOA估计中的声速测量更具有实用性。
技术实现思路
本专利技术的目的在于克服现有方法存在的各种不足,提供一种用于二维DOA估计的声速测量方法,实现对二维DOA估计中的探测路径声速大小进行精确的测量。另外,本专利技术另一目的在于提供一种用于二维DOA估计的声速测量装置。本专利技术方法的原理是基于从目标反射回来的声波在两个正交的均匀线阵上形成的方向角相加或相减等于的关系,得出两个均匀线阵上的旋转算子与二维DOA估计中的探测路径声速间的关系式,最后通过计算两均匀线阵上的旋转算子来求得最终探测路径声速大小的。为了实现上述目的,本专利技术至少通过以下技术方案实现。一种用于二维波达方向估计的声速测量方法,其包括如下步骤:步骤一:确定水平均匀线阵和竖直均匀线阵的接收信号模型;步骤二:计算水平均匀线阵和竖直均匀线阵的旋转算子;步骤三:建立声波在正交均匀线阵上两个方向角之间的关系;步骤四:根据两方向角之间的关系求出探测路径上声速v的大小。进一步地,步骤一具体包括:正交均匀线阵的水平和竖直均匀线阵都有M个接收阵元,M为正整数,且阵元间间距为d,窄带目标声源为S,中心频率为f,目标对应于水平线阵的方向角为θx,对应于竖直线阵的方向角为θy;其中水平均匀线阵的M个阵元为x1,x2,…,xM,竖直均匀线阵的M个阵元为y1,y2,…,yM;将水平线阵中的M个阵元分为平移矢量为d的两个子阵列Xh和Yh;子阵列Xh由水平线阵中的第一到第M-1个阵元组成,即有:xh1(t)=x1(t),xh2(t)=x2(t),…,xh(M-1)(t)=xM-1(t)其中,xh1(t),xh2(t),…,xh(M-1)(t)分别是子阵列Xh上第一个阵元到第M-1个阵元接收到的信号;子阵列Yh由水平线阵的第二到第M个阵元组成,即有:yh1(t)=x2(t),yh2(t)=x3(t),…,yh(M-1)(t)=xM(t)其中,yh1(t),yh2(t),…,yh(M-1)(t)分别是子阵列Yh上第一个阵元到第M-1个阵元接收到的信号;x1(t),x2(t),…,xM(t)是水平线阵中第一个阵元到第M个阵元的接收信号;子阵列Xh接收信号中,以第一个阵元xh1为参考点,则第一个阵元接收的信号为:xh1(t)=s(t)+nhx1(t)其中s(t)表示目标信号,nhx1(t)表示子阵列Xh第一个阵元上的噪声;接收信号满足窄带条件,即当信号延迟远小于带宽倒数时,延迟作用相当于使基带信号产生一个相移;那么子阵列Xh上第m个阵元在同一时刻接收到的信号为:xhm(t)=s(t)am(θx)+nhxm(t),m=1,2,…,M-1其中am(θx)中v表示声波在探测路径上的速度,nhxm(t)表示子阵列Xh上第m个阵元上的噪声;由于子阵列Yh和子阵列Xh的相对平移矢量为d,那么子阵列Yh上第m个阵元在同一时刻接收到的信号为:其中nhym(t)表示子阵列Xh上第m个阵元上的噪声;将子阵列Xh和子阵列Yh上的各阵元的接收信号排列成列向量形式,则上面两式可写成矢量形式:Xh(t)=Axs(t)+Nhx(t)(公式1)Yh(t)=AxΦxs(t)+Nhy(t)(公式2)其中,为(M-1)×1的子阵列Xh导向矢量矩阵,Xh(t)=[xh1(t),xh2(t),…,xh(M-1)(t)]T为(M-1)×1的子阵列Xh接收信号矩阵,Yh(t)=[yh1(t),yh2(t),…,yh(M-1)(t)]T为(M-1)×1的子阵列Yh接收信号矩阵,Nhx(t)=[nhx1(t),nhx2(t),…,nhx(M-1)(t)]T为(M-1)×1的子阵列Xh噪声矩阵,Nhy(t)=[nhy1(t),nhy2(t),…,nhy(M-1)(t)]T为(M-1)×1的子阵列Yh噪声矩阵;而Φx是把子阵Xh和Yh的输出联系起来的一个因子,称旋转算子,它包含了目标反射回来的信号波前在任意水平均匀线阵中一个阵元偶之间的相位延迟信息,表示为:对于竖直均匀线阵,同样将其分为两个平移矢量为d的子阵列Xv和Yv;子阵列Xv由竖直阵列的第一到第M-1个阵元组成,则有:xv1(t)=y1(t),xv2(t)=y2(t),…,xv(M-1)(t)=yM-1(t)其中,xv1(t),xv2(t),…,xv(M-1)(t)分别是子阵列Xv上第一个阵元到第M-1个阵元接收到的信号;子阵列Yv由竖直阵列的第二到第M个阵元组成,则有:yv1(t)=y2(t),yv2(t)=y3(t),…,yv(M-1)(t)=yM(t)其中,yv1(t),yv2(t),…,yv(M-1)(t)分别是子阵列Yv上第一个阵元到第M-1个阵元接收到的信号;y1(t),y2(t),…,yM(t)是竖直线阵中第一个阵元到第M个阵元的接收信号;以第一个阵元xv1为参考阵元,那么子阵列Xv和子阵列Yv中第m个阵元的接收信号分别为:xvm(t)=s(t)am(θy)+nvxm(t),m=1,2,…,M-1其本文档来自技高网
...
一种用于二维波达方向估计的声速测量方法及装置

【技术保护点】
一种用于二维波达方向估计的声速测量方法,其特征在于包括如下步骤:步骤一:确定水平均匀线阵和竖直均匀线阵的接收信号模型;步骤二:计算水平均匀线阵和竖直均匀线阵的旋转算子;步骤三:建立声波在正交均匀线阵上两个方向角之间的关系;步骤四:根据两方向角之间的关系求出探测路径上声速v的大小。

【技术特征摘要】
1.一种用于二维波达方向估计的声速测量方法,其特征在于包括如下步骤:步骤一:确定水平均匀线阵和竖直均匀线阵的接收信号模型;步骤二:计算水平均匀线阵和竖直均匀线阵的旋转算子;步骤三:建立声波在正交均匀线阵上两个方向角之间的关系;步骤四:根据两方向角之间的关系求出探测路径上声速v的大小。2.根据权利要求1所述的一种用于二维波达方向估计的声速测量方法,其特征在于步骤一具体包括:正交均匀线阵的水平和竖直均匀线阵都有M个接收阵元,M为正整数,且阵元间间距为d,窄带目标声源为S,中心频率为f,目标对应于水平线阵的方向角为θx,对应于竖直线阵的方向角为θy;其中水平均匀线阵的M个阵元为x1,x2,…,xM,竖直均匀线阵的M个阵元为y1,y2,…,yM;将水平线阵中的M个阵元分为平移矢量为d的两个子阵列Xh和Yh;子阵列Xh由水平线阵中的第一到第M-1个阵元组成,即有:xh1(t)=x1(t),xh2(t)=x2(t),…,xh(M-1)(t)=xM-1(t)其中,xh1(t),xh2(t),…,xh(M-1)(t)分别是子阵列Xh上第一个阵元到第M-1个阵元接收到的信号;子阵列Yh由水平线阵的第二到第M个阵元组成,即有:yh1(t)=x2(t),yh2(t)=x3(t),…,yh(M-1)(t)=xM(t)其中,yh1(t),yh2(t),…,yh(M-1)(t)分别是子阵列Yh上第一个阵元到第M-1个阵元接收到的信号;x1(t),x2(t),…,xM(t)是水平线阵中第一个阵元到第M个阵元的接收信号;子阵列Xh接收信号中,以第一个阵元xh1为参考点,则第一个阵元接收的信号为:xh1(t)=s(t)+nhx1(t)其中s(t)表示目标信号,nhx1(t)表示子阵列Xh第一个阵元上的噪声;接收信号满足窄带条件,即当信号延迟远小于带宽倒数时,延迟作用相当于使基带信号产生一个相移;那么子阵列Xh上第m个阵元在同一时刻接收到的信号为:xhm(t)=s(t)am(θx)+nhxm(t),m=1,2,…,M-1其中am(θx)中v表示声波在探测路径上的速度,nhxm(t)表示子阵列Xh上第m个阵元上的噪声;由于子阵列Yh和子阵列Xh的相对平移矢量为d,那么子阵列Yh上第m个阵元在同一时刻接收到的信号为:yhm(t)=s(t)j2πfvdsinθxam(θx)+nhym(t),m=1,2,...,M-1]]>其中nhym(t)表示子阵列Xh上第m个阵元上的噪声;将子阵列Xh和子阵列Yh上的各阵元的接收信号排列成列向量形式,则上面两式可写成矢量形式:Xh(t)=Axs(t)+Nhx(t)(公式1)Yh(t)=AxΦxs(t)+Nhy(t)(公式2)其中,为(M-1)×1的子阵列Xh导向矢量矩阵,Xh(t)=[xh1(t),xh2(t),…,xh(M-1)(t)]T为(M-1)×1的子阵列Xh接收信号矩阵,Yh(t)=[yh1(t),yh2(t),…,yh(M-1)(t)]T为(M-1)×1的子阵列Yh接收信号矩阵,Nhx(t)=[nhx1(t),nh...

【专利技术属性】
技术研发人员:宁更新王波文宁秋燕杨萃张军冯义志季飞
申请(专利权)人:华南理工大学
类型:发明
国别省市:广东;44

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1