CN111830460A - 基于序贯music的doa估计方法 - Google Patents
基于序贯music的doa估计方法 Download PDFInfo
- Publication number
- CN111830460A CN111830460A CN202010717462.0A CN202010717462A CN111830460A CN 111830460 A CN111830460 A CN 111830460A CN 202010717462 A CN202010717462 A CN 202010717462A CN 111830460 A CN111830460 A CN 111830460A
- Authority
- CN
- China
- Prior art keywords
- array
- sub
- vector
- music
- estimation method
- 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
- 238000000034 method Methods 0.000 title claims abstract description 25
- 239000013598 vector Substances 0.000 claims abstract description 27
- 238000001228 spectrum Methods 0.000 claims abstract description 13
- 239000011159 matrix material Substances 0.000 claims description 30
- 238000003491 array Methods 0.000 claims description 11
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 abstract description 10
- 238000007635 classification algorithm Methods 0.000 abstract description 2
- 238000004422 calculation algorithm Methods 0.000 description 20
- 238000004088 simulation Methods 0.000 description 7
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种基于序贯MUSIC的DOA估计方法,其解决了现有多重信号分类算法对均匀线性阵列做信号参数估计时,计算过程的复杂度高、成本高,不利于工程实现的技术问题,其首先对均匀线性阵列建立信号模型,其次将阵列划分为若干子阵,然后通过对接收和A/D采样通道的时分复用,分别计算每个子阵的协方差矩阵,并对其进行特征值分解,得到噪声子空间;最后利用各子阵导向矢量与噪声子空间的正交性,构造基于子阵的联合谱函数。本发明广泛用于阵列信号测向技术领域。
Description
技术领域
本发明涉及阵列信号测向技术领域,具体而言,涉及一种基于序贯MUSIC的DOA估计方法。
背景技术
波达方向(direction of arrival,DOA)估计广泛应用于雷达、声呐、无线通信、无源定位、导航、地震探测等领域有着广泛应用。
均匀线性阵列是一种较常见的阵列形式,其结构简单,可以估计电磁波的1维DOA参数。在均匀线阵基础上,结合多重信号分类(Multiple Signal Classfication,MUSIC)算法可以实现对信号源的超分辨、高精度测向。因此,其在民用通信、军用侦察等技术领域具有广阔的应用前景。
现有MUSIC算法的处理流程为:(1)估计阵列协方差矩阵;(2)对阵列协方差矩阵进行特征值分解,并得到噪声子空间;(3)利用导向矢量与噪声子空间的正交性,构造谱函数,在一维角度域进行搜索,实现对信号源波达方向的估计。由于测向算法中一般要求阵元与接收通道及A/D采样通道一一对应,因此当阵列中阵元数量较多时,测向系统需要与阵元数相同的接收通道和A/D采样通道,系统复杂度、体积、功耗、成本均较高,不利于工程实现。同时,MUSIC算法包含的协方差矩阵计算、特征值分解、以及谱函数计算等步骤的计算量均与阵元个数有关,较大的阵元数,必然带来较大的计算复杂度,为MUSIC算法的应用带来困难。
发明内容
本发明就是为了解决现有多重信号分类算法对均匀线性阵列做信号参数估计时,计算过程的复杂度高、成本高,不利于工程实现的技术问题,提供了一种降低复杂度、成本,有利于工程实现的基于序贯MUSIC的DOA估计方法。
本发明提供一种基于序贯MUSIC的DOA估计方法,包括以下步骤:
第一步:对均匀线性阵列建立信号模型;
第二步:将阵列划分为若干子阵;
第三步:通过对接收和A/D采样通道的时分复用,分别计算每个子阵的协方差矩阵,并对其进行特征值分解,得到噪声子空间;
第四步:利用各子阵导向矢量与噪声子空间的正交性,构造基于子阵的联合谱函数。
优选地,第一步的过程是:
针对N个信号入射到由M个阵元构成的均匀阵列,建立信号模型为:
X(t)=A(θ)S(t)+N(t) (1)
公式(1)中:
X(t)为数据接收矢量,X(t)=[x1(t),x2(t),…,xM(t)]T;
N(t)为噪声矢量,N(t)=[n1(t),n2(t),…,nM(t)];
S(t)为空间信号源矢量,S(t)=[s1(t),s2(t),…,sN(t)];
A(θ)为阵列流型矩阵,A(θ)=[a(θ1) a(θ2) … a(θN)];
其中导向矢量a(θi)可表示为:
导向矢量a(θi)的公式中,d为阵元间距,θi为第i个信源入射角度,λ为波长。
优选地,第二步的过程是:
选取第1阵元为参考阵元,将其余阵元等间隔划分为K个子阵,每个子阵含L个阵元,记参考阵元及第k个子阵的阵元编号为ek=[1,k1,…kL],则tk时刻,第k个子阵接收信号可表示为:Xk(tk)=X(ek,tk),与其对应的导向矢量可表示为:
优选地,第三步的过程是:
公式(3)中,Λk,S是大特征值组成的对角阵,Λk,N是小特征值组成的对角阵,Uk,S是第k个子阵对应的信号子空间,Uk,N是第k个子阵对应的噪声子空间。
优选地,第四步的过程是:
利用子阵导向矢量与各自噪声子空间正交的性质,构造基于各子阵联合谱的MUSIC表达式为:
通过谱峰搜索,找到极大值点对应的角度就是要估计的DOA信息。
本发明的有益效果是:本发明通过将阵列划分为若干子阵,通过对接收和A/D采样通道的时分复用,大大地降低测向系统的复杂度和成本;同时由于各子阵阵元数减少,其协方差矩阵及特征值分解部分的计算量亦将有效减少。
本发明进一步的特征,将在以下具体实施方式的描述中,得以清楚地记载。
附图说明
图1是均匀线阵结构图,阵元采用均匀线性排布方式,信号入射方向定义如图中所示;
图2是本发明所述方法的接收通道分时复用示意图;
图3是本发明中基于序贯MUSIC的DOA估计谱图;
图4是本发明中序贯MUSIC算法与现有技术的MUSIC算法精度对比。
具体实施方式
以下参照附图,以具体实施例对本发明作进一步详细说明。
本发明通过对均匀线阵进行信号建模,以及构造子阵的方式,实现对各子阵协方差矩阵的估计,在此基础上通过对各子阵协方差矩阵的特征值分解得到与其对应的噪声子空间,最后利用子阵的导向矢量与其噪声子空间正交,实现基于序贯MUSIC算法的DOA估计。
基于序贯MUSIC的DOA估计方法包括以下步骤:
步骤1,针对N个信号入射到由M个阵元构成的均匀线性阵列,建立信号模型为:
X(t)=A(θ)S(t)+N(t) (1)
公式(1)中:
X(t)为数据接收矢量,X(t)=[x1(t),x2(t),…,xM(t)]T;
N(t)为噪声矢量,N(t)=[n1(t),n2(t),…,nM(t)];
S(t)为空间信号源矢量,S(t)=[s1(t),s2(t),…,sN(t)];
A(θ)为阵列流型矩阵,A(θ)=[a(θ1) a(θ2) … a(θN)];
其中导向矢量a(θi)可表示为:
导向矢量a(θi)的公式中,d为阵元间距,θi为第i个信源入射角度,λ为波长。
步骤2,选取第1阵元为参考阵元,将其余阵元等间隔划分为K(K≥2)个子阵,每个子阵含L个阵元(L≥N),记参考阵元及第k个子阵的阵元编号为ek=[1,k1,…kL],则tk时刻,第k个子阵接收信号可表示为:Xk(tk)=X(ek,tk),与其对应的导向矢量可表示为:
公式(3)中,Λk,S是大特征值组成的对角阵,Λk,N是小特征值组成的对角阵,Uk,S是第k个子阵对应的信号子空间,Uk,N是第k个子阵对应的噪声子空间。
步骤4,利用子阵导向矢量与各自噪声子空间正交的性质,构造基于各子阵联合谱的MUSIC表达式为:
通过谱峰搜索,找到极大值点对应的角度就是要估计的DOA信息。
下面介绍相关仿真实验结果:
仿真实验1:考虑由11个阵元构成的均匀线阵,阵元间距d=λ/2,信号源入射方向θ=15°,采样快拍数2028。将整个阵列划分为2个子阵,其中子阵1的阵元序号为:1,3,5,7,9,11;子阵2的阵元序号为:1,2,4,6,8,10;子阵1采样的快拍数为1:1024;子阵2采样的快拍数为1025:2048。取SNR从0dB开始,以2dB为间隔,变化到20dB。在每个SNR进行100次蒙特卡洛仿真,分别统计经典MUSIC算法和序贯MUSIC算法的参数估计精度,仿真结果如图4所示。由仿真结果可知,序贯MUSIC算法在精度略有下降的情况下,有效地减少了接收通道数,降低了系统的复杂度。
仿真实验2:考虑由22个阵元构成的均匀线阵,阵元间距d=λ/2,信号源入射方向θ=15°,采样快拍数300。将整个阵列划分为3个子阵,其中子阵1的阵元序号为:1,4,7,10,13,16,19,22;子阵2的阵元序号为:1,2,5,8,11,14,17,20;子阵3的阵元序号为:1,3,6,9,12,15,18,21;子阵1采样的快拍数为1:100;子阵2采样的快拍数为101:200;子阵3采样的快拍数为201:300。在Intel i5-3470CPU,4GBRAM的计算机上,对协方差矩阵计算,特征值分解,谱函数计算各个步骤的运行时间进行统计,结果如表1所示。仿真结果表明,序贯MUSIC算法的在各个环节的计算时间均小于经典MUSIC算法。
表一:序贯MUSIC算法与现有技术的MUSIC算法计算时间对比(s)
对比项 | 现有技术MUSIC算法 | 序贯MUSIC算法 |
协方差矩阵计算 | 7.2900e<sup>-05</sup> | 1.8999e<sup>-05</sup> |
特征值分解 | 8.7755e<sup>-05</sup> | 4.7320e<sup>-05</sup> |
谱函数计算 | 3.6939e<sup>-05</sup> | 2.1397e<sup>-05</sup> |
综上,本发明公开了一种基于接收通道时分复用的序贯MUSIC方法。该方法首先进行阵列信号建模,然后将阵列划分为K个子阵,分别估计每个子阵的协方差矩阵和噪声子空间;最后利用子阵导向矢量与其对应的噪声子空间正交,构造基于K个子阵的序贯MUSIC表达式,通过对序贯MUSIC谱峰极大值的搜索,实现对DOA参数的估计。仿真结果表明:该算法在参数估计精度略有下降的情况下,极大地减少了系统的通道数,减少了系统的体积、功耗、成本,降低了侧向系统的复杂度和算法的计算复杂度,提高了系统的可实现性。
由于本发明采用接收通道时分复用技术,大大减少了接收通道数量,降低了计算复杂度,为MUSIC算法的工程应用奠定了基础。
以上所述仅对本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。
Claims (5)
1.一种基于序贯MUSIC的DOA估计方法,其特征在于,包括以下步骤:
第一步:对均匀线性阵列建立信号模型;
第二步:将阵列划分为若干子阵;
第三步:通过对接收和A/D采样通道的时分复用,分别计算每个子阵的协方差矩阵,并对其进行特征值分解,得到噪声子空间;
第四步:利用各子阵导向矢量与噪声子空间的正交性,构造基于子阵的联合谱函数。
2.根据权利要求1所述的基于序贯MUSIC的DOA估计方法,其特征在于:
所述第一步的过程是:
针对N个信号入射到由M个阵元构成的均匀阵列,建立信号模型为:
X(t)=A(θ)S(t)+N(t) (1)
公式(1)中:
X(t)为数据接收矢量,X(t)=[x1(t),x2(t),…,xM(t)]T;
N(t)为噪声矢量,N(t)=[n1(t),n2(t),…,nM(t)];
S(t)为空间信号源矢量,S(t)=[s1(t),s2(t),…,sN(t)];
A(θ)为阵列流型矩阵,A(θ)=[a(θ1) a(θ2) … a(θN)];
其中导向矢量a(θi)可表示为:
导向矢量a(θi)的公式中,d为阵元间距,θi为第i个信源入射角度,λ为波长。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010717462.0A CN111830460B (zh) | 2020-07-23 | 2020-07-23 | 基于序贯music的doa估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010717462.0A CN111830460B (zh) | 2020-07-23 | 2020-07-23 | 基于序贯music的doa估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111830460A true CN111830460A (zh) | 2020-10-27 |
CN111830460B CN111830460B (zh) | 2022-11-01 |
Family
ID=72925136
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010717462.0A Expired - Fee Related CN111830460B (zh) | 2020-07-23 | 2020-07-23 | 基于序贯music的doa估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111830460B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1761171A (zh) * | 2004-10-15 | 2006-04-19 | 北京大学 | 基于达波信号方向估计的信号处理系统和方法 |
CN101150345A (zh) * | 2006-09-18 | 2008-03-26 | 中国人民解放军空军装备研究院雷达与电子对抗研究所 | 适用于非平稳噪声背景下相干信源的测向方法 |
CN102175990A (zh) * | 2011-01-27 | 2011-09-07 | 西安交通大学 | 基于龙伯格观测器与子空间更新的波达方向跟踪方法及装置 |
CN102279387A (zh) * | 2011-07-18 | 2011-12-14 | 西安电子科技大学 | Mimo雷达的目标到达角估计方法 |
CN107340512A (zh) * | 2017-06-29 | 2017-11-10 | 电子科技大学 | 一种基于子阵划分的近远场混合源被动定位方法 |
CN108430090A (zh) * | 2017-02-15 | 2018-08-21 | 华为技术有限公司 | 确定和用于确定doa的方法以及接入网设备和终端 |
CN110031794A (zh) * | 2019-04-16 | 2019-07-19 | 中国人民解放军国防科技大学 | 一种基于差分共性阵重构的相干信源doa估计方法 |
US20200103515A1 (en) * | 2018-09-28 | 2020-04-02 | Panasonic Intellectual Property Management Co., Ltd. | Radar apparatus and radar method |
-
2020
- 2020-07-23 CN CN202010717462.0A patent/CN111830460B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1761171A (zh) * | 2004-10-15 | 2006-04-19 | 北京大学 | 基于达波信号方向估计的信号处理系统和方法 |
CN101150345A (zh) * | 2006-09-18 | 2008-03-26 | 中国人民解放军空军装备研究院雷达与电子对抗研究所 | 适用于非平稳噪声背景下相干信源的测向方法 |
CN102175990A (zh) * | 2011-01-27 | 2011-09-07 | 西安交通大学 | 基于龙伯格观测器与子空间更新的波达方向跟踪方法及装置 |
CN102279387A (zh) * | 2011-07-18 | 2011-12-14 | 西安电子科技大学 | Mimo雷达的目标到达角估计方法 |
CN108430090A (zh) * | 2017-02-15 | 2018-08-21 | 华为技术有限公司 | 确定和用于确定doa的方法以及接入网设备和终端 |
CN107340512A (zh) * | 2017-06-29 | 2017-11-10 | 电子科技大学 | 一种基于子阵划分的近远场混合源被动定位方法 |
US20200103515A1 (en) * | 2018-09-28 | 2020-04-02 | Panasonic Intellectual Property Management Co., Ltd. | Radar apparatus and radar method |
CN110031794A (zh) * | 2019-04-16 | 2019-07-19 | 中国人民解放军国防科技大学 | 一种基于差分共性阵重构的相干信源doa估计方法 |
Non-Patent Citations (4)
Title |
---|
XIANG-TIAN MENG ET AL.: "Real-Valued DOA Estimation for Non-Circular Sources via Reduced-Order Polynomial Rooting", 《IEEE ACCESS》 * |
刘帅 等: "基于四元数MUSIC的锥面共形阵列极化-DOA联合估计", 《系统工程与电子技术》 * |
张琛: "基于极化敏感阵列的多维参数估计算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
闫锋刚 等: "高效超分辨波达方向估计算法综述", 《系统工程电子与技术》 * |
Also Published As
Publication number | Publication date |
---|---|
CN111830460B (zh) | 2022-11-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103901417B (zh) | L型阵列mimo雷达低复杂度空间目标二维角度估计方法 | |
CN108896954B (zh) | 互质阵中一种基于联合实值子空间的波达角估计方法 | |
CN109239657A (zh) | 装载嵌套阵无人机平台下的辐射源高精度定位方法 | |
CN109490820A (zh) | 一种基于平行嵌套阵的二维doa估计方法 | |
CN110244273A (zh) | 一种基于均匀分布式阵列的目标角度估计方法 | |
CN113671439B (zh) | 基于非均匀智能超表面阵列的无人机集群测向系统及方法 | |
CN112462363B (zh) | 非均匀稀疏极化阵列相干目标参数估计方法 | |
CN108120953A (zh) | 一种基于波达方向估计的无线电定位方法 | |
CN111273269B (zh) | 基于ipso-bp的频率分集阵列的雷达目标定位方法 | |
CN111965591A (zh) | 一种基于四阶累积量矢量化dft的测向估计方法 | |
CN108872930A (zh) | 扩展孔径二维联合对角化doa估计方法 | |
CN113759303B (zh) | 一种基于粒子群算法的无网格波达角估计方法 | |
CN109946663B (zh) | 一种线性复杂度的Massive MIMO目标空间方位估计方法和装置 | |
CN105572629B (zh) | 一种适用于任意阵列结构的低运算复杂度的二维测向方法 | |
CN111368256B (zh) | 一种基于均匀圆阵的单快拍测向方法 | |
CN108594165B (zh) | 一种基于期望最大化算法的窄带信号波达方向估计方法 | |
CN105068090B (zh) | 一种利用单快拍数据进行gnss天线阵干扰抑制的方法 | |
CN110208736B (zh) | 基于四阶累量的非圆信号均匀阵列波达方向角估计方法 | |
CN111830460B (zh) | 基于序贯music的doa估计方法 | |
CN116699511A (zh) | 一种多频点信号波达方向估计方法、系统、设备及介质 | |
CN113341371B (zh) | 一种基于l阵和二维esprit算法的doa估计方法 | |
CN111366891B (zh) | 一种基于伪协方差矩阵的均匀圆阵单快拍测向方法 | |
CN114609580A (zh) | 一种基于非圆信号的无孔互质阵列设计方法 | |
CN114265004A (zh) | 一种基于子空间对消的干扰下的目标角度估计方法 | |
CN114879136B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20221101 |
|
CF01 | Termination of patent right due to non-payment of annual fee |