CN109633522B - 基于改进的music算法的波达方向估计方法 - Google Patents

基于改进的music算法的波达方向估计方法 Download PDF

Info

Publication number
CN109633522B
CN109633522B CN201811599589.6A CN201811599589A CN109633522B CN 109633522 B CN109633522 B CN 109633522B CN 201811599589 A CN201811599589 A CN 201811599589A CN 109633522 B CN109633522 B CN 109633522B
Authority
CN
China
Prior art keywords
matrix
array
complex
real
antenna 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
Application number
CN201811599589.6A
Other languages
English (en)
Other versions
CN109633522A (zh
Inventor
李晓婷
胡山峰
王君方
张宝强
张凡
李�瑞
赵微
张忠波
刘洋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shaanxi Fenghuo Communication Group Co Ltd
Original Assignee
Shaanxi Fenghuo Communication Group Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shaanxi Fenghuo Communication Group Co Ltd filed Critical Shaanxi Fenghuo Communication Group Co Ltd
Priority to CN201811599589.6A priority Critical patent/CN109633522B/zh
Publication of CN109633522A publication Critical patent/CN109633522A/zh
Application granted granted Critical
Publication of CN109633522B publication Critical patent/CN109633522B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Electromagnetism (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本发明属于通信技术领域，公开了基于改进的MUSIC算法的波达方向估计方法。该方法包括：建立对称阵列模型，对各天线阵元接收的射频信号进行采样，正交下变频和滤波得到复信号，对复信号构成的阵列输出矩阵做线性变换得到复矩阵，进而处理得到实对称矩阵，对实对称矩阵进行特征值分解，构建噪声子空间矩阵，构造阵列方向矩阵，利用噪声子空间和阵列方向矩阵构造空间谱函数，进而进行谱峰局部极值搜索，得到函数极大值对应的角度值即为波达方向。本发明能够在不损失算法性能的情况下将特征值分解由复矩阵运算转换为实矩阵运算，并利用实方向矩阵进行峰值搜索，从而能有效地估计出信号的波达方向，提高了波达方向估计的准确性并且降低了计算的复杂度。

Description

基于改进的MUSIC算法的波达方向估计方法
技术领域
本发明涉及通信技术领域，尤其涉及基于改进的MUSIC算法的波达方向估计方法。
背景技术
在无线电通信、导航、雷达、声纳、电子侦察等领域，获得信号的来波方向极其重要。在军事无线通信领域，获得高分辨率的来波方向不仅有助于对我方目标的导航、跟踪，还可以有助于我方对敌方目标进行准确定位以及打击。在民用无线通信领域，准确的来波方向估计能够提升系统容量和减少多址干扰等。
空间谱表示信号在各个方向上的能量分布，若能够估计出空间谱就可以得到信号源的波达方向，所以空间谱估计也常称为波达方向(英文全称：Direction Of Arrival，英文缩写：DOA)估计，DOA估计技术的关键在于利用空间不同位置的天线阵列，接收来自不同方向的多个信号源发出的信号，运用现代信号处理技术估计出信号源的来波方向。在多种空间谱估计算法中，多重信号分类(英文全称：Multiple Signal Classification，英文缩写：MUSIC)算法是一种基于子空间分解的算法，其特点为测向误差小、精度高、分辨率高、灵敏度高，可以对入射信号数目、到达方向及波形的强度做出无偏估计。因此研究MUSIC算法对于促进更广泛的应用有着十分重要的意义。
但是，传统的MUSIC算法在协方差矩阵计算及特征分解时涉及大量的复数运算，导致算法的复杂度高，不利于硬件实现。因此，如何在降低计算复杂度的同时又不损失算法性能就显得尤为重要。
发明内容
本发明的实施例提供基于改进的MUSIC算法的波达方向估计方法，能够在不损失算法性能的情况下将特征值分解由复矩阵运算转换为实矩阵运算，并利用实方向矩阵进行峰值搜索，从而能有效地估计出信号的波达方向。此发明提高了波达方向估计的准确性并且降低了计算的复杂度。
为达到上述目的，本发明的实施例采用如下技术方案：
步骤1、建立对称阵列模型，天线阵列为均匀线阵或者平面均匀圆阵，阵元个数为偶数；当天线阵元为均匀线阵时，以阵列的对称轴为基准，向右将阵列对称轴右侧天线阵元依次编号为
Figure BDA0001922122910000021
向左将阵列对称轴左侧天线阵元依次编号为
Figure BDA0001922122910000022
当天线阵元为平面均匀圆阵时，以其中任一阵元为起点，按顺时针方向依次编号为1，2...M，阵元间隔为d；
步骤2、对各天线阵元接收的射频模拟信号进行K次采样，其中，第m个阵元天线的第k次采样得到的采样信号为sm(k)，m＝1，2，…M，k＝1，2，…K；K为快拍数，K＞2fs/Δf，fs为采样频率，Δf为需要分辨的最小频率，M为天线阵元个数，M为偶数；
步骤3、对各天线阵元的K个采样信号做I、Q正交下变频和滤波，得到各天线阵元的K个复信号，并对第一阵元的K个复信号做快速傅里叶变换，判断快速傅里叶变换的最大值是否超过预设的门限值，如果是，则执行步骤4；
步骤4、利用各天线阵元的K个复信号构成阵列输出矩阵X，对阵列输出矩阵X做线性变换得到复矩阵Y，利用复矩阵Y进一步计算得到复共轭对称矩阵R，利用复共轭对称矩阵R进一步计算得到实对称矩阵RX；
步骤5、对实对称矩阵RX进行特征值分解，得到实对称矩阵RX的特征值和对应的特征向量，并将特征值按照从大到小的顺序进行排序；
步骤6、基于最小描述长度MDL准则，利用排序后的实对称矩阵RX的特征值估计信号源个数，进而根据信号源个数，利用实对称矩阵RX的特征向量得到噪声子空间矩阵；
步骤7、基于对称阵列模型构造阵列方向矩阵，利用噪声子空间和阵列方向矩阵构造空间谱函数，进而利用空间谱函数进行谱峰局部极值搜索，得到空间谱函数极大值对应的角度值，角度值即为信号源的入射方向。
本发明的实施例提供基于改进的MUSIC算法的波达方向估计方法，首先建立天线的对称阵列模型，通过对中心对称阵列结构的数学分析，利用采样数据及其共轭信息得到了一种线性变换方法，能够在不损失算法性能的情况下将特征值分解由复矩阵运算转换为实矩阵运算，并利用实方向矩阵进行峰值搜索，从而能有效地估计出信号的波达方向。此发明提高了波达方向估计的准确性并且降低了计算的复杂度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的基于改进的MUSIC算法的波达方向估计方法的流程示意图；
图2为四阵元均匀直线阵的阵列结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1为本发明实施例提供的基于改进的MUSIC算法的波达方向估计方法的流程示意图，如图1所示，
步骤1、建立对称阵列模型，天线阵列为均匀线阵或者平面均匀圆阵，阵元个数为偶数；当天线阵元为均匀线阵时，以阵列的对称轴为基准，向右将阵列对称轴右侧天线阵元依次编号为
Figure BDA0001922122910000041
向左将阵列对称轴左侧天线阵元依次编号为
Figure BDA0001922122910000042
当天线阵元为平面均匀圆阵时，以其中任一阵元为起点，按顺时针方向依次编号为1，2...M，阵元间隔为d；
步骤2、对各天线阵元接收的射频模拟信号进行K次采样，其中，第m个阵元天线的第k次采样得到的采样信号为sm(k)，m＝1，2，…M，k＝1，2，…K；K为快拍数，K＞2fs/Δf，fs为采样频率，Δf为需要分辨的最小频率，M为天线阵元个数，M为偶数；
步骤3、对各天线阵元的K个采样信号做I、Q正交下变频和滤波，得到各天线阵元的K个复信号，并对第一阵元的K个复信号做快速傅里叶变换，判断快速傅里叶变换的最大值是否超过预设的门限值，如果是，则执行步骤4；
进一步的，利用公式
Figure BDA0001922122910000043
对第一阵元的K个复信号做快速傅里叶变换；
其中，第m个天线阵元第k次采样对应的复信号为xm(k)＝Im(k)+jQm(k)，m＝1，2，…M，k＝1，2，…K；N为快速傅里叶变换的点数，
Figure BDA0001922122910000051
步骤4、利用各天线阵元的K个复信号构成阵列输出矩阵X，对阵列输出矩阵X做线性变换得到复矩阵Y，利用复矩阵Y进一步计算得到复共轭对称矩阵R，利用复共轭对称矩阵R进一步计算得到实对称矩阵RX；
进一步的，步骤4具体包括：
(4.1)利用各天线阵元的K个复信号构成阵列输出矩阵X，X＝[X1(k)，X2(k)，...，Xm(k)，...，XM(k)]T，Xm(k)＝[xm(1)，xm(2)，...xm(k)...，xm(K)]，m＝1，2，…M；
(4.2)构造线性变换矩阵
Figure BDA0001922122910000052
根据线性变换矩阵T和阵列输出矩阵X，计算得到复矩阵Y＝TX；
其中，
Figure BDA0001922122910000053
M为表示阵元个数，M为偶数；
(4.3)利用复矩阵Y，计算得到复共轭对称矩阵R＝Y×YH，其中，上标H表示矩阵的共轭转置；
(4.4)取复共轭对称矩阵R的实部，得到实对称矩阵RX＝Re{R}，其中，Re{·}表示取矩阵的实部。
步骤5、对实对称矩阵RX进行特征值分解，得到实对称矩阵RX的特征值和对应的特征向量，并将特征值按照从大到小的顺序进行排序；
进一步的，步骤5具体包括：
(5.1)利用预设公式：RX＝UAUH，对实对称矩阵RX进行特征值分解；
其中，Λ为对角矩阵，其对角线上的元素为实对称矩阵RX的特征值，U是RX的特征矢量矩阵；上标H表示矩阵的共轭转置；
(5.2)将RX的特征值按大小顺序排列，得到λ1≥λ2≥…≥λM。
步骤6、基于最小描述长度MDL准则，利用排序后的实对称矩阵RX的特征值估计信号源个数，进而根据信号源个数，利用实对称矩阵RX的特征向量得到噪声子空间矩阵。
进一步的，步骤6具体包括：
(6.1)基于最小描述长度MDL准则，利用特征值构造信号源数目估计模型：
Figure BDA0001922122910000061
其中，D的取值是0，1，…，M-1中的一个。
(6.2)利用信号源数目估计模型，计算得到空间信号源数目的估计值
Figure BDA0001922122910000062
Figure BDA0001922122910000063
其中，arg{min MDL(D)}表示MDL(D)取最小值时，D的取值；
(6.3)利用实对称矩阵RX的
Figure BDA0001922122910000064
个较小特征值
Figure BDA0001922122910000065
对应的特征向量构成噪声子空间矩阵
Figure BDA0001922122910000066
噪声子空间矩阵Un的维度为
Figure BDA0001922122910000067
步骤7、基于对称阵列模型构造阵列方向矩阵，利用噪声子空间和阵列方向矩阵构造空间谱函数，进而利用空间谱函数进行谱峰局部极值搜索，得到空间谱函数极大值对应的角度值，角度值即为信号源的入射方向。
优选的，步骤7具体包括：
(7.1)基于对称阵列模型构造阵列方向矩阵
Figure BDA0001922122910000071
第n个列矢量为
Figure BDA0001922122910000072
，进一步计算得到阵列的实方向矩阵
Figure BDA0001922122910000073
(7.2)利用实方向矩阵的各个列矢量
Figure BDA0001922122910000074
和噪声子空间矩阵Un构造空间谱函数：
Figure BDA0001922122910000075
其中，θ为信号的入射角度；
(7.3)计算PMUSIC(θ)的
Figure BDA0001922122910000076
个极值点，进一步得到
Figure BDA0001922122910000077
个极值点对应的角度值
Figure BDA0001922122910000078
角度值即为信号源的入射方向。
参照图2，以四阵元均匀直线阵(M＝4)为例，对线性变换的计算过程进行原理性说明：
为了符合阵列在结构上的对称性，将阵元1和阵元2组成子阵1，阵元3和阵元4组成子阵2，信号入射角度用θ表示，则子阵1的方向矩阵A1和子阵2的方向矩阵A2可以分别表示为
Figure BDA0001922122910000079
Figure BDA00019221229100000710
Figure BDA0001922122910000081
Figure BDA0001922122910000082
其中，阵元间距
Figure BDA0001922122910000083
λ表示入射信号波长。
进一步得到子阵1的输出矩阵X1和子阵2的输出矩阵X2：
Figure BDA0001922122910000084
Figure BDA0001922122910000085
其中，S(k)为时域连续K个采样快拍的入射信号构成的矩阵，N1(k)为子阵1的噪声矩阵，N2(k)和子阵2噪声矩阵。
由于两个子阵的方向矩阵A1和A2是互为共轭对称的。于是构造一个线性变换矩阵T：
Figure BDA0001922122910000086
Figure BDA0001922122910000087
对X(k)作线性变换：
Figure BDA0001922122910000088
则有，
Figure BDA0001922122910000089
为一个实矩阵，其中：
Figure BDA00019221229100000810
N(k)为阵列的噪声矩阵。
本发明的实施例提供基于改进的MUSIC算法的波达方向估计方法，首先建立天线的对称阵列模型，通过对中心对称阵列结构的数学分析，利用采样数据及其共轭信息得到了一种线性变换方法，能够在不损失算法性能的情况下将特征值分解由复矩阵运算转换为实矩阵运算，并利用实方向矩阵进行峰值搜索，从而能有效地估计出信号的波达方向。此发明提高了波达方向估计的准确性并且降低了计算的复杂度。
本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成，前述的程序可以存储于一计算机可读取存储介质中，该程序在执行时，执行包括上述方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
以上，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。

Claims (2)

1.一种基于改进的MUSIC算法的波达方向估计方法，其特征在于，包括以下步骤：
步骤1、建立对称阵列模型，天线阵列为均匀线阵或者平面均匀圆阵，阵元个数为偶数；当天线阵元为均匀线阵时，以阵列的对称轴为基准，向右将阵列对称轴右侧天线阵元依次编号为
Figure FDA0003972638610000011
向左将阵列对称轴左侧天线阵元依次编号为
Figure FDA0003972638610000012
当天线阵元为平面均匀圆阵时，以其中任一阵元为起点，按顺时针方向依次编号为1，2...M，阵元间隔为d；
步骤2、对各天线阵元接收的射频模拟信号进行K次采样，其中，第m个阵元天线的第k次采样得到的采样信号为sm(k)，m＝1，2，…M，k＝1，2，…K；K为快拍数，K＞2fs/Δf，fs为采样频率，Δf为需要分辨的最小频率，M为天线阵元个数，M为偶数；
步骤3、对所述各天线阵元的K个采样信号做I、Q正交下变频和滤波，得到所述各天线阵元的K个复信号，并对第一阵元的K个复信号做快速傅里叶变换，判断所述快速傅里叶变换的最大值是否超过预设的门限值，如果是，则执行步骤4；
在步骤3中，利用公式
Figure FDA0003972638610000013
对第一阵元的K个复信号做快速傅里叶变换；
其中，第m个天线阵元第k次采样对应的复信号为xm(k)＝Im(k)+jQm(k)，m＝1，2，…M，k＝1，2，…K；N为快速傅里叶变换的点数，
Figure FDA0003972638610000014
步骤4、利用所述各天线阵元的K个复信号构成阵列输出矩阵X，对阵列输出矩阵X做线性变换得到复矩阵Y，利用复矩阵Y进一步计算得到复共轭对称矩阵R，利用复共轭对称矩阵R进一步计算得到实对称矩阵RX；
所述步骤4包括：
(4.1)利用所述各天线阵元的K个复信号构成阵列输出矩阵X，X＝[X1(k)，X2(k)，...，Xm(k)，...，XM(k)]T，Xm(k)＝[xm(1)，xm(2)，...xm(k)...，xm(K)]，m＝1，2，…M；
(4.2)构造线性变换矩阵
Figure FDA0003972638610000021
根据线性变换矩阵T和阵列输出矩阵X，计算得到复矩阵Y＝TX；
其中，
Figure FDA0003972638610000022
M为表示阵元个数，M为偶数；
(4.3)利用复矩阵Y，计算得到复共轭对称矩阵R-Y×YH，其中，上标H表示矩阵的共轭转置；
(4.4)取所述复共轭对称矩阵R的实部，得到实对称矩阵RX＝Re{R}，其中，Re{·}表示取矩阵的实部；
步骤5、对实对称矩阵RX进行特征值分解，得到实对称矩阵RX的特征值和对应的特征向量，并将所述特征值按照从大到小的顺序进行排序；
步骤6、基于最小描述长度MDL准则，利用排序后的实对称矩阵RX的特征值估计信号源个数，进而根据信号源个数，利用实对称矩阵RX的特征向量得到噪声子空间矩阵；
所述步骤6具体为：
(6.1)基于最小描述长度MDL准则，利用特征值构造信号源数目估计模型：
Figure FDA0003972638610000031
其中，D的取值是0，1，…，M-1中的一个；
(6.2)利用信号源数目估计模型，计算得到空间信号源数目的估计值
Figure FDA0003972638610000032
Figure FDA0003972638610000033
其中，arg{minMDL(D)}表示MDL(D)取最小值时，D的取值；
(6.3)利用实对称矩阵RX的
Figure FDA0003972638610000034
个较小特征值
Figure FDA0003972638610000035
对应的特征向量构成噪声子空间矩阵
Figure FDA0003972638610000036
噪声子空间矩阵Un的维度为
Figure FDA0003972638610000037
步骤7、基于对称阵列模型构造阵列方向矩阵，利用噪声子空间和阵列方向矩阵构造空间谱函数，进而利用空间谱函数进行谱峰局部极值搜索，得到所述空间谱函数极大值对应的角度值，所述角度值即为所述信号源的入射方向；
所述步骤7具体包括：
(7.1)基于对称阵列模型构造阵列方向矩阵
Figure FDA0003972638610000038
第n个列矢量为
Figure FDA0003972638610000039
，进一步计算得到阵列的实方向矩阵
Figure FDA00039726386100000310
(7.2)利用实方向矩阵的各个列矢量
Figure FDA0003972638610000041
和噪声子空间矩阵Un构造空间谱函数：
Figure FDA0003972638610000042
其中，θ为信号的入射角度；
(7.3)计算PMUSIC(θ)的
Figure FDA0003972638610000043
个极值点，进一步得到
Figure FDA0003972638610000044
个极值点对应的角度值
Figure FDA0003972638610000045
所述角度值即为所述信号源的入射方向。
2.根据权利要求1所述的方法，其特征在于，所述步骤5中具体包括：
(5.1)利用预设公式：RX＝UΛUH，对实对称矩阵RX进行特征值分解；
其中，人为对角矩阵，其对角线上的元素为实对称矩阵RX的特征值，U是RX的特征矢量矩阵；上标H表示矩阵的共轭转置；
(5.2)将RX的特征值按大小顺序排列，得到λ1≥λ2≥...≥λM。
CN201811599589.6A 2018-12-26 2018-12-26 基于改进的music算法的波达方向估计方法 Active CN109633522B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811599589.6A CN109633522B (zh) 2018-12-26 2018-12-26 基于改进的music算法的波达方向估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811599589.6A CN109633522B (zh) 2018-12-26 2018-12-26 基于改进的music算法的波达方向估计方法

Publications (2)

Publication Number Publication Date
CN109633522A CN109633522A (zh) 2019-04-16
CN109633522B true CN109633522B (zh) 2023-03-17

Family

ID=66077749

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811599589.6A Active CN109633522B (zh) 2018-12-26 2018-12-26 基于改进的music算法的波达方向估计方法

Country Status (1)

Country Link
CN (1) CN109633522B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110320492B (zh) * 2019-07-16 2023-01-03 哈尔滨工程大学 一种基于导向矢量对称特性的低复杂度doa估计方法
CN112666558B (zh) * 2019-10-16 2024-05-14 深圳开阳电子股份有限公司 一种适用于汽车fmcw雷达的低复杂度music测向方法及装置
CN111175693A (zh) * 2020-01-19 2020-05-19 河北科技大学 一种波达方向估计方法及波达方向估计装置
CN113138367A (zh) * 2020-01-20 2021-07-20 中国科学院上海微系统与信息技术研究所 一种目标定位方法、装置、电子设备及存储介质
CN111551924B (zh) * 2020-06-10 2022-11-04 重庆圭研科技有限公司 一种数字信号处理方法
CN111965593B (zh) * 2020-07-10 2024-03-08 北京邮电大学 一种信号的波达方向确定方法及装置
CN112526446B (zh) * 2020-11-27 2022-08-09 石家庄铁道大学 用于评价电磁干扰源doa估计算法分辨率的方法
CN112666513B (zh) * 2020-12-11 2024-05-07 中国人民解放军63892部队 一种改进式music波达方向估计方法
CN113009410A (zh) * 2021-02-18 2021-06-22 西北工业大学 一种浅海多径环境下目标doa估计联合处理方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002048853A (ja) * 2000-08-02 2002-02-15 Matsushita Electric Ind Co Ltd 電波到来方向推定装置及び指向性可変送受信装置
CN1523372A (zh) * 2003-02-21 2004-08-25 重庆邮电学院 一种基于td-scdma无线定位来波方向的估计方法
JP2014119343A (ja) * 2012-12-17 2014-06-30 Toshiba Corp 到来方向推定装置及び方法
CN105548953A (zh) * 2016-01-11 2016-05-04 电子科技大学 利用多个天线同步发射不同信号的单天线接收机测向方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002048853A (ja) * 2000-08-02 2002-02-15 Matsushita Electric Ind Co Ltd 電波到来方向推定装置及び指向性可変送受信装置
CN1523372A (zh) * 2003-02-21 2004-08-25 重庆邮电学院 一种基于td-scdma无线定位来波方向的估计方法
JP2014119343A (ja) * 2012-12-17 2014-06-30 Toshiba Corp 到来方向推定装置及び方法
CN105548953A (zh) * 2016-01-11 2016-05-04 电子科技大学 利用多个天线同步发射不同信号的单天线接收机测向方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
DOA measurements on indoor channel based on music and MDL processing;Kurdi Y 等;《2008 IEEE Antennas and Propagation Society International Symposium and USNC/URSI National Radio Science Meeting》;20080909;全文 *
Millimeter-Wave Channel Estimation Based on 2-D Beamspace MUSIC Method;Ziyu Guo 等;《IEEE Transactions on Wireless Communications》;20170602;全文 *
Robust DOA estimation for coherent signals based on modified MUSIC algorithm;Yongmei Sun 等;《2012 Third International Conference on Intelligent Control and Information Processing》;20121224;全文 *
The PARAFAC-MUSIC Algorithm for DOA Estimation with Doppler Frequency in a MIMO Radar System;Nan Wang 等;《International Journal of Antennas and Propagation》;20140623;全文 *
基于最小模MUSIC算法的电力谐波/间谐波检测方法;卢洪瑞 等;《江苏科技大学学报》;20180615;第32卷(第03期);全文 *
矢量线阵二维波达方位估计的方法;赵微 等;《声学技术》;20081015;全文 *

Also Published As

Publication number Publication date
CN109633522A (zh) 2019-04-16

Similar Documents

Publication Publication Date Title
CN109633522B (zh) 基于改进的music算法的波达方向估计方法
CN105912791B (zh) 虚拟互质阵列中基于局部搜索的doa估计方法
CN109061554B (zh) 一种基于空间离散网格动态更新的目标到达角度估计方法
CN109375154B (zh) 一种冲击噪声环境下基于均匀圆阵的相干信号参数估计方法
CN108120967B (zh) 一种平面阵列doa估计方法及设备
CN111123192B (zh) 一种基于圆形阵列和虚拟扩展的二维doa定位方法
CN109375152B (zh) 电磁矢量嵌套l阵下低复杂度的doa与极化联合估计方法
CN112051540B (zh) 一种快速高精度测向方法
CN111965591B (zh) 一种基于四阶累积量矢量化dft的测向估计方法
CN112130111A (zh) 一种大规模均匀十字阵列中单快拍二维doa估计方法
CN108363048B (zh) 一种基于块稀疏的极化mimo雷达的角度估计方法
Qi et al. Time-frequency DOA estimation of chirp signals based on multi-subarray
Yilmazer et al. DOA estimation using matrix pencil and ESPRIT methods using single and multiple snapshots
CN111693947A (zh) 基于互质阵列doa估计的改进music方法
CN113296050B (zh) 基于各向异性阵列的极化和角度参数联合估计方法
Gupta et al. Angle of arrival detection by ESPRIT method
CN113032721B (zh) 一种低计算复杂度的远场和近场混合信号源参数估计方法
CN113253192B (zh) 一种用于非圆信号的互质线阵级联doa估计方法
CN112763972B (zh) 基于稀疏表示的双平行线阵二维doa估计方法及计算设备
CN112363108B (zh) 信号子空间加权超分辨的波达方向检测方法及系统
CN113341371B (zh) 一种基于l阵和二维esprit算法的doa估计方法
Biswas et al. New high resolution direction of arrival estimation using Compressive Sensing
CN115248413A (zh) 一种适用于非均匀线阵的离格信号波达方向估计方法
Ghobadzadeh et al. Low-complexity 2D root-MUSIC pairing for an L-shaped array
yun Kim et al. A cascade AOA estimation technique with rectangular array antenna

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