CN113219398B - 远场窄带无线电信号波达方向估计方法 - Google Patents

远场窄带无线电信号波达方向估计方法 Download PDF

Info

Publication number
CN113219398B
CN113219398B CN202010573169.1A CN202010573169A CN113219398B CN 113219398 B CN113219398 B CN 113219398B CN 202010573169 A CN202010573169 A CN 202010573169A CN 113219398 B CN113219398 B CN 113219398B
Authority
CN
China
Prior art keywords
root
music
polynomial
gcd
algorithm
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
CN202010573169.1A
Other languages
English (en)
Other versions
CN113219398A (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.)
Harbin Institute of Technology Weihai
Original Assignee
Harbin Institute of Technology Weihai
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 Harbin Institute of Technology Weihai filed Critical Harbin Institute of Technology Weihai
Priority to CN202010573169.1A priority Critical patent/CN113219398B/zh
Publication of CN113219398A publication Critical patent/CN113219398A/zh
Application granted granted Critical
Publication of CN113219398B publication Critical patent/CN113219398B/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
    • 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

Abstract

本发明属于阵列信号处理技术领域，具体的说是一种能够显著降低运算复杂度的基于最大公因式的广义多项式降阶求根的远场窄带无线电信号波达方向估计方法，其特征在于,接收辐射源信号后，获取辐射源信号的求根多项式，然后通过根与系数的关系，构造近似求根多项式及其导数；将所述近似求根多项式及其导数组成病态方程组后，提取方程组的最大公因式，在得到含有真实DOA信息的最大公因式后，直接通过求根运算得到信号的波达方向。

Description

远场窄带无线电信号波达方向估计方法
技术领域：
本发明属于阵列信号处理技术领域，具体的说是一种能够显著降低运算复杂度的基于最大公因式的广义多项式降阶求根的远场窄带无线电信号波达方向估计方法。
背景技术：
为了在复杂电磁环境下对空间目标进行识别、定位、跟踪，波达方向(directionof arrival，DOA)估计需要利用空间平面内的天线阵列对不同信号的来波方向进行准确地估计，在雷达、声呐、导航等领域中具有重要的研究价值。以多重信号分类(MUSIC)为典型代表的搜索类算法的提出，使得传统测向理论正式迈入超分辨的时代。但在超分辨测向由理论研究向装备研发的工程化转型过程中，具有庞大计算复杂度的MUSIC算法阻碍了工程化的推进。为此，求根类波达方向估计技术的诞生则开启了新的篇章。
root-MUSIC算法是使用最为广泛的求根类波达方向估计技术。假设天线单元个数为M，为了避免MUSIC的谱峰搜索，root-MUSIC 算法通过对一个包含目标信号DOA信息的2(M-1)阶求根多项式进行求根运算来获取DOA，其计算复杂度与(2(M-1))3成正比。现如今，为了追求超分辨算法具有良好角度分辨力的优势，接收端通常采用大型天线阵列，比如：相控阵雷达、超视距雷达和MIMO通信系统等。此时，更高的2(M-1)阶求根多项式将会给测向系统带来计算冗余，同时降低了测向系统对接收数据的实时处理能力。
由于root-MUSIC算法在特征值分解(Eigenvalue Decomposition, EVD)和求根运算两个过程都进行复值运算，酉root-MUSIC算法 (U-root-MUSIC)利用前后向平滑和酉变换等数学方法，在EVD阶段实现了实值计算。受U-root-MUSIC算法的启发，近年来的全实值root-MUSIC算法(RV-root-MUSIC)通过对协方差矩阵的实部进行实值EVD，实现了EVD和求根运算的双重实值计算过程，极大地减少了计算复杂度。
但是，对于众多实值求根类波达方向估计技术而言，例如， U-root-MUSIC算法和RV-root-MUSIC算法等，虽然实值计算可以降低部分计算复杂度，但是求根多项式的阶数依旧为2(M-1)。在面对大型天线阵列的时候，高阶多项式所带来的计算负担显然已经成为了制约测向理论迈向装备研发的掣肘。为此，假设目标信号个数为L(通常L＜＜M)，若能构建一个阶数仅与L相关且仅包含目标信号DOA信息的降阶求根多项式，那么将大幅度地降低求根运算过程中的计算复杂度。
发明内容：
本发明针对现有技术中复求根多项式阶数过高的问题，提出一种基于最大公因式的广义多项式降阶求根波达方向估计新方法，通过求根多项式中系数与根的内在关系，构造了近似的求根多项式；另一方面，将近似多项式与其导数联立并组成病态方程组，通过提取方程组中的最大公因式，实现了高效的波达方向估计；同时，最大公因式的阶数仅与目标信号的个数有关，且最大公因式中仅包含目标信号的波达方向，理论分析和试验结果表明，该方法相对于求根类波达方向估计技术，显著地降低了求根过程中的计算复杂度。
本发明通过以下措施达到：
一种远场窄带无线电信号波达方向估计方法，其特征在于,接收辐射源信号后，获取辐射源信号的求根多项式，然后通过根与系数的关系，构造近似求根多项式及其导数；将所述近似求根多项式及其导数组成病态方程组后，提取方程组的最大公因式，在得到含有真实 DOA信息的最大公因式后，直接通过求根运算得到信号的波达方向。
本发明所述辐射源信号可以通过天线阵列接收获得，假设天线阵列中设有M个相互独立的阵元，以d等间距组成均匀线阵(ULA)，考虑空间中存在L个远场窄带信号从方向θ＝[θ1,…,θL]入射到阵列，其中，假设L先验已知，阵元间距d满足d≤λ/2以避免相位模糊，λ为窄带信号的波长，则天线阵列接收辐射源信号为：
Figure BDA0002550386060000031
其中，θk为空间中第k∈[1,L]个信号的来波方向，A(θ)为M×L维的阵列流型矩阵，s(t)为L×1维的入射信号矢量，n(t)为M×1维的加性高斯白噪声矢量，a(θ)为A(θ)的列向量，表示为：
a(z)＝[1,z,z2,…,zM-1]T,
其中，(·)T为转置运算，z＝jφ，φ＝(2π/λ)dsinθ；
M×M维阵列协方差矩阵为：
Figure BDA0002550386060000032
其中，(·)H为共轭转置运算，Rss＝E[s(t)sH(t)]为L×L的信号协方差矩阵，
Figure BDA0002550386060000033
为噪声功率，
Figure BDA0002550386060000036
为单位矩阵。理论上Rxx不可知，但可用N个快拍数据对其估计：
Figure BDA0002550386060000034
其复值特征值分解可以表示为：
Figure BDA0002550386060000035
本发明所述获取辐射源信号的求根多项式是指通过现有方法，如 root-MUSIC算法、U-root-MUSIC算法或RV-root-MUSIC算法，获得求根类波达方向估计技术的求根多项式,具体分别如下：
root-MUSIC算法的求根多项式为：
Figure BDA0002550386060000041
其中，(·)-1为矩阵的逆运算，ξk为多项式froot-MUSIC(z)的系数。 U-root-MUSIC算法利用前后向平滑技术构造了半实值求根多项式，可表示为：
Figure BDA0002550386060000042
其中，ζk为多项式fU-root-MUSIC(z)的系数，
Figure BDA0002550386060000043
U是由单位矩阵
Figure BDA0002550386060000044
和反对角单位矩阵
Figure BDA0002550386060000045
构成的酉矩阵，且定义为：
Figure BDA0002550386060000046
RV-root-MUSIC算法通过对
Figure BDA0002550386060000047
的实部进行特征值分解，构造了实值求根多项式，可表示为：
Figure BDA0002550386060000048
其中，ηk是多项式fRV-root-MUSIC(z)的系数；假设，span(·)为空间张成运算，(·)为共轭运算，那么
Figure BDA0002550386060000049
本发明所述通过根与系数的关系，构造近似求根多项式及其导数,
其中以root-MUSIC算法为例，利用根与系数的内在关系将 froot-MUSIC(z)因式分解为：
Figure BDA0002550386060000051
其中，zi为含有真实DOA的根，zj为余下不含有DOA信息的根, 由于zi和
Figure BDA0002550386060000052
分别位于单位圆两侧且彼此的距离非常近，即
Figure BDA0002550386060000053
因此，froot-MUSIC(z)可近似为：
Figure BDA0002550386060000054
其中
Figure BDA0002550386060000055
froot-MUSIC(z)的导数可进一步表示为：
Figure BDA0002550386060000056
其中
Figure BDA0002550386060000057
由于Qroot-MUSIC,1(z)和Γ(z)无法被Qroot-MUSIC,2(z)整除，因此 Qroot-MUSIC,2(z)分别是froot-MUSIC(z)和
Figure BDA0002550386060000058
的一重和二重不可约多项式。
本发明中组成病态方程组，具体为:
将froot-MUSIC(z)和
Figure BDA0002550386060000059
联立，组成病态方程组：
Figure BDA0002550386060000061
并且得出结论，二重不可约多项式Qroot-MUSIC,2(z)是病态方程组 IAEroot-MUSIC中的最大公因式,这意味着在root-MUSIC算法中，原始的求根多项式和其导数的最大公因式包含且仅包含真实的DOA信息,因此，对froot-MUSIC(z)进行求根运算获取DOA信息的传统处理过程可转换为对如下最大公因式(greatest common divisor，GCD)进行求根运算：
Figure BDA0002550386060000062
可以清晰地看出，GCDroot-MUSIC阶数仅为L。
本发明所述快速提取最大公因式，包括：
(1)首先将IAEroot-MUSIC中两个多项式转化为具有相同阶数的首一多项式：
Figure BDA0002550386060000063
其中，ξM-1和ψM-2分别为froot-MUSIC(z)和
Figure BDA0002550386060000064
中具有最高阶数的项的系数,由于GCDroot-MUSIC是由多个相似的根所构成，因此，按照下式进行递归运算：
gK+1(z)＝gK-1(z)-gK(z)，K∈[3,2+(2(M-1)-2-L)×2]，
最终，可以得到含有真实DOA信息的GCDroot-MUSIC
GCDroot-MUSIC＝gK(z)；
(2)fU-root-MUSIC(z)具有和froot-MUSIC(z)相同的根与系数的内在关系,因此，在fU-root-MUSIC(z)中，根的近似性和froot-MUSIC(z)相同，即
Figure BDA0002550386060000071
相似地，GCDU-root-MUSIC表示为：
Figure BDA0002550386060000072
而对于fRV-root-MUSIC(z)而言，它的根呈共轭和共轭对称分布,因此，在fRV-root-MUSIC(z)中，根的近似性应为
Figure BDA0002550386060000073
同理， GCDRV-root-MUSIC可以表示为：
Figure BDA0002550386060000074
可以看出，GCDU-root-MUSIC和GCDRV-root-MUSIC的阶数分别仅为L和2L。
本发明在得到含有真实DOA信息的最大公因式后，直接通过求根运算可以得到所有DOAs的估计值：
Figure BDA0002550386060000075
其中deg(GCD)代表GCD的阶数。
本发明针对求根类波达方向估计技术中求根多项式阶数过高的问题，提出一种基于最大公因式的广义多项式降阶求根波达方向估计新方法，通过求根多项式中系数与根的内在关系，对多项式进行近似并求其导数；同时将两者组成病态方程组后，快速提取方程组中的最大公因式，实现了求根多项式的降阶，显著地降低了计算复杂度，为波达方向估计的工程化推进提供了重要的理论支撑。
附图说明
图1是本发明的流程图。
图2是实施例1中本发明与root-MUSIC算法、U-root-MUSIC算法和 RV-root-MUSIC算法中根的分布图，其中M＝12，SNR＝5dB，N＝100， L＝2，θ1＝10°，θ2＝30°。
图3是是实施例1中本发明算法根的分布图，其中M＝12，SNR＝5dB, N＝100，L＝2，θ1＝10°，θ2＝30°。
图4是是实施例1中本发明与不同算法的RMSE随输入信噪比的变化情况，其中M＝12，N＝100，L＝2，θ1＝10°，θ2＝30°。
图5是是实施例1中本发明与不同算法的RMSE随快拍数的变化情况，其中M＝12，SNR＝5dB,L＝2,θ1＝10°,θ2＝30°。
图6是是实施例1中本发明与不同算法的计算效率随阵元数的变化情况，其中SNR＝5dB，N＝100，L＝2，θ1＝10°，θ2＝30°。
图7是不同算法中计算复杂度的对比。
具体实施方式：
下面结合附图和实施例，对本发明作进一步的说明：
本发明针对求根类波达方向估计技术中求根多项式阶数过高的问题，提出一种基于最大公因式的广义多项式降阶求根波达方向估计新方法，通过求根多项式中系数与根的内在关系，对多项式进行近似并求其导数；同时将两者组成病态方程组后，快速提取方程组中的最大公因式，实现了求根多项式的降阶，显著地降低了计算复杂度，为波达方向估计的工程化推进提供了重要的理论支撑。
如附图1所示,本发明通过以下步骤实现：
第一步，利用天线阵列接收辐射源信号，所述第一步包括以下步骤：
(1)假设天线阵列中设有M个相互独立的阵元，以d等间距组成均匀线阵(ULA)，考虑空间中存在L个远场窄带信号从方向θ＝[θ1,…,θL]入射到阵列，其中，假设L先验已知，阵元间距d满足 d≤λ/2以避免相位模糊，λ为窄带信号的波长，则天线阵列接收辐射源信号为：
Figure BDA0002550386060000091
其中，θk为空间中第k∈[1,L]个信号的来波方向，A(θ)为M×L维的阵列流型矩阵，s(t)为L×1维的入射信号矢量，n(t)为M×1维的加性高斯白噪声矢量，a(θ)为A(θ)的列向量，可以表示为：
a(z)＝[1,z,z2,…,zM-1]T,
其中，(·)T为转置运算，z＝jφ，φ＝(2π/λ)dsinθ；
(2)M×M维阵列协方差矩阵为：
Figure BDA0002550386060000092
其中，(·)H为共轭转置运算，Rss＝E[s(t)sH(t)]为L×L的信号协方差矩阵，
Figure BDA0002550386060000093
为噪声功率，
Figure BDA0002550386060000094
为单位矩阵，理论上Rxx不可知，但可以用N快拍数据对其估计：
Figure BDA0002550386060000095
其复值特征值分解可以表示为：
Figure BDA0002550386060000096
第二步：以root-MUSIC算法、U-root-MUSIC算法和 RV-root-MUSIC算法为例，获得求根类波达方向估计技术的求根多项式，所述第二步包括以下步骤：
(1)root-MUSIC算法的求根多项式为：
Figure BDA0002550386060000101
其中，(·)-1为矩阵的逆运算，ξk为多项式froot-MUSIC(z)的系数。
(2)U-root-MUSIC算法利用前后向平滑技术构造了半实值求根多项式，可表示为：
Figure BDA0002550386060000102
其中，ζk为多项式fU-root-MUSIC(z)的系数，
Figure BDA0002550386060000103
U是由单位矩阵
Figure BDA0002550386060000104
和反对角单位矩阵
Figure BDA0002550386060000105
构成的酉矩阵，且定义为：
Figure BDA0002550386060000106
(3)RV-root-MUSIC算法通过对
Figure BDA0002550386060000107
的实部进行特征值分解，构造了实值求根多项式，可表示为：
Figure BDA0002550386060000108
其中，ηk是多项式fRV-root-MUSIC(z)的系数；假设，span(·)为空间张成运算，(·)为共轭运算，那么
Figure BDA0002550386060000109
第三步：通过多项式中系数与根的内在关系，对求根多项式进行近似并求其导数，所述第三步包括以下步骤：
(1)因为噪声投影矩阵
Figure BDA0002550386060000111
和
Figure BDA0002550386060000112
是Hermitian矩阵，所以复系数ξk和ζk均是中心共轭对称，即ξk＝ξ-k，ζk＝ζ-k。因此， froot-MUSIC(z)和fU-root-MUSIC(z)中的根均以共轭对称对的形式呈现，即如果z0是froot-MUSIC(z)或fU-root-MUSIC(z)的根，那么
Figure BDA0002550386060000113
同样是对应的根。
然而，由于噪声投影矩阵
Figure BDA0002550386060000114
是对称矩阵，所以实系数ηk是对称相等，即ηk＝η-k。因此，fRV-root-MUSIC(z)中的根以共轭和共轭对称对的形式呈现，即
Figure BDA0002550386060000115
均是fRV-root-MUSIC(z)的根。
(2)为了简化而不失一般性，首先以root-MUSIC算法为例，利用前文根与系数的内在关系将froot-MUSIC(z)因式分解为：
Figure BDA0002550386060000116
其中，zi为含有真实DOA的根，zj为余下不含有DOA信息的根。值得注意的是，zi和
Figure BDA0002550386060000117
分别位于单位圆两侧且彼此的距离非常近，即
Figure BDA0002550386060000118
因此，froot-MUSIC(z)可近似为：
Figure BDA0002550386060000119
其中
Figure BDA00025503860600001110
froot-MUSIC(z)的导数可以进一步表示为：
Figure BDA0002550386060000121
其中
Figure BDA0002550386060000122
由于Qroot-MUSIC,1(z)和Γ(z)无法被Qroot-MUSIC,2(z)整除，因此 Qroot-MUSIC,2(z)分别是froot-MUSIC(z)和
Figure BDA0002550386060000123
的一重和二重不可约多项式。
第四步：组成病态方程组，所述第四步包括以下步骤：
将froot-MUSIC(z)和
Figure BDA0002550386060000124
联立，组成病态方程组：
Figure BDA0002550386060000125
并且得出结论，二重不可约多项式Qroot-MUSIC,2(z)是病态方程组 IAEroot-MUSIC中的最大公因式。这意味着在root-MUSIC算法中，原始的求根多项式和其导数的最大公因式包含且仅包含真实的DOA信息。因此，对froot-MUSIC(z)进行求根运算获取DOA信息的传统处理过程可以转换为对如下最大公因式(greatest common divisor，GCD)进行求根运算：
Figure BDA0002550386060000126
可以清晰地看出，GCDroot-MUSIC阶数仅为L。
第五步：快速提取最大公因式，所述第五步包括以下步骤：
(1)为了高效地提取出GCDroot-MUSIC从而实现求根多项式阶数的降低，首先将IAEroot-MUSIC中两个多项式转化为具有相同阶数的首一多项式：
Figure BDA0002550386060000131
其中，ξM-1和ψM-2分别为froot-MUSIC(z)和
Figure BDA0002550386060000132
中具有最高阶数的项的系数。由于GCDroot-MUSIC是由多个相似的根所构成，因此，可以按照下式进行递归运算：
gK+1(z)＝gK-1(z)-gK(z)，K∈[3,2+(2(M-1)-2-L)×2]，
最终，可以得到含有真实DOA信息的GCDroot-MUSIC
GCDroot-MUSIC＝gK(z)。
(2)根据前文可知，fU-root-MUSIC(z)具有和froot-MUSIC(z)相同的根与系数的内在关系。因此，在fU-root-MUSIC(z)中，根的近似性和 froot-MUSIC(z)相同，即
Figure BDA0002550386060000133
相似地，GCDU-root-MUSIC可以表示为：
Figure BDA0002550386060000134
而对于fRV-root-MUSIC(z)而言，它的根呈共轭和共轭对称分布。因此，在fRV-root-MUSIC(z)中，根的近似性应为
Figure BDA0002550386060000135
同理， GCDRV-root-MUSIC可以表示为：
Figure BDA0002550386060000136
Figure BDA0002550386060000142
可以看出，GCDU-root-MUSIC和GCDRV-root-MUSIC的阶数分别仅为L和2L。
值得注意的是，对于具有中心共轭对称系数的求根类波达方向估计技术而言，可以参考froot-MUSIC(z)和fU-root-MUSIC(z)来通过提取最大公因式实现求根多项式的降阶；同时，最大公因式的阶数仅为L。而对于具有对称相等系数的求根类波达方向估计技术而言，可以参考 fRV-root-MUSIC(z)来实现求根多项式的降阶，最大公因式的阶数为2L。
第六步：获得信号的波达方向，所述第六步包括以下步骤：
在得到含有真实DOA信息的最大公因式后，直接通过求根运算可以得到所有DOAs的估计值：
Figure BDA0002550386060000141
其中deg(GCD)代表GCD的阶数。
如图7所示,比较本发明与root-MUSIC算法、U-root-MUSIC算法和RV-root-MUSIC算法的计算复杂度，其中o(·)表示实值计算的计算复杂度。考虑到一次复值计算需要进行四次实值计算，因此，在EVD和求根过程中复值计算的计算复杂度为实值计算的四倍。另一方面，本发明只降低了求根多项式的阶数，所以EVD的计算复杂度保持不变。考虑到root-MUSIC算法和U-root-MUSIC算法中最大公因式GCDroot-MUSIC和GCDU-root-MUSIC的阶数仅为L，因此，两者在求根过程中仅需要4×o(L3)的计算复杂度。而GCDRV-root-MUSIC的阶数为2L，因此，对其进行求根运算的计算复杂度仅为o((2L)3)。由图7可知，本发明相比于传统求根类波达方向估计技术而言，显著地降低了求根过程中的计算复杂度。
实施例1：
本例提出了一种基于最大公因式的广义多项式降阶求根波达方向估计新方法：
第一步，利用天线阵列接收辐射源信号，所述第一步包括以下步骤：
(1)假设天线阵列中设有M个相互独立的阵元，以d等间距组成均匀线阵(ULA)，考虑空间中存在L个远场窄带信号从方向θ＝[θ1,…,θL]入射到阵列，其中，假设L先验已知，阵元间距d满足 d≤λ2以避免相位模糊，λ为窄带信号的波长，则天线阵列接收辐射源信号为：
Figure BDA0002550386060000151
其中，θk为空间中第k∈[1,L]个信号的来波方向，A(θ)为M×L维的阵列流型矩阵，s(t)为L×1维的入射信号矢量，n(t)为M×1维的加性高斯白噪声矢量，a(θ)为A(θ)的列向量，可以表示为：
a(z)＝[1,z,z2,…,zM-1]T,
其中，(·)T为转置运算，z＝jφ，φ＝(2π/λ)dsinθ；
(2)M×M维阵列协方差矩阵为：
Figure BDA0002550386060000152
其中，(·)H为共轭转置运算，Rss＝E[s(t)sH(t)]为L×L的信号协方差矩阵，
Figure BDA0002550386060000153
为噪声功率，
Figure BDA0002550386060000154
为单位矩阵。理论上Rxx不可知，但可以用N快拍数据对其估计：
Figure BDA0002550386060000155
其复值特征值分解可以表示为：
Figure BDA0002550386060000161
第二步：本例以root-MUSIC算法、U-root-MUSIC算法和 RV-root-MUSIC算法为例，获得求根类波达方向估计技术的求根多项式，具体包括以下步骤：
(1)以root-MUSIC算法的求根多项式为：
Figure BDA0002550386060000162
其中，(·)-1为矩阵的逆运算，ξk为多项式froot-MUSIC(z)的系数。
(2)U-root-MUSIC算法利用前后向平滑技术构造了半实值求根多项式，可表示为：
Figure BDA0002550386060000163
其中，ζk为多项式fU-root-MUSIC(z)的系数，
Figure BDA0002550386060000164
U是由单位矩阵
Figure BDA0002550386060000165
和反对角单位矩阵
Figure BDA0002550386060000166
构成的酉矩阵，且定义为：
Figure BDA0002550386060000167
(3)RV-root-MUSIC算法通过对
Figure BDA0002550386060000168
的实部进行特征值分解，构造了实值求根多项式，可表示为：
Figure BDA0002550386060000169
其中，ηk为多项式fRV-root-MUSIC(z)的系数；假设，span(·)为空间张成运算，(·)为共轭运算，那么
Figure BDA0002550386060000171
第三步：通过多项式中系数与根的内在关系，对求根多项式进行近似并求其导数，所述第三步包括以下步骤：
(1)因为噪声投影矩阵
Figure BDA0002550386060000172
和
Figure BDA0002550386060000173
是Hermitian矩阵，所以复系数ξk和ζk均是中心共轭对称，即ξk＝ξ-k，ζk＝ζ-k；因此， froot-MUSIC(z)和fU-root-MUSIC(z)中的根均以共轭对称对的形式呈现，即如果z0是froot-MUSIC(z)或fU-root-MUSIC(z)的根，那么
Figure BDA0002550386060000174
同样是对应的根。
然而，由于噪声投影矩阵
Figure BDA0002550386060000175
是对称矩阵，所以实系数ηk是对称相等，即ηk＝η-k，因此，fRV-root-MUSIC(z)中的根以共轭和共轭对称对的形式呈现，即
Figure BDA0002550386060000176
均是fRV-root-MUSIC(z)的根。
(2)为了简化而不失一般性，首先以root-MUSIC算法为例，利用前文中根与系数的内在关系将froot-MUSIC(z)进行因式分解为：
Figure BDA0002550386060000177
其中，zi为含有真实DOA的根，zj为余下不含有DOA信息的根。值得注意的是，zi和
Figure BDA0002550386060000178
分别位于单位圆两侧且彼此的距离非常近，即
Figure BDA0002550386060000179
因此，froot-MUSIC(z)可近似为：
Figure BDA00025503860600001710
其中
Figure BDA0002550386060000181
froot-MUSIC(z)的导数可以进一步表示为：
Figure BDA0002550386060000182
其中
Figure BDA0002550386060000183
由于Qroot-MUSIC,1(z)和Γ(z)无法被Qroot-MUSIC,2(z)整除，因此 Qroot-MUSIC,2(z)分别是froot-MUSIC(z)和
Figure BDA0002550386060000184
的一重和二重不可约多项式。
第四步：组成病态方程组，所述第四步包括以下步骤：
将froot-MUSIC(z)和
Figure BDA0002550386060000185
联立，构成病态方程组：
Figure BDA0002550386060000186
并且得出结论，二重不可约多项式Qroot-MUSIC,2(z)是病态方程组 IAEroot-MUSIC中的最大公因式。这意味着在root-MUSIC算法中，原始的求根多项式和其导数的最大公因式包含且仅包含真实的DOA信息。因此，对froot-MUSIC(z)进行求根运算获取DOA信息的传统处理过程可以转换为对如下最大公因式(greatest common divisor，GCD)进行求根运算：
Figure BDA0002550386060000187
可以清晰地看出，GCDroot-MUSIC阶数仅为L。
第五步：快速提取最大公因式，所述第五步包括以下步骤：
(1)为了高效地提取出GCDroot-MUSIC从而实现求根多项式阶数的降低，首先将IAEroot-MUSIC中两个多项式转化为具有相同阶数的首一多项式：
Figure BDA0002550386060000191
其中，ξM-1和ψM-2分别为froot-MUSIC(z)和
Figure BDA0002550386060000192
中具有最高阶数的项的系数。由于GCDroot-MUSIC是由多个相似的根所构成，因此，可以按照下式进行递归运算：
gK+1(z)＝gK-1(z)-gK(z)，K∈[3,2+(2(M-1)-2-L)×2]，
最终，可以得到含有真实DOA信息的GCDroot-MUSIC
GCDroot-MUSIC＝gK(z)。
(2)根据前文可知，fU-root-MUSIC(z)具有和froot-MUSIC(z)相同的根与系数的内在关系。因此，在fU-root-MUSIC(z)中，根的近似性和 froot-MUSIC(z)相同，即
Figure BDA0002550386060000193
相似地，GCDU-root-MUSIC可以表示为：
Figure BDA0002550386060000194
而对于fRV-root-MUSIC(z)而言，它的根呈共轭和共轭对称分布。因此，在fRV-root-MUSIC(z)中，根的近似性应为
Figure BDA0002550386060000195
同理， GCDRV-root-MUSIC可以表示为：
Figure BDA0002550386060000201
，可以看出，GCDU-root-MUSIC和GCDRV-root-MUSIC的阶数分别仅为L和2L。
值得注意的是，对于具有中心共轭对称系数的求根类波达方向估计技术，可以参考froot-MUSIC(z)和fU-root-MUSIC(z)来通过提取最大公因式实现求根多项式的降阶；同时，最大公因式的阶数仅为L。而对于具有对称相等系数的求根类波达方向估计技术，可以参考fRV-root-MUSIC(z)来实现求根多项式的降阶，最大公因式的阶数为2L。
第六步：获得信号的波达方向，所述第六步包括以下步骤：
在得到含有真实DOA信息的最大公因式后，直接通过求根运算可以得到所有DOAs的估计值：
Figure BDA0002550386060000202
其中deg(GCD)代表GCD的阶数。
本例所记载技术方案的性能可通过以下仿真说明：
仿真条件如下:假设采用12阵元的阵元间距为d＝λ2的ULA阵型，两个入射信号的方向为θ1＝10°和θ2＝30°。为了进一步评价本发明的性能，设置蒙特卡洛实验次数为500，将均方根误差(Root Mean Square Error,RMSE)作为评价指标，并引入克拉美罗界(Cramér-Rao Lower Bound,CRLB)作为对比。
仿真1，设置阵元数M＝12，信源数L＝2，比较本例所记载技术方案与root-MUSIC算法、U-root-MUSIC算法和RV-root-MUSIC算法中根的分布，其结果如图2和图3所示。
由图2和图3可知，root-MUSIC算法、U-root-MUSIC算法和 RV-root-MUSIC算法需要计算2(M-1)＝22个根，其中root-MUSIC算法、U-root-MUSIC算法各有两对根贴近单位圆，RV-root-MUSIC算法有四对根贴近单位圆。而针对root-MUSIC算法和U-root-MUSIC算法，本发明仅需要计算L＝2个根；针对RV-root-MUSIC算法，本例仅需要计算2L＝4个根。这个结果证实了此前分析的正确性。
仿真2，设置快拍数N＝100，比较本发明与root-MUSIC算法、 U-root-MUSIC算法和RV-root-MUSIC算法的RMSE随输入信噪比 (Signal-to-Noise Ratio,SNR)的变化情况，其结果如图4所示。
由图4可知，root-MUSIC算法、U-root-MUSIC算法和 RV-root-MUSIC算法在SNR≤0dB时性能略优于本发明，但是随着SNR 增加，本发明和其他三种传统算法性能十分接近。
仿真3，设置信噪比SNR＝5dB，比较本发明与root-MUSIC算法、 U-root-MUSIC算法和RV-root-MUSIC算法的RMSE随快拍数的变化情况，其结果如图5所示。
由图5可知，root-MUSIC算法、U-root-MUSIC算法和 RV-root-MUSIC算法在快拍数N≤80时性能略优于本发明，随着N的增加，本发明的RMSE和其他三种传统算法性能十分接近。
仿真4，比较本发明与不同算法的计算效率，其中通过在Intel(R) Core(TM)i5-9400 2.90GHz CPU处理器和16GB RAM内存的同一PC 环境中运行MATLAB代码，从CPU时间的角度对计算效率进行等效评估。
由图6可知，相比另外三种算法，本发明的计算效率更高；结合上述仿真可知，本发明同时可以提供优良的角度估计精度。
综上,可以说明本发明相对于现有技术具有显著的优越性。

Claims (1)

1.一种远场窄带无线电信号波达方向估计方法，其特征在于,接收辐射源信号后，获取辐射源信号的求根多项式，然后通过根与系数的关系，构造近似求根多项式及其导数；将所述近似求根多项式及其导数组成病态方程组后，提取方程组的最大公因式，在得到含有真实DOA信息的最大公因式后，直接通过求根运算得到信号的波达方向；其中快速提取最大公因式，包括：
(1)首先将IAEroot-MUSIC中两个多项式转化为具有相同阶数的首一多项式：
Figure FDA0003699924130000011
其中，ξM-1和ψM-2分别为froot-MUSIC(z)和
Figure FDA0003699924130000012
中具有最高阶数的项的系数,由于GCDroot-MUSIC是由多个相似的根所构成，因此，按照下式进行递归运算：
gK+1(z)＝gK-1(z)-gK(z)，K∈[3,2+(2(M-1)-2-L)×2]，
最终，可得到含有真实DOA信息的GCDroot-MUSIC
GCDroot-MUSIC＝gK(z)；
(2)fU-root-MUSIC(z)具有和froot-MUSIC(z)相同的根与系数的内在关系,因此，在fU-root-MUSIC(z)中，根的近似性和froot-MUSIC(z)相同，即
Figure FDA0003699924130000013
相似地，GCDU-root-MUSIC表示为：
Figure FDA0003699924130000014
而fRV-root-MUSIC(z)的根呈共轭和共轭对称分布,因此，在fRV-root-MUSIC(z)中，根的近似性应为
Figure FDA0003699924130000015
GCDRV-root-MUSIC表示为：
Figure FDA0003699924130000021
在得到含有真实DOA信息的最大公因式后，直接通过求根运算得到所有DOAs的估计值：
Figure FDA0003699924130000022
其中deg(GCD)代表GCD的阶数；
所述辐射源信号通过天线阵列接收获得，假设天线阵列中设有M个相互独立的阵元，以d等间距组成均匀线阵(ULA)，考虑空间中存在L个远场窄带信号从方向θ＝[θ1,…,θL]入射到阵列，其中，假设L先验已知，阵元间距d满足d≤λ/2以避免相位模糊，λ为窄带信号的波长，则天线阵列接收辐射源信号为：
Figure FDA0003699924130000023
其中，θk为空间中第k∈[1,L]个信号的来波方向，A(θ)为M×L维的阵列流型矩阵，s(t)为L×1维的入射信号矢量，n(t)为M×1维的加性高斯白噪声矢量，a(θ)为A(θ)的列向量，表示为：
a(z)＝[1,z,z2,…,zM-1]T,
其中，(·)T为转置运算，z＝jφ，φ＝(2π/λ)dsinθ；
M×M维阵列协方差矩阵为：
Figure FDA0003699924130000024
其中，(·)H为共轭转置运算，Rss＝E[s(t)sH(t)]为L×L的信号协方差矩阵，
Figure FDA0003699924130000025
为噪声功率，
Figure FDA0003699924130000028
为单位矩阵，理论上Rxx不可知，但用N个快拍数据对其估计：
Figure FDA0003699924130000026
其复值特征值分解表示为：
Figure FDA0003699924130000027
所述获取辐射源信号的求根多项式是指通过root-MUSIC算法、U-root-MUSIC算法或RV-root-MUSIC算法，获得求根类波达方向估计技术的求根多项式,具体分别如下：
root-MUSIC算法的求根多项式为：
Figure FDA0003699924130000031
其中，(·)-1为矩阵的逆运算，ξk为多项式froot-MUSIC(z)的系数；
U-root-MUSIC算法利用前后向平滑技术构造了半实值求根多项式，表示为：
Figure FDA0003699924130000032
其中，ζk为多项式fU-root-MUSIC(z)的系数，
Figure FDA0003699924130000033
U是由单位矩阵
Figure FDA0003699924130000034
和反对角单位矩阵
Figure FDA0003699924130000035
构成的酉矩阵，且定义为：
Figure FDA0003699924130000036
RV-root-MUSIC算法通过对
Figure FDA0003699924130000037
的实部进行特征值分解，构造了实值求根多项式，表示为：
Figure FDA0003699924130000038
其中，ηk是多项式fRV-root-MUSIC(z)的系数；假设，span(·)为空间张成运算，(·)为共轭运算，那么
Figure FDA0003699924130000039
通过根与系数的关系，构造近似求根多项式及其导数,其中root-MUSIC算法的情况下,利用根与系数的内在关系将froot-MUSIC(z)因式分解为：
Figure FDA00036999241300000310
其中，zi为含有真实DOA的根，zj为余下不含有DOA信息的根,由于zi和
Figure FDA0003699924130000041
分别位于单位圆两侧且彼此的距离非常近，即
Figure FDA0003699924130000042
因此，froot-MUSIC(z)近似为：
Figure FDA0003699924130000043
其中
Figure 1
；froot-MUSIC(z)的导数进一步表示为：
Figure FDA0003699924130000045
其中
Figure FDA0003699924130000046
由于Qroot-MUSIC,1(z)和Γ(z)无法被Qroot-MUSIC,2(z)整除，因此Qroot-MUSIC,2(z)分别是froot-MUSIC(z)和
Figure FDA0003699924130000047
的一重和二重不可约多项式；
所述组成病态方程组具体为:
将froot-MUSIC(z)和
Figure FDA0003699924130000048
联立，组成病态方程组：
Figure FDA0003699924130000049
并且得出结论二重不可约多项式Qroot-MUSIC,2(z)是病态方程组IAEroot-MUSIC中的最大公因式,这意味着对froot-MUSIC(z)进行求根运算获取DOA信息的传统处理过程转换为对如下最大公因式进行求根运算：
Figure FDA0003699924130000051
CN202010573169.1A 2020-06-22 2020-06-22 远场窄带无线电信号波达方向估计方法 Active CN113219398B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010573169.1A CN113219398B (zh) 2020-06-22 2020-06-22 远场窄带无线电信号波达方向估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010573169.1A CN113219398B (zh) 2020-06-22 2020-06-22 远场窄带无线电信号波达方向估计方法

Publications (2)

Publication Number Publication Date
CN113219398A CN113219398A (zh) 2021-08-06
CN113219398B true CN113219398B (zh) 2022-09-13

Family

ID=77085751

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010573169.1A Active CN113219398B (zh) 2020-06-22 2020-06-22 远场窄带无线电信号波达方向估计方法

Country Status (1)

Country Link
CN (1) CN113219398B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102707264A (zh) * 2012-06-13 2012-10-03 西安电子科技大学 基于圆形阵列双基地mimo雷达的波达方向估计方法
CN103323827A (zh) * 2013-05-27 2013-09-25 杭州电子科技大学 基于快速傅里叶变换的mimo雷达系统角度估计方法
CN104123462A (zh) * 2014-07-22 2014-10-29 西安电子科技大学 实多项式求根实现均匀线阵的谱music方法
CN108398659A (zh) * 2018-02-11 2018-08-14 华南理工大学 一种矩阵束与求根music结合的波达方向估计方法
CN108896954A (zh) * 2018-06-07 2018-11-27 河海大学 互质阵中一种基于联合实值子空间的波达角估计方法
CN110197112A (zh) * 2019-03-19 2019-09-03 上海大学 一种基于协方差修正的波束域Root-MUSIC方法
CN110749858A (zh) * 2019-09-26 2020-02-04 南京航空航天大学 一种基于多项式求根的展开互质阵测向估计方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102707264A (zh) * 2012-06-13 2012-10-03 西安电子科技大学 基于圆形阵列双基地mimo雷达的波达方向估计方法
CN103323827A (zh) * 2013-05-27 2013-09-25 杭州电子科技大学 基于快速傅里叶变换的mimo雷达系统角度估计方法
CN104123462A (zh) * 2014-07-22 2014-10-29 西安电子科技大学 实多项式求根实现均匀线阵的谱music方法
CN108398659A (zh) * 2018-02-11 2018-08-14 华南理工大学 一种矩阵束与求根music结合的波达方向估计方法
CN108896954A (zh) * 2018-06-07 2018-11-27 河海大学 互质阵中一种基于联合实值子空间的波达角估计方法
CN110197112A (zh) * 2019-03-19 2019-09-03 上海大学 一种基于协方差修正的波束域Root-MUSIC方法
CN110749858A (zh) * 2019-09-26 2020-02-04 南京航空航天大学 一种基于多项式求根的展开互质阵测向估计方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Real-Valued DOA Estimation for Non-Circular Sources via Reduced-Order Polynomial Rooting;Xiangtian Meng et al.;《IEEE Access》;20191030;第7卷;第158892-158903页 *
基于实值分解技术的循环root-MUSIC算法;刘志刚 等;《系统仿真学报》;20060930;第18卷(第9期);第2438-2441页 *
基于谱分解的降阶求根MUSIC算法;闫锋刚 等;《电子与信息学报》;20171031;第39卷(第10期);第2421-2427页 *

Also Published As

Publication number Publication date
CN113219398A (zh) 2021-08-06

Similar Documents

Publication Publication Date Title
WO2021139208A1 (zh) 一种基于特定频率组合信号的一维doa估计方法
CN106980106B (zh) 阵元互耦下的稀疏doa估计方法
CN108896954B (zh) 互质阵中一种基于联合实值子空间的波达角估计方法
CN107589399B (zh) 基于多采样虚拟信号奇异值分解的互质阵列波达方向估计方法
CN110109050B (zh) 嵌套阵列下基于稀疏贝叶斯的未知互耦的doa估计方法
CN109375152B (zh) 电磁矢量嵌套l阵下低复杂度的doa与极化联合估计方法
Yang et al. A unified array geometry composed of multiple identical subarrays with hole-free difference coarrays for underdetermined DOA estimation
CN112379327A (zh) 一种基于秩损估计的二维doa估计与互耦校正方法
CN109557504B (zh) 一种近场窄带信号源的定位方法
CN111965591A (zh) 一种基于四阶累积量矢量化dft的测向估计方法
CN112731278A (zh) 一种部分极化信号的角度与极化参数欠定联合估计方法
Tan et al. Covariance matrix reconstruction for direction finding with nested arrays using iterative reweighted nuclear norm minimization
CN111366893A (zh) 一种均匀圆阵未知互耦条件下的非圆信号方位角估计方法
CN111352063A (zh) 一种均匀面阵中基于多项式求根的二维测向估计方法
Akkar et al. Reactance domains unitary MUSIC algorithms based on real-valued orthogonal decomposition for electronically steerable parasitic array radiator antennas
CN116299150B (zh) 一种均匀面阵中降维传播算子的二维doa估计方法
CN113219398B (zh) 远场窄带无线电信号波达方向估计方法
CN115329261A (zh) 一种基于空间平滑稀疏重构的mimo雷达低仰角估计方法
CN113219400B (zh) 基于无空洞互质阵列结构的欠定波达方向估计方法
CN114325568A (zh) 脉冲噪声环境下基于bnc的嵌套阵列非圆信号doa估计方法
CN114397619A (zh) 基于非均匀稀疏阵列二维定位算法
WO2021068494A1 (zh) 基于平面互质阵列虚拟域张量空间谱搜索的高分辨精确二维波达方向估计方法
Chen et al. A novel ULA-difference-coarray-based DOA estimation method for general coherent signals
CN112763972A (zh) 基于稀疏表示的双平行线阵二维doa估计方法及计算设备
CN112698263A (zh) 一种基于正交传播算子的单基地互质mimo阵列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