CN105138743A - 一种基于单矢量水听器的相干信号源方位估计方法 - Google Patents
一种基于单矢量水听器的相干信号源方位估计方法 Download PDFInfo
- Publication number
- CN105138743A CN105138743A CN201510477827.6A CN201510477827A CN105138743A CN 105138743 A CN105138743 A CN 105138743A CN 201510477827 A CN201510477827 A CN 201510477827A CN 105138743 A CN105138743 A CN 105138743A
- Authority
- CN
- China
- Prior art keywords
- theta
- coherent
- max
- objective function
- coherent signal
- 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
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于单矢量水听器的相干信号源方位估计方法。包括以下步骤:确定三个通道中两个相干信号的加权矩阵,根据加权矩阵将接收信号分解成两个子信号;确定子信号的MUSIC谱为目标函数;根据目标函数进行最优搜索,确定两个相干信号源的方位角。本发明通过确定三个通道中两个相干信号的加权矩阵,实现两个相干信号源的分离,解决了单只矢量水听器两个相干信号源条件下,信号方位合并为一个的问题,实现了单只矢量水听器两个相干信号源的方位估计。
Description
技术领域
本发明属于基阵信号处理领域。尤其涉及一种基于单矢量水听器的相干信号源方位估计方法。
背景技术
与声压水听器相比,矢量水听器可以同时、共点测量声场的声压与振速,获得更多的声场信息,依靠单个矢量水听器就可以完成水下目标的方位估计,且不存在多值模糊的问题。因此在许多安装空间有限的应用场所,采取单只矢量水听器作为声呐的接收器,是目前水声领域比较流行而且非常实用的工程应用方式。
相干信号源的估计问题始终是空间谱估计中的研究热点问题。对声压阵或声矢量阵而言,当信号源完全相干时,阵列接收的数据协方差矩阵的秩降为1,这会导致某些相干源的导向矢量与噪声子空间不完全正交,从而导致常规的MUSIC算法无法正确估计信号源方位;而对单只矢量水听器而言,当信号源完全相干时,两个相干信号的方向矢量会合并为一个,常规的空间谱估计方法估计出的信号方位为合成的信号方位,不能正确估计两个相干信号源的方位。因此,在相干信号源情况下正确估计信号方向的关键问题,一方面是如何通过预处理使得信号协方差矩阵的秩得到有效恢复,从而正确估计信号源的方向;另一方面是如何通过预处理将相干信号源分开,再对其进行方位估计。目前关于解相干的处理主要是通过恢复协方差矩阵的秩实现的,基本有两大类:一类是降维处理,如基于空间平滑、基于矩阵重构两类算法,这类算法的解相干性能是通过降低自由度换取的,不适合单只矢量水听器相干信号源方位估计;另一类是非降维处理,主要有频域平滑算法、Toeplitz方法、虚拟阵列变换法等。这类算法与降维算法相比最大的优点在于阵列的自由度不变,但这类算法往往针对的是特定环境,也不适合单只矢量水听器相干信号源方位估计。
发明内容
本发明的目的是提供一种操作方法简单、效率高的,基于单矢量水听器的相干信号源方位估计方法。
一种基于单矢量水听器的相干信号源方位估计方法,包括以下步骤,
步骤一:确定声压P通道、振速Vx通道和振速Vy通道中两个相干信号的加权矩阵T,根据加权矩阵T将接收数据X(t)分成两个子信号X1(t)和X2(t),其中X1(t)=TX(t);
步骤二:确定子信号的MUSIC谱为目标函数,子信号的MUSIC谱为:
其中: 或 是第i个水声信号源在单只矢量水听器上的方向矢量,i≤2,和为搜索角度,是搜索到的两个相干信号的幅度比,是对原始接收数据协方差矩阵进行特征分解得到噪声子空间;
步骤三:根据目标函数进行最优搜索,找到目标函数最大值对应的搜索角度和即为两个相干信号源的入射角度θ1和θ2。
本发明一种基于单矢量水听器的相干信号源方位估计方法,还可以包括:
1、目标函数进行最优搜索的过程为:
(1)确定搜索范围并设置初始参数根据初始参数对接收数据进行分解得到子信号X1(t),进而求出初始参数条件下的目标函数值,定义该目标函数值为最大值,并将初始参数值定义为估计到的参数值,即
(2)重置参数求出在新的参数条件下的目标函数
(3)比较fmax和fnew的大小,将其中大的一个定义为新的目标函数的最大值,并将其对应的参数定义为新的估计到的参数值,即:
fmax=max(fmax,fnew)
<θ1,θ2,k>=<θ1,θ2,k>|max(fmax,fnew)
(4)重复步骤(2)~(3),直到将搜索范围内的所有参数全部搜索完毕。
有益效果:
本发明不同于传统的通过恢复信号协方差矩阵的秩来实现解相干的算法,本发明方法提出了单只矢量水听器三个通道中两个相干信号的加权矩阵,根据加权矩阵实现接收数据中两个相干信号的分解,解决了单只矢量水听器两个相干信号源条件下,信号方位合并为一个的问题,实现了单只矢量水听器两个相干信号源的方位估计。
附图说明
图1为本发明的方法流程图。
具体实施方式
下面将结合附图对本发明作进一步详细说明。
本发明方法从另一方面考虑,不同于传统的通过恢复信号协方差矩阵的秩来实现解相干的算法,而是从接收数据入手,先确定三个通道中两个相干信号的加权矩阵,根据加权矩阵将接收信号分解成两个子信号,再对其中一个子信号进行信号协方差矩阵的构造和方位估计。由于两个相干信号的加权矩阵包含了两个信号的方位信息,因此可实现两个相干信号源的方位估计。
本发明的目的是这样实现的:
步骤S1:确定三个通道中两个相干信号的加权矩阵,根据加权矩阵将接收信号分解成两个子信号;
步骤S2:确定子信号的MUSIC谱为目标函数;
步骤S3:根据目标函数进行最优搜索,确定两个相干信号源的方位角;
本发明的核心技术在于提出了三个通道中两个相干信号的加权矩阵,并根据加权矩阵将接收信号分解成两个子信号。两个相干信号的加权矩阵包含了两个相干信号的方位角和两个相干信号的幅度比,当搜索的角度和分别等于两个相干信号源的方位角θ1和θ2,且搜索到的两个相干信号的幅度比等于真实的两个相干信号的幅度比k时,则实现对接收数据中两个相干信号源的正确分解,进而子信号的MUSIC谱会产生一个最大值。因此定义子信号的MUSIC谱为目标函数,找出目标函数最大值点对应的两个搜索角度即为两个相干信号源的方位角。
下面将结合附图1对本发明加以详细说明:
步骤S1:两个相干信号源s1(t)、s2(t),分别从不同角度θ1、θ2入射到单只矢量水听器上时,由于s1(t)、s2(t)为相干信号,满足s2(t)=ks1(t),其中k为常数,因此两个信号的方向矢量在单只矢量水听器接收数据模型中将会合并为一个,也就是说,两个信号的方位合并为一个。此时若用常规的MUSIC算法估计其方位,得到的是两个信号的合成方位,不能正确估计出两个相干信号源的方位。
本发明方法通过提出的三个通道中两个相干信号的加权矩阵,实现两个相干信号源的分离。具体的说就是先确定声压P通道、振速Vx通道、振速Vy通道中两个相干信号的加权矩阵T,根据加权矩阵T将接收数据X(t)分成两个子信号X1(t)和X2(t),如X1(t)表达式可记为:
X1(t)=TX(t)
T是相干信号s1(t)在各个通道中的加权矩阵,包含了两个相干信号的方位角和两个相干信号的幅度比,当搜索的角度和分别等于两个相干信号源的方位角θ1和θ2,且搜索到的两个相干信号的幅度比等于真实的两个相干信号的幅度比k时,则实现对接收数据中两个相干信号源的正确分解。
步骤S2:确定子信号的MUSIC谱为目标函数。具体的说就是我们定义目标函数为分解后得到的子信号的MUSIC谱:
式中: 或 是第i(i≤2)个水声信号源在单只矢量水听器上的方向矢量,表达式为:
a(θi)=[1cosθisinθi]T
是对原始接收数据协方差矩阵进行特征分解得到噪声子空间。
当搜索的角度和分别等于两个相干信号源的方位角θ1和θ2,且搜索到的两个相干信号的幅度比等于真实的两个相干信号的幅度比k时,则实现对接收数据中两个相干信号源的正确分解,目标函数将会取最大值。
步骤S3:根据目标函数进行最优搜索,确定两个相干信号源的方位角。具体地说就是通过最优搜索,找到目标函数最大值对应的搜索角度和即为两个相干信号源的入射角度θ1和θ2。具体步骤如下:
(1)确定搜索范围并设置初始参数根据初始参数对接收数据进行分解得到子信号X1(t),进而求出初始参数条件下的目标函数值,定义该目标函数值为最大值,并将初始参数值定义为估计到的参数值,即
(2)重置参数求出在新的参数条件下的目标函数
(3)比较fmax和fnew的大小,将其中大的一个定义为新的目标函数的最大值,并将其对应的参数定义为新的估计到的参数值,即:
fmax=max(fmax,fnew)
<θ1,θ2,k>=<θ1,θ2,k>|max(fmax,fnew)
(4)重复步骤(2)、(3),直到将搜索范围内的所有参数全部搜索完毕,此时得到的θ1、θ2即为估计到的两个相干信号源的方位。
本发明两个相干信号的加权矩阵包含了两个相干信号的方位角和两个相干信号的幅度比,当搜索的角度和分别等于两个相干信号源的方位角θ1和θ2,且搜索到的两个相干信号的幅度比等于真实的两个相干信号的幅度比k时,则实现对接收数据中两个相干信号源的正确分解。定义子信号的MUSIC谱为目标函数。当MUSIC谱取最大值时,则实现对接收数据中两个相干信号源的正确分解,此时对应的两个搜索角度即为两个相干信号源的方位角。
Claims (2)
1.一种基于单矢量水听器的相干信号源方位估计方法,其特征在于:包括以下步骤,
步骤一:确定声压P通道、振速Vx通道和振速Vy通道中两个相干信号的加权矩阵T,根据加权矩阵T将接收数据X(t)分成两个子信号X1(t)和X2(t),其中X1(t)=TX(t);
步骤二:确定子信号的MUSIC谱为目标函数,子信号的MUSIC谱为:
其中: 或 是第i个水声信号源在单只矢量水听器上的方向矢量,i≤2,和为搜索角度,是搜索到的两个相干信号的幅度比,是对原始接收数据协方差矩阵进行特征分解得到噪声子空间;
步骤三:根据目标函数进行最优搜索,找到目标函数最大值对应的搜索角度和即为两个相干信号源的入射角度θ1和θ2。
2.根据权利要求1所述的一种基于单矢量水听器的相干信号源方位估计方法,其特征在于:所述的目标函数进行最优搜索的过程为:
(1)确定搜索范围并设置初始参数根据初始参数对接收数据进行分解得到子信号X1(t),进而求出初始参数条件下的目标函数值,定义该目标函数值为最大值,并将初始参数值定义为估计到的参数值,即
(2)重置参数求出在新的参数条件下的目标函数
(3)比较fmax和fmax的大小,将其中大的一个定义为新的目标函数的最大值,并将其对应的参数定义为新的估计到的参数值,即:
fmax=max(fmax,fnew)
<θ1,θ2,k>=<θ1,θ2,k>|max(fmax,fnew)
(4)重复步骤(2)~(3),直到将搜索范围内的所有参数全部搜索完毕。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510477827.6A CN105138743B (zh) | 2015-08-06 | 2015-08-06 | 一种基于单矢量水听器的相干信号源方位估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510477827.6A CN105138743B (zh) | 2015-08-06 | 2015-08-06 | 一种基于单矢量水听器的相干信号源方位估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105138743A true CN105138743A (zh) | 2015-12-09 |
CN105138743B CN105138743B (zh) | 2018-12-07 |
Family
ID=54724090
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510477827.6A Active CN105138743B (zh) | 2015-08-06 | 2015-08-06 | 一种基于单矢量水听器的相干信号源方位估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105138743B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106202892A (zh) * | 2016-06-30 | 2016-12-07 | 哈尔滨工业大学(威海) | 一种基于噪声子空间单一矢量的快速doa估计算法 |
CN106646350A (zh) * | 2016-09-08 | 2017-05-10 | 哈尔滨工程大学 | 一种单只矢量水听器各通道幅度增益不一致时的修正方法 |
CN111596262A (zh) * | 2020-05-07 | 2020-08-28 | 武汉大学 | 矢量水听器及基于该矢量水听器的多目标方位估计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0616705A4 (en) * | 1991-12-12 | 1994-10-26 | Essex Corp | RADIATION WITH DELAYED TIME. |
CN102255671A (zh) * | 2011-07-06 | 2011-11-23 | 哈尔滨工程大学 | 一种单矢量传感器水声多址通信方法 |
CN102981146A (zh) * | 2012-11-19 | 2013-03-20 | 哈尔滨工程大学 | 一种单矢量水听器被动定位方法 |
-
2015
- 2015-08-06 CN CN201510477827.6A patent/CN105138743B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0616705A4 (en) * | 1991-12-12 | 1994-10-26 | Essex Corp | RADIATION WITH DELAYED TIME. |
CN102255671A (zh) * | 2011-07-06 | 2011-11-23 | 哈尔滨工程大学 | 一种单矢量传感器水声多址通信方法 |
CN102981146A (zh) * | 2012-11-19 | 2013-03-20 | 哈尔滨工程大学 | 一种单矢量水听器被动定位方法 |
Non-Patent Citations (3)
Title |
---|
PALANISAMY P 等: "Two-dimensional DOA estimation of coherent signals using acoustic vector sensor array", 《SIGNAL PROCESSING》 * |
孟春霞 等: "单矢量水听器多目标方位估计的算法研究", 《舰船科学技术》 * |
王绪虎 等: "矢量水听器阵波束域 MUSIC 算法研究", 《计算机工程与应用》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106202892A (zh) * | 2016-06-30 | 2016-12-07 | 哈尔滨工业大学(威海) | 一种基于噪声子空间单一矢量的快速doa估计算法 |
CN106646350A (zh) * | 2016-09-08 | 2017-05-10 | 哈尔滨工程大学 | 一种单只矢量水听器各通道幅度增益不一致时的修正方法 |
CN106646350B (zh) * | 2016-09-08 | 2019-06-14 | 哈尔滨工程大学 | 一种单只矢量水听器各通道幅度增益不一致时的修正方法 |
CN111596262A (zh) * | 2020-05-07 | 2020-08-28 | 武汉大学 | 矢量水听器及基于该矢量水听器的多目标方位估计方法 |
CN111596262B (zh) * | 2020-05-07 | 2023-03-10 | 武汉敏声新技术有限公司 | 矢量水听器及基于该矢量水听器的多目标方位估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105138743B (zh) | 2018-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wen | Computationally efficient DOA estimation algorithm for MIMO radar with imperfect waveforms | |
CN107329108B (zh) | 基于内插虚拟阵列协方差矩阵Toeplitz化重建的互质阵列波达方向估计方法 | |
CN107102291B (zh) | 基于虚拟阵列内插的无网格化互质阵列波达方向估计方法 | |
CN109815942B (zh) | 基于海洋环境噪声信号的简正波特征提取方法 | |
CN105589056B (zh) | 一种多目标远近场混合源定位方法 | |
Saeed et al. | Underwater optical sensor networks localization with limited connectivity | |
CN106972882B (zh) | 基于虚拟域空间功率谱估计的互质阵列自适应波束成形方法 | |
CN103901395B (zh) | 一种冲击噪声环境下相干信号波达方向动态跟踪方法 | |
CN104408278A (zh) | 一种基于干扰噪声协方差矩阵估计的稳健波束形成方法 | |
CN104898113A (zh) | 一种多输入多输出雷达波形设计方法 | |
CN104698433A (zh) | 基于单快拍数据的相干信号doa估计方法 | |
CN105376673A (zh) | 基于空间分析的麦克风阵列处理器 | |
CN111123192A (zh) | 一种基于圆形阵列和虚拟扩展的二维doa定位方法 | |
WO2008020870A3 (en) | M-algorithm multiuser detector with correlation based pruning | |
CN105138743A (zh) | 一种基于单矢量水听器的相干信号源方位估计方法 | |
CN106093866A (zh) | 一种适用于空心球阵列的声源定位方法 | |
CN105188133A (zh) | 一种基于准平稳信号局部协方差匹配的kr子空间doa估计方法 | |
CN103323810A (zh) | 一种l阵方位角和俯仰角配对的信号处理方法 | |
CN101666876B (zh) | 一种双线阵声纳左右舷识别方法及相应系统 | |
CN103513249B (zh) | 一种宽带相干模基信号处理方法及系统 | |
Huang et al. | Toeplitz matrix completion for direction finding using a modified nested linear array | |
CN116366170A (zh) | 一种相关噪声下矢量水听器声压振速信号的合并方法 | |
CN105372644A (zh) | 一种基于动态重修正的自适应波束形成方法及系统 | |
CN112152959B (zh) | 水声通信探测信号处理方法、装置、设备及存储介质 | |
Qin et al. | Cross-correlation matrix Root-MUSIC algorithm for bistatic multiple-input multiple-output radar |
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 |