CN113341371B - 一种基于l阵和二维esprit算法的doa估计方法 - Google Patents
一种基于l阵和二维esprit算法的doa估计方法 Download PDFInfo
- Publication number
- CN113341371B CN113341371B CN202110596959.6A CN202110596959A CN113341371B CN 113341371 B CN113341371 B CN 113341371B CN 202110596959 A CN202110596959 A CN 202110596959A CN 113341371 B CN113341371 B CN 113341371B
- Authority
- CN
- China
- Prior art keywords
- matrix
- array
- subarray
- dimensional
- signal
- 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
本发明属于雷达信号处理技术领域,具体涉及一种基于L阵和二维ESPRIT算法的DOA估计方法。本发明提出了一种基于L阵和二维ESPRIT算法的DOA估计方法,在阵列接收到辐射源信号后,利用L阵的阵列特点对阵列数据进行处理,构造用于二维ESPRIT算法运算的数据矩阵,得到估计角度,完成信号源的方向测量估计。本发明的有益效果为,本发明利用L阵构造大矩阵,得到的矩阵形式简单,在后续的测向过程中保持了较低的计算复杂度,同时有效去除了统计独立噪声源的影响,改善了在高斯噪声环境下的分辨能力和统计性能。
Description
技术领域
本发明属于雷达信号处理技术领域,具体涉及一种基于L阵和二维ESPRIT算法的DOA估计方法。
背景技术
在雷达信号处理中,基于L阵实现对目标雷达辐射源信号的波达方向(directionof arrival,DOA)估计一直是研究的重点,然而随着现代电磁频谱的日益密集、脉冲密度越来越高以及大功率电子设备的应用,电磁环境越来越复杂,传统的基于MUSIC算法的DOA估计方法面临着测向精度不高,计算复杂度大的问题,同时,这些方法大多数研究仅适用于一维DOA估计,针对实际空间下二维角度的测向能力有所欠缺,在L阵上实现二维DOA估计仍然是一个比较大的挑战。所以寻求新的技术和手段是当务之急,迫切需要提升雷达接收阵列的信号处理能力、测向精确性能。
ESPRIT算法其原理在于利用数据协方差矩阵信号子空间的旋转不变特性求解信号的入射角等信息,由于不需要进行谱峰搜索,其复杂度和计算量与MUSIC算法相比大幅降低,在DOA估计中得到了广泛的应用。针对基于L阵的二维测向,此前也有一些二维ESPRIT方法提出,但这些方法都存在构造形式复杂,计算量高的问题,且不能很好的抑制噪声影响,在进行雷达信号处理时不能很好地平衡性能与计算成本之间的关系,阵列测向结果始终会存在较大误差。
发明内容
本发明的目的,是针对上述问题,提出了一种基于L阵和二维ESPRIT算法的DOA估计方法,在阵列接收到辐射源信号后,利用L阵的阵列特点对阵列数据进行处理,构造用于二维ESPRIT算法运算的数据矩阵,得到估计角度,完成信号源的方向测量估计。
本发明的技术方案为:
将L阵划分为四个子阵,通过对空间中目标信号接收数据进行处理,计算L阵中三个无共用阵元的子阵之间的互协方差矩阵,利用三个互协方差矩阵组成一个用于L阵子空间分解的大矩阵,根据旋转不变性得到两个旋转矩阵,对两个矩阵基于二维ESPRIT算法进行处理,通过特征分解得到两个旋转矩阵的特征值并进行特征值配对,最后求解得到信号的方位角和仰角,完成二维DOA估计。
本发明基于L阵的信号接收模型如图1所示,信号si(t)的波达方向(DOA)用矢量表示:其中θ为方位角,范围为θ∈(-π,π),定义为波达方向在xoy平面的投影与x正半轴的夹角,为仰角,范围为定义为波达方向在yoz平面的投影与z正半轴的夹角。相邻阵元之间的间距为d。
为了简化分析过程和保证所构建的数学模型的合理性,本发明中基于L阵的信号接收模型基于如下假设:(1)不考虑接收天线阵列的阵元通道幅相误差;(2)阵列接收信号均为远场点源发出的窄带信号,信号的中心频率相同且已知,信号波长λ大于阵元间距d的两倍,即λ≥2d;(3)接收信号的相位是随机的,不考虑严格非圆信号等已知信号部分先验信息下的特殊情况。
如图2所示,均匀L阵可划分为四个子阵:x轴上的前M个阵元组成子阵X1,后M个阵元组成子阵X2;y轴上的前M个阵元组成子阵Y1,后M个阵元组成子阵Y2。
子阵X1、X2和Y1、Y2接收的数据矩阵分别为:
其中,xl、yl分别表示阵列X与阵列Y标号为l(l=0,1,2,…,M,M+1)的阵元接收的数据。
基于L阵和二维ESPRIT的DOA估计方法,主要包括以下步骤:
S1、利用子阵X1、X2和Y1、Y2接收的数据矩阵构造得到三个互协方差矩阵:
由于子阵X1和Y2、子阵X2和Y1、子阵X2和Y2之间均没有共用阵元,且各阵元接收的噪声相互独立,因此得到的三个互协方差矩阵中都不存在噪声项。
S2、将步骤S1中三个互协方差矩阵堆叠成一个大矩阵CL,矩阵的构造形式为
S3、对步骤S2中构造的大矩阵CL进行奇异值分解,将前N个大特征值对应的特征向量组成的矩阵作为信号子空间的估计得到信号子空间的估计ES。
其中,E1、E2和E3为组成矩阵ES的三个M×M维块矩阵;
S4、计算两个旋转矩阵ΨX和ΨY,计算公式为
对这两个旋转矩阵进行特征分解,得到
其中,λ1,…,λN和γ1,…,γN分别是旋转矩阵ΨX和ΨY的特征值,ΦX和ΦY是特征值矩阵,T1和T2是由旋转矩阵ΨX和ΨY的特征向量组成的正交矩阵。
S5、基于二维ESPRIT算法进行二维角度参数配对,确定两个旋转矩阵ΨX和ΨY的特征值之间的对应关系,具体步骤为:
S52、取出矩阵的对角元素u1,…,uN,并对这些元素取复数相角,并按照相角大小,从大到小对对角元素进行排序,得到排序后的对角元素取矩阵ΨY特征值γ1,…,γN的相角,并按照相角的大小按从大到小对γ1,…,γN进行排序,得到排序后的特征值顺序按照排序后的序号得到配对关系
S53、根据步骤S52中对角元素与矩阵T1中特征向量的对应关系,调整ΨX的特征值顺序(若在矩阵中的行序为j,那么其对应特征向量为矩阵T1的第j个列向量,对应的ΨX的特征值为λj),得到调整之后矩阵ΨX的特征值以及配对关系为
S54、根据步骤S52、S53得到的两组配对关系,得到矩阵ΨX和ΨY的特征值之间的配对关系为
其中,
本发明的有益效果为,本发明利用L阵构造大矩阵,得到的矩阵形式简单,在后续的测向过程中保持了较低的计算复杂度,同时有效去除了统计独立噪声源的影响,改善了在高斯噪声环境下的分辨能力和统计性能。
附图说明
图1为L阵信号接收示意图;
图2为L阵子阵划分示意图;
图3为RMSE随SNR变化曲线;
图4为RMSE随子阵数阵元数变化曲线;
图5为仿真运行时间随子阵阵元数变化曲线。
具体实施方式
下面将结合附图和仿真,对本发明的性能进行说明。
利用Matlab对本发明所提方法进行仿真验证。计算机仿真环境:MicrosoftWindows 10操作***,Matlab 2020a软件,AMD R7-4800U处理器(支持AVX指令集),16GBDDR4-3200内存。
(一)本发明方法的DOA测向估计性能
仿真1:方法的方向精度随信号信噪比变化情况
空间中的L阵阵元总数为13,四个子阵的阵元数均为7,仿真信号数据长度为500个快拍,信号的来波方向为(60°,60°)。不同信噪比下均进行1000次Monte-Carlo实验,得到本发明方法的角度估计RMSE随SNR变化曲线如图3所示。
结论分析:从仿真结果来看,本发明L阵ESPRIT新方法的具有较高的波达方向估计精度,性能上的损失较少,测向效果颇佳。
仿真2:方法的方向精度随接收阵列子阵阵元数变化情况
空间中的L形阵列,四个子阵的阵元数均为M,仿真数据长度为200个快拍,信噪比为-10dB,信号的来波方向为(60°,60°),不同子阵数M下进行500次Monte-Carlo实验得到本发明ESPRIT算法的RMSE随子阵阵元数M的变化曲线如图4所示。
结论分析:图4的结果表明,本发明方法的统计性能随子阵阵元数增多而增强。
(二)本发明方法的运行时间
仿真1:方法的仿真时间
空间中的L形阵列,四个子阵的阵元数均为M,信号的来波方向为(30°,30°),仿真数据长度为500个快拍,信噪比为0dB。对本发明方法进行仿真,得到不同子阵阵元数下该算法进行1000次Monte-Carlo实验的仿真时长。
将仿真时长数据绘制成曲线图得到该算法仿真时长随子阵阵元数变化曲线如图5所示。
结论分析:从图5的结果来看,随着子阵阵元数的增长,本发明算法的时长也随之增长,但增长幅度较小,这说明本发明计算成本较小,计算复杂度低,有着计算上的优势。
综合以上仿真结果,本发明所提方法在计算成本和估计性能之间取得了更好的平衡。
Claims (1)
1.一种基于L阵和二维ESPRIT算法的DOA估计方法,建立基于L阵的信号接收模型,信号si(t)的波达方向(DOA)用矢量表示,其中θ为方位角,范围为θ∈(-π,π),定义为波达方向在xoy平面的投影与x正半轴的夹角,为仰角,范围为定义为波达方向在yoz平面的投影与z正半轴的夹角,相邻阵元之间的间距为d;均匀L阵划分为四个子阵:x轴上的前M个阵元组成子阵X1,后M个阵元组成子阵X2;y轴上的前M个阵元组成子阵Y1,后M个阵元组成子阵Y2;子阵X1、X2和Y1、Y2接收的数据矩阵分别为:
其中xl、yl分别表示阵列X与阵列Y标号为l的阵元接收的数据,l=0,1,2,…,M,M+1;
其特征在于,DOA估计方法包括以下步骤:
S1、利用子阵X1、X2和Y1、Y2接收的数据矩阵构造得到三个互协方差矩阵:
由于子阵X1和Y2、子阵X2和Y1、子阵X2和Y2之间均没有共用阵元,且各阵元接收的噪声相互独立,因此得到的三个互协方差矩阵中都不存在噪声项;
S2、将步骤S1中三个互协方差矩阵堆叠成一个矩阵CL,矩阵的构造形式为
S3、对步骤S2中构造的矩阵CL进行奇异值分解,将前N个特征值对应的特征向量组成的矩阵作为信号子空间的估计得到信号子空间的估计ES:
其中,E1、E2和E3为组成矩阵ES的三个M×M维块矩阵;
S4、计算两个旋转矩阵ΨX和ΨY,计算公式为
对这两个旋转矩阵进行特征分解,得到
其中,λ1,…,λN和γ1,…,γN分别是旋转矩阵ΨX和ΨY的特征值,ΦX和ΦY是特征值矩阵,T1和T2是由旋转矩阵ΨX和ΨY的特征向量组成的正交矩阵;
S5、基于二维ESPRIT算法进行二维角度参数配对,确定两个旋转矩阵ΨX和ΨY的特征值之间的对应关系,具体步骤为:
S52、取出矩阵的对角元素u1,…,uN,并对这些元素取复数相角,并按照相角大小,从大到小对对角元素进行排序,得到排序后的对角元素取矩阵ΨY特征值γ1,…,γN的相角,并按照相角的大小按从大到小对γ1,…,γN进行排序,得到排序后的特征值顺序按照排序后的序号得到配对关系
S53、根据步骤S52中对角元素与矩阵T1中特征向量的对应关系,调整ΨX的特征值顺序,调整方式为,若在矩阵中的行序为j,那么其对应特征向量为矩阵T1的第j个列向量,对应的ΨX的特征值为λj,得到调整之后矩阵ΨX的特征值以及配对关系为
S54、根据步骤S52、S53得到的两组配对关系,得到矩阵ΨX和ΨY的特征值之间的配对关系为
其中,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110596959.6A CN113341371B (zh) | 2021-05-31 | 2021-05-31 | 一种基于l阵和二维esprit算法的doa估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110596959.6A CN113341371B (zh) | 2021-05-31 | 2021-05-31 | 一种基于l阵和二维esprit算法的doa估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113341371A CN113341371A (zh) | 2021-09-03 |
CN113341371B true CN113341371B (zh) | 2022-03-08 |
Family
ID=77472091
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110596959.6A Active CN113341371B (zh) | 2021-05-31 | 2021-05-31 | 一种基于l阵和二维esprit算法的doa估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113341371B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115308680A (zh) * | 2022-08-04 | 2022-11-08 | 中国科学院微小卫星创新研究院 | 二维doa估计方法、***及计算机可读介质 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6278406B1 (en) * | 1998-03-24 | 2001-08-21 | Nec Corporation | Direction finder and device for processing measurement results for the same |
CN102253363A (zh) * | 2011-03-29 | 2011-11-23 | 西安交通大学 | 基于l型阵列的相干信号二维波达方向估计装置及其方法 |
CN102279381A (zh) * | 2011-03-29 | 2011-12-14 | 西安交通大学 | 一种基于l型阵列二维波达方向估计的装置及其方法 |
CN106054123A (zh) * | 2016-06-06 | 2016-10-26 | 电子科技大学 | 一种稀疏l阵及其二维doa估计方法 |
CN106872935A (zh) * | 2017-03-20 | 2017-06-20 | 北京理工大学 | 一种基于四元数的电磁矢量传感器阵列波达方向估计方法 |
CN108663653A (zh) * | 2018-05-17 | 2018-10-16 | 西安电子科技大学 | 基于l形电磁矢量传感器阵列的波达方向估计方法 |
CN110286351A (zh) * | 2019-07-12 | 2019-09-27 | 电子科技大学 | 一种基于l型嵌套阵的二维doa估计方法及装置 |
CN110286350A (zh) * | 2019-07-12 | 2019-09-27 | 电子科技大学 | 一种l型稀疏阵doa估计的精确配对方法及装置 |
CN110954859A (zh) * | 2019-11-22 | 2020-04-03 | 宁波大学 | 基于l型阵列的二维非相干分布式非圆信号参数估计方法 |
CN111175690A (zh) * | 2019-10-29 | 2020-05-19 | 宁波大学 | 联合对角化l型mimo雷达圆与非圆混合测向方法 |
CN111983554A (zh) * | 2020-08-28 | 2020-11-24 | 西安电子科技大学 | 非均匀l阵下的高精度二维doa估计 |
CN112763972A (zh) * | 2020-12-30 | 2021-05-07 | 长沙航空职业技术学院 | 基于稀疏表示的双平行线阵二维doa估计方法及计算设备 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8446253B2 (en) * | 2009-03-11 | 2013-05-21 | Checkpoint Systems, Inc. | Localization using virtual antenna arrays in modulated backscatter RFID systems |
CN114185042A (zh) * | 2015-09-17 | 2022-03-15 | 松下电器产业株式会社 | 雷达装置 |
US10871543B2 (en) * | 2018-06-12 | 2020-12-22 | Kaam Llc | Direction of arrival estimation of acoustic-signals from acoustic source using sub-array selection |
-
2021
- 2021-05-31 CN CN202110596959.6A patent/CN113341371B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6278406B1 (en) * | 1998-03-24 | 2001-08-21 | Nec Corporation | Direction finder and device for processing measurement results for the same |
CN102253363A (zh) * | 2011-03-29 | 2011-11-23 | 西安交通大学 | 基于l型阵列的相干信号二维波达方向估计装置及其方法 |
CN102279381A (zh) * | 2011-03-29 | 2011-12-14 | 西安交通大学 | 一种基于l型阵列二维波达方向估计的装置及其方法 |
CN106054123A (zh) * | 2016-06-06 | 2016-10-26 | 电子科技大学 | 一种稀疏l阵及其二维doa估计方法 |
CN106872935A (zh) * | 2017-03-20 | 2017-06-20 | 北京理工大学 | 一种基于四元数的电磁矢量传感器阵列波达方向估计方法 |
CN108663653A (zh) * | 2018-05-17 | 2018-10-16 | 西安电子科技大学 | 基于l形电磁矢量传感器阵列的波达方向估计方法 |
CN110286351A (zh) * | 2019-07-12 | 2019-09-27 | 电子科技大学 | 一种基于l型嵌套阵的二维doa估计方法及装置 |
CN110286350A (zh) * | 2019-07-12 | 2019-09-27 | 电子科技大学 | 一种l型稀疏阵doa估计的精确配对方法及装置 |
CN111175690A (zh) * | 2019-10-29 | 2020-05-19 | 宁波大学 | 联合对角化l型mimo雷达圆与非圆混合测向方法 |
CN110954859A (zh) * | 2019-11-22 | 2020-04-03 | 宁波大学 | 基于l型阵列的二维非相干分布式非圆信号参数估计方法 |
CN111983554A (zh) * | 2020-08-28 | 2020-11-24 | 西安电子科技大学 | 非均匀l阵下的高精度二维doa估计 |
CN112763972A (zh) * | 2020-12-30 | 2021-05-07 | 长沙航空职业技术学院 | 基于稀疏表示的双平行线阵二维doa估计方法及计算设备 |
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 |
---|---|
CN113341371A (zh) | 2021-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109633525B (zh) | 一种快速的均匀圆阵空间谱测向实现方法 | |
CN107092004B (zh) | 基于信号子空间旋转不变性的互质阵列波达方向估计方法 | |
CN108761381B (zh) | 一种基于阵列天线的目标信源探测方法 | |
CN109633522B (zh) | 基于改进的music算法的波达方向估计方法 | |
CN108663653B (zh) | 基于l形电磁矢量传感器阵列的波达方向估计方法 | |
CN110244272B (zh) | 基于秩一去噪模型的波达方向估计方法 | |
CN109633520B (zh) | 一种均匀圆阵超分辨率空间谱估计方法 | |
CN111239678A (zh) | 一种基于l型阵列的二维doa估计方法 | |
Roy et al. | Comparative performance of ESPRIT and MUSIC for direction-of-arrival estimation | |
CN109696657B (zh) | 一种基于矢量水听器的相干声源定位方法 | |
CN106526531A (zh) | 基于三维天线阵列的改进传播算子二维doa估计算法 | |
CN112130111A (zh) | 一种大规模均匀十字阵列中单快拍二维doa估计方法 | |
CN108120953A (zh) | 一种基于波达方向估计的无线电定位方法 | |
CN111983554A (zh) | 非均匀l阵下的高精度二维doa估计 | |
CN111693947A (zh) | 基于互质阵列doa估计的改进music方法 | |
CN112255629A (zh) | 基于联合uca阵列的序贯esprit二维不相干分布源参数估计方法 | |
CN113341371B (zh) | 一种基于l阵和二维esprit算法的doa估计方法 | |
Hou et al. | A robust capon beamforming approach for sparse array based on importance resampling compressive covariance sensing | |
CN113759303B (zh) | 一种基于粒子群算法的无网格波达角估计方法 | |
CN113296050B (zh) | 基于各向异性阵列的极化和角度参数联合估计方法 | |
CN112327292B (zh) | 一种二维稀疏阵列doa估计方法 | |
CN112363108B (zh) | 信号子空间加权超分辨的波达方向检测方法及*** | |
Khudhair et al. | Estimation of direction of arrival for antenna array based on esprit and multiple signal classification algorithms | |
WO2010066306A1 (en) | Apparatus and method for constructing a sensor array used for direction of arrival (doa) estimation | |
CN113126021B (zh) | 一种基于三平行线阵的单快拍二维doa估计方法 |
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 |