CN109814065A - 基于相位因子加权的波束形成方法 - Google Patents
基于相位因子加权的波束形成方法 Download PDFInfo
- Publication number
- CN109814065A CN109814065A CN201910030601.XA CN201910030601A CN109814065A CN 109814065 A CN109814065 A CN 109814065A CN 201910030601 A CN201910030601 A CN 201910030601A CN 109814065 A CN109814065 A CN 109814065A
- Authority
- CN
- China
- Prior art keywords
- phase factor
- phase
- signal
- wave beam
- beam output
- 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 29
- 238000003032 molecular docking Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 9
- 238000004364 calculation method Methods 0.000 abstract description 5
- 230000015572 biosynthetic process Effects 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000004304 visual acuity Effects 0.000 description 1
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于相位因子加权的波束形成方法,所述方法将阵列接收信号变换到波束域,然后计算波束序列信号频点处相位与幅度,求取相位因子及其方差,利用方差倒数加权计算波束输出功率,实现水下目标探测和DOA估计。本发明适用于窄带和宽带信号以及任意阵型,具有较好的稳健性,与常规波束形成方法相比,在不明显增加计算量情况下,能显著降低旁瓣,提高水下弱目标的探测性能。
Description
技术领域
本发明属于水声信号处理,涉及线阵或圆阵等任意形状的阵列,具体涉及一种基于相位因子加权的波束形成方法。
背景技术
常规波束形成方法在实际应用中,由于其不需要目标先验信息、稳健性好、计算量小、实时性好等优点,从而在声纳、通信、雷达等领域得到广泛运用,但常规波束形成方法的旁瓣高、分辨力受瑞利限限制,导致对弱目标的探测性能下降、空间分辨率低。改进的常规波束形成算法,利用波束序列相位因子的稳定性,在只增加少量计算量的同时,显著降低常规波束形成的旁瓣,减小主瓣宽度,提高弱目标的检测性能。该方法不受阵型限制,适用于窄带和宽带信号。
发明内容
发明目的:针对现有技术无法实现对弱目标信号的有效探测的问题,本发明提供基于相位因子加权的波束形成方法,降低常规波束形成的旁瓣,提高常规波束形成算法对弱目标信号的探测性能。
技术方案:基于相位因子加权的波束形成方法,所述方法基于波束序列相位因子的稳定性,将阵列接收信号变换到波束域,对波束序列计算相位因子,利用相位因子的方差得到加权后的波束输出功率;包括如下步骤:
(1)对接收信号分段作FFT,并计算波束输出;
(2)对每段数据的波束输出分别求取幅度值和相位值;
(3)对波束输出求取相位因子;
(4)计算各段相位因子的方差;
(5)利用步骤(4)中的方差倒数对波束输出功率进行加权,得到新的功率输出。
进一步的,步骤(1)所述接收信号表达式如下:
所述波束输出为:
其中ym(t)为第m个阵元的接收信号,Yi,k(θ)为扫描角θ处的波束输出,K为信源个数,Ak、fk、分别为第k个信号的幅度、频率和初相,τm为源信号到第m个阵元的延时,M为阵元数,θ为扫描角度,θ0为信号方位,为频域离散化后的信号频率,fs为采样频率,d为阵元间距,c为声速,N为每段数据的长度。
进一步的,步骤(2)中对于第i段频点为k的信号的波束输出的幅度值和相位值分别如下:
其中,k为fk对应的频点,f′k为频点k对应的频率,fs为采样频率,N为每段数据的长度,d为阵元间距,c为声速,θ为扫描角度,θ0为信号方位角。
进一步的,步骤(3)中所求的相位因子为:
其中分别为第k个信号,第i段、第i+1段和第i+2段的相位。
进一步的,步骤(4)所述方差值计算表达式如下:
δk=Var(ψi,k(θ))
其中Var(.)表示计算方差。
进一步的,步骤(5)所得到的新的波束功率输出为:
其中I为信号段数。
本发明将阵列信号变换到波束域后,计算波束序列信号频点处的幅度和相位,得到相位因子,计算相位因子方差加权后的波束输出功率,之后进行水下弱信号探测和DOA估计。
有益效果:本发明与现有技术相比,其显著的效果在于:本发明不仅保持了现有的波束形成方法不需要目标先验信息、稳健性好、计算量低、实时性好的特点;并且,该方法适用于线阵或圆阵等任意形状的阵列,能够显著降低常规波束形成的旁瓣,减小主瓣宽度,提高水下弱目标的探测性能。
附图说明
图1是本发明的流程图;
图2是本发明与常规波束形成的波束输出功率图。
具体实施方式
为了详细的说明本发明所公开的技术方案,下面结合说明书附图及具体实施方式对本发明做进一步的阐述。
本发明公开的是一种基于相位因子加权的波束形成方法,该方法步骤流程图如图1所示。
所述方法基于波束序列相位因子的稳定性,对波束序列求取相位因子,并利用相位因子的方差重新计算波束输出功率,包括如下步骤:
(1)对接收信号分段作FFT(傅里叶变换),并计算波束输出;
(2)对每段数据的波束输出分别求取幅度值和相位值;
(3)对波束输出求取相位因子;
(4)计算各段相位因子的方差;
(5)利用步骤(4)中的方差倒数对波束输出功率进行加权,求得波束输出功率。
下面对上述步骤做具体的说明,以便本领域技术人员在现有技术的基础上理解该技术方案,且能够清楚的实施,并得到本发明所述的技术效果。
线列阵的水听器数目为M=32,水听器间距d=7.1429m,采样频率fs=1000Hz,声速c=1500m/s,信号频率为fo=70Hz,快拍数据共I=10段,每段N=1024点数据,信噪比为SNR=-20dB。
将接收数据分为I段,每段N点,得到接收阵列数据为:
其中,nm为均值为0,方差为σ2的高斯白噪声,各个阵元的噪声相互独立;τm为第m个阵元的延时。
对接收数据分段进行FFT,得到第i段数据的频域接收阵列数据为:
对每段数据作频域波束形成得:
其中k为信号频率对应的频点,Xm(k)为第m个阵元、第k个频点的谱值。
对每段数据的波束输出计算幅度和相位:
Ai(θ)=sqrt{Im2[Yi(θ)]+Re2[Yi(θ)]}
利用各段波束输出的相位值计算相位因子:
计算各段相位因子余弦值的方差:
δk=Var(ψi,k(θ))
计算相位因子方差加权后的波束输出功率:
为了进一步的说明本发明的技术效果,图2给出了-20dB信噪比条件下基于相位因子加权的波束形成方法和常规波束形成方法的波束输出功率图,两种方法都能准确估计目标的方位,但基于相位因子加权的波束形成方法将波束图的旁瓣降低了约25dB,更有利于水下弱目标的探测,同时减小了主瓣宽度,提高了空间分辨率。
Claims (6)
1.基于相位因子加权的波束形成方法,其特征在于:所述方法基于波束序列相位因子的稳定性,将阵列接收信号变换到波束域,对波束序列计算相位因子,利用相位因子的方差得到加权后的波束输出功率;包括如下步骤:
(1)对接收信号分段作FFT,并计算波束输出;
(2)对每段数据的波束输出分别求取幅度值和相位值;
(3)对波束输出求取相位因子;
(4)计算各段相位因子的方差;
(5)利用相位因子的方差进行加权计算得到波束输出功率。
2.根据权利要求1所述的基于相位因子加权的波束形成方法,其特征在于:步骤(1)所述接收信号表达式如下:
所述波束输出为:
其中ym(t)为第m个阵元的接收信号,Yi,k(θ)为扫描角θ处的波束输出,K为信源个数,Ak、fk、分别为第k个信号的幅度、频率和初相,τm为源信号到第m个阵元的延时,M为阵元数,θ为扫描角度,θ0为信号方位,为频域离散化后的信号频率,fs为采样频率,d为阵元间距,c为声速,N为每段数据的长度。
3.根据权利要求1所述的基于相位因子加权的波束形成方法,其特征在于:步骤(2)中对于第i段频点为k的信号的波束输出的幅度值和相位值分别如下:
式中:k为fk对应的频点,f′k为k对应的频率,fs为采样频率,N为每段数据的长度,d为阵元间距,c为声速,θ为扫描角度,θ0为信号方位角。
4.根据权利要求1所述的基于相位因子加权的波束形成方法,其特征在于:步骤(3)中所求的相位因子为:
其中分别为第k个信号,第i段、第i+1段和第i+2段的相位。
5.根据权利要求1所述的基于相位因子加权的波束形成方法,其特征在于:步骤(4)所述相位因子的方差如下:
δk=Var(ψi,k(θ))。
6.根据权利要求1所述的基于相位因子加权的波束形成方法,其特征在于:步骤(5)相位因子方差加权后的波束输出功率为:
其中I为信号段数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910030601.XA CN109814065B (zh) | 2019-01-14 | 2019-01-14 | 基于相位因子加权的波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910030601.XA CN109814065B (zh) | 2019-01-14 | 2019-01-14 | 基于相位因子加权的波束形成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109814065A true CN109814065A (zh) | 2019-05-28 |
CN109814065B CN109814065B (zh) | 2020-07-31 |
Family
ID=66604243
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910030601.XA Expired - Fee Related CN109814065B (zh) | 2019-01-14 | 2019-01-14 | 基于相位因子加权的波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109814065B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110391820A (zh) * | 2019-06-11 | 2019-10-29 | 东南大学 | 一种基于dft的规避同频干扰的新型通信接收方法 |
CN110609271A (zh) * | 2019-10-29 | 2019-12-24 | 海鹰企业集团有限责任公司 | 一种基于空间变迹的波束旁瓣抑制方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120326925A1 (en) * | 2007-05-21 | 2012-12-27 | Spatial Digital Systems | Digital beam-forming apparatus and technique for a multi-beam global positioning system (gps) receiver |
CN105137437A (zh) * | 2015-07-20 | 2015-12-09 | 中国科学院声学研究所 | 一种基于空域相位方差加权的目标检测方法 |
-
2019
- 2019-01-14 CN CN201910030601.XA patent/CN109814065B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120326925A1 (en) * | 2007-05-21 | 2012-12-27 | Spatial Digital Systems | Digital beam-forming apparatus and technique for a multi-beam global positioning system (gps) receiver |
CN105137437A (zh) * | 2015-07-20 | 2015-12-09 | 中国科学院声学研究所 | 一种基于空域相位方差加权的目标检测方法 |
Non-Patent Citations (4)
Title |
---|
王思秀 等: "基于频域相位方差加权的线性调频信号检测方法", 《计算机应用》 * |
程娜: "超声成像波束合成理论与算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
郑恩明 等: "一种基于瞬时相位方差加权的线谱检测器", 《电子与信息学报》 * |
郑恩明 等: "基于波束域相位稳定性的目标检测方法研究", 《电子学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110391820A (zh) * | 2019-06-11 | 2019-10-29 | 东南大学 | 一种基于dft的规避同频干扰的新型通信接收方法 |
CN110391820B (zh) * | 2019-06-11 | 2021-06-11 | 东南大学 | 一种基于dft的规避同频干扰的通信接收方法 |
CN110609271A (zh) * | 2019-10-29 | 2019-12-24 | 海鹰企业集团有限责任公司 | 一种基于空间变迹的波束旁瓣抑制方法 |
CN110609271B (zh) * | 2019-10-29 | 2022-12-13 | 海鹰企业集团有限责任公司 | 一种基于空间变迹的波束旁瓣抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109814065B (zh) | 2020-07-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108375763B (zh) | 一种应用于多声源环境的分频定位方法 | |
CN108828522B (zh) | 一种利用垂直阵lcmv波束形成的水下目标辐射噪声测量方法 | |
CN107179535A (zh) | 一种基于畸变拖曳阵的保真增强波束形成的方法 | |
CN103076594B (zh) | 一种基于互相关的水声脉冲信号双阵元定位的方法 | |
CN111123192B (zh) | 一种基于圆形阵列和虚拟扩展的二维doa定位方法 | |
CN109799495B (zh) | 一种用于高保真阵列处理的宽带时延估计方法 | |
CN103064077B (zh) | 一种基于幅值加权的目标检测方法和设备 | |
CN101813772B (zh) | 一种快速宽带频域扩展拖曳阵波束形成方法 | |
Zeng et al. | Eigenspace-based minimum variance beamformer combined with Wiener postfilter for medical ultrasound imaging | |
CN108549052B (zh) | 一种时频-空域联合加权的圆谐域伪声强声源定位方法 | |
CN109581388B (zh) | 一种实时三维成像声纳的近场宽视角波束形成方法 | |
CN111025273B (zh) | 一种畸变拖曳阵线谱特征增强方法及*** | |
CN109725285B (zh) | 一种基于mvdr协方差矩阵元素自适应相角转换的doa估计方法 | |
CN109541548A (zh) | 一种基于匹配场的空气声呐定位方法 | |
CN108761394A (zh) | 一种基于时空处理的高分辨低旁瓣去卷积谱估计方法 | |
CN109814065A (zh) | 基于相位因子加权的波束形成方法 | |
CN105334435A (zh) | 一种基于任意阵形的自适应局部放电超声监测方法 | |
CN109541526A (zh) | 一种利用矩阵变换的圆环阵方位估计方法 | |
CN108957389A (zh) | 一种实数域多通道信号目标方位估计方法 | |
CN111427045B (zh) | 基于分布式多输入-多输出声纳的水下目标逆散射成像方法 | |
CN109541572B (zh) | 一种基于线性环境噪声模型的子空间方位估计方法 | |
Zhao et al. | Passive estimation of target velocity based on cross-spectrum histogram | |
CN111722178B (zh) | 一种基于指向性模型数值求解的远场窄带信号来波方向估计方法 | |
CN102183755A (zh) | 一种基于柯西-高斯模型的新型高分辨方位估计方法 | |
CN113075645B (zh) | 一种基于主成分分析-密度聚类的畸变阵形线谱增强方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200731 |