CN111337893B - 一种基于实值稀疏贝叶斯学习的离格doa估计方法 - Google Patents
一种基于实值稀疏贝叶斯学习的离格doa估计方法 Download PDFInfo
- Publication number
- CN111337893B CN111337893B CN201911317378.3A CN201911317378A CN111337893B CN 111337893 B CN111337893 B CN 111337893B CN 201911317378 A CN201911317378 A CN 201911317378A CN 111337893 B CN111337893 B CN 111337893B
- Authority
- CN
- China
- Prior art keywords
- matrix
- real
- grid
- representing
- data
- 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
Images
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- 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/02—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 radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于实值稀疏贝叶斯学习的离格DOA估计方法,属于雷达信号处理领域,通过构造一个粗糙的离格网格,利用酉矩阵将复数流型矩阵转换为实数矩阵,将DOA估计问题转化到实数域,极大地降低了计算复杂度,同时采用固定的步长来更新网格点的位置,有效的提高了DOA估计的精度,并结合奇异值分解(Singular Value Decomposition,SVD)来降低矩阵维数以减少计算量,同时将采样网格视为可调参数,使用固定的步长来更新网格位置,并利用期望最大化(Expectation Maximization,EM)算法迭代地更新采样网格,从而有效地提高估计精度。
Description
技术领域
本发明属于雷达信号处理领域,涉及雷达信号的角度估计,具体地说是一种基于实值稀疏贝叶斯学习的适用于雷达信号波达方向估计的方法。
背景技术
近几十年来,波达方向(Direction of Arrival,DOA)估计一直是雷达信号处理的一个重要研究内容,针对雷达信号波达方向估计的问题,人们提出了大量行之有效的方法。稀疏贝叶斯学习(Sparse Bayesian Learning,SBL)是近些年来稀疏信号恢复领域内比较流行的一种方法,例如在文献:J.Dai,X.Bao,W.Xu,and C.Chang,“Root sparse Bayesianlearning for off-grid DOA estimation,”IEEE Signal Processing Letters,vol.24,no.1, pp.46-50,2017.中,提出了一种求根SBL的方法,将网格作为可调参数,并通过求解多项式的根来更新网格点。然而大多数现有的基于SBL的DOA估计方法都是在复数域进行求解,运算量比较大,计算复杂度也比较高,如何在降低计算复杂度的同时保持较高的估计精度是一个亟待解决的问题。
发明内容
针对现有方法的不足,本发明提出了一种基于实值稀疏贝叶斯学习的离格(off-grid) DOA估计方法,通过构造一个粗糙的离格网格,利用酉矩阵将复数流型矩阵转换为实数矩阵,并结合奇异值分解(Singular Value Decomposition,SVD)来降低矩阵维数以减少计算量,同时将采样网格视为可调参数,使用固定的步长来更新网格位置,并利用期望最大化(Expectation Maximization,EM)算法迭代地更新采样网格,从而有效地提高估计精度。
用于实现本发明的技术解决方案包括如下步骤:
步骤1:接收***接收到的雷达信号经过匹配滤波后,在接收机处得到在t时刻包含K组DOA信息的数据向量y(t)。
步骤5:设置迭代次数计数变量i=1,初始化背景噪声的精度α、信号方差向量δ以及角度偏移值ε。
步骤6:利用期望最大化准则,更新背景噪声精度α和信号方差向量δ。
步骤7:更新角度偏移值ε。
步骤8:利用步骤(7)中求出的角度偏移值ε更新网格点。
步骤10:对方差向量δ进行谱峰搜索,得到K个极大值点对应的角度,即为目标角度的最终估计值。
本发明的有益效果:
本发明提出了一种基于实值稀疏贝叶斯学习的离格DOA估计方法,通过构造一个粗糙的离格网格,利用期望最大化算法迭代地进行网格细化,不断地更新网格点的位置,从而消除由离格模型引起的误差,并通过酉变换,将DOA估计问题转化到实数域,极大地降低了计算复杂度,同时采用固定的步长来更新网格点的位置,有效的提高了DOA 估计的精度。
附图说明
图1是本发明实施流程图;
图2是500次蒙特卡洛实验条件下,快拍数T=30,网格间隔r=3的情况下,检测2个目标时本发明与求根SBL方法估计DOA的均方根误差(RMSE)比较图。
图3是500次蒙特卡洛实验条件下,信噪比SNR为10dB,快拍数T=30的情况下,检测2个目标时本发明与求根SBL方法估计DOA的运算时间比较图。
具体实施方式
下面结合附图对本发明作进一步说明。
如图1所示,本发明提出的一种基于实值稀疏贝叶斯学习的离格DOA估计方法,其实施包括如下步骤:
(1)接收***接收到的雷达信号经过匹配滤波后,在接收机处得到在t时刻包含K组DOA信息的数据向量y(t),t=1,2,...,T,其中T表示快拍数。
M表示接收端的天线根数,(·)T表示矩阵的转置,d表示相邻阵元之间的距离,λ表示信号的波长,diag(·)表示取对角运算,εl表示网格点上的角度偏移值,s(t)表示t时刻接收信号在上的向量表示,n(t)表示t时刻的零均值高斯白噪声。进而,将数据向量y(t)扩展到T快拍上,可以得到新的数据矩阵:
(3)定义一个酉矩阵QM:
其中:的列表示K个最大奇异值所对应的奇异向量,是由K个最大奇异值作为对角线元素的对角矩阵,Σc与Σs类似,是奇异值分解的特定结构,且奇异值分解的目的是利用分解得到的矩阵Vs来降低矩阵的维度,继而用Vs右乘矩阵得到降维的数据模型:
(6)利用期望最大化(Expectation Maximization,EM)准则,更新背景噪声精度α和信号方差向量δ:
其中:a,b为极小的正常数(例如:a=b=0.0001),rt表示数据矩阵R中的第t列向量,Δ=diag(δ),(·)T表示矩阵转置,tr(·)表示矩阵的迹,Ξt=μtμt T+Σ,[·]ll表示矩阵的第l个对角线元素。
(7)更新角度偏移值ε:
ε=P-1v,
(10)对方差向量δ进行谱峰搜索,得到K个极大值点对应的角度,即为目标角度的最终估计值。
下面结合仿真实验对本发明的效果做进一步说明。
为了评估本方法的性能,考虑一均匀线性阵列***,阵元间距为电磁波的半波长,假设远场有两个相互独立的目标,分别随机取自范围[-30°,-20°],[10°,20°]。在所有的试验中,背景噪声均假设为高斯白噪声,蒙特卡洛实验次数为500。
实验条件
实验1,采用本发明与求根SBL方法,在阵元个数M=6,快拍数为T=30,网格间隔r=3的情况下,信噪比SNR从-10dB到10dB变化时,估计DOA的均方根误差(RMSE) 的比较实验,仿真结果如图2所示。
实验2,采用本发明与求根SBL方法,在阵元个数M=8,信噪比SNR为10dB,快拍数为T=30的情况下,网格间隔r分别取1°,2°,4°,6°,8°,10°时,估计DOA的运算时间的比较实验,仿真结果如图3所示。
实验分析
从图2可以看出,随着信噪比的增加,所有方法估计目标角度的RMSE都显著降低,但本发明相较于求根SBL方法有更好的估计性能。
从图3可以看出,随着网格间隔变大,所有方法估计目标角度的运算时间都显著降低,但本发明的运算时间明显低于求根SBL方法。
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明,它们并非用以限制本发明的保护范围,凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。
Claims (1)
1.一种基于实值稀疏贝叶斯学习的离格DOA估计方法,其特征在于,包括如下步骤:
步骤1:接收***接收到的雷达信号经过匹配滤波后,在接收机处得到在t时刻包含K组DOA信息的数据向量y(t);
步骤5:设置迭代次数计数变量i=1,初始化背景噪声的精度α、信号方差向量δ以及角度偏移值ε;
步骤6:利用期望最大化准则,更新背景噪声精度α和信号方差向量δ;
步骤7:更新角度偏移值ε;
步骤8:利用步骤7中求出的角度偏移值ε更新网格点;
步骤10:对方差向量δ进行谱峰搜索,得到K个极大值点对应的角度,即为目标角度的最终估计值;
所述步骤2中y(t)一阶泰勒展开的近似表达式为:
M表示接收端的天线根数,(·)T表示矩阵的转置,d表示相邻阵元之间的距离,λ表示信号的波长,diag(·)表示取对角运算,εl表示网格点上的角度偏移值,s(t)表示t时刻接收信号在上的向量表示,n(t)表示t时刻的零均值高斯白噪声;
所述步骤2中新的数据矩阵Y的表达式为:
所述步骤3中酉矩阵QM的表达式为:
其中:IM表示一个M×M维的单位矩阵,JM表示一个M×M维的逆向单位矩阵;
所述步骤6中更新背景噪声精度α和信号方差向量δ的方法如下:
所述步骤7中更新角度偏移值ε的方法如下:
ε=P-1v,
所述步骤8中更新网格点的方法如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911317378.3A CN111337893B (zh) | 2019-12-19 | 2019-12-19 | 一种基于实值稀疏贝叶斯学习的离格doa估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911317378.3A CN111337893B (zh) | 2019-12-19 | 2019-12-19 | 一种基于实值稀疏贝叶斯学习的离格doa估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111337893A CN111337893A (zh) | 2020-06-26 |
CN111337893B true CN111337893B (zh) | 2022-09-16 |
Family
ID=71183319
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911317378.3A Active CN111337893B (zh) | 2019-12-19 | 2019-12-19 | 一种基于实值稀疏贝叶斯学习的离格doa估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111337893B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111880143B (zh) * | 2020-07-07 | 2023-05-23 | 西安交通大学 | 改进稀疏贝叶斯学习的高精度定位方法、存储介质及设备 |
CN112132760B (zh) * | 2020-09-14 | 2024-02-27 | 北京大学 | 基于可学习可微分矩阵逆及矩阵分解的图像恢复方法 |
CN112924947B (zh) * | 2020-11-28 | 2023-10-27 | 海南大学 | 一种基于实值稀疏贝叶斯学习的mimo雷达稳健doa估计方法 |
CN112948606B (zh) * | 2020-12-14 | 2022-10-21 | 西南交通大学 | 一种基于自适应网格的信号估计方法及装置 |
CN113609450B (zh) * | 2021-02-22 | 2023-11-03 | 天津大学 | 基于LCGAMP网络和1-Bit量化的大规模稀疏阵列DoA估计方法 |
CN113281795B (zh) * | 2021-05-06 | 2023-08-01 | 国网福建省电力有限公司 | 基于北斗和超宽带***中联合稀疏恢复室内定位方法 |
CN113534040B (zh) * | 2021-05-31 | 2023-08-11 | 河海大学 | 基于加权二阶稀疏贝叶斯的相干信源离栅doa估计方法 |
CN113447815B (zh) * | 2021-07-08 | 2022-08-23 | 江苏大学 | 一种基于实值esprit的电机故障在线检测方法及*** |
CN113673317B (zh) * | 2021-07-12 | 2023-04-07 | 电子科技大学 | 基于原子范数最小化可降维的二维离格doa估计方法 |
CN115267654B (zh) * | 2022-07-29 | 2024-05-10 | 青岛理工大学 | 一种基于酉变换的求根稀疏渐近最小方差离格方位估计方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104537249A (zh) * | 2015-01-04 | 2015-04-22 | 西安电子科技大学 | 基于稀疏贝叶斯学习的波达方向角估计方法 |
CN108337199A (zh) * | 2018-01-17 | 2018-07-27 | 江苏大学 | 一种基于稀疏贝叶斯学习的大规模mimo通信***的下行链路信道估计方法 |
CN108459296A (zh) * | 2018-01-17 | 2018-08-28 | 江苏大学 | 一种基于离格稀疏贝叶斯学习的nested阵列波达方向估计方法 |
CN109490819A (zh) * | 2018-11-16 | 2019-03-19 | 南京邮电大学 | 一种基于稀疏贝叶斯学习的离格波达方向估计方法 |
CN109639604A (zh) * | 2019-01-10 | 2019-04-16 | 江苏大学 | 一种基于实值稀疏贝叶斯学习的大规模mimo通信***的下行链路信道估计方法 |
-
2019
- 2019-12-19 CN CN201911317378.3A patent/CN111337893B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104537249A (zh) * | 2015-01-04 | 2015-04-22 | 西安电子科技大学 | 基于稀疏贝叶斯学习的波达方向角估计方法 |
CN108337199A (zh) * | 2018-01-17 | 2018-07-27 | 江苏大学 | 一种基于稀疏贝叶斯学习的大规模mimo通信***的下行链路信道估计方法 |
CN108459296A (zh) * | 2018-01-17 | 2018-08-28 | 江苏大学 | 一种基于离格稀疏贝叶斯学习的nested阵列波达方向估计方法 |
CN109490819A (zh) * | 2018-11-16 | 2019-03-19 | 南京邮电大学 | 一种基于稀疏贝叶斯学习的离格波达方向估计方法 |
CN109639604A (zh) * | 2019-01-10 | 2019-04-16 | 江苏大学 | 一种基于实值稀疏贝叶斯学习的大规模mimo通信***的下行链路信道估计方法 |
Non-Patent Citations (3)
Title |
---|
Real-Valued Sparse Bayesian Learning Approach for Massive MIMO Channel Estimation;Lei Zhou 等;《IEEE WIRELESS COMMUNICATIONS LETTERS》;20191113;全文 * |
基于稀疏贝叶斯学习的高效DOA估计方法;孙磊 等;《电子与信息学报》;20130531;第35卷(第5期);全文 * |
基于酉变换和稀疏贝叶斯学习的离格DOA估计;高阳 等;《通信学报》;20170630;第38卷(第6期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111337893A (zh) | 2020-06-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111337893B (zh) | 一种基于实值稀疏贝叶斯学习的离格doa估计方法 | |
CN109490819B (zh) | 一种基于稀疏贝叶斯学习的离格波达方向估计方法 | |
CN109655799B (zh) | 基于iaa的协方差矩阵向量化的非均匀稀疏阵列测向方法 | |
WO2018094565A1 (zh) | 脉冲噪声下的波束成形方法及装置 | |
CN110109050B (zh) | 嵌套阵列下基于稀疏贝叶斯的未知互耦的doa估计方法 | |
CN109061554B (zh) | 一种基于空间离散网格动态更新的目标到达角度估计方法 | |
CN108663653B (zh) | 基于l形电磁矢量传感器阵列的波达方向估计方法 | |
CN110244272B (zh) | 基于秩一去噪模型的波达方向估计方法 | |
CN109375154B (zh) | 一种冲击噪声环境下基于均匀圆阵的相干信号参数估计方法 | |
CN111337873B (zh) | 一种基于稀疏阵的doa估计方法 | |
CN113376569B (zh) | 基于最大似然的嵌套阵稀疏表示波达方向估计方法 | |
CN108957390B (zh) | 一种存在互耦时基于稀疏贝叶斯理论的到达角估计方法 | |
CN107544051A (zh) | 嵌套阵列基于k‑r子空间的波达方向估计方法 | |
CN113032721B (zh) | 一种低计算复杂度的远场和近场混合信号源参数估计方法 | |
CN113835063B (zh) | 一种无人机阵列幅相误差与信号doa联合估计方法 | |
CN108398659B (zh) | 一种矩阵束与求根music结合的波达方向估计方法 | |
CN112731273A (zh) | 一种基于稀疏贝叶斯的低复杂度信号波达方向估计方法 | |
CN106125039B (zh) | 基于局域联合处理的改进空时自适应单脉冲测角方法 | |
CN111273269A (zh) | 基于ipso-bp的频率分集阵列的雷达目标定位方法 | |
CN108594165B (zh) | 一种基于期望最大化算法的窄带信号波达方向估计方法 | |
CN113671439A (zh) | 基于非均匀智能超表面阵列的无人机集群测向***及方法 | |
CN109917330B (zh) | 一种存在相位误差时基于稀疏正交匹配追踪理论的到达角估计方法 | |
CN117092585A (zh) | 单比特量化DoA估计方法、***和智能终端 | |
CN116699511A (zh) | 一种多频点信号波达方向估计方法、***、设备及介质 | |
CN110927663A (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 |