CN110954860A - 一种doa和极化参数估计方法 - Google Patents
一种doa和极化参数估计方法 Download PDFInfo
- Publication number
- CN110954860A CN110954860A CN201911305688.3A CN201911305688A CN110954860A CN 110954860 A CN110954860 A CN 110954860A CN 201911305688 A CN201911305688 A CN 201911305688A CN 110954860 A CN110954860 A CN 110954860A
- Authority
- CN
- China
- Prior art keywords
- updating
- doa
- sparse
- vector
- polarization
- 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
- 230000010287 polarization Effects 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 41
- 239000013598 vector Substances 0.000 claims abstract description 47
- 230000014509 gene expression Effects 0.000 claims abstract description 28
- 230000005855 radiation Effects 0.000 claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000006870 function Effects 0.000 claims description 6
- 230000003595 spectral effect Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- 239000000654 additive Substances 0.000 claims description 3
- 230000000996 additive effect Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000013398 bayesian method Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
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
- 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
- G01S3/143—Systems for determining direction or deviation from predetermined direction by vectorial combination of signals derived from differently oriented antennae
-
- 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/78—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 electromagnetic waves other than radio waves
- G01S3/782—Systems for determining direction or deviation from predetermined direction
-
- 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/8027—By vectorial composition of signals received by plural, differently-oriented transducers
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于非网格块稀疏贝叶斯学习的DOA和极化参数估计方法,包括利用矢量传感器阵列接收信号的导向矢量构建非网格信号模型;在稀疏贝叶斯学习框架下构建块稀疏向量;对块稀疏向量施加二阶稀疏分层先验;计算隐含变量和超参数的更新表达式;求解隐含变量和超参数更新结果;对源信号进行稀疏重构,求得目标辐射源的DOA和极化参数估计。本发明通过构建块稀疏向量并对其施加二阶分层先验促进块间稀疏和块内稀疏,提高重构精度,进而提高估计性能,解决了现有技术在非理想环境下测向精度差的问题。
Description
技术领域
本发明属于阵列信号处理技术领域,具体涉及一种DOA和极化参数估计方法。
背景技术
DOA估计是阵列信号处理领域的研究热点,广泛应用于雷达、声纳、无线通信等实际应用系统。相比于传统的标量阵列,矢量传感器阵列能够充分利用入射信号的空域信息和极化信息,有助于实现高精度DOA估计。
基于矢量传感器阵列的测向方法主要可以分为:子空间类和稀疏重构类。子空间类代表方法有:极化-MUSIC方法、极化-ESPRIT方法和四阶累积量类方法,该类方法在低信噪比、小快怕数等非理想情况下测向性能不理想;基于矢量传感器阵列的稀疏重构类方法目前研究较少,代表方法有:信号重构、加权“组-lasso”和稀疏贝叶斯方法。
现有稀疏重构类方法假设目标辐射源恰巧落在划分好的网格上,然而对于实际测向系统来说上述假设是不合理的,并且以上方法在实现过程中没有考虑块稀疏结构。
发明内容
本发明所要解决的技术问题是针对上述现有技术的不足,提供一种基于非网格分层块稀疏贝叶斯理论的DOA和极化参数估计方法,在非理想情况下仍能具有良好的估计性能。
为实现上述技术目的,本发明采取的技术方案为:
一种DOA和极化参数估计方法,包括:
步骤1:基于矢量传感器阵列,根据源信号导向矢量的一阶泰勒展开式构建非网格信号模型;
步骤2:基于步骤1构建的非网格信号模型,在稀疏贝叶斯学习框架下构建块稀疏向量;
步骤3:对步骤2构建的块稀疏向量施加二阶稀疏分层先验;
步骤4:计算隐含变量和超参数的更新表达式;
步骤5:基于步骤4的更新表达式,求解隐含变量和超参数更新结果;
步骤6:根据步骤5的更新结果,对源信号进行稀疏重构,求得目标辐射源的DOA和极化参数估计。
为优化上述技术方案,采取的具体措施还包括:
上述的步骤1包括:
步骤1.1:获取信号空域采样数据:
设M为双极化矢量传感器阵列的阵元数,K为信源数;
对于极化方向d,天线阵列接收信号向量为:
步骤1.2:构建非网格信号模型:
上述的步骤2包括:
基于步骤1构建的非网格信号模型,对X[d]进行向量化处理:
上述的步骤3包括:
对步骤2构建的块稀疏向量施加二阶稀疏分层先验:
第二层为两个服从Gamma分布的超先验:
上述的步骤4为:基于变分贝叶斯理论将后验分布的概率密度函数进行变分近似,计算各隐含变量和超参数的更新表达式:
步骤4.4:更新ν[d]:
q(ν[d])服从Gamma分布,ν[d]的更新表达式为:
步骤4.6:更新Δθ:
通过最小化似然函数,Δθ的更新表达式:
上述的步骤5为:
根据步骤4.1-4.6,基于KL散度收敛原则,交替迭代更新各隐含变量和超参数直至求得更新结果。
上述的步骤6包括:
步骤6.3:根据DOA估计结果,估计极化参数,极化辅助角和极化相位差的估计结果分别为:
本发明具有以下有益效果:
不同于传统的子空间类方法和现有的基于网格化的稀疏重构类方法,本发明的一种基于非网格块稀疏贝叶斯学习的DOA和极化参数估计方法,构建块稀疏向量并施加二阶分层稀疏先验,能够促进块间稀疏度和内部稀疏度,同时减少重构误差;在低信噪比、小快拍数下仍具有良好的估计精度。
附图说明
图1是本发明的流程示意图;
图2是相同条件下本发明方法、稀疏重构方法(DPE-SR)、长矢量MUSIC(LV-MUSIC)方法的测向性能以及CRB下界对比图。
具体实施方式
以下结合附图对本发明的实施例作进一步详细描述。
参照图1,本发明的一种DOA和极化参数估计方法,包括:
步骤1:基于矢量传感器阵列,根据源信号导向矢量的一阶泰勒展开式构建非网格信号模型:
步骤1.1:获取信号空域采样数据:
设M为双极化矢量传感器阵列的阵元数,K为信源数;
对于极化方向d,天线阵列接收信号向量为:
步骤1.2:构建非网格信号模型:
步骤2:基于步骤1构建的非网格信号模型,在稀疏贝叶斯学习框架下构建块稀疏向量,包括:
基于步骤1构建的非网格信号模型,对X[d]进行向量化处理:
步骤3:对步骤2构建的块稀疏向量施加二阶稀疏分层先验:
第二层为两个服从Gamma分布的超先验:
步骤4:计算隐含变量和超参数的更新表达式:
基于变分贝叶斯理论将后验分布的概率密度函数进行变分近似,计算各隐含变量和超参数的更新表达式:
步骤4.4:更新ν[d]:
q(ν[d])服从Gamma分布,ν[d]的更新表达式为:
步骤4.6:更新Δθ:
通过最小化似然函数,Δθ的更新表达式:
步骤5:基于步骤4的更新表达式,求解隐含变量和超参数更新结果:
根据步骤4.1-4.6,基于KL散度收敛原则,交替迭代更新各隐含变量和超参数直至求得更新结果。
步骤6:根据步骤5的更新结果,对源信号进行稀疏重构,求得目标辐射源的DOA和极化参数估计,包括:
步骤6.1:重构源信号分量:
步骤6.2:DOA估计:
步骤6.3:极化参数估计:
根据DOA估计结果,极化辅助角和极化相位差的估计结果分别为:
图2是相同条件下本发明提出的非网格分块稀疏贝叶斯方法、稀疏重构方法(DPE-SR)、长矢量MUSIC(LV-MUSIC)方法的测向性能以及CRB下界对比图。由图2可知,在其他条件相同的情况下,相比于DPE-SR和LV-MUSIC方法,本发明方法具有更好的估计精度,特别在低信噪比(0dB)时,优势更为明显。
综上,本发明公开了一种基于非网格块稀疏贝叶斯学习的DOA和极化参数估计方法,通过构建块稀疏向量并对其施加二阶分层先验促进块间稀疏和块内稀疏,提高重构精度,进而提高估计性能,解决了现有技术在非理想环境下测向精度差的问题。
以上仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,应视为本发明的保护范围。
Claims (7)
1.一种DOA和极化参数估计方法,其特征在于,包括:
步骤1:基于矢量传感器阵列,根据源信号导向矢量的一阶泰勒展开式构建非网格信号模型;
步骤2:基于步骤1构建的非网格信号模型,在稀疏贝叶斯学习框架下构建块稀疏向量;
步骤3:对步骤2构建的块稀疏向量施加二阶稀疏分层先验;
步骤4:计算隐含变量和超参数的更新表达式;
步骤5:基于步骤4的更新表达式,求解隐含变量和超参数更新结果;
步骤6:根据步骤5的更新结果,对源信号进行稀疏重构,求得目标辐射源的DOA和极化参数估计。
2.根据权利要求1所述的一种DOA和极化参数估计方法,其特征在于,所述步骤1包括:
步骤1.1:获取信号空域采样数据:
设M为双极化矢量传感器阵列的阵元数,K为信源数;
对于极化方向d,天线阵列接收信号向量为:
步骤1.2:构建非网格信号模型:
6.根据权利要求5所述的一种DOA和极化参数估计方法,其特征在于,所述步骤5为:
根据步骤4.1-4.6,基于KL散度收敛原则,交替迭代更新各隐含变量和超参数直至求得更新结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911305688.3A CN110954860B (zh) | 2019-12-18 | 2019-12-18 | 一种doa和极化参数估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911305688.3A CN110954860B (zh) | 2019-12-18 | 2019-12-18 | 一种doa和极化参数估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110954860A true CN110954860A (zh) | 2020-04-03 |
CN110954860B CN110954860B (zh) | 2021-06-29 |
Family
ID=69982378
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911305688.3A Active CN110954860B (zh) | 2019-12-18 | 2019-12-18 | 一种doa和极化参数估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110954860B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113325364A (zh) * | 2021-07-15 | 2021-08-31 | 金陵科技学院 | 一种基于数据压缩的空时联合测向方法 |
CN113625219A (zh) * | 2021-06-23 | 2021-11-09 | 沈阳航空航天大学 | 基于稀疏贝叶斯学习的矢量共形阵列doa-极化参数联合估计方法 |
CN114428225A (zh) * | 2022-04-06 | 2022-05-03 | 中国电子科技集团公司第二十九研究所 | 一种基于准匹配滤波的多辐射源到达角度测量方法及装置 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7274622B1 (en) * | 2005-05-23 | 2007-09-25 | The United States Of America Represented By The Secretary Of The Navy | Nonlinear techniques for pressure vector acoustic sensor array synthesis |
EP2293094A1 (en) * | 2009-09-01 | 2011-03-09 | Fujitsu Limited | Method of estimating direction of arrival and apparatus thereof |
CN103941220A (zh) * | 2014-04-25 | 2014-07-23 | 电子科技大学 | 一种基于稀疏重构的网格外目标波达方向估计方法 |
CN107436421A (zh) * | 2017-07-24 | 2017-12-05 | 哈尔滨工程大学 | 一种稀疏贝叶斯学习框架下混合信号doa估计方法 |
CN107817465A (zh) * | 2017-10-12 | 2018-03-20 | 中国人民解放军陆军工程大学 | 超高斯噪声背景下的基于无网格压缩感知的doa估计方法 |
CN108445462A (zh) * | 2018-02-05 | 2018-08-24 | 江苏大学 | 一种基于稀疏贝叶斯学习的双基地mimo雷达的dod和doa估计方法 |
CN109061551A (zh) * | 2018-08-02 | 2018-12-21 | 西北工业大学 | 一种基于多项式求根的无网格稀疏谱估计方法 |
CN109307855A (zh) * | 2018-08-02 | 2019-02-05 | 西北工业大学 | 基于网格误差模型的无网格稀疏近似最小方差doa估计方法 |
CN109444810A (zh) * | 2018-12-24 | 2019-03-08 | 哈尔滨工程大学 | 一种非负稀疏贝叶斯学习框架下的互质阵列非网格doa估计方法 |
CN110031796A (zh) * | 2019-02-28 | 2019-07-19 | 重庆工业职业技术学院 | 一种三维多快拍无网格压缩波束形成声源识别方法 |
-
2019
- 2019-12-18 CN CN201911305688.3A patent/CN110954860B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7274622B1 (en) * | 2005-05-23 | 2007-09-25 | The United States Of America Represented By The Secretary Of The Navy | Nonlinear techniques for pressure vector acoustic sensor array synthesis |
EP2293094A1 (en) * | 2009-09-01 | 2011-03-09 | Fujitsu Limited | Method of estimating direction of arrival and apparatus thereof |
CN103941220A (zh) * | 2014-04-25 | 2014-07-23 | 电子科技大学 | 一种基于稀疏重构的网格外目标波达方向估计方法 |
CN107436421A (zh) * | 2017-07-24 | 2017-12-05 | 哈尔滨工程大学 | 一种稀疏贝叶斯学习框架下混合信号doa估计方法 |
CN107817465A (zh) * | 2017-10-12 | 2018-03-20 | 中国人民解放军陆军工程大学 | 超高斯噪声背景下的基于无网格压缩感知的doa估计方法 |
CN108445462A (zh) * | 2018-02-05 | 2018-08-24 | 江苏大学 | 一种基于稀疏贝叶斯学习的双基地mimo雷达的dod和doa估计方法 |
CN109061551A (zh) * | 2018-08-02 | 2018-12-21 | 西北工业大学 | 一种基于多项式求根的无网格稀疏谱估计方法 |
CN109307855A (zh) * | 2018-08-02 | 2019-02-05 | 西北工业大学 | 基于网格误差模型的无网格稀疏近似最小方差doa估计方法 |
CN109444810A (zh) * | 2018-12-24 | 2019-03-08 | 哈尔滨工程大学 | 一种非负稀疏贝叶斯学习框架下的互质阵列非网格doa估计方法 |
CN110031796A (zh) * | 2019-02-28 | 2019-07-19 | 重庆工业职业技术学院 | 一种三维多快拍无网格压缩波束形成声源识别方法 |
Non-Patent Citations (3)
Title |
---|
ANUP DAS: "Theoretical and Experimental Comparison of Off-Grid Sparse Bayesian Direction-of-Arrival Estimation Algorithms", 《IEEE ACCESS》 * |
冯晓敏等: "基于变分稀疏贝叶斯的off-grid DOA 估计", 《海南大学学报自然科学版》 * |
赵嫔姣: "混合信号的波达方向估计算法研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113625219A (zh) * | 2021-06-23 | 2021-11-09 | 沈阳航空航天大学 | 基于稀疏贝叶斯学习的矢量共形阵列doa-极化参数联合估计方法 |
CN113625219B (zh) * | 2021-06-23 | 2023-10-13 | 沈阳航空航天大学 | 基于稀疏贝叶斯学习的矢量共形阵列doa-极化参数联合估计方法 |
CN113325364A (zh) * | 2021-07-15 | 2021-08-31 | 金陵科技学院 | 一种基于数据压缩的空时联合测向方法 |
CN114428225A (zh) * | 2022-04-06 | 2022-05-03 | 中国电子科技集团公司第二十九研究所 | 一种基于准匹配滤波的多辐射源到达角度测量方法及装置 |
CN114428225B (zh) * | 2022-04-06 | 2022-06-14 | 中国电子科技集团公司第二十九研究所 | 一种基于准匹配滤波的多辐射源到达角度测量方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN110954860B (zh) | 2021-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110954860B (zh) | 一种doa和极化参数估计方法 | |
CN107436421B (zh) | 一种稀疏贝叶斯学习框架下混合信号doa估计方法 | |
CN109298383B (zh) | 一种基于变分贝叶斯推断的互质阵波达方向角估计方法 | |
CN111337893B (zh) | 一种基于实值稀疏贝叶斯学习的离格doa估计方法 | |
CN103941220B (zh) | 一种基于稀疏重构的网格外目标波达方向估计方法 | |
CN110109050B (zh) | 嵌套阵列下基于稀疏贝叶斯的未知互耦的doa估计方法 | |
CN110197112B (zh) | 一种基于协方差修正的波束域Root-MUSIC方法 | |
CN111337873B (zh) | 一种基于稀疏阵的doa估计方法 | |
CN109061554A (zh) | 一种基于空间离散网格动态更新的目标到达角度估计方法 | |
CN113835063B (zh) | 一种无人机阵列幅相误差与信号doa联合估计方法 | |
CN109239649A (zh) | 一种阵列误差条件下的互质阵列doa估计新方法 | |
CN111880143B (zh) | 改进稀疏贝叶斯学习的高精度定位方法、存储介质及设备 | |
CN113671439A (zh) | 基于非均匀智能超表面阵列的无人机集群测向系统及方法 | |
CN114415110B (zh) | 一种非负稀疏贝叶斯学习的直接定位方法 | |
CN114462454B (zh) | 基于深度学习的无网格信号源doa估计方法 | |
CN115356729A (zh) | 一种近场非均匀采样的直接稀疏成像方法 | |
CN109696651B (zh) | 一种基于m估计的低快拍数下波达方向估计方法 | |
CN117523344B (zh) | 一种基于相位质量加权卷积神经网络的干涉相位解缠方法 | |
CN112946564B (zh) | 基于dnn的波束空间的doa估计方法、装置及计算机存储介质 | |
CN110364829B (zh) | 一种稀疏直线阵天线 | |
CN117195482A (zh) | 一种大规模阵列无网格doa估计方法及系统 | |
CN114415106B (zh) | 基于改进lamp网络的互耦阵列doa估计方法 | |
CN116027264A (zh) | 一种字典矩阵迭代优化的somp离网格直接定位方法 | |
CN109298381A (zh) | 一种基于变分贝叶斯推断的互质阵相干信号方位角估计方法 | |
CN114019445A (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 |