CN105607033A - 基于正交均匀线阵的水下波达方向估计方法及*** - Google Patents
基于正交均匀线阵的水下波达方向估计方法及*** Download PDFInfo
- Publication number
- CN105607033A CN105607033A CN201610126961.6A CN201610126961A CN105607033A CN 105607033 A CN105607033 A CN 105607033A CN 201610126961 A CN201610126961 A CN 201610126961A CN 105607033 A CN105607033 A CN 105607033A
- Authority
- CN
- China
- Prior art keywords
- array
- even linear
- linear array
- module
- matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 76
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 15
- 238000003491 array Methods 0.000 claims abstract description 11
- 239000000523 sample Substances 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 91
- 230000008569 process Effects 0.000 claims description 46
- 238000000354 decomposition reaction Methods 0.000 claims description 28
- 239000013598 vector Substances 0.000 claims description 20
- 230000003993 interaction Effects 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 238000010304 firing Methods 0.000 claims description 2
- 238000002604 ultrasonography Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 abstract 1
- 230000008901 benefit Effects 0.000 description 3
- 230000004807 localization Effects 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000004304 visual acuity Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/80—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
- G01S3/802—Systems for determining direction or deviation from predetermined direction
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于正交均匀线阵的水下波达方向估计方法及***,所述估计***包括数据采集处理与控制模块、发射模块、接收模块、输出模块和电源模块,其中,所述接收模块中超声波接收探头阵列为正交的均匀直线阵列,其中每个阵列包括多个超声接收探头,并且数量相同。本发明实现声波在水中传播速度不确定的情况下对目标进行精确定位,对ESPRIT算法进行了改进,在算法中消去了声速,使得最后的估计结果和声速无关,从而排除了声速的不确定对波达方向估计的影响,使得其实用性更高,并且可在不同水声环境下进行测量,水下目标定位的精确性也更高。
Description
技术领域
本发明涉及目标定位的技术领域,特别涉及一种基于正交均匀线阵的水下波达方向估计方法及***。
背景技术
空间信号的波达方向(简称:DOA)估计问题是阵列信号处理所研究的一类重要问题。所谓阵列信号处理就是指在空间的不同位置放置多个传感器组成阵列,利用阵列对空间信号进行接收和处理,目的是抑制无用的干扰和噪声,提取有用的信号特征或信息。DOA估计便是用这种处理方法来估计信号的入射角。
现有的DOA估计方法主要有波束形成法、高分辨率谱估计法、子空间拟合类算法和子空间分解类算法。第一,波束形成法将传统时域傅立叶谱估计扩展成相应的空域形式,即通过空域各阵列的接收数据代替传统时域处理中的时域数据。但扩展至空域后,阵列的角度分辨力受到空域“傅立叶限”的限制。空域“傅立叶限”就是阵列的物理孔径限,即对位于一个波束宽度内的空间目标不可分辨。所以,提高空域处理精度的有效方法就是增大天线孔径(等效于减小波束宽度),但对于许多实际应用环境而言,增大天线孔径是不现实的。第二,高分辨率谱估计法是一种空域谱估计方法,它由许多时域非线性谱估计方法推广而成。它假定信号源在空间中是连续分布的,信号是空间平稳的随机过程。但这样的假设在大多数空间谱估计中是不成立的,因而这些方法具有局限性;并且这种方法不能有效利用加性噪声的统计特性,因而其分辨性能较差。第三,子空间拟合类算法将各种不同的方位估计方法用统一的算法结构联系起来,使其协方差矩阵的估计误差达到最小。这种算法精确度高,分辨率强,但是一种迭代算法,运算量非常大。第四,子空间分解类算法通过对阵列接收数据的数学分解(如奇异值分解、特征分解及QR分解等),将接收数据划分为两个相互正交的子空间:一个与信号源的阵列流型空间一致的信号子空间和一个与信号子空间正交的噪声子空间。它利用两个子空间的正交性可以大大提高算法的分辨力。典型代表算法为ESPRIT算法。ESPRIT算法建立在子空间旋转不变技术的基础上,不需要全空间搜索,减少了运算量。因此,ESPRIT算法具有现实可行、分辨率高和运算量较小的特点,从而克服了前三种方法存在的不足,更具优越性。但是ESPRIT算法假定信号在介质中的传播速度是确定不变的,没有考虑速度对算法结果的影响,在速度变化较大、不确定的情况下容易产生较大的误差,大大降低了目标定位的精确度。
发明内容
本发明的目的在于克服现有技术的缺点与不足,提供一种基于正交均匀线阵的水下波达方向估计方法及***,实现声波在水中传播速度不确定的情况下对目标进行精确定位,从而提高水下目标定位的精确性,且该方法分辨力高,运算量适中。
根据公开的实施例,本发明第一方面提出了一种基于正交均匀线阵的水下波达方向估计方法,所述估计方法包括下列步骤:
S1、建立具有正交均匀线阵的超声波接收探头阵列,所述正交均匀线阵为水平均匀线阵和竖直均匀线阵;
S2、推导所述水平均匀线阵的接收信号模型;
S3、计算所述水水平均匀线阵和所述水竖直均匀线阵的接收信号表达式;
S4、分别计算所述水平均匀线阵的两个协方差矩阵Rhxx和Rhxy以及所述竖直均匀线阵的两个协方差矩阵Rvxx和Rvxy,对上述协方差矩阵进行特征值分解,利用求出的最小特征值得出Chxx和Chxy、Cvxx和Cvxy,分别计算矩阵束{Chxx,Chxy}和{Cvxx,Cvxy}的广义特征值分解,得到特征值λx1,λx2,…,λxK和λy1,λy2,…,λyK;
S5、建立声波在所述水平均匀线阵以及所述竖直均匀线阵上的阵列方向角θxi和θyi之间的关系式;
S6、根据阵列方向角θxi和θyi之间的关系式,对所述水平均匀线阵的第一旋转算子以及所述竖直均匀线阵的第二旋转算子中的对角元素进行配对;
S7、根据配对结果,并联系对角元素与θxi之间的关系式,计算出θxi的表达式。
进一步地,所述步骤S1具体为:
在水中放置两个正交的所述水平均匀线阵和所述竖直均匀线阵,其中,每个均匀线阵都有M个接收阵元,且阵元间间距为d,K个窄带目标声源分别为S1,S2,…,SK,中心频率为f,上述目标声源对应于所述水平均匀线阵的方向角分别为θx1,θx2,…,θxK,对应于所述竖直均匀线阵的方向角分别为θy1,θy2,…,θyK。
进一步地,所述步骤S2具体为:
S21、以第一个阵元为参考点,则第一个阵元接收的信号为:
其中si(t)表示第i个源信号,n1(t)表示第一个阵元上的噪声;
S22、第m个阵元在同一时刻接收到的信号为:
其中λi表示第i个目标源反射回来的声波波长,nm(t)表示第m个阵元上的噪声;
S23、将各阵元的接收信号排列成列向量形式,则所述水平均匀线阵接收的信号可用以下矢量式子表示:
X(t)=AS(t)+N(t)(公式1)
其中,为M×K的导向矢量矩阵,
X(t)=[x1(t),x2(t),…,xM(t)]T为M×1的接收信号矩阵,
S(t)=[s1(t),s2(t),…,sK(t)]T为K×1的源信号矩阵,
N(t)=[n1(t),n2(t),…,nM(t)]T为M×1的噪声矩阵。
进一步地,所述步骤S3中计算所述水平均匀线阵的接收信号表达式过程如下:
S301、将所述水平均匀线阵中的M个阵元分为两个平移矢量为d的子阵列Zhx和Zhy,
所述子阵列Zhx由所述水平均匀线阵的第一到第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)分别是所述子阵列Zhx上第一个阵元到第M-1个阵元接收到的信号;x1(t),x2(t),…,xM(t)分别是水平阵列上第一个阵元到第M个阵元。
所述子阵列Zhy由所述水平均匀线阵的第二到第M个阵元组成,则有:
yh1(t)=x2(t),yh2(t)=x3(t),…,yh(M-1)(t)=xM(t)
其中,yh1(t),yh2(t),…,yh(M-1)(t)分别是所述子阵列Zhy上第一个阵元到第M-1个阵元接收到的信号;
S302、将所述子阵列Zhx和Zhy中第m个阵元的接收信号分别表示为:
其中nhxm(t)和nhym(t)分别为子阵Zhx和Zhy上第m个阵元的加性噪声;
S303、将上式改写成矢量形式:
Xh(t)=AxS(t)+Nhx(t)
Yh(t)=AxΦxS(t)+Nhy(t)
其中
矩阵Φx为K×K的对角矩阵,又称第一旋转算子,其对角元素包含K个信号的波前在任意一个阵元偶之间的相位延迟信息,表示为:
所述步骤S3中计算所述竖直均匀线阵的接收信号表达式过程如下:
S3X1、将所述竖直均匀线阵中的M个阵元分为两个平移矢量为d的子阵列Zvx和Zvy,
所述子阵列Zvx由所述竖直均匀线阵的第一到第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)分别是所述子阵列Zvx上第一个到第M-1个阵元接收到的信号;y1(t),y2(t),…,yM(t)分别是竖直阵列上第一个阵元到第M个阵元。
所述子阵列Zvy由所述竖直均匀线阵的第二到第M个阵元组成,则有:
yv1(t)=y2(t),yv2(t)=y3(t),…,yv(M-1)(t)=yM(t)
其中,yv1(t),yv2(t),…,yv(M-1)(t)分别是所述子阵列Zvy上第一个阵元到第M-1个阵元接收到的信号;
S3X2、将所述子阵列Zvx和Zvy中第m个阵元的接收信号分别为:
其中nvxm(t)和nvym(t)分别为子阵Zvx和Zvy上第m个阵元的加性噪声;
S3X3、将上式写成矢量形式:
Xv(t)=AyS(t)+Nvx(t)
Yv(t)=AyΦyS(t)+Nvy(t)
其中
矩阵Φy为K×K的对角矩阵,又称第二旋转算子,其对角元素包含了K个信号的波前在任意一个阵元偶之间的相位延迟信息,表示为:
进一步地,所述步骤S4中计算矩阵束{Chxx,Chxy}的广义特征值分解过程如下:
S401、Xh(t)的协方差矩阵表示为:
Rhxx=E[Xh(t)Xh H(t)]=AxRssAH+σx 2I
其中Rss=E{S(t)SH(t)},为信源部分协方差矩阵;
S402、Xh(t)和Yh(t)的互协方差矩阵表示为:
Rhxy=E{Xh(t)Yh H(t)}=AxRssΦx HAx H+σx 2Z;
S403、对上述协方差矩阵以及互协方差矩阵进行特征值分解得到最小特征值为σ2,利用σ2可以得到矩阵束其中Chxy=Rhxy-σx 2Z=AxRssΦx HAx H;
S404、计算所述矩阵束{Chxx,Chxy}的广义特征值分解,得到非零特征值λx1,λx2,…,λxK,因此由公式(2)可记:
其中φxi为矩阵Φx上的对角元素,且φxi∈{λx1,λx2,…,λxK},i=1,2,…,K;
所述步骤S4中计算矩阵束{Cvxx,Cvxy}的广义特征值分解过程如下:
S4X1、Xv(t)的协方差矩阵表示为:
Rvxx=E[Xv(t)Xv H(t)]=AyRssAy H+σy 2I
其中Rss=E{S(t)SH(t)},为信源部分协方差矩阵;
S4X2、Xv(t)和Yv(t)的互协方差矩阵表示为:
Rvxy=E{Xv(t)Yv H(t)}=AyRssΦy HAy H+σy 2Z;
S4X3、对上述协方差矩阵以及互协方差矩阵进行特征值分解得到最小特征值σy 2,利用σy 2得到矩阵束y其中Cvxy=Rvxy-σy 2Z=AyRssΦy HAy H;
S4X4、计算所述矩阵束{Cvxx,Cvxy}的广义特征值分解,得到非零特征值λy1,λy2,…,λyK,因此由公式(3)可记:
其中φyi为矩阵Φy上的对角元素,且φyi∈{λy1,λy2,…,λyK},i=1,2,…,K。
进一步地,所述步骤S5具体为
步骤S5、建立声波在所述水平均匀线阵以及所述竖直均匀线阵上的阵列方向角θxi和θyi之间的关系式;
S51、分别以所述水平均匀线阵以及所述竖直均匀线阵所在的方向建立一个直角坐标系,其中,θv1和θh1分别为第一象限中声波与两个均匀线阵法线的夹角,θv2和θh2分别为第二象限中声波与两个均匀线阵法线的夹角;
S52、当声波从第一象限中入射时,此时有θxi=θh1,θyi=θv1,则θxi和θyi的关系为:
S53、当声波从第二象限中入射时,此时有θxi=-θh2,θyi=θv2,则θxi和θyi的关系为:
S54、根据公式(6)和公式(7)得到:
进一步地,所述步骤S6具体为
S61、将公式(8)带入公式(5),则有:
S62、根据公式(4)和公式(9)可知,若配对成功,则有以下式子成立:
其中,λi是第i个目标源反射回来的声波波长,且有λi=vi/f;
S63、将arg(λx1),arg(λx2),…,arg(λxK)按照各自的平方大小顺序从大到小排列得到序列Η;将arg(λy1),arg(λy2),…,arg(λyK)按照各自的平方大小顺序从小到大排列得到序列V,于是有:
其中hi为序列Η中的第i个元素;vi为序列V中的第i个元素。
进一步地,所述θxi的表达式;
当arg(φxi)大于0时,θxi取正值,即目标声源在x轴负方向;当arg(φ)xi小于0时,θxi取负值,即目标声源在x轴正方向。
根据公开的实施例,本发明第二方面提出了一种基于正交均匀线阵的水下波达方向估计***,所述估计***包括数据采集处理与控制模块、发射模块、接收模块、输出模块和电源模块,
其中,所述数据采集处理与控制模块包括处理器、D/A转换器和一个A/D转换器,并且分别与所述发射模块、所述接收模块、所述输出模块和所述电源模块相连,能够通过设定参数来控制所述发射模块,使所述发射模块发射指定的声波信号;还能够对所述接收模块传送回来的数据进行采集与处理,计算出精确的波达方向角;该模块还能够将处理完毕的数据传送给所述输出模块,然后通过USB接口输出到外部装置或者在显示器上显示出来。
所述发射模块包括至少一个超声波发送探头,并且与所述数据采集处理与控制模块和所述电源模块相连,能够根据所述数据采集处理与控制模块的指令,获取指定的信号并对其进行超声发射;
所述接收模块包括两个以相互垂直方式放置的超声波接收探头阵列,并与所述数据采集处理与控制模块和所述电源模块相连,能够接收目标声源发射回来的信号,并将其传送给所述数据采集处理与控制模块;
所述输出模块包括至少一个USB接口和至少一个显示器组成,并且与所述数据采集处理与控制模块模块和所述电源模块相连,能够提供人机交互,将所述数据采集处理与控制模块模块中处理好的数据通过USB接口输出到外部装置或者在显示器上显示出来;
所述电源模块与上述各模块相连,为上述各模块供电。
进一步地,所述数据采集处理与控制模块中处理器为TMS320VC5509A型号的DSP芯片,所述接收模块中超声波接收探头阵列为正交的均匀直线阵列,其中每个阵列包括多个超声接收探头,并且数量相同。
本发明相对于现有技术具有如下的优点及效果:
1、本发明相对于利用传统ESPRIT算法来进行水下目标定位的方法更具有实用性,并且精确度也更高。在复杂的水声环境中,声波的速度并不确定,即每个目标声源反射回来的声波速度大小不相同,而传统ESPRIT算法假定声速是不变且相同的,没有考虑声速的不确定性对估计算法的影响,这使得在水声环境下传统ESPRIT算法不具有实用性。本发明采用了两个正交的均匀线阵,对ESPRIT算法进行了改进,在算法中消去了声速,使得最后的估计结果和声速无关,从而排除了声速的不确定对波达方向估计的影响,使得其实用性更高,并且可在不同水声环境下进行测量,目标定位的精度也更高。
2、本发明对传统ESPRIT算法进行了改进,从而保留了ESPRIT算法分辨率高的优点,并且改进后算法的运算量和复杂度也没有增加,保证了算法的可行性。
3、本发明的估计***可行性强,安装简单。除此之外,现代处理器计算处理能力的不断提高,这使得本发明所使用的处理器等芯片的集成度高,并且计算能力强,从而保证了本发明的可行性。
附图说明
图1是本发明估计***的硬件结构模块图;
图2是本发明估计***所用的正交均匀线阵模型;
图3是本发明中水平均匀线阵的接收信号模型;
图4是本发明公开的基于正交均匀线阵的水下波达方向估计方法的流程图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚、明确,以下参照附图并举实施例对本发明进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
实施例一
本实施例提供的基于正交均匀线阵的水下波达方向估计方法,该估计方法的流程步骤图参照附图4所示,包含以下几个具体步骤:
步骤S1、建立正交均匀线阵的场景模型;
在水中放置两个正交的均匀线阵,其模型场景如图2所示。每个均匀线阵都有M个接收阵元,且阵元间间距为d,K个窄带目标声源分别为S1,S2,…,SK,中心频率为f,这些目标声源对应于水平线阵的方向角分别为θx1,θx2,…,θxK,对应于竖直线阵的方向角分别为θy1,θy2,…,θyK。
步骤S2、推导水平均匀线阵的接收信号模型;
水平均匀线阵的模型场景如图3所示。以第一个阵元为参考点,则第一个阵元接收的信号为:
其中si(t)表示第i个源信号,n1(t)表示第一个阵元上的噪声。
接收信号满足窄带条件,即当信号延迟远小于带宽倒数时,延迟作用相当于使基带信号产生一个相移。那么第m个阵元在同一时刻接收到的信号为:
其中λi表示第i个目标源反射回来的声波波长,nm(t)表示第m个阵元上的噪声。将各阵元的接收信号排列成列向量形式,则整个水平线阵接收的信号可用以下矢量式子表示:
X(t)=AS(t)+N(t)(公式1)
其中,为M×K的导向矢量矩阵,
X(t)=[x1(t),x2(t),…,xM(t)]T为M×1的接收信号矩阵,
S(t)=[s1(t),s2(t),…,sK(t)]T为K×1的源信号矩阵,
N(t)=[n1(t),n2(t),…,nM(t)]T为M×1的噪声矩阵。
步骤S3、计算水平阵列和竖直阵列的接收信号表达式;
将水平均匀线阵中的M个阵元分为两个平移矢量为d的子阵列Zhx和Zhy。子阵列Zhx由水平均匀线阵的第一到第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)分别是子阵列Zhx上第一个阵元到第M-1个阵元接收到的信号;x1(t),x2(t),…,xM(t)分别是水平阵列上第一个阵元到第M个阵元。
子阵列Zhy由水平阵列的第二到第M个阵元组成,则有:
yh1(t)=x2(t),yh2(t)=x3(t),…,yh(M-1)(t)=xM(t)
其中,yh1(t),yh2(t),…,yh(M-1)(t)分别是子阵列Zhy上第一个阵元到第M-1个阵元接收到的信号。
那么两个子阵列中第m个阵元的接收信号分别为:
其中nhxm(t)和nhym(t)分别为子阵Zhx和Zhy上第m个阵元的加性噪声。将上式写成矢量形式:
Xh(t)=AxS(t)+Nhx(t)
Yh(t)=AxΦxS(t)+Nhy(t)
其中
矩阵Φx为K×K的对角矩阵,它是把子阵Zhx和Zhy的输出联系起来的酉阵,也称第一旋转算子,其对角元素包含了K个信号的波前在任意一个阵元偶之间的相位延迟信息,表示为:
同样,将竖直均匀线阵中的M个阵元也分为两个平移矢量为d的子阵列Zvx和Zvy。子阵列Zvx由竖直均匀线阵的第一到第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)分别是子阵列Zvx上第一个阵元到第M-1个阵元接收到的信号。
子阵列Zvy由竖直阵列的第二到第M个阵元组成,则有:
yv1(t)=y2(t),yv2(t)=y3(t),…,yv(M-1)(t)=yM(t)
其中,yv1(t),yv2(t),…,yv(M-1)(t)分别是子阵列Zvy上第一个阵元到第M-1个阵元接收到的信号;y1(t),y2(t),…,yM(t)分别是竖直阵列上第一个阵元到第M个阵元。
那么两个子阵列中第m个阵元的接收信号分别为:
其中nvxm(t)和nvym(t)分别为子阵Zvx和Zvy上第m个阵元的加性噪声。将上式写成矢量形式:
Xv(t)=AyS(t)+Nvx(t)
Yv(t)=AyΦyS(t)+Nvy(t)
其中矩阵Φy为K×K的对角矩阵,它是把子阵Zvx和Zvy的输出联系起来的酉阵,也称第二旋转算子,其对角元素包含了K个信号的波前在任意一个阵元偶之间的相位延迟信息,表示为:
步骤S4、分别计算水平均匀线阵的两个协方差矩阵Rhxx和Rhxy以及竖直均匀线阵的两个协方差矩阵Rvxx和Rvxy,对上述矩阵协方差矩阵进行特征值分解,利用求出的最小特征值得出Chxx和Chxy、Cvxx和Cvxy,分别计算矩阵束{Chxx,Chxy}和{Cvxx,Cvxy}的广义特征值分解,得到特征值。
求解第一旋转算子Φx和第二旋转算子Φy。Xh(t)的协方差矩阵可以表示为:
Rhxx=E[Xh(t)Xh H(t)]=AxRssAx H+σx 2I
其中Rss=E{S(t)SH(t)},为信源部分协方差矩阵。
Xh(t)和Yh(t)的互协方差矩阵为:
Rhxy=E{Xh(t)Yh H(t)}=AxRssΦx HAx H+σx 2Z
对矩阵协方差矩阵进行特征值分解得到最小特征值为σ2,利用σ2可以得到矩阵束{Chxx,Chxy},其中Chxx=Rhxx-σx 2I=AxRssAx H,Chxy=Rhxy-σx 2Z=AxRssΦx HAx H。计算矩阵束{Chxx,Chxy}的广义特征值分解,得到非零特征值λx1,λx2,…,λxK,它们一一对应着矩阵Φx对角线上的元素,但对应关系并不确定,因此由公式(2)可记:
其中φxi为矩阵Φx上的对角元素,且φxi∈{λx1,λx2,…,λxK},i=1,2,…,K。
根据以上步骤,同理可以求得竖直均匀线阵的两个协方差矩阵Rvxx和Rvxy,Xv(t)的协方差矩阵可以表示为:
Rvxx=E[Xv(t)Xv H(t)]=AyRssAy H+σy 2I
其中Rss=E{S(t)SH(t)},为信源部分协方差矩阵。
Xv(t)和Yv(t)的互协方差矩阵为:
Rvxy=E{Xv(t)Yv H(t)}=AyRssΦy HAy H+σy 2Z
其中Rss=E{S(t)SH(t)},为信源部分协方差矩阵。Xv(t)和Yv(t)的互协方差矩阵为:对矩阵协方差矩阵进行特征值分解得到最小特征值σy 2,利用σy 2得到矩阵束{Cvxx,Cvxy},其中Cvxx=Rvxx-σy 2I=AyRssAy H,Cvxy=Rvxy-σy 2Z=AyRssΦy HAy H。
计算矩阵束{Cvxx,Cvxy}的广义特征值分解,得到非零特征值λy1,λy2,…,λyK,它们一一对应着矩阵Φy对角线上的元素,但对应关系并不确定,因此由公式(3)可记:
其中φyi为矩阵Φy上的对角元素,且φyi∈{λy1,λy2,…,λyK},i=1,2,…,K。
步骤S5、建立声波在水平均匀线阵以及竖直均匀线阵的阵列方向角θxi和θyi之间的关系式;
在图2中分别以这两个正交均匀线阵所在的方向可以建立一个直角坐标系。因为只是在水面的位置对水面下方的目标进行定位,因此可只考虑两个象限的声波入射情况。其中,θv1和θh1分别为第一象限①中声波与两个均匀线阵法线的夹角;θv2和θh2分别为第二象限②中声波与两个均匀线阵法线的夹角。
当声波从第一象限①中入射时,此时有θxi=θh1,θyi=θv1,则θxi和θyi的关系为:
当声波从第二象限②中入射时,由于处在x轴上的阵列信号是以处在x轴最负方向的阵元为参考阵元的,并且子阵Zhx也在子阵Zhy的负x轴方向。因此当声波从象限②中入射时,参考阵元是最晚接收到信号的,子阵Zhx中的阵元也比子阵Zhy中对应的阵元晚接收到信号,从而可以得到时延参数τ小于0,又因为所以此时有θxi=-θh2,θyi=θv2,则θxi和θyi的关系为:
根据公式(6)和公式(7)可以得到:
将公式(8)带入公式(5),则有:
步骤S6、根据阵列方向角θxi和θyi之间的关系式,对两个阵列的旋转算子中的对角元素进行配对;
对矩阵Φx和矩阵Φy上的对角元素φxi与φyi进行配对。
将公式(8)带入公式(5),则有:
根据公式(4)和公式(9)可知,若配对成功,则有以下式子成立:
其中,λi是第i个目标源反射回来的声波波长,且有由于每个目标源反射回来的声波速度vi差异不大,故对于每个目标源来说,arg2(φxi)与arg2(φyi)的和都是一定的。
将arg(λx1),arg(λx2),…,arg(λxK)按照各自的平方大小顺序从大到小排列得到序列Η;将arg(λy1),arg(λy2),…,arg(λyK)按照各自的平方大小顺序从小到大排列得到序列V。于是有:
其中hi为序列Η中的第i个元素;vi为序列V中的第i个元素。
步骤S7、根据配对结果,并联系对角元素与θxi之间的关系式,计算出θxi的表达式;
根据配对结果求出θxi的大小。
根据公式(4)和公式(9)可以约去λi得:
当arg(φxi)大于0时,θxi取正值,即目标声源在x轴负方向;当arg(φ)xi小于0时,θxi取负值,即目标声源在x轴正方向。
根据以上算法流程可知,本发明提出的改进算法不需要知道声速的大小就可以对θxi进行精确的估计,即可以在声速不确定的情况下估计出波达方向角θxi的值,克服了传统ESPRIT算法的缺点。
实施例二
本实施例提供的基于正交均匀线阵的水下波达方向估计***,其硬件结构模块组成如图1所示,包括数据采集处理与控制模块、发射模块、接收模块、输出模块和电源模块。
数据采集处理与控制模块由一个处理器、一个D/A转换器和一个A/D转换器组成,并且与发射模块、接收模块、输出模块和电源模块相连。它能够通过设定参数来控制发射模块,使发射模块发射指定的声波信号;还能够对接收模块传送回来的数据进行采集与处理,运用本发明算法计算出精确的波达方向角。除此之外,它还能够将处理好的数据传送给输出模块,使它们通过USB接口输出到外部装置或者在显示器上显示出来。
发射模块由一个超声波发送探头组成,并且与数据采集处理与控制模块和电源模块相连。它能够根据数据采集处理与控制模块的指令,获取指定的信号并对其进行超声发射;
接收模块由两个以相互垂直方式放置的超声波接收探头阵列组成,并与数据采集处理与控制模块和电源模块相连。它能够接收目标声源发射回来的信号,并将其传送给数据采集处理与控制模块。
输出模块由一个USB接口和一个显示器组成,并且与数据采集处理与控制模块模块和电源模块相连。它能够提供人机交互,将数据采集处理与控制模块模块中处理好的数据通过USB接口输出到外部装置或者在显示器上显示出来。
电源模块由一个电源组成,并且与数据采集处理与控制模块、发射模块、接收模块和输出模块相连。它能够为这些模块供电。
本发明装置的主要工作流程如下:操作者通过数据采集处理与控制模块输入对应的参数,使该模块产生所需要的数字信号,该模块中的D/A转换器将其转换成模拟信号后,传送给发射模块,超声波发射探头将对其进行发射。接收模块中的垂直阵列接收到信号后将其传送给数据采集处理与控制模块,信号经过该模块中的A/D转换器数字信号后变为数字信号,然后模块中的处理器利用本发明算法对其进行计算。最后数据采集处理与控制模块将计算得出的结果传给输出模块,使结果通过USB接口输出给外部装置或者在显示屏上显示出来。工作的全程所有模块都由电源模块供电。
实施例三
本实施例提供的基于正交均匀线阵的水下波达方向估计***,包括数据采集处理与控制模块、发射模块、接收模块、输出模块和电源模块。数据采集处理与控制模块采用DSP芯片实现(如:TI公司的TMS320VC5509A型号的DSP芯片),此DSP芯片可实现A/D转换和D/A转换的功能,并能够实现均匀线阵的旋转算子和最终波达方向的计算;发射模块使用一个超声波发射探头;接收模块使用两个正交的均匀直线阵列,其中每个阵列包括多个超声接收探头,并且数量相同;输出模块使用一个USB接口和一个LCD显示屏。图1即为本实施例所述装置的硬件结构模块图。
本发明的主要工作步骤具体如下:
步骤1:在水下放置4个目标声源,与水平阵列的方向角分别为35°,60°,-30°和-70°。设定阵列接收模块和发射模块中的各项参数。将接收模块中的每个均匀线阵中的阵元个数M定为9,则每个均匀线阵中前8个阵元为子阵Zx,后8个阵元为子阵Zy。两个子阵列之间的距离d设为5mm。通过数据采集处理与控制模块控制发射模块使待发射的超声频率为100kHz,脉冲长度5ms。
步骤2:对超声接收探头接收到的目标声源信号进行采样;水平方向均匀线阵接收到的信号分别为xx1(t),xx2(t),…,xx8(t)和yx1(t),yx2(t),…,yx8(t),竖直方向均匀线阵接收的信号分别为xy1(t),xy2(t),…,xy8(t)和yy1(t),yy2(t),…,yy8(t)。共采样接收200次,并将接收到的信号传递给数据采集处理与控制模块模块做运算处理。
步骤3:信号在数据采集处理与控制模块中的处理步骤具体如下:
31)将处在水平方向上的均匀线阵接收到的信号排成矢量形式Xh(t)和Yh(t),计算Xh(t)的协方差矩阵Rhxx=E[Xh(t)Xh H(t)],Xh(t)和Yh(t)之间的互协方差矩阵Rhxy=E{Xh(t)Yh H(t)}。同时对竖直方向上的均匀线阵接收到的信号也进行相同处理,得到Rvxx=E[Xv(t)Xv H(t)]和Rvxy=E{Xv(t)Yv H(t)}。
32)对水平阵列中的两个协方差矩阵Rhxx和Rhxy进行特征值分解,得到最小的特征值σ2,从而有Chxx=Rhxx-σ2I=ARssAH和Chxy=Rhxy-σ2Z=ARssΦHAH。同时对竖直阵列中的两个协方差矩阵进行相同的处理,得到Cvxx和Cvxy。
33)分别计算矩阵束{Chxx,Chxy}和{Cvxx,Cvxy}的广义特征值分解,得到λx1,λx2,…,λxK和λy1,λy2,…,λyK。
34)将arg(λx1),arg(λx2),…,arg(λxK)按照各自的平方大小顺序从大到小排列得到序列Η,将arg(λy1),arg(λy2),…,arg(λyK)按照各自的平方大小顺序从小到大排列得到序列V。然后把H中第i个元素hi的值赋给arg(φxi),V中第i个元素vi的值赋给arg(φyi)。
34)根据匹配得出的arg(φxi)和arg(φyi),联立式和式对θxi进行求解,最终求得
步骤4:将计算出的方向角信息存储下来,并传送给输出模块,使其通过USB接口输出给外部装置或者显示在LCD显示屏上。根据本发明算法估计出的方向角分别为35°,60°,-30°和-70°,与实际角度相同,说明估计结果正确,本发明方法及***可行。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (10)
1.一种基于正交均匀线阵的水下波达方向估计方法,其特征在于,所述估计方法包括下列步骤:
S1、建立具有正交均匀线阵的超声波接收探头阵列,所述正交均匀线阵为水平均匀线阵和竖直均匀线阵;
S2、推导所述水平均匀线阵的接收信号模型;
S3、计算所述水平均匀线阵和所述竖直均匀线阵的接收信号表达式;
S4、分别计算所述水平均匀线阵的两个协方差矩阵Rhxx和Rhxy以及所述竖直均匀线阵的两个协方差矩阵Rvxx和Rvxy,对上述协方差矩阵进行特征值分解,利用求出的最小特征值得出Chxx和Chxy、Cvxx和Cvxy,分别计算矩阵束{Chxx,Chxy}和{Cvxx,Cvxy}的广义特征值分解,得到特征值λx1,λx2,…,λxK和λy1,λy2,…,λyK;
S5、建立声波在所述水平均匀线阵以及所述竖直均匀线阵上的阵列方向角θxi和θyi之间的关系式;
S6、根据阵列方向角θxi和θyi之间的关系式,对所述水平均匀线阵的第一旋转算子以及所述竖直均匀线阵的第二旋转算子中的对角元素进行配对;
S7、根据配对结果,并联系对角元素与θxi之间的关系式,计算出θxi的表达式。
2.根据权利要求1所述的一种基于正交均匀线阵的水下波达方向估计方法,其特征在于,所述步骤S1具体为:
在水中放置两个正交的所述水平均匀线阵和所述竖直均匀线阵,其中,每个均匀线阵都有M个接收阵元,且阵元间间距为d,K个窄带目标声源分别为S1,S2,…,SK,中心频率为f,上述目标声源对应于所述水平均匀线阵的方向角分别为θx1,θx2,…,θxK,对应于所述竖直均匀线阵的方向角分别为θy1,θy2,…,θyK。
3.根据权利要求2所述的一种基于正交均匀线阵的水下波达方向估计方法,其特征在于,所述步骤S2具体为:
S21、以第一个阵元为参考点,则第一个阵元接收的信号为:
其中si(t)表示第i个源信号,n1(t)表示第一个阵元上的噪声;
S22、第m个阵元在同一时刻接收到的信号为:
其中λi表示第i个目标源反射回来的声波波长,nm(t)表示第m个阵元上的噪声;
S23、将各阵元的接收信号排列成列向量形式,则所述水平均匀线阵接收的信号可用以下矢量式子表示:
X(t)=AS(t)+N(t)(公式1)
其中,为M×K的导向矢量矩阵,
X(t)=[x1(t),x2(t),…,xM(t)]T为M×1的接收信号矩阵,
S(t)=[s1(t),s2(t),…,sK(t)]T为K×1的源信号矩阵,
N(t)=[n1(t),n2(t),…,nM(t)]T为M×1的噪声矩阵。
4.根据权利要求3所述的一种基于正交均匀线阵的水下波达方向估计方法,其特征在于,所述步骤S3中计算所述水平均匀线阵的接收信号表达式过程如下:
S301、将所述水平均匀线阵中的M个阵元分为两个平移矢量为d的子阵列Zhx和Zhy,
所述子阵列Zhx由所述水平均匀线阵的第一到第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)分别是所述子阵列Zhx上第一个阵元到第M-1个阵元接收到的信号;x1(t),x2(t),…,xM(t)分别是水平阵列上第一个阵元到第M个阵元;
所述子阵列Zhy由所述水平均匀线阵的第二到第M个阵元组成,则有:
yh1(t)=x2(t),yh2(t)=x3(t),…,yh(M-1)(t)=xM(t)
其中,yh1(t),yh2(t),…,yh(M-1)(t)分别是所述子阵列Zhy上第一个阵元到第M-1个阵元接收到的信号;
S302、将所述子阵列Zhx和Zhy中第m个阵元的接收信号分别表示为:
其中nhxm(t)和nhym(t)分别为子阵Zhx和Zhy上第m个阵元的加性噪声;
S303、将上式改写成矢量形式:
Xh(t)=AxS(t)+Nhx(t)
Yh(t)=AxΦxS(t)+Nhy(t)
其中
矩阵Φx为K×K的对角矩阵,又称第一旋转算子,其对角元素包含K个信号的波前在任意一个阵元偶之间的相位延迟信息,表示为:
所述步骤S3中计算所述竖直均匀线阵的接收信号表达式过程如下:
S3X1、将所述竖直均匀线阵中的M个阵元分为两个平移矢量为d的子阵列Zvx和Zvy,
所述子阵列Zvx由所述竖直均匀线阵的第一到第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)分别是所述子阵列Zvx上第一个到第M-1个阵元接收到的信号;y1(t),y2(t),…,yM(t)分别是竖直阵列上第一个阵元到第M个阵元;
所述子阵列Zvy由所述竖直均匀线阵的第二到第M个阵元组成,则有:
yv1(t)=y2(t),yv2(t)=y3(t),…,yv(M-1)(t)=yM(t)
其中,yv1(t),yv2(t),…,yv(M-1)(t)分别是所述子阵列Zvy上第一个阵元到第M-1个阵元接收到的信号;
S3X2、将所述子阵列Zvx和Zvy中第m个阵元的接收信号分别为:
其中nvxm(t)和nvym(t)分别为子阵Zvx和Zvy上第m个阵元的加性噪声;
S3X3、将上式写成矢量形式:
Xv(t)=AyS(t)+Nvx(t)
Yv(t)=AyΦyS(t)+Nvy(t)
其中
矩阵Φy为K×K的对角矩阵,又称第二旋转算子,其对角元素包含了K个信号的波前在任意一个阵元偶之间的相位延迟信息,表示为:
5.根据权利要求4所述的一种基于正交均匀线阵的水下波达方向估计方法,其特征在于,所述步骤S4中计算矩阵束{Chxx,Chxy}的广义特征值分解过程如下:
S401、Xh(t)的协方差矩阵表示为:
Rhxx=E[Xh(t)Xh H(t)]=AxRssAx H+σx 2I
其中Rss=E{S(t)SH(t)},为信源部分协方差矩阵;
S402、Xh(t)和Yh(t)的互协方差矩阵表示为:
Rhxy=E{Xh(t)Yh H(t)}=AxRssΦx HAx H+σx 2Z;
S403、对上述协方差矩阵以及互协方差矩阵进行特征值分解得到最小特征值为σx 2,利用σx 2可以得到矩阵束{Chxx,Chxy},其中Chxy=Rhxy-σx 2Z=AxRssΦx HAx H;
S404、计算所述矩阵束{Chxx,Chxy}的广义特征值分解,得到非零特征值λx1,λx2,…,λxK,因此由公式(2)可记:
其中φxi为矩阵Φx上的对角元素,且φxi∈{λx1,λx2,…,λxK},i=1,2,…,K;
所述步骤S4中计算矩阵束{Cvxx,Cvxy}的广义特征值分解过程如下:
S4X1、Xv(t)的协方差矩阵表示为:
Rvxx=E[Xv(t)Xv H(t)]=AyRssAy H+σy 2I
其中Rss=E{S(t)SH(t)},为信源部分协方差矩阵;
S4X2、Xv(t)和Yv(t)的互协方差矩阵表示为:
Rvxy=E{Xv(t)Yv H(t)}=AyRssΦy HAy H+σy 2Z;
S4X3、对上述协方差矩阵以及互协方差矩阵进行特征值分解得到最小特征值σy 2,利用σy 2得到矩阵束其中Cvxy=Rvxy-σy 2Z=AyRssΦy HAy H;
S4X4、计算所述矩阵束{Cvxx,Cvxy}的广义特征值分解,得到非零特征值λy1,λy2,…,λyK,因此由公式(3)可记:
其中φyi为矩阵Φy上的对角元素,且φyi∈{λy1,λy2,…,λyK},i=1,2,…,K。
6.根据权利要求5所述的一种基于正交均匀线阵的水下波达方向估计方法,其特征在于,所述步骤S5具体为:
S51、分别以所述水平均匀线阵以及所述竖直均匀线阵所在的方向建立一个直角坐标系,其中,θv1和θh1分别为第一象限中声波与两个均匀线阵法线的夹角,θv2和θh2分别为第二象限中声波与两个均匀线阵法线的夹角;
S52、当声波从第一象限中入射时,此时有θxi=θh1,θyi=θv1,则θxi和θyi的关系为:
S53、当声波从第二象限中入射时,此时有θxi=-θh2,θyi=θv2,则θxi和θyi的关系为:
S54、根据公式(6)和公式(7)得到:
7.根据权利要求6所述的一种基于正交均匀线阵的水下波达方向估计方法,其特征在于,所述步骤S6具体为:
S61、将公式(8)带入公式(5),则有:
S62、根据公式(4)和公式(9)可知,若配对成功,则有以下式子成立:
其中,λi是第i个目标源反射回来的声波波长,且有λi=vi/f;
S63、将arg(λx1),arg(λx2),…,arg(λxK)按照各自的平方大小顺序从大到小排列得到序列Η;将arg(λy1),arg(λy2),…,arg(λyK)按照各自的平方大小顺序从小到大排列得到序列V,于是有:
其中hi为序列Η中的第i个元素;vi为序列V中的第i个元素。
8.根据权利要求7所述的一种基于正交均匀线阵的水下波达方向估计方法,其特征在于,所述θxi的表达式;
当arg(φxi)大于0时,θxi取正值,即目标声源在x轴负方向;当arg(φx)i小于0时,θxi取负值,即目标声源在x轴正方向。
9.一种基于正交均匀线阵的水下波达方向估计***,其特征在于,所述估计***包括数据采集处理与控制模块、发射模块、接收模块、输出模块和电源模块,
其中,所述数据采集处理与控制模块包括处理器、D/A转换器和一个A/D转换器,并且分别与所述发射模块、所述接收模块、所述输出模块和所述电源模块相连,能够通过设定参数来控制所述发射模块,使所述发射模块发射指定的声波信号;还能够对所述接收模块传送回来的数据进行采集与处理,计算出精确的波达方向角;该模块还能够将处理完毕的数据传送给所述输出模块,然后通过USB接口输出到外部装置或者在显示器上显示出来;
所述发射模块包括至少一个超声波发送探头,并且与所述数据采集处理与控制模块和所述电源模块相连,能够根据所述数据采集处理与控制模块的指令,获取指定的信号并对其进行超声发射;
所述接收模块包括两个以相互垂直方式放置的超声波接收探头阵列,并与所述数据采集处理与控制模块和所述电源模块相连,能够接收目标声源发射回来的信号,并将其传送给所述数据采集处理与控制模块;
所述输出模块包括至少一个USB接口和至少一个显示器组成,并且与所述数据采集处理与控制模块模块和所述电源模块相连,能够提供人机交互,将所述数据采集处理与控制模块模块中处理好的数据通过USB接口输出到外部装置或者在显示器上显示出来;
所述电源模块与上述各模块相连,为上述各模块供电。
10.根据权利要求9所述的一种基于正交均匀线阵的水下波达方向估计***,其特征在于,所述数据采集处理与控制模块中处理器为TMS320VC5509A型号的DSP芯片,所述接收模块中超声波接收探头阵列为正交的均匀直线阵列,其中每个阵列包括多个超声接收探头,并且数量相同。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610126961.6A CN105607033B (zh) | 2016-03-07 | 2016-03-07 | 基于正交均匀线阵的水下波达方向估计方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610126961.6A CN105607033B (zh) | 2016-03-07 | 2016-03-07 | 基于正交均匀线阵的水下波达方向估计方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105607033A true CN105607033A (zh) | 2016-05-25 |
CN105607033B CN105607033B (zh) | 2019-01-18 |
Family
ID=55987111
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610126961.6A Expired - Fee Related CN105607033B (zh) | 2016-03-07 | 2016-03-07 | 基于正交均匀线阵的水下波达方向估计方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105607033B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106441543A (zh) * | 2016-12-09 | 2017-02-22 | 华南理工大学 | 基于三维正交阵的水下探测路径声速测量方法及装置 |
CN106500820A (zh) * | 2016-10-13 | 2017-03-15 | 华南理工大学 | 一种用于二维波达方向估计的声速测量方法及装置 |
CN107414831A (zh) * | 2017-08-01 | 2017-12-01 | 广东工业大学 | 一种机械臂关节坐标计算***及方法 |
CN107942284A (zh) * | 2017-10-31 | 2018-04-20 | 华南理工大学 | 基于二维正交非均匀线阵的水下波达方向估计方法与装置 |
CN108008348A (zh) * | 2017-11-16 | 2018-05-08 | 华南理工大学 | 基于可调夹角均匀线阵的水下波达方向估计方法及装置 |
CN108398659A (zh) * | 2018-02-11 | 2018-08-14 | 华南理工大学 | 一种矩阵束与求根music结合的波达方向估计方法 |
CN108414967A (zh) * | 2018-04-11 | 2018-08-17 | 华南理工大学 | 基于夹角可调双l阵的水下二维波达方向估计方法与装置 |
CN108519576A (zh) * | 2018-03-21 | 2018-09-11 | 华南理工大学 | 基于夹角可调非均匀线阵的水下波达方向估计方法与装置 |
CN108802735A (zh) * | 2018-06-15 | 2018-11-13 | 华南理工大学 | 一种用于未知声速环境的水下目标定位及测速方法和装置 |
CN109407048A (zh) * | 2018-11-12 | 2019-03-01 | 华南理工大学 | 基于非圆信号和夹角可调阵的水下doa估计方法与装置 |
CN109521392A (zh) * | 2018-10-24 | 2019-03-26 | 华南理工大学 | 基于非圆信号和l型线阵的水下一维doa估计方法和装置 |
CN109581274A (zh) * | 2018-12-13 | 2019-04-05 | 华南理工大学 | 基于夹角可调三维阵的非圆信号水下doa估计方法和装置 |
CN109581275A (zh) * | 2018-12-13 | 2019-04-05 | 华南理工大学 | 基于非圆信号和三维正交阵的二维水下doa估计方法和装置 |
CN109884580A (zh) * | 2019-02-22 | 2019-06-14 | 华南理工大学 | 水下一维doa估计方法和装置 |
CN110133594A (zh) * | 2018-02-09 | 2019-08-16 | 北京搜狗科技发展有限公司 | 一种声源定位方法、装置和用于声源定位的装置 |
CN111427044A (zh) * | 2020-04-19 | 2020-07-17 | 中国电子科技集团公司电子科学研究院 | 水域目标的定位***及定位方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140269200A1 (en) * | 2013-03-12 | 2014-09-18 | Teledyne Instruments, Inc. D/B/A Teledyne Benthos | Passive acoustic bearing estimation via ultra short baseline wideband methods |
CN104133217A (zh) * | 2014-07-17 | 2014-11-05 | 华南理工大学 | 一种水下运动目标与水流的三维速度联合测定方法及装置 |
CN104931929A (zh) * | 2015-06-11 | 2015-09-23 | 华南理工大学 | 基于线阵综合声速补偿的近场波达方向估计方法及装置 |
CN104977576A (zh) * | 2015-06-15 | 2015-10-14 | 西安交通大学 | 基于十字阵雷达***的分布式紧凑感知矩阵追踪方法 |
CN105004413A (zh) * | 2015-06-11 | 2015-10-28 | 华南理工大学 | 用于水下目标定位的声传播路径综合速度测定方法与装置 |
-
2016
- 2016-03-07 CN CN201610126961.6A patent/CN105607033B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140269200A1 (en) * | 2013-03-12 | 2014-09-18 | Teledyne Instruments, Inc. D/B/A Teledyne Benthos | Passive acoustic bearing estimation via ultra short baseline wideband methods |
CN104133217A (zh) * | 2014-07-17 | 2014-11-05 | 华南理工大学 | 一种水下运动目标与水流的三维速度联合测定方法及装置 |
CN104931929A (zh) * | 2015-06-11 | 2015-09-23 | 华南理工大学 | 基于线阵综合声速补偿的近场波达方向估计方法及装置 |
CN105004413A (zh) * | 2015-06-11 | 2015-10-28 | 华南理工大学 | 用于水下目标定位的声传播路径综合速度测定方法与装置 |
CN104977576A (zh) * | 2015-06-15 | 2015-10-14 | 西安交通大学 | 基于十字阵雷达***的分布式紧凑感知矩阵追踪方法 |
Non-Patent Citations (1)
Title |
---|
甄佳奇等: ""基于正交直线阵列的二维相干源测向方法"", 《吉林大学学报(工学版)》 * |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106500820A (zh) * | 2016-10-13 | 2017-03-15 | 华南理工大学 | 一种用于二维波达方向估计的声速测量方法及装置 |
CN106441543A (zh) * | 2016-12-09 | 2017-02-22 | 华南理工大学 | 基于三维正交阵的水下探测路径声速测量方法及装置 |
CN107414831A (zh) * | 2017-08-01 | 2017-12-01 | 广东工业大学 | 一种机械臂关节坐标计算***及方法 |
CN107942284A (zh) * | 2017-10-31 | 2018-04-20 | 华南理工大学 | 基于二维正交非均匀线阵的水下波达方向估计方法与装置 |
CN107942284B (zh) * | 2017-10-31 | 2020-02-18 | 华南理工大学 | 基于二维正交非均匀线阵的水下波达方向估计方法与装置 |
CN108008348A (zh) * | 2017-11-16 | 2018-05-08 | 华南理工大学 | 基于可调夹角均匀线阵的水下波达方向估计方法及装置 |
WO2019095912A1 (zh) * | 2017-11-16 | 2019-05-23 | 华南理工大学 | 基于可调夹角均匀线阵的水下波达方向估计方法及装置 |
CN110133594A (zh) * | 2018-02-09 | 2019-08-16 | 北京搜狗科技发展有限公司 | 一种声源定位方法、装置和用于声源定位的装置 |
CN108398659A (zh) * | 2018-02-11 | 2018-08-14 | 华南理工大学 | 一种矩阵束与求根music结合的波达方向估计方法 |
CN108398659B (zh) * | 2018-02-11 | 2020-06-19 | 华南理工大学 | 一种矩阵束与求根music结合的波达方向估计方法 |
CN108519576A (zh) * | 2018-03-21 | 2018-09-11 | 华南理工大学 | 基于夹角可调非均匀线阵的水下波达方向估计方法与装置 |
CN108519576B (zh) * | 2018-03-21 | 2021-07-20 | 华南理工大学 | 基于夹角可调非均匀线阵的水下波达方向估计方法与装置 |
CN108414967A (zh) * | 2018-04-11 | 2018-08-17 | 华南理工大学 | 基于夹角可调双l阵的水下二维波达方向估计方法与装置 |
CN108802735A (zh) * | 2018-06-15 | 2018-11-13 | 华南理工大学 | 一种用于未知声速环境的水下目标定位及测速方法和装置 |
CN108802735B (zh) * | 2018-06-15 | 2020-08-18 | 华南理工大学 | 一种用于未知声速环境的水下目标定位及测速方法和装置 |
CN109521392A (zh) * | 2018-10-24 | 2019-03-26 | 华南理工大学 | 基于非圆信号和l型线阵的水下一维doa估计方法和装置 |
CN109521392B (zh) * | 2018-10-24 | 2022-12-16 | 华南理工大学 | 基于非圆信号和l型线阵的水下一维doa估计方法和装置 |
CN109407048A (zh) * | 2018-11-12 | 2019-03-01 | 华南理工大学 | 基于非圆信号和夹角可调阵的水下doa估计方法与装置 |
CN109407048B (zh) * | 2018-11-12 | 2023-09-29 | 华南理工大学 | 基于非圆信号和夹角可调阵的水下doa估计方法与装置 |
CN109581275A (zh) * | 2018-12-13 | 2019-04-05 | 华南理工大学 | 基于非圆信号和三维正交阵的二维水下doa估计方法和装置 |
CN109581275B (zh) * | 2018-12-13 | 2023-09-29 | 华南理工大学 | 基于非圆信号和三维正交阵的二维水下doa估计方法和装置 |
CN109581274A (zh) * | 2018-12-13 | 2019-04-05 | 华南理工大学 | 基于夹角可调三维阵的非圆信号水下doa估计方法和装置 |
CN109581274B (zh) * | 2018-12-13 | 2023-09-29 | 华南理工大学 | 基于夹角可调三维阵的非圆信号水下doa估计方法和装置 |
CN109884580A (zh) * | 2019-02-22 | 2019-06-14 | 华南理工大学 | 水下一维doa估计方法和装置 |
CN111427044A (zh) * | 2020-04-19 | 2020-07-17 | 中国电子科技集团公司电子科学研究院 | 水域目标的定位***及定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105607033B (zh) | 2019-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105607033A (zh) | 基于正交均匀线阵的水下波达方向估计方法及*** | |
CN108008348B (zh) | 基于可调夹角均匀线阵的水下波达方向估计方法及装置 | |
CA2941477C (en) | Adaptive beamformer for sonar imaging | |
Hawkes et al. | Effects of sensor placement on acoustic vector-sensor array performance | |
CN106500820B (zh) | 一种用于二维波达方向估计的声速测量方法及装置 | |
CN107942284B (zh) | 基于二维正交非均匀线阵的水下波达方向估计方法与装置 | |
CN101644765B (zh) | 一种用于水声换能器线阵的幅相误差校正方法 | |
CN109581275B (zh) | 基于非圆信号和三维正交阵的二维水下doa估计方法和装置 | |
CN108535682B (zh) | 一种基于旋转非均匀双l阵的水下二维doa估计方法及装置 | |
CN113126028B (zh) | 一种基于多个麦克风阵列的噪声源定位方法 | |
CN106249225B (zh) | 稀疏圆形声矢量传感器阵列四元数esprit参数估计方法 | |
Najeem et al. | Open lake experiment for direction of arrival estimation using acoustic vector sensor array | |
CN106249196B (zh) | 三分量声矢量传感器稀疏阵列四元数解模糊方法 | |
Wang et al. | Array calibration of angularly dependent gain and phase uncertainties with carry-on instrumental sensors | |
CN109725285B (zh) | 一种基于mvdr协方差矩阵元素自适应相角转换的doa估计方法 | |
CN108845325A (zh) | 拖曳线列阵声纳子阵误差失配估计方法 | |
CN104931929B (zh) | 基于线阵综合声速补偿的近场波达方向估计方法及装置 | |
CN110109053B (zh) | 一种未知声速环境下快速doa估计方法 | |
CN106441543A (zh) | 基于三维正交阵的水下探测路径声速测量方法及装置 | |
Sun et al. | Array geometry calibration for underwater compact arrays | |
CN109521392B (zh) | 基于非圆信号和l型线阵的水下一维doa估计方法和装置 | |
CN109581274B (zh) | 基于夹角可调三维阵的非圆信号水下doa估计方法和装置 | |
CN109407048B (zh) | 基于非圆信号和夹角可调阵的水下doa估计方法与装置 | |
CN106054122B (zh) | 基于数字信号处理机的时域宽带信号频域闭环测向方法 | |
JPH07504505A (ja) | 受動音響ブイより成るアンテナによって雑音源の位置を決定する方法および装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190118 |
|
CF01 | Termination of patent right due to non-payment of annual fee |