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个阵元在同一时刻接收到的信号为:
其中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
其中nvxm(t)和nvym(t)分别为子阵Xv和Yv上第m个阵元的加性噪声;将上式写成矢量形式有:
Xv(t)=Ays(t)+Nvx(t) (公式4)
Yv(t)=AyΦys(t)+Nvy(t) (公式5)
其中,为(M-1)×1的子阵列Xv导向矢量矩阵,
Xv(t)=[xv1(t),xv2(t),…,xv(M-1)(t)]T为(M-1)×1的子阵列Xv接收信号矩阵,
Yv(t)=[yv1(t),yv2(t),…,yv(M-1)(t)]T为(M-1)×1的子阵列Yv接收信号矩阵,
Nvx(t)=[nvx1(t),nvx2(t),…,nvx(M-1)(t)]T为(M-1)×1的子阵列Xv噪声矩阵,
Nvy(t)=[nvy1(t),nvy2(t),…,nvy(M-1)(t)]T为(M-1)×1的子阵列Yv噪声矩阵;
同样,Φy是把子阵Xv和Yv的输出联系起来的一个因子,也称旋转算子,它包含了目标反射回来的信号波前在任意竖直均匀线阵中一个阵元偶之间的相位延迟信息,表示为:
3.根据权利要求2所述的一种用于二维波达方向估计的声速测量方法,其特征在于步骤二具体包括:
对于水平均匀线阵,计算Xh(t)的协方差矩阵Rhxx与Xh(t)和Yh(t)的互协方差矩阵Rhxy,然后对其进行特征值分解得到最小特征值为σx2,利用σx2计算Chxx=Rhxx-σx2I和Chxy=Rhxy-σx2Z;其中I是(M-1)×(M-1)阶的单位矩阵,Z也是(M-1)×(M-1)阶的矩阵,并且
最后计算矩阵束{Chxx,Chxy}的广义特征值分解,得到非零特征值λx,此特征值就是旋转算子的值,即:
对于竖直均匀线阵,同样计算Xv(t)的协方差矩阵Rvxx与Xv(t)和Yv(t)的互协方差矩阵Rvxy,然后对其进行特征值分解得到最小特征值为σy2,利用σy2计算Cvxx=Rvxx-σy2I和Cvxy=Rvxy-σy2Z。最后计算矩阵束{Cvxx,Cvxy}的广义特征值分解,得到非零特征值λy,此特征值就是旋转算子的值,即:
4.根据权利要求3所述的一种用于二维波达方向估计的声速测量方法,其特征在于步骤三具体包括:
在水面位置对水下的探测路径声速进行测量,只考虑从水平均匀线阵下方反射的声波,反射回来的声波分为两种,一种是声波①在水平均匀线阵上到达的第一个阵元不是参考阵元;另一种是声波②在水平均匀线阵上到达的第一个阵元是参考阵元,,设θx1和θy1分别为声波①与两个均匀线阵法线的夹角,它们相加之和为θx2和θy2分别为声波②与两个均匀线阵法线的夹角,它们的相加之和也为
当声波在水平均匀线阵上到达的第一个阵元不是参考阵元,即属于声波①的情形时,参考阵元x1是最晚接收到信号的,因此时延参数τ是小于0的,又因为所以此时声波在水平线阵上的方向角是负的,则有θx=-θx1,θy=θy1,由此可得θx和θy的关系为:
当声波在水平均匀线阵上到达的第一个阵元是参考阵元,即属于声波②的情形时,有θx=θx2,θy=θy2,由此可得θx和θy的关系为:
根据公式(9)和公式(10)可以得到:
5.根据权利要求4所述的一种用于二维波达方向估计的声速测量方法,其特征在于步骤四具体包括:
将公式(11)带入公式(8)中,则有:
则联立公式(7)和公式(12)消去θx可得:
整理可得速度v的最终表达为:
因此,利用正交均匀线阵接收到的声波信号求出两个旋转算子的值即λx和λy,然后将其带入公式(14)便可计算出二维DOA估计中探测路径上声速v的大小。
6.根据权利要求5所述的一种用于二维波达方向估计的声速测量方法,其特征在于:当有多个目标即多条探测路径时,两个均匀线阵上的旋转算子是一个对角矩阵,每一个对角元素都对应于一个目标,求得的特征值也有多个,它们一一对应着对角矩阵上的对角元素,最后利用这些特征值求得二维DOA估计中每条探测路径上的声速。
7.实现权利要求1~6任一项所述测量方法的装置,其特征在于包括数据采集处理与控制模块、发射模块、接收模块、输出模块和电源模块;数据采集处理与控制模块由一个处理器、一个D/A转换器和一个A/D转换器组成,并且与发射模块、接收模块、输出模块和电源模块相连;处理器能够通过设定参数来控制发射模块,使发射模块发射指定的声波信号;还能够对接收模块传送回来的数据进行采集与处理,然后运计算出二维DOA估计中精准的水下探测路径声速大小;处理器还将处理好的数据传送给输出模块,再通过USB接口输出到外部装置或者在显示器上显示出来。
8.根据权利要求7所述的装置,其特征在于发射模块由一个阻抗匹配电路和一个超声波发送探头组成,并且与数据采集处理与控制模块和电源模块相连,能够根据数据采集处理与控制模块的指令,获取指定的信号并进行超声发射;
接收模块由两个相互正交的超声波接收探头阵列和阻抗匹配电路组成,与数据采集处理与控制模块和电源模块相连,能够接收从目标发射回来的信号,并传送给数据采集处理与控制模块;
输出模块由一个USB接口和一个显示器组成,并且与数据采集处理与控制模块模块和电源模块相连,输出模块能够提供人机交互,将数据采集处理与控制模块模块中处理好的数据通过USB接口输出到外部装置或者在显示器上显示出来。