CN113126021B - 一种基于三平行线阵的单快拍二维doa估计方法 - Google Patents
一种基于三平行线阵的单快拍二维doa估计方法 Download PDFInfo
- Publication number
- CN113126021B CN113126021B CN202110416816.2A CN202110416816A CN113126021B CN 113126021 B CN113126021 B CN 113126021B CN 202110416816 A CN202110416816 A CN 202110416816A CN 113126021 B CN113126021 B CN 113126021B
- Authority
- CN
- China
- Prior art keywords
- matrix
- dimensional
- signal
- snapshot
- array
- 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
- 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
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于电子对抗技术领域,具体涉及一种基于三平行线阵的单快拍二维DOA估计方法。本发明采用的技术方案是:将三平行线阵接收的单快拍数据与参考阵元数据相乘来构造三个满秩的伪矩阵,基于所构造的三个伪矩阵之间存在的旋转不变性,利用二维ESPRIT方法求解得到两个旋转矩阵,通过对旋转矩阵的特征分解实现二维角度参数配对并求解入射信号的方位角和仰角,完成二维DOA估计。本发明的有益效果为,本发明可以利用单个快拍的阵列接收数据对信号的DOA进行二维估计,算法拥有较低的自由度损失,在一定信噪比和阵列接收条件下拥有较好的分辨能力和统计性能。
Description
技术领域
本发明属于电子对抗技术领域,具体涉及一种基于三平行线阵的单快拍二维DOA估计方法。
背景技术
在电子对抗中,阵列能够截获到用于DOA估计的有效信息往往是不足的,极限情况下只有一个快拍。例如,随着跳频通信和短时通信技术在军事通信中的广泛使用,具有反侦察属性的跳频(F-H)信号、突发(Burst)信号由于跳变方式未知或持续时间短,阵列只能截获到极少个快拍。此外,机载、舰载等高机动平台对测向系统的实时性要求很高,无法通过长时间采样来积累快拍数据,只能利用单个快拍的数据。随着空间电磁对抗形势日趋激烈,这类极限单快拍下的DOA估计问题越来越普遍。
MUSIC、ESPRIT等传统的DOA估计方法均基于多快拍信号假设,在单快拍下由于数据不足以支撑对子空间进行准确划分,估计性能极差甚至失效。而目前大部分基于空间平滑或者构造伪协方差矩阵的单快拍DOA估计方法都无可避免地存在较大的阵列孔径和自由度损失,一些文献提出了基于严格非圆信号特征的伪协方差矩阵构造方法,该方法虽无自由度损失,但是要求所有信号具有零初始相位,这个约束条件在实践中是几乎无法实现的。此外,单快拍下由于需要对阵列数据进行复杂的数学重构,大多数单快拍DOA估计研究基于一维线性阵列,只能实现一维DOA估计,在较为复杂的二维平面阵上实现单快拍二维DOA估计仍然是一个比较大的挑战。
发明内容
针对上述问题,本发明提出了一种基于三平行线阵的单快拍二维DOA估计方法。
本发明采用的技术方案是:将三平行线阵接收的单快拍数据与参考阵元数据相乘来构造三个满秩的伪矩阵,基于所构造的三个伪矩阵之间存在的旋转不变性,利用二维ESPRIT方法求解得到两个旋转矩阵,通过对旋转矩阵的特征分解实现二维角度参数配对并求解入射信号的方位角和仰角,完成二维DOA估计。
本发明基于的三平行线阵信号接收模型示意图见图1,三平行线阵由三条平行的线阵组成,位于中间的线阵阵元数为M+1,其余两线阵的阵元数为M。信号的波达方向(DOA)用矢量表示:其中θ为方位角,范围为θ∈(-π,π),定义为波达方向在xoy平面的投影与x正半轴的夹角,为仰角,范围为定义为波达方向在yoz平面的投影与z正半轴的夹角。
为了简化分析过程和保证所构建的数学模型的合理性,本发明中的阵列信号接收模型基于如下假设:(1)不考虑接收天线阵列的阵元通道辐相误差;(2)阵列接收信号均为远场点源发出的窄带信号,信号的中心频率相同且已知,信号波长λ大于阵元间距d的两倍,即λ≥2d;(3)接收信号的相位是随机的,不考虑严格非圆信号等已知信号部分先验信息下的特殊情况。
其中
X=[x0,x1,…,xM-1]T
Y=[y0,y1,…,yM]T
Z=[z0,z1,…,zM-1]T
基于三平行线阵的单快拍二维DOA估计方法的具体步骤如下:
S1、利用三平行线阵各阵元接收的单快拍数据构造三个满秩的伪矩阵,三个伪矩阵的构造形式分别为:
S4、利用步骤S3中得到的信号子空间的估计ES,计算两个旋转矩阵Ψ1和Ψ2,计算公式为
其中
E1、E2和E3为组成矩阵ES的三个M×N维块矩阵;
S5、对步骤S4中得到的两个旋转矩阵分别进行特征分解,得到相应的特征向量矩阵和特征值对角矩阵,特征分解过程如下
其中,T1表示矩阵Ψ1特征分解后得到的特征向量矩阵,δ1,…,δN表示矩阵Ψ1的特征值,T2表示矩阵Ψ2特征分解后得到的特征向量矩阵,γ1,…,γN表示矩阵Ψ2的特征值;
S6、进行二维角度参数配对,确定两个旋转矩阵Ψ1和Ψ2的特征值之间的对应关系,具体步骤为:
S62、取出矩阵的对角元素u1,…,uN,并对所有对角元素取复数相角,并按照相角大小,从大到小对对角元素进行排序,得到排序后的对角元素取矩阵Ψ2特征值γ1,…,γN的相角,并按照相角的大小按从大到小对γ1,…,γN进行排序,得到排序后的特征值顺序按照排序后的序号得到配对关系
S63、根据步骤S61中对角元素与矩阵T1中特征向量的对应关系,调整Ψ1的特征值顺序(若在矩阵中的行序为j,那么其对应特征向量为矩阵T1的第j个列向量,对应的Ψ1的特征值为δj),得到调整之后矩阵Ψ1的特征值以及配对关系为
S64、根据步骤S62、S63得到的两组配对关系,得到矩阵Ψ1和Ψ2的特征值之间的配对关系为
其中
本发明的有益效果为,本发明可以利用单个快拍的阵列接收数据对信号的波达方向(DOA)进行二维估计,算法拥有较低的自由度损失,在一定信噪比和阵列接收条件下拥有较好的分辨能力和统计性能。
附图说明
图1为三平行线阵信号接收模型示意图;
图2为信噪比为10dB时二维DOA估计散布图;
图3为信噪比为20dB时二维DOA估计散布图;
图4为二维DOA估计均方误差随信噪比变化曲线,其中线(a)为等刻度纵坐标;(b)为对数刻度纵坐标。
具体实施方式
结合计算机仿真实验对本发明提出的方法的实施过程和效果进行描述。
仿真中阵列接收的噪声类型为零均值高斯白噪声,且信号与噪声之间、不同阵元的噪声之间均是统计独立的;且仿真的单快拍信号的初始相位是随机的。
(一)本发明方法的目标分辨能力
仿真场景一:空间中的三平行线阵阵元总数为37,其中中间一条线阵的阵元数为13,其余两条线阵阵元数为12,利用本发明方法进行DOA估计,信噪比为10dB时进行2000次Monte-Carlo实验得到二维DOA估计散布图如图2所示,
仿真场景二:接收阵列结构同仿真场景二,信噪比为20dB时进行2000次Monte-Carlo实验得到二维DOA估计散布图如图3所示。
从图2、图3的两幅散布图来看,多次独立实验的估计值均能够围绕着真实值的位置形成三个点“簇”,“簇”与“簇”之间可明显区分,本发明方法能够对接收到的三独立信源单快拍信号在DOA上进行区分。仿真结果验证了本发明提出的基于三平行线阵的单快拍测向算法的有效性。
(二)本发明方法的统计性能
通过二维DOA估计(单位为度)角度值的均方误差(RMSE)对算法的统计性能进行评估,二维角度估计值均方误差定义为
仿真场景三:空间中的三平行线阵阵元总数为31,其中中间一条线阵的阵元数为11,其余两条线阵阵元数为10,信号的二维DOA为利用本发明方法进行DOA估计,在5~40dB范围内以5dB为步长改变信噪比,不同信噪比下进行1000次Monte-Carlo实验,得到本发明方法的二维DOA估计均方误差随信噪比变化曲线如图4所示,其中图4(a)为等刻度纵坐标下的曲线,图4(b)为对数刻度纵坐标下的曲线。
从图4的曲线来看,本发明方法的二维角度估计误差随信噪比的提高而减小,本发明提出的统计性能与信噪比呈正相关。
Claims (1)
1.一种基于三平行线阵的单快拍二维DOA估计方法,三平行线阵中位于中间的线阵阵元数为M+1,其余两线阵的阵元数为M,假设三平行线阵接收到N个波达方向不同的单快拍信号为S=[s1,s2,…,sN]T,各信号的方位角为θ=[θ1,θ2,…,θN],仰角为三平行线阵接收的单快拍数据矩阵为
其中
X=[x0,x1,…,xM-1]T
Y=[y0,y1,…,yM]T
Z=[z0,z1,…,zM-1]T
其特征在于,所述DOA估计方法包括以下步骤:
S1、利用三平行线阵各阵元接收的单快拍数据构造三个满秩的伪矩阵,三个伪矩阵的构造形式分别为:
E1、E2和E3为组成矩阵ES的三个M×N维块矩阵;
S4、利用步骤S3中得到的信号子空间的估计ES,计算两个旋转矩阵Ψ1和Ψ2,计算公式为:
S5、对步骤S4中得到的两个旋转矩阵分别进行特征分解,得到相应的特征向量矩阵和特征值对角矩阵,特征分解过程如下
其中,T1表示矩阵Ψ1特征分解后得到的特征向量矩阵,Φ1表示矩阵Ψ1特征分解后得到的特征值对角矩阵,δ1,…,δN表示矩阵Ψ1的特征值,T2表示矩阵Ψ2特征分解后得到的特征向量矩阵,Φ2表示矩阵Ψ2特征分解后得到的特征值对角矩阵,γ1,…,γN表示矩阵Ψ2的特征值;
S6、进行二维角度参数配对,确定两个旋转矩阵Ψ1和Ψ2的特征值之间的对应关系,具体步骤为:
S62、取出矩阵的对角元素u1,…,uN,并对所有对角元素取复数相角,按照相角大小,从大到小对对角元素进行排序,得到排序后的对角元素取矩阵Ψ2特征值γ1,…,γN的相角,并按照相角的大小按从大到小对γ1,…,γN进行排序,得到排序后的特征值顺序按照排序后的序号得到配对关系
S63、根据步骤S61中对角元素与矩阵T1中特征向量的对应关系,调整Ψ1的特征值顺序,调制原则是,若在矩阵中的行序为j,那么其对应特征向量为矩阵T1的第j个列向量,对应的Ψ1的特征值为δj,从而得到调整之后矩阵Ψ1的特征值以及配对关系为
S64、根据步骤S62、S63得到的两组配对关系,得到矩阵Ψ1和Ψ2的特征值之间的配对关系为
其中
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110416816.2A CN113126021B (zh) | 2021-04-19 | 2021-04-19 | 一种基于三平行线阵的单快拍二维doa估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110416816.2A CN113126021B (zh) | 2021-04-19 | 2021-04-19 | 一种基于三平行线阵的单快拍二维doa估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113126021A CN113126021A (zh) | 2021-07-16 |
CN113126021B true CN113126021B (zh) | 2022-03-29 |
Family
ID=76777493
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110416816.2A Active CN113126021B (zh) | 2021-04-19 | 2021-04-19 | 一种基于三平行线阵的单快拍二维doa估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113126021B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114879138B (zh) * | 2022-07-13 | 2022-09-27 | 广东大湾区空天信息研究院 | 一种毫米波雷达二维角度计算方法、装置及相关设备 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101325807A (zh) * | 2008-07-24 | 2008-12-17 | 中国人民解放军理工大学 | 信号波达方向估计方法 |
CN104931918A (zh) * | 2015-05-06 | 2015-09-23 | 集美大学 | 基于复数盲源分离的船舶doa方位估计方法 |
CN108037493A (zh) * | 2017-11-27 | 2018-05-15 | 西安电子科技大学 | 基于多基地雷达的杂波子空间下双门限恒虚警检测方法 |
CN108318877A (zh) * | 2017-11-27 | 2018-07-24 | 西安电子科技大学 | 基于多基地雷达的杂波子空间下目标检测方法 |
CN108957391A (zh) * | 2018-07-24 | 2018-12-07 | 北京理工大学 | 一种基于嵌套阵列的l型天线阵的二维波达方向估计方法 |
CN110244258A (zh) * | 2019-06-12 | 2019-09-17 | 南京航空航天大学 | 用于双平行阵二维测向中扩展doa矩阵方法 |
CN111368256A (zh) * | 2020-03-23 | 2020-07-03 | 电子科技大学 | 一种基于均匀圆阵的单快拍测向方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5320792B2 (ja) * | 2008-03-28 | 2013-10-23 | 富士通株式会社 | 到来方向推定装置、到来方向推定方法および到来方向推定プログラム |
JP5686920B1 (ja) * | 2014-07-02 | 2015-03-18 | 眞吉 西本 | アレイアンテナビーム幅内の量子化多重・狭ビーム形成方法、アレイアンテナビーム幅内の量子化多重・狭ビーム形成装置およびレーダシステム |
-
2021
- 2021-04-19 CN CN202110416816.2A patent/CN113126021B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101325807A (zh) * | 2008-07-24 | 2008-12-17 | 中国人民解放军理工大学 | 信号波达方向估计方法 |
CN104931918A (zh) * | 2015-05-06 | 2015-09-23 | 集美大学 | 基于复数盲源分离的船舶doa方位估计方法 |
CN108037493A (zh) * | 2017-11-27 | 2018-05-15 | 西安电子科技大学 | 基于多基地雷达的杂波子空间下双门限恒虚警检测方法 |
CN108318877A (zh) * | 2017-11-27 | 2018-07-24 | 西安电子科技大学 | 基于多基地雷达的杂波子空间下目标检测方法 |
CN108957391A (zh) * | 2018-07-24 | 2018-12-07 | 北京理工大学 | 一种基于嵌套阵列的l型天线阵的二维波达方向估计方法 |
CN110244258A (zh) * | 2019-06-12 | 2019-09-17 | 南京航空航天大学 | 用于双平行阵二维测向中扩展doa矩阵方法 |
CN111368256A (zh) * | 2020-03-23 | 2020-07-03 | 电子科技大学 | 一种基于均匀圆阵的单快拍测向方法 |
Non-Patent Citations (4)
Title |
---|
"2-D angle of arrival estimation with two parallel uniform linear arrays for coherent signals";T.Q.Xia等;《2007 IEEE radar conference》;20071231;第244-247页 * |
"一种基于双平行线阵相干源二维波达方向估计的新方法";曾操等;《雷达科学与技术》;20030831;第1卷(第2期);第104-108页 * |
"二维波达方向估计方法研究";夏铁骑;《中国优秀博硕士学位论文全文数据库(博士)信息科技辑》;20090415;第I136-34页 * |
"卫星干扰源定位及干扰抑制技术研究";王爱莹;《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》;20160315;第I136-1730页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113126021A (zh) | 2021-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wong et al. | Closed-form direction finding and polarization estimation with arbitrarily spaced electromagnetic vector-sensors at unknown locations | |
CN110824415B (zh) | 一种基于多发多收阵列的稀疏波达方向角度估计方法 | |
CN109490820B (zh) | 一种基于平行嵌套阵的二维doa估计方法 | |
CN109655799A (zh) | 基于iaa的协方差矩阵向量化的非均匀稀疏阵列测向方法 | |
Wang et al. | 2-D spatial smoothing for multipath coherent signal separation | |
CN107064926B (zh) | 空域色噪声背景下的双基地mimo雷达角度估算方法 | |
CN106526531A (zh) | 基于三维天线阵列的改进传播算子二维doa估计算法 | |
CN112130111A (zh) | 一种大规模均匀十字阵列中单快拍二维doa估计方法 | |
CN113189592A (zh) | 考虑幅相互耦误差的车载毫米波mimo雷达测角方法 | |
CN113126021B (zh) | 一种基于三平行线阵的单快拍二维doa估计方法 | |
CN113075649B (zh) | 一种适用于分布式网络化雷达的信号级直接定位方法 | |
Grundmann et al. | Investigation of direction of arrival estimation using characteristic modes | |
CN112711000B (zh) | 基于最小化准则的电磁矢量互质面阵张量功率谱估计方法 | |
Khudhair et al. | Estimation of direction of arrival for antenna array based on esprit and multiple signal classification algorithms | |
Ma et al. | DOD/DOA and polarization estimation in MIMO systems with spatially spread dipole quints | |
Zhang et al. | Estimation of fading coefficients in the presence of multipath propagation | |
CN113341371B (zh) | 一种基于l阵和二维esprit算法的doa估计方法 | |
Schurwanz et al. | Compressive sensing for direction-of-arrival estimation using an electronically steered multiple-input multiple-output array | |
Hong et al. | Comparison of MUSIC and ESPRIT for direction of arrival estimation of jamming signal | |
Nechaev et al. | The Research of the Digital Beamforming Algorithm for Optimal Noise Reduction in a Cylindrical Antenna Array with Directive Radiators | |
Kuroda et al. | DOA estimation and pairing method in 2D‐ESPRIT using triangular antenna array | |
CN113917389A (zh) | 一种相控阵协同探测系统和差波束角度估计方法 | |
Liu et al. | A semi-real-valued capon algorithm for 2-D DOA estimation of coherent signals with L-shape array | |
Taillefer et al. | ESPRIT algorithm for a seven-element regular-hexagonal shaped ESPAR antenna | |
CN113740797B (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 |