CN106093845B - 一种基于伪空间谱搜索的快速doa估计方法 - Google Patents
一种基于伪空间谱搜索的快速doa估计方法 Download PDFInfo
- Publication number
- CN106093845B CN106093845B CN201610365394.XA CN201610365394A CN106093845B CN 106093845 B CN106093845 B CN 106093845B CN 201610365394 A CN201610365394 A CN 201610365394A CN 106093845 B CN106093845 B CN 106093845B
- Authority
- CN
- China
- Prior art keywords
- matrix
- doa
- represent
- music
- quick
- 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.)
- Expired - Fee Related
Links
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/04—Details
- G01S3/12—Means for determining sense of direction, e.g. by combining signals from directional antenna or goniometer search coil with those from non-directional antenna
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估计方法,包括如下步骤:1:***接收到的雷达信号经过匹配滤波后,在接收机处得到的输出表示为x(tl),并构造一个新的数据矩阵X。2:计算X的协方差矩阵RX,并对其进行特征值分解,求得其噪声子空间UXN。3:定义矩阵J,构造一个新的数据矩阵Y=X+J*X;4:计算Y的协方差矩阵RY,并对其进行特征值分解,求得其噪声子空间UYN。5:构造基于矩阵Y的DOA的空间谱函数Pmusic(θ),并对其在进行谱峰搜索,找出N个极大值点对应的N个DOA角度的绝对值。6:构造基于矩阵X的DOA的空间谱函数P′music(θ),并将分别代入函数P′music(θ),然后得到并比较和的值,若较大则若较大则得到范围内的N个DOA角度值。本发明提出的方法,能使搜索范围减半,减小运算量,节省运算时间。
Description
技术领域
本发明属于雷达信号处理领域,涉及采用空间谱搜索方法的雷达信号的角度估计,具体地说是一种基于伪空间谱搜索的快速波达方向估计的方法。
背景技术
近几十年来,波达方向(Direction of Arrival,DOA)估计一直是空间谱估计的一个重要内容,在空间谱中利用谱峰搜索的方法进行目标角度估计时,其分辨率明显优于传统的常规波束形成(CBF)法。针对利用谱峰搜索方法的信号DOA估计问题,人们提出了大量行之有效的方法。例如在文献:Schmidt R O.Multiple emitter location and signalparameter estimation.IEEE Trans.on AP.1986.34(3):276~280中,提出了多重信号分类(MUSIC)算法;例如在文献:Stoica P,Nehorai A.MUSIC,maximum likelihood,andCramer-Rao bound:futher results and comparisons.In Proc.ICASSP,1989:2605~2608中,提出了加权MUSIC算法;这些算法具有很高的分辨力、估计精度和稳定性,然而这些方法在进行谱峰搜索时,其搜索区间为因此在运算过程中运算量较为大,搜索时间较长。
发明内容
针对现有方法的不足,本发明提出了一种新型的基于伪空间谱搜索的快速DOA估计方法,该方法在现有的MUSIC算法上进行改进,通过构造雷达信号的伪空间谱使得谱峰搜索范围减半,只需搜索就能得到范围内的DOA值,减小运算量,缩短运算时间。
用于实现本发明的技术解决方案包括如下步骤:
一种基于伪空间谱搜索的快速DOA估计方法,包括如下步骤:
步骤1:接收***接收到的雷达信号经过匹配滤波后,在接收机处得到的输出可以表示为x(tl),并构造一个新的数据矩阵X=[x(t1),x(t2),...,x(tL)];其中l=1,2,...,L;L表示快拍数;
步骤2:计算X的协方差矩阵RX=XXH,并对RX进行特征值分解,求得X的噪声子空间UXN;
步骤3:定义矩阵J,构造一个新的数据矩阵Y:Y=X+J*X;
步骤4:计算Y的协方差矩阵RY=YYH,并对RY进行特征值分解,求得其噪声子空间UYN;
步骤5:构造基于矩阵Y的DOA的空间谱函数Pmusic(θ),并对其在进行谱峰搜索,找出N个极大值点对应的角度,即为N个DOA角度估计值的绝对值,记为
步骤6:构造基于矩阵X的DOA的空间谱函数P′music(θ),并将分别代入函数P′music(θ),把正角度的函数值记为把负角度的函数值记为其中j=1,2,...,N;然后比较和的值,若较大则若较大则得到范围内的N个DOA角度值
进一步地,步骤1中x(tl)=As(tl)+n(tl);A=[a1,a2,...,an]表示M×N维的阵列流型矩阵,其中M表示发射阵列的阵元个数,N表示接收阵列的阵元个数;
导向矢量an=a(θn)=[1,exp(j2πdsinθn/λ),...,exp(j2π(M-1)dsinθn/λ)]T,其中θn表示第n组的真实DOA,n=1,2,...,N;d表示阵元间距,λ表示电磁波的波长,(·)T表示转置;
s(tl)表示tl时刻一个N维的发射信号向量,n(tl)表示tl时刻一个M维的零均值高斯白噪声。
进一步地,步骤2中所述的对RX进行特征值分解如下:
其中,ΣXS表示协方差矩阵RX的前N个最大特征值构成的对角矩阵,UXS表示与其对应的特征向量,ΣXN表示其余的M-N个特征值构成的对角矩阵,与其对应的特征向量为UXN,UXN为矩阵X的噪声子空间。
进一步地,步骤3中所述定义的矩阵J具体为:
进一步地,步骤4中所述对RY进行特征值分解如下:
其中,ΣYS表示协方差矩阵RY的N个最大特征值构成的对角矩阵,UYS表示与其对应的特征向量,ΣYN表示其余的M-N个特征值构成的对角矩阵,与其对应的特征向量为UYN,UYN为矩阵Y的噪声子空间。
进一步地,步骤5中所述构造的空间谱函数Pmusic(θ)的表达式如下:
其中:θ取值范围为
b(θ)=[1+exp(j2π(M-1)dsinθ/λ),exp(j2πdsinθ/λ)+exp(j2π(M-2)dsinθ/λ),...,exp(j2π(M-1)dsinθ/λ)+1]T。
进一步地,步骤6中所述构造基于矩阵X的函数P′music(θ)的表达式如下:
本发明的有益效果:
与现有方法相比,本发明提出了一种新型的基于伪空间谱搜索的快速DOA估计方法,使搜索范围减半,减小运算量,节省运算时间。
附图说明
图1是本发明实施流程图。
图2是在信噪比为10,快拍数为100次的情况下,本发明的空间谱函数Pmusic(θ)图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明。
如图1所示,本发明提出的快速DOA估计方法包括如下步骤:
(1)***接收到的雷达信号经过匹配滤波后,在接收机得到的输出表示为x(tl)=As(tl)+n(tl),l=1,2,...,L式中:
L表示快拍数;
A=[a1,a2,...,an]表示M×N维的阵列流型矩阵,其中M表示发射阵列的阵元个数,N表示接收阵列的阵元个数,其中
导向矢量an=a(θn)=[1,exp(j2πdsinθn/λ),...,exp(j2π(M-1)dsinθn/λ)]Tθn,n=1,2,...,N分别表示第n组的真实DOA,d表示阵元间距,λ表示电磁波的波长,(·)T表示转置;
s(tl)表示tl时刻一个N维的发射信号向量,n(tl)表示tl时刻一个M维的零均值高斯白噪声。
并构造一个M×L维的数据矩阵X:X=[x(t1),x(t2),...,x(tL)]。
(2)计算X的协方差矩阵RX=XXH,其中(·)H表示共轭转置。对RX进行特征值分解:
其中,ΣXS表示协方差矩阵RX的前N个最大特征值构成的对角矩阵,UXS表示与其对应的特征向量,ΣXN表示其余的M-N个特征值构成的对角矩阵,与其对应的特征向量为UXN,这里称为矩阵X的噪声子空间。
(3)定义矩阵J:
构造一个新的数据矩阵Y:Y=X+J*X。
(4)计算Y的协方差矩阵RY=YYH,其中(·)H表示共轭转置。对RY进行特征值分解:
其中,ΣYS表示协方差矩阵RY的N个最大特征值构成的对角矩阵,UYS表示与其对应的特征向量,ΣYN表示其余的M-N个特征值构成的对角矩阵,与其对应的特征向量为UYN,这里称为矩阵Y的噪声子空间。
(5)构造基于矩阵Y的DOA的空间谱函数Pmusic(θ):
其中:θ取值范围为
并对其在进行谱峰搜索,找出N个极大值点对应的角度,即为N个DOA角度的估计值的绝对值,将其记为
(6)构造基于矩阵X的函数P′music(θ):
其中θ取值范围为并将分别代入P′music(θ),把正角度的函数值记为把负角度的函数值记为其中j=1,2,...,N;然后比较和的值,若较大则若较大则得到范围内的N个DOA角度值。
实施例
下面结合仿真实验对本发明的效果做进一步说明。
为了评估本方法的性能,考虑***,发射阵列与接收列均为阵元,间距为电磁波的半波长的均匀线阵,发射阵列的阵元个数M=11,接收阵列的阵元个数N=10,假设远场有三个相互独立的目标,分别位于θ1=-32°,θ2=7°,θ3=23°。在所有的试验中,背景噪声假设为高斯白噪声,快拍数L=100。
实验条件
采用本发明在信噪比(SNR)为10dB时对目标角度进行角度估计,仿真结果如图2和表1所示。
表1
实验分析
从图2可以看出,本发明可以准确地估计出目标角度的绝对值j=1,2,3。
从表1可以看出,本发明中和的值具有较大的差异,可以通过比较和值精确地判别估计目标角度的真实值(正负号)。
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明,它们并非用以限制本发明的保护范围,凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于伪空间谱搜索的快速DOA估计方法,其特征在于,包括如下步骤:
步骤1:接收***接收到的雷达信号经过匹配滤波后,在接收机处得到的输出可以表示为x(tl),并构造一个新的数据矩阵X=[x(t1),x(t2),...,x(tL)];其中l=1,2,...,L;L表示快拍数;
步骤2:计算X的协方差矩阵RX=XXH,并对RX进行特征值分解,求得X的噪声子空间UXN;
步骤3:定义矩阵J,构造一个新的数据矩阵Y:Y=X+J*X;所述定义的矩阵J具体为:
步骤4:计算Y的协方差矩阵RY=YYH,并对RY进行特征值分解,求得其噪声子空间UYN;
步骤5:构造基于矩阵Y的DOA的空间谱函数Pmusic(θ),并对其在进行谱峰搜索,找出N个极大值点对应的角度,即为N个DOA角度估计值的绝对值,记为
步骤6:构造基于矩阵X的DOA的空间谱函数P′music(θ),并将分别代入函数P′music(θ),把正角度的函数值记为把负角度的函数值记为其中j=1,2,...,N;然后比较和的值,若较大则若较大则得到范围内的N个DOA角度值
2.根据权利要求1所述的一种基于伪空间谱搜索的快速DOA估计方法,其特征在于,步骤1中x(tl)的表达式为:x(tl)=As(tl)+n(tl);A=[a1,a2,...,an]表示M×N维的阵列流型矩阵,其中M表示发射阵列的阵元个数,N表示接收阵列的阵元个数;
导向矢量an=a(θn)=[1,exp(j2πdsinθn/λ),...,exp(j2π(M-1)dsinθn/λ)]T,其中θn表示第n组的真实DOA,n=1,2,...,N;d表示阵元间距,λ表示电磁波的波长,(·)T表示转置;
s(tl)表示tl时刻一个N维的发射信号向量,n(tl)表示tl时刻一个M维的零均值高斯白噪声。
3.根据权利要求1所述的一种基于伪空间谱搜索的快速DOA估计方法,其特征在于,步骤2中所述的对RX进行特征值分解如下:
其中,ΣXS表示协方差矩阵RX的前N个最大特征值构成的对角矩阵,与其对应的特征向量为UXS,ΣXN表示其余的M-N个特征值构成的对角矩阵,与其对应的特征向量为UXN,UXN为矩阵X的噪声子空间。
4.根据权利要求1所述的一种基于伪空间谱搜索的快速DOA估计方法,其特征在于,步骤4中所述对RY进行特征值分解如下:
其中,ΣYS表示协方差矩阵RY的N个最大特征值构成的对角矩阵,与其对应的特征向量为UYS,ΣYN表示其余的M-N个特征值构成的对角矩阵,与其对应的特征向量为UYN,UYN为矩阵Y的噪声子空间。
5.根据权利要求1所述的一种基于伪空间谱搜索的快速DOA估计方法,其特征在于,步骤5中所述构造的空间谱函数Pmusic(θ)的表达式如下:
其中:
d表示阵元间距。
6.根据权利要求2所述的一种基于伪空间谱搜索的快速DOA估计方法,其特征在于,步骤6中所述构造基于矩阵X的函数P′music(θ)的表达式如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610365394.XA CN106093845B (zh) | 2016-05-27 | 2016-05-27 | 一种基于伪空间谱搜索的快速doa估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610365394.XA CN106093845B (zh) | 2016-05-27 | 2016-05-27 | 一种基于伪空间谱搜索的快速doa估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106093845A CN106093845A (zh) | 2016-11-09 |
CN106093845B true CN106093845B (zh) | 2018-06-26 |
Family
ID=57230227
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610365394.XA Expired - Fee Related CN106093845B (zh) | 2016-05-27 | 2016-05-27 | 一种基于伪空间谱搜索的快速doa估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106093845B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109270484B (zh) * | 2018-07-24 | 2023-01-06 | 南京航空航天大学 | 一种基于运动合成阵列的多信源doa估计方法 |
CN109407047B (zh) * | 2018-09-19 | 2022-08-23 | 江苏大学 | 一种基于秩损求根的幅相误差校准和波达方向估计方法 |
CN109781791A (zh) * | 2019-02-22 | 2019-05-21 | 广西大学 | 基于空间谱估计算法的电阻抗成像方法 |
CN110082741B (zh) * | 2019-03-14 | 2022-07-15 | 哈尔滨工程大学 | 一种基于伪数据重构的超分辨波达角估计算法 |
CN110133588B (zh) * | 2019-05-14 | 2021-08-06 | 普联技术有限公司 | 一种天线定位方法、装置及设备 |
CN113391257B (zh) * | 2020-03-13 | 2023-04-18 | 光宝科技新加坡私人有限公司 | 用于物件的角度估测的运算装置及物件的角度估测方法 |
CN116973835B (zh) * | 2023-07-29 | 2024-01-30 | 同方工业有限公司 | 一种多信号混叠的伪空间谱累积测向方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104698433A (zh) * | 2015-03-16 | 2015-06-10 | 电子科技大学 | 基于单快拍数据的相干信号doa估计方法 |
-
2016
- 2016-05-27 CN CN201610365394.XA patent/CN106093845B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104698433A (zh) * | 2015-03-16 | 2015-06-10 | 电子科技大学 | 基于单快拍数据的相干信号doa估计方法 |
Non-Patent Citations (5)
Title |
---|
A NOVEL METHOD IN COMPENSATING RANDOM PERTURBATIONS IN EIGEN-BASED SUBSPACE METHOD;K.M.Cheung 等;《Proceedings of IEEE International Symposium on Circuits and Systems》;19941231;37-40 * |
An improved fast Root-MUSIC algorithm for DOA estimation;Xiong Jing 等;《2012 International Conference on Image Analysis and Signal Processing》;20121231;1-3 * |
Fast Algorithm for Estimating 2-D DOA in Coherent Signals Case;Jian-Feng Gu 等;《IEEE》;20061231;1-4 * |
一种快速波达方向估计算法;吴建新 等;《西安电子科技大学学报(自然科学版)》;20090430;第36卷(第2期);263-268 * |
基于解相干的改进MUSIC算法DOA估计;周小军 等;《无线电工程》;20141231;第44卷(第12期);18-21 * |
Also Published As
Publication number | Publication date |
---|---|
CN106093845A (zh) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106093845B (zh) | 一种基于伪空间谱搜索的快速doa估计方法 | |
CN109188344B (zh) | 脉冲噪声环境下基于互循环相关music算法信源个数与来波方向角估计方法 | |
CN107064892B (zh) | 基于张量子空间和旋转不变的mimo雷达角度估计算法 | |
CN106610483B (zh) | 基于张量子空间和谱峰搜索的mimo雷达角度估计算法 | |
CN107015213B (zh) | 基于music算法的双基地mimo雷达角度估算方法 | |
CN111337893B (zh) | 一种基于实值稀疏贝叶斯学习的离格doa估计方法 | |
CN103353588B (zh) | 基于天线均匀平面阵的二维波达方向角估计方法 | |
Zhang et al. | Application of artificial bee colony algorithm to maximum likelihood DOA estimation | |
CN109143197B (zh) | 一种基于辅助阵元的极化mimo雷达的2d-doa和极化参数估计方法 | |
Abou Chaaya et al. | Localization of spatially distributed near-field sources with unknown angular spread shape | |
CN105334488B (zh) | 基于信源数估计的栅格偏移优化目标到达角估计方法 | |
CN107340512A (zh) | 一种基于子阵划分的近远场混合源被动定位方法 | |
CN113253192B (zh) | 一种用于非圆信号的互质线阵级联doa估计方法 | |
CN107907853A (zh) | 一种基于均匀圆阵差分相位的单分布源doa估计方法 | |
CN106154217B (zh) | Ula和uca中互耦未知时基于空间谱伪峰消除的自校准方法 | |
CN112180339A (zh) | 一种基于稀疏处理的雷达回波信号精确测向方法 | |
Bai et al. | Association of DOA estimation from two ULAs | |
Reaz et al. | A comprehensive analysis and performance evaluation of different direction of arrival estimation algorithms | |
CN104502885A (zh) | 基于变换矩阵的特征值差分信源数估计方法 | |
WO2020140658A1 (zh) | 波达方向估计方法及装置、雷达、可读存储介质 | |
CN105388450B (zh) | 一种基于目标角度分离的双基地mimo***dod和doa估计方法 | |
Ning et al. | A velocity independent MUSIC algorithm for DOA estimation | |
CN109407047A (zh) | 一种基于秩损求根的幅相误差校准和波达方向估计方法 | |
CN103792509B (zh) | 电磁信号的二维波达方向角估计方法 | |
CN112363106B (zh) | 基于量子粒子群的信号子空间波达方向检测方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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: 20180626 Termination date: 20200527 |