【技术实现步骤摘要】
本专利技术涉及声速测量的
,尤其是指一种用于二维波达方向估计的声速测量方法及装置。
技术介绍
空间信号的波达方向(简称:DOA)估计是阵列信号处理领域的一项重要技术,它利用阵列对空间信号进行接收和处理,提取有用的信号特征或信息,以此来估计信号的入射角,并常常用于水下目标定位中。二维DOA估计可以对目标的方向角和距离进行估计,但最终的估计结果受探测路径上的声速大小影响很大,如果探测路径声速测量得不精准,就会对二维DOA估计的结果产生影响。而在进行水下目标定位时,二维DOA估计通常是利用一片海水区域内的声速来代替每条探测路径上的声速,但是海水环境复杂,且受水流因素影响较大,每一条探测路径上的声速其实是不一样的,即二维DOA估计中对探测路径声速的确定是不准确的,这将导致二维DOA估计的估计误差较大。因此精准的水下探测路径声速测量对二维DOA估计至关重要。目前,海水声速测量方法主要分为两种:一种是直接法,另外一种是间接法。直接法是用水声设备直接在海洋现场测量海水声速,它可分为时差法、共振声谱法、驻波干涉法与相位比较法等方法。而间接法则是根据海水中的深度、温度和盐度三个主要影响海水声速大小的因素,经过成百上千次的分析和实践,提出一种经验公式,利用此经验公式来测量海水声速。根据经验公式的不同,间接法可分为DelGrosso声速算法、Wilson声速算法和Chen-Millero-Li声速算法三种声速算法。直接法和间接法虽然原理简单且操作较方便,但它们都需要在目标位置已知的情况下才能对某一条探测路径上的声速进行测量,而在二维DOA估计中,目标位置都未知的且是用来 ...
【技术保护点】
一种用于二维波达方向估计的声速测量方法,其特征在于包括如下步骤:步骤一:确定水平均匀线阵和竖直均匀线阵的接收信号模型;步骤二:计算水平均匀线阵和竖直均匀线阵的旋转算子;步骤三:建立声波在正交均匀线阵上两个方向角之间的关系;步骤四:根据两方向角之间的关系求出探测路径上声速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
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。