CN106646350B - 一种单只矢量水听器各通道幅度增益不一致时的修正方法 - Google Patents
一种单只矢量水听器各通道幅度增益不一致时的修正方法 Download PDFInfo
- Publication number
- CN106646350B CN106646350B CN201610810833.3A CN201610810833A CN106646350B CN 106646350 B CN106646350 B CN 106646350B CN 201610810833 A CN201610810833 A CN 201610810833A CN 106646350 B CN106646350 B CN 106646350B
- Authority
- CN
- China
- Prior art keywords
- modifying factor
- channel
- single vector
- vector hydrophone
- music
- 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.)
- Active
Links
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
- G01S3/803—Systems for determining direction or deviation from predetermined direction using amplitude comparison of signals derived from receiving transducers or transducer systems having differently-oriented directivity characteristics
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
本发明具体涉及一种单只矢量水听器各通道幅度增益不一致时的修正方法。本发明包括:用常规的MVDR算法及MUSIC算法对单只矢量水听器接收到的数据进行方位估计;确定Vx通道及Vy通道修正因子的搜索范围;设定修正因子,并用其对接收数据进行修正,以修正后信号归一化MUSIC谱的取值范围为目标函数进行最优搜索,确定最优修正因子;用最优修正因子对接收数据进行修正并进行方位估计。本发明在单只矢量水听器各个通道幅度增益不一致的条件下,可通过修正因子对接收数据进行修正,解决了单只矢量水听器各通道幅度增益不一致条件下,信号方位估计精度低、存在“伪峰”等问题,实现了单只矢量水听器各通道幅度增益不一致时的方位估计。
Description
技术领域
本发明具体涉及一种单只矢量水听器各通道幅度增益不一致时的修正方法。
背景技术
与声压水听器相比,矢量水听器具有三个通道,分别为声压p通道、振速vx通道、振速vy通道,因此矢量水听器可以同时、共点测量声场的声压与振速,获得更多的声场信息,依靠单个矢量水听器就可以完成水下目标的方位估计,且不存在多值模糊的问题。但是,单只矢量水听器方位估计的性能严重依赖于良好的通道特性和外部环境。当单只矢量水听器的声压和振速通道的幅度和相位特性基本一致时,单只矢量水听器可以正确反映声压和振速的信息,进而进行高分辨方位估计;但是,在实际的工程应用中,由于传感器敏感元件的差异、硬件放大电路构造的差异、安装的差异以及传输信道的不同等原因,单只矢量水听器的声压和振速通道的幅度和相位特性很难做到严格一致。因此研究单只矢量水听器通道的幅相误差给方位估计带来的影响十分必要。
为了解决单只矢量水听器通道幅相特性不一致导致空间谱的谱峰变宽,甚至出现“伪峰”的问题,本发明提供了一种单只矢量水听器各通道幅度增益不一致时的修正方法。本发明方法通过最优搜索,实现对单只矢量水听器各通道幅度增益一致性的修正,解决了单只矢量水听器各通道幅度增益不一致条件下,信号方位估计精度低、存在“伪峰”等问题,实现了单只矢量水听器各通道幅度增益不一致时的方位估计。
发明内容
本发明的目的是提供一种单只矢量水听器各通道幅度增益不一致时的修正方法。
本发明的目的是这样实现的:
步骤S1:用常规的MVDR方法及MUSIC方法对单矢量水听器接收到的数据进行方位估计;
步骤S2:确定Vx通道及Vy通道修正因子的搜索范围;
步骤S3:设定修正因子,并用其对接收数据进行修正,以修正后信号归一化MUSIC谱的取值范围为目标函数进行最优搜索,确定最优修正因子;
步骤S4:用最优修正因子对接收数据进行修正并进行方位估计;
本发明的核心技术在于根据常规MVDR及MUSIC方法方位估计得到的结果,对比分析可判断各通道幅度增益相对大小,从而确定Vx通道及Vy通道修正因子的搜索范围;并以修正因子修正后的数据进行方位估计得到的归一化MUSIC谱的取值范围为目标函数进行最优搜索,当目标函数取最大值时,对应的修正因子即为最优修正因子,由最优修正因子修正后的数据,可实现单只矢量水听器各通道幅度增益不一致时的方位估计。
本发明的有益效果在于:与目前存在的常规的单只矢量水听器目标方位估计方法相比,本发明在单只矢量水听器各个通道幅度增益不一致的条件下,可通过修正因子对接收数据进行修正,解决了单只矢量水听器各通道幅度增益不一致条件下,信号方位估计精度低、存在“伪峰”等问题,实现了单只矢量水听器各通道幅度增益不一致时的方位估计。
附图说明
图1为本发明的方法流程图。
图2为通道幅度增益不一致时常规MVDR及MUSIC方法估计结果。
图3为常规方法与本文方法对单目标方位估计结果
具体实施方式
下面结合附图对本发明做进一步描述:
本发明提供的是一种单只矢量水听器各通道幅度增益不一致时的修正算法。首先,用MVDR算法(即最小方差无失真响应方位估计器)及MUSIC算法(即多重信号分类算法)对单矢量水听器接收到的数据进行方位估计;其次,根据常规MVDR算法及MUSIC算法方位估计得到的结果判断各通道幅度增益相对大小,确定Vx通道及Vy通道修正因子的搜索范围;再次,设定修正因子,并用其对接收数据进行修正,以修正后信号归一化MUSIC谱的取值范围为目标函数进行最优搜索,目标函数取最大值时对应的修正因子即为最优修正因子;最后,用最优修正因子对接收数据进行修正并进行方位估计。本发明方法通过最优搜索,实现对单只矢量水听器各通道幅度增益一致性的修正,解决了单只矢量水听器各通道幅度增益不一致条件下,信号方位估计精度低、存在“伪峰”等问题,实现了单只矢量水听器各通道幅度增益不一致时的方位估计。
步骤S1:用常规的MVDR方法及MUSIC方法对单矢量水听器接收到的数据进行方位估计。当单只矢量水听器各通道幅度增益不一致时,通道幅度误差可以通过在单只矢量水听器接收数据模型中引入与方位无关的通道幅度增益对角阵来描述,其数据模型如下:
X(t)=Γ(A(θ)x(t)+N(t))
式中:Γ是通道的幅度增益对角阵。θ为入射声波的水平方位角;A(θ)是单只矢量水听器阵列流型;x(t)为传感器接收的声压波形;N(t)为噪声数据矢量。假定Γ满足:
Γ=diag(ηp,ηx,ηy)
式中:ηx,ηy,ηp分别是3个通道的幅度增益系数。由于此处幅度增益为实数,所以此时的接收数据模型为:
求其协方差矩阵为:
R1=E{XXH}=ΓRΓH
式中:H表示共轭转置,E{·}表示集合平均,R为理想接收数据模型的协方差矩阵,注意到R1与R是酉相似矩阵,R1的逆矩阵为:
R1'=(ΓRΓH)'=(ΓH)'R'Γ'
对R1进行特征分解有:
式中:VS=ΓUS;VN=ΓUN,US是由大的特征值对应的特征矢量张成的子空间,即信号子空间,而UN是由小的特征值对应的特征矢量张成的子空间,即噪声子空间。则此时得到单只矢量水听器各通道幅度增益不一致时的MVDR谱与MUSIC谱分别为:
对上面的两个式子进行峰值搜索,其峰值对应的方位即为估计到的目标方位。由上面两式可知,当单只矢量水听器各通道幅度增益不一致时,将会对单只矢量水听器目标方位估计结果产生严重影响。图2给出了ηx,ηy,ηp不同取值情况下,对MVDR、MUSIC方法方位估计结果的影响。其中,信号频率为f=3000Hz,采样频率为fs=30000Hz,入射角度为θ=100°,所加噪声为零均值的高斯白噪声,信噪比为SNR=20dB。
步骤S2:根据步骤S1得到的常规MVDR及MUSIC目标方位估计结果,以理论分析为基础,通过对比分析确定Vx通道及Vy通道修正因子的搜索范围
步骤S3:确定修正后信号归一化MUSIC谱的取值范围为目标函数,并进行最优搜索。具体的说就是我们定义目标函数为修正后信号归一化MUSIC谱的取值范围:
式中:为由修正因子修正后的数据进行MUSIC谱估计得到的归一化MUSIC谱。
然后,根据目标函数进行最优搜索,从而确定最优修正因子的取值。具体地说就是通过最优搜索,找到目标函数最大值对应的修正因子和具体步骤如下:
(1)设置初始修正因子根据初始修正因子修正后的数据进行MUSIC谱估计,进而求出初始修正因子条件下的目标函数值,定义该目标函数值为最大值,并将初始修正因子值定义为最优修正因子,即
(2)重置修正因子求出在新的修正因子条件下的目标函数
(3)比较fmax和fnew的大小,将其中大的一个定义为新的目标函数的最大值,并将其对应的修正因子定义为新的最优修正因子,即:
fmax=max(fmax、fnew)
<kx、ky>=<kx、ky>|max(fmax、fnew)
(4)重复步骤(2)、(3),直到将搜索范围内的所有参数全部搜索完毕,此时得到的kx、ky即为最优修正因子。
步骤S4:用步骤S3中所得最优修正因子kx、ky对单只矢量水听器接收数据进行修正,并对修正后的数据进行方位估计,最终得到单只矢量水听器各通道幅度增益不一致时的目标方位估计。
本发明性能设计实例:
单目标方位估计的仿真实例
单只矢量水听器方位估计,入射信号频率为f=3000Hz,采样频率为fs=30000Hz,入射角度为θ=100°,所加噪声为零均值的高斯白噪声,信噪比为SNR=20dB。当单只矢量水听器各通道幅度增益分别为:ηp=0.1,ηx=ηy=1,的特殊情况时,常规的MVDR、MUSIC谱估计方法与本发明所提出的修正的MVDR、MUSIC方法得到的方位估计结果如图3所示。
如图3所示的仿真结果可以看出,当单只矢量水听器各通道幅度增益不一致时,常规的MVDR方法和MUSIC方法目标分辨性能恶化,且有时产生伪峰。修正后的MVDR方法和MUSIC方法仍具有较高的目标分辨性能,且不存在伪峰。
Claims (3)
1.一种单只矢量水听器各通道幅度增益不一致时的修正方法,其特征在于,包括如下步骤:
步骤S1:用常规的MVDR算法及MUSIC算法对单只矢量水听器接收到的数据进行方位估计;
步骤S2:确定Vx通道及Vy通道修正因子的搜索范围;
步骤S3:设定修正因子,并用其对接收数据进行修正,以修正后信号归一化MUSIC谱的取值范围为目标函数进行最优搜索,确定最优修正因子;
确定修正后信号归一化MUSIC谱的取值范围为目标函数,并进行最优搜索;具体的说就是我们定义目标函数为修正后信号归一化MUSIC谱的取值范围:
式中:为由修正因子修正后的数据进行MUSIC谱估计得到的归一化MUSIC谱;
然后,根据目标函数进行最优搜索,从而确定最优修正因子的取值;具体地说就是通过最优搜索,找到目标函数最大值对应的修正因子和具体步骤如下:
(1)设置初始修正因子根据初始修正因子修正后的数据进行MUSIC谱估计,进而求出初始修正因子条件下的目标函数值,定义该目标函数值为最大值,并将初始修正因子值定义为最优修正因子,即
(2)重置修正因子求出在新的修正因子条件下的目标函数
(3)比较fmax和fnew的大小,将其中大的一个定义为新的目标函数的最大值,并将其对应的修正因子定义为新的最优修正因子,即:
fmax=max(fmax、fnew)
<kx、ky>=<kx、ky>|max(fmax、fnew)
(4)重复步骤(2)、(3),直到将搜索范围内的所有参数全部搜索完毕,此时得到的kx、ky即为最优修正因子;
步骤S4:用最优修正因子对接收数据进行修正并进行方位估计。
2.根据权利要求1所述的一种单只矢量水听器各通道幅度增益不一致时的修正方法,其特征在于:所述的根据常规MVDR及MUSIC算法方位估计,根据得到的结果,对比分析可判断各通道幅度增益相对大小,从而确定Vx通道及Vy通道修正因子的搜索范围。
3.根据权利要求1所述的一种单只矢量水听器各通道幅度增益不一致时的修正方法,其特征在于:定义由所述的修正因子修正后的数据进行方位估计得到的归一化MUSIC谱的取值范围为目标函数,当目标函数取最大值时,对应的修正因子即为最优修正因子,由最优修正因子修正后的数据,实现单只矢量水听器各通道幅度增益不一致时的方位估计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610810833.3A CN106646350B (zh) | 2016-09-08 | 2016-09-08 | 一种单只矢量水听器各通道幅度增益不一致时的修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610810833.3A CN106646350B (zh) | 2016-09-08 | 2016-09-08 | 一种单只矢量水听器各通道幅度增益不一致时的修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106646350A CN106646350A (zh) | 2017-05-10 |
CN106646350B true CN106646350B (zh) | 2019-06-14 |
Family
ID=58851634
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610810833.3A Active CN106646350B (zh) | 2016-09-08 | 2016-09-08 | 一种单只矢量水听器各通道幅度增益不一致时的修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106646350B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111273265A (zh) * | 2020-03-02 | 2020-06-12 | 中国人民解放军国防科技大学 | 基于demon谱的单矢量水听器自主测距装置及方法 |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108469599B (zh) * | 2018-02-28 | 2021-11-23 | 哈尔滨工程大学 | 一种声矢量传感器幅度加权music测向方法 |
CN108845287B (zh) * | 2018-04-09 | 2022-08-02 | 哈尔滨工程大学 | 小生境遗传算法的单矢量水听器相干源快速方位估计方法 |
CN109239652B (zh) * | 2018-09-12 | 2023-07-14 | 中国船舶重工集团公司第七0五研究所 | 一种水下声源水平方位检测方法 |
CN110361690B (zh) * | 2019-06-14 | 2022-12-13 | 哈尔滨工程大学 | 一种平面障板条件下单矢量水听器测向方法 |
CN113093098B (zh) * | 2021-04-09 | 2023-05-16 | 河南理工大学 | 基于lp范数补偿的轴向不一致矢量水听器阵列测向方法 |
CN114047499B (zh) * | 2022-01-12 | 2022-05-17 | 宁波博海深衡科技有限公司武汉分公司 | 传感器阵列信号幅度不一致时的信号归一化方法及设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102353937A (zh) * | 2011-07-06 | 2012-02-15 | 哈尔滨工程大学 | 单矢量有源平均声强器 |
CN103267966A (zh) * | 2013-05-20 | 2013-08-28 | 西北工业大学 | 一种单个二维压差式矢量水听器方位估计的方法 |
WO2014088724A1 (en) * | 2012-12-03 | 2014-06-12 | Qualcomm Incorporated | Direction of arrival estimation using linear array |
CN105138743A (zh) * | 2015-08-06 | 2015-12-09 | 哈尔滨工程大学 | 一种基于单矢量水听器的相干信号源方位估计方法 |
-
2016
- 2016-09-08 CN CN201610810833.3A patent/CN106646350B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102353937A (zh) * | 2011-07-06 | 2012-02-15 | 哈尔滨工程大学 | 单矢量有源平均声强器 |
WO2014088724A1 (en) * | 2012-12-03 | 2014-06-12 | Qualcomm Incorporated | Direction of arrival estimation using linear array |
CN103267966A (zh) * | 2013-05-20 | 2013-08-28 | 西北工业大学 | 一种单个二维压差式矢量水听器方位估计的方法 |
CN105138743A (zh) * | 2015-08-06 | 2015-12-09 | 哈尔滨工程大学 | 一种基于单矢量水听器的相干信号源方位估计方法 |
Non-Patent Citations (2)
Title |
---|
Modified MVDR algorithm for DOA estimation using acoustic vector hydrophone;Guo-long Liang et al.;《2011 IEEE International Conference on Computer Science and Automation Engineering》;20110612;全文 |
单矢量水听器多目标方位估计的算法研究;孟春霞 等;《船舰科学技术》;20070430;第29卷(第2期);第123-126页 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111273265A (zh) * | 2020-03-02 | 2020-06-12 | 中国人民解放军国防科技大学 | 基于demon谱的单矢量水听器自主测距装置及方法 |
CN111273265B (zh) * | 2020-03-02 | 2020-09-29 | 中国人民解放军国防科技大学 | 基于demon谱的单矢量水听器自主测距装置及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106646350A (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106646350B (zh) | 一种单只矢量水听器各通道幅度增益不一致时的修正方法 | |
US10602267B2 (en) | Sound signal processing apparatus and method for enhancing a sound signal | |
CN108375763B (zh) | 一种应用于多声源环境的分频定位方法 | |
US7583808B2 (en) | Locating and tracking acoustic sources with microphone arrays | |
CN111474521B (zh) | 多径环境中基于麦克风阵列的声源定位方法 | |
CN108710103B (zh) | 基于稀疏阵列的强弱多目标超分辨测向与信源数估计方法 | |
See et al. | Parametric sensor array calibration using measured steering vectors of uncertain locations | |
CN110544490A (zh) | 一种基于高斯混合模型和空间功率谱特征的声源定位方法 | |
de Paolo et al. | Skill assessment of resolving ocean surface current structure using compact-antenna-style HF radar and the MUSIC direction-finding algorithm | |
Huang et al. | Two-step spherical harmonics ESPRIT-type algorithms and performance analysis | |
CN110888105A (zh) | 一种基于卷积神经网络和接收信号强度的doa估计方法 | |
CN111580042B (zh) | 一种基于相位优化的深度学习测向方法 | |
CN110196407B (zh) | 一种基于频率预估的单矢量水听器信号来波方向估计方法 | |
Hu et al. | Decoupled direction-of-arrival estimations using relative harmonic coefficients | |
Nesta et al. | Enhanced multidimensional spatial functions for unambiguous localization of multiple sparse acoustic sources | |
Okane et al. | Resolution improvement of wideband direction-of-arrival estimation" Squared-TOPS" | |
CN108957389A (zh) | 一种实数域多通道信号目标方位估计方法 | |
Lee et al. | Deep learning-based ultra short baseline underwater positioning | |
Wang et al. | Adaptive frequency response calibration method for microphone arrays | |
Suo et al. | Single acoustic vector sensor DOA enhanced by unsaturated bistable stochastic resonance with linear amplitude response constrained | |
KR20200061550A (ko) | 렌즈 안테나 기반 2단계 도래각 추정 방법 | |
CN113030847B (zh) | 一种用于双通道测向系统的深度学习数据集生成方法 | |
CN111722178B (zh) | 一种基于指向性模型数值求解的远场窄带信号来波方向估计方法 | |
Bo et al. | A novel wideband DOA estimation method using direction-free focusing matrix | |
JP2005077205A (ja) | 音源方向推定装置、信号の時間遅延推定装置及びコンピュータプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |