CN112881971B - 电磁方向性互耦效应下相干干扰源测向方法 - Google Patents

电磁方向性互耦效应下相干干扰源测向方法 Download PDF

Info

Publication number
CN112881971B
CN112881971B CN202110046094.6A CN202110046094A CN112881971B CN 112881971 B CN112881971 B CN 112881971B CN 202110046094 A CN202110046094 A CN 202110046094A CN 112881971 B CN112881971 B CN 112881971B
Authority
CN
China
Prior art keywords
coherent
array
angle
matrix
mutual coupling
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
CN202110046094.6A
Other languages
English (en)
Other versions
CN112881971A (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN202110046094.6A priority Critical patent/CN112881971B/zh
Publication of CN112881971A publication Critical patent/CN112881971A/zh
Application granted granted Critical
Publication of CN112881971B publication Critical patent/CN112881971B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Direction-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/02Direction-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/14Systems for determining direction or deviation from predetermined direction
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE 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/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing 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)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供了一种电磁方向性互耦效应下相干干扰源测向方法,计算阵列得到的协方差矩阵并进行特征值分解,利用均匀线阵阵列互耦矩阵的结构特征,分解阵列实际导向矢量,并构造方向性电磁耦合效应下估计相干信号到达角的代价函数,对代价函数进行求解并进行多维谱峰搜索,估计到达角,使用估计出的相干信源的到达角估计角度依赖的阵列互耦向量,使用MVP算法估计出相干信源的到达角。本发明解决接收阵列存在角度相关互耦条件下的相干信源的到达角估计问题,本发明能够实现阵列互耦的自校正,并且在低信噪比条件下,能够实现相干信源高分辨率的到达角估计。

Description

电磁方向性互耦效应下相干干扰源测向方法
技术领域
本发明涉及信号处理领域,是一种相干信号源的到达角估计方法,特别涉及一种相干干扰源测向方法。
背景技术
信号空间谱估计技术广泛应用于雷达、通信、声纳、勘测、射电天文及生物医学工程等众多军事及国民经济领域。最具代表性的空间谱估计算法有子空间分解类算法和子空间拟合类算法,这些算法在理想的情况下都具有高分辨率的空间谱估计性能。在实际工程应用中,超分辨方法对阵元位置误差,通道不一致以及阵元间互耦等导致的阵列流型误差非常敏感。在实际环境中,由于多径传播,欺骗式干扰等因素使得空间中存在大量的相干信号,在这种情况下,MUSIC、ESPRIT等子空间分解类算法无法有效分辨或测向。
目前空间谱估计方法和阵元互耦校正方法可分为两类:有源校正法和自校正法,有源校正法需要预先知道校正源的精确位置信息,但已知精确位置的校正源较难得到。常用的自校正方法有辅助阵元法和秩损法。这两种方法,都存在盲角,且在互耦效应较强时,将会受到空间谱中伪峰的影响,导致估计性能大幅度下降。
对于相干信源的处理方法,目前有空间平滑类算法、矩阵重构类算法和非降维处理算法,另一方面基于子空间拟合类的空间谱估计算法对于相干信源具有良好的估计性能。
发明内容
为了克服现有技术的不足,本发明提供一种电磁方向性互耦效应下相干干扰源测向方法。本发明的目的是解决在同时受相干信源和方向性电磁耦合效应影响,超分辨算法无法有效分辨或测向问题。而提出了一种方向性电磁耦合效应下相干信源空域信息联合估计技术。可显著提高非理想条件下,方向估计的成功率及分辨率,且实现了阵元互耦角度依赖情况下的自校准。
本发明解决其技术问题所采用的技术方案包括以下步骤:
步骤1:计算阵列得到的接收信号x(n)的协方差矩阵Rx,并对Rx进行特征值分解;
步骤2:利用均匀线阵阵列互耦矩阵的结构特征,分解阵列实际导向矢量;
对于任意到达角θ,
Figure BDA0002897205170000021
分解为:
Figure BDA0002897205170000022
其中:
Figure BDA0002897205170000023
[T1(θ)]p,r表示矩阵T1(θ)中第p行第r列的元素,[c(θ)]p=[C(θ)]1,p,T(θ)和c(θ)的维度分别为M×m和m×1,m为阵元互耦长度。
步骤3:使用子空间拟合算法,利用新分解的阵列实际导向矢量T(θ),构造方向性电磁耦合效应下估计相干信号到达角的代价函数;
步骤4:对代价函数进行求解并进行多维谱峰搜索,估计到达角;
对于式(7),固定
Figure BDA0002897205170000024
求出Γ的最小二乘解:
Figure BDA0002897205170000025
将式(9)代入式(7)得到:
Figure BDA0002897205170000026
或:
Figure BDA0002897205170000027
其中,
Figure BDA0002897205170000028
为矩阵/>
Figure BDA0002897205170000029
零空间上的正交投影矩阵,/>
Figure BDA00028972051700000210
I为与/>
Figure BDA00028972051700000211
维度相同的单位阵,令θ∈[-90°,90°],以步长为θstep划分网格,求每个网格点上/>
Figure BDA00028972051700000212
的值,其最大值所对应的网格点的多维坐标即为估计出的方向性电磁耦合效应下相干信源的到达角
Figure BDA00028972051700000213
步骤5:使用估计出的相干信源的到达角
Figure BDA00028972051700000214
估计角度依赖的阵列互耦向量;
此时,
Figure BDA00028972051700000215
已知,则将/>
Figure BDA00028972051700000216
代入式(5)求得/>
Figure BDA00028972051700000217
得到/>
Figure BDA00028972051700000218
将/>
Figure BDA00028972051700000219
代入式(9),得到矩阵/>
Figure BDA00028972051700000220
从式(8)可知,矩阵/>
Figure BDA00028972051700000221
由与角度相关的互耦向量c(θq)与辅助参量Z组成,根据/>
Figure BDA00028972051700000222
组成特性,互耦向量由式(12)计算获得:
Figure BDA0002897205170000031
步骤6:得到角度相关的互耦矢量c(θq)后,使用MVP(modified variableprojection)算法估计出相干信源的到达角。
所述步骤1中,对于M个阵元的均匀线阵,空间中有Q个相干的期望信号,则M元均匀线阵在存在方向性电磁耦合效应下阵列的接收模型表示为:
Figure BDA0002897205170000032
其中,θq为第q个相干信源的到达角,q=1,2,...,Q,N为快拍数,C(θq)为与第q个相干信号相关的互耦矩阵,相关的互耦a(θq)为与第q个相干信号相关的阵列导向矢量,sq(n)表示第q个相干信号的信号向量,n(n)为独立同分布加性高斯白噪声矢量,定义
Figure BDA0002897205170000033
为方向性电磁耦合效应下与第q个信号相关的阵列实际导向矢量;
采用有限次样本计算的阵列数据协方差矩阵由下式计算:
Figure BDA0002897205170000034
对数据协方差矩阵进行特征值分解得到:
Figure BDA0002897205170000035
其中,Us和Un分别表示信号子空间和噪声子空间,∑s为协方差矩阵Q个较大的特征值组成的对角阵,∑n为(M-Q)个较小的特征值组成的对角阵。
所述步骤3中,加权子空间拟合问题,表示为:
Figure BDA0002897205170000036
其中,W为最有权矩阵,Z为辅助参量,是一个维度为Q×Q的满秩矩阵,
Figure BDA0002897205170000037
由于
Figure BDA0002897205170000038
中包含未知参量c(θ),所以将式(6)转换成如下形式,得到代价函数:
Figure BDA0002897205170000039
其中,
Figure BDA00028972051700000310
且:
Figure BDA0002897205170000041
所述MVP算法获得相干信源的到达角的步骤为:
代价函数写为:
V=tr(PAKKH),K=UsW1/2 (13)
代价函数的一阶导数为:
Figure BDA0002897205170000042
式中,矩阵D为A的一阶导数;
在解的过程中需要确定搜索的方向:
θk+1=θk-Δθ,Δθ=H-1V′ (15)
其中,θk=[θk,1,θk,2,...,θk,Q]T为第k次迭代的初始角度,k=1,2,...,K,θk,1表示第k次迭代第q个相干信源的来向,Δθ=[Δθ1,Δθ2,...,ΔθQ]为第k次迭代时Q个信源的信号搜索方向;
Figure BDA0002897205170000043
利用式(14)和(16)确定矩阵V′和H,从而由式(15)得到角度的搜索方向,进而得到更新的角度θk+1,当完成K次迭代后停止迭代,得到相干信源的到达角。
本发明的有益效果在于本发明提供的一种方向性电磁耦合效应下相干信源空域信息联合估计方法,可以解决接收阵列存在角度相关互耦条件下的相干信源的到达角估计问题,本发明能够实现阵列互耦的自校正,并且在低信噪比条件下,能够实现相干信源高分辨率的到达角估计。
附图说明
图1为本发明的相干干扰源测向流程示意图。
图2为本发明DOA估计结果。
图3为本发明到达角估计均方根误差与信噪比关系示意图。
图4为本发明互耦矢量估计均方根误差与信噪比关系示意图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
为了加深对本发明的理解,下面将结合实施例对本发明作进一步详述,该实施例仅用于解释本发明,并不构成对本发明保护范围的限定。
步骤1:计算阵列得到的接收信号x(n)的协方差矩阵Rx,并对Rx进行特征值分解;
对于M个阵元的均匀线阵,空间中有Q个相干的期望信号,则M元均匀线阵在存在方向性电磁耦合效应下阵列的接收模型表示为:
Figure BDA0002897205170000051
其中,θq为第q个相干信源的到达角,q=1,2,...,Q,N为快拍数,C(θq)为与第q个相干信号相关的互耦矩阵,相关的互耦a(θq)为与第q个相干信号相关的阵列导向矢量,sq(n)表示第q个相干信号的信号向量,n(n)为独立同分布加性高斯白噪声矢量,定义
Figure BDA0002897205170000052
为方向性电磁耦合效应下与第q个信号相关的阵列实际导向矢量;
在实际情况下,采用有限次样本计算的阵列数据协方差矩阵由下式计算:
Figure BDA0002897205170000053
对数据协方差矩阵进行特征值分解得到:
Figure BDA0002897205170000054
其中,Us和Un分别表示信号子空间和噪声子空间,∑s为协方差矩阵Q个较大的特征值组成的对角阵,∑n为(M-Q)个较小的特征值组成的对角阵;
步骤2:利用均匀线阵阵列互耦矩阵的结构特征,分解阵列实际导向矢量;
对于任意到达角θ,
Figure BDA0002897205170000055
分解为:
Figure BDA0002897205170000056
其中:
Figure BDA0002897205170000057
[T1(θ)]p,r表示矩阵T1(θ)中第p行第r列的元素,[c(θ)]p=[C(θ)]1,p,T(θ)和c(θ)的维度分别为M×m和m×1,m为阵元互耦长度。
步骤3:使用子空间拟合算法,利用新分解的阵列实际导向矢量T(θ),构造方向性电磁耦合效应下估计相干信号到达角的代价函数;
加权子空间拟合问题,表示为:
Figure BDA0002897205170000061
其中,W为最有权矩阵,Z为辅助参量,是一个维度为Q×Q的满秩矩阵,
Figure BDA0002897205170000062
由于
Figure BDA0002897205170000063
中包含未知参量c(θ),所以将式(6)转换成如下形式,得到代价函数:
Figure BDA0002897205170000064
其中,
Figure BDA0002897205170000065
且:
Figure BDA0002897205170000066
步骤4:对代价函数进行求解并进行多维谱峰搜索,估计到达角;
对于式(7),固定
Figure BDA0002897205170000067
求出Γ的最小二乘解:
Figure BDA0002897205170000068
将式(9)代入式(7)得到:
Figure BDA0002897205170000069
或:
Figure BDA00028972051700000610
其中,
Figure BDA00028972051700000611
为矩阵/>
Figure BDA00028972051700000612
零空间上的正交投影矩阵,/>
Figure BDA00028972051700000613
I为与/>
Figure BDA00028972051700000614
维度相同的单位阵,令θ∈[-90°,90°],以步长为θstep划分网格,求每个网格点上/>
Figure BDA00028972051700000615
的值,其最大值所对应的网格点的多维坐标即为估计出的方向性电磁耦合效应下相干信源的到达角
Figure BDA00028972051700000616
步骤5:使用估计出的相干信源的到达角
Figure BDA00028972051700000617
估计角度依赖的阵列互耦向量;
此时,
Figure BDA0002897205170000071
已知,则将/>
Figure BDA0002897205170000072
代入式(5)求得/>
Figure BDA0002897205170000073
得到/>
Figure BDA0002897205170000074
将/>
Figure BDA0002897205170000075
代入式(9),得到矩阵/>
Figure BDA0002897205170000076
从式(8)可知,矩阵/>
Figure BDA0002897205170000077
由与角度相关的互耦向量c(θq)与辅助参量Z组成,根据/>
Figure BDA0002897205170000078
组成特性,互耦向量由式(12)计算获得:
Figure BDA0002897205170000079
步骤6:得到角度相关的互耦矢量c(θq)后,使用MVP(modified variableprojection)算法进一步估计出相干信源的到达角;
当角度相关的互耦矢量c(θq)后,直接将c(θq)代入式(6),获得辅助参量Z的最小二乘解
Figure BDA00028972051700000710
并将/>
Figure BDA00028972051700000711
代入式(6),得到/>
Figure BDA00028972051700000712
通过网格划分找到/>
Figure BDA00028972051700000713
的最大值估计出/>
Figure BDA00028972051700000714
但这种方法需要再次进行多维搜索,算法的计算量过大。因此,接下来使用MVP算法获得相干信源的到达角;
代价函数写为:
V=tr(PAKKH),K=UsW1/2 (13)
代价函数的一阶导数为:
Figure BDA00028972051700000715
式中,矩阵D为A的一阶导数;
在解的过程中需要确定搜索的方向:
θk+1=θk-Δθ,Δθ=H-1V′ (15)
其中,θk=[θk,1,θk,2,...,θk,Q]T为第k次迭代的初始角度,k=1,2,...,K,θk,1表示第k次迭代第q个相干信源的来向,Δθ=[Δθ1,Δθ2,...,ΔθQ]为第k次迭代时Q个信源的信号搜索方向;
Figure BDA00028972051700000716
利用式(14)和(16)确定矩阵V′和H,从而由式(15)得到角度的搜索方向,进而得到更新的角度θk+1,当完成K次迭代后停止迭代,得到相干信源的到达角。
如图1-图4所示,本实施例提供一种方向性电磁耦合效应下相干信源空域信息联合估计技术,如图1所示,过程如下。
步骤1:根据均匀线阵接收信号的数学模型,获得阵列得到的接收信号的数据矩阵x(n),并计算接收信号的协方差矩阵Rx,并对Rx进行特征值分解。
步骤2:根据均匀线阵导向矢量的空间域模型与均匀线阵阵列互耦矩阵的Toeplitz结构,将存在互耦的均匀线阵实际导向矢量进行模拟等效,并重新分解阵列实际导向矢量。
步骤3:使用子空间拟合算法,利用新分解的阵列实际导向矢量,构造方向性电磁耦合效应下估计相干信号到达角的代价函数。
步骤4:对代价函数进行求解,在阵列可接收来向范围内,对可能来向角度进行网格划分,在每一个网格点上求代价函数解的具体值,并搜索得到最大值所在的网格点,即为初次估计出的到达角。
步骤5:分析步骤4中解出的辅助矩阵,并使用初次估计出的到达角估计出估计角度依赖的阵列互耦矢量。
步骤6:得到角度相关的互耦矢量后,为降低算法的计算复杂度,使用MVP算法迭代估计出更高精度的相干信源到达角,MVP算法的初始值使用初次估计出的到达角。
具体实施例:
初始化参数设置:接收阵列为阵元数M=8的均匀线阵,接收Q=2个相干信号,相干信号的来向分别为-30°和-35°,互耦长度为3,互耦向量为c1=[1 0.7361-0.4250i 0.2000+0.3464i]T,c2=[1 0.4500+0.7794i 0.2474-0.2474i]T。信源的信噪比为10dB,采样快拍数为200。
进行500次蒙特卡洛实验,并对估计出的到达角和互耦矢量求解其随信噪比变化的均方根误差(RMSE),其均方根误差公式分别为:
Figure BDA0002897205170000081
Figure BDA0002897205170000082
采用本发明方法可以估计出相干信源的到达角及互耦矢量,图2为多次实验的DOA估计结果,圆圈代表实际信源到达角,*为多次估计出的到达角,如图2所示,本发明可准确估计出互耦角度依赖情况下的相干信源的到达角,图3为到达角估计均方根误差与信噪比关系示意图,图3中,使用常用的MUSIC算法完全失效,无法估计互耦角度依赖情况下的相干信源的到达角,所提算法估计性能接近克拉美罗下界(CRLB),图4为互耦矢量估计均方根误差与信噪比关系示意图。由图2图-4可见,本实施例方法具有更好的到达角和互耦矢量估计性能,本发明方法具有明显的优越性。

Claims (4)

1.一种电磁方向性互耦效应下相干干扰源测向方法,其特征在于包括下述步骤:
步骤1:计算阵列得到的接收信号x(n)的协方差矩阵Rx,并对Rx进行特征值分解;
步骤2:利用均匀线阵阵列互耦矩阵的结构特征,分解阵列实际导向矢量;
对于任意到达角θ,
Figure FDA0002897205160000011
分解为:
Figure FDA0002897205160000012
其中:
Figure FDA0002897205160000013
[T1(θ)]p,r表示矩阵T1(θ)中第p行第r列的元素,[c(θ)]p=[C(θ)]1,p,T(θ)和c(θ)的维度分别为M×m和m×1,m为阵元互耦长度;
步骤3:使用子空间拟合算法,利用新分解的阵列实际导向矢量T(θ),构造方向性电磁耦合效应下估计相干信号到达角的代价函数;
步骤4:对代价函数进行求解并进行多维谱峰搜索,估计到达角;
对于式(7),固定
Figure FDA0002897205160000014
求出Γ的最小二乘解:
Figure FDA0002897205160000015
将式(9)代入式(7)得到:
Figure FDA0002897205160000016
或:
Figure FDA0002897205160000017
其中,
Figure FDA0002897205160000018
为矩阵/>
Figure FDA0002897205160000019
零空间上的正交投影矩阵,/>
Figure FDA00028972051600000110
I为与/>
Figure FDA00028972051600000111
维度相同的单位阵,令θ∈[-90°,90°],以步长为θstep划分网格,求每个网格点上/>
Figure FDA00028972051600000112
的值,其最大值所对应的网格点的多维坐标即为估计出的方向性电磁耦合效应下相干信源的到达角
Figure FDA00028972051600000113
步骤5:使用估计出的相干信源的到达角
Figure FDA00028972051600000114
估计角度依赖的阵列互耦向量;
此时,
Figure FDA00028972051600000115
已知,则将/>
Figure FDA00028972051600000116
代入式(5)求得/>
Figure FDA00028972051600000117
得到/>
Figure FDA00028972051600000118
将/>
Figure FDA0002897205160000021
代入式(9),得到矩阵/>
Figure FDA0002897205160000022
从式(8)可知,矩阵/>
Figure FDA0002897205160000023
由与角度相关的互耦向量c(θq)与辅助参量Z组成,根据/>
Figure FDA0002897205160000024
组成特性,互耦向量由式(12)计算获得:
Figure FDA0002897205160000025
步骤6:得到角度相关的互耦矢量c(θq)后,使用MVP(modified variable projection)算法估计出相干信源的到达角。
2.根据权利要求1所述的电磁方向性互耦效应下相干干扰源测向方法,其特征在于:
所述步骤1中,对于M个阵元的均匀线阵,空间中有Q个相干的期望信号,则M元均匀线阵在存在方向性电磁耦合效应下阵列的接收模型表示为:
Figure FDA0002897205160000026
其中,θq为第q个相干信源的到达角,q=1,2,…,Q,N为快拍数,C(θq)为与第q个相干信号相关的互耦矩阵,相关的互耦a(θq)为与第q个相干信号相关的阵列导向矢量,sq(n)表示第q个相干信号的信号向量,n(n)为独立同分布加性高斯白噪声矢量,定义
Figure FDA0002897205160000027
为方向性电磁耦合效应下与第q个信号相关的阵列实际导向矢量;
采用有限次样本计算的阵列数据协方差矩阵由下式计算:
Figure FDA0002897205160000028
对数据协方差矩阵进行特征值分解得到:
Figure FDA0002897205160000029
其中,Us和Un分别表示信号子空间和噪声子空间,∑s为协方差矩阵Q个较大的特征值组成的对角阵,∑n为(M-Q)个较小的特征值组成的对角阵。
3.根据权利要求1所述的电磁方向性互耦效应下相干干扰源测向方法,其特征在于:
所述步骤3中,加权子空间拟合问题,表示为:
Figure FDA00028972051600000210
其中,W为最有权矩阵,Z为辅助参量,是一个维度为Q×Q的满秩矩阵,
Figure FDA00028972051600000211
由于
Figure FDA00028972051600000212
中包含未知参量c(θ),所以将式(6)转换成如下形式,得到代价函数:
Figure FDA0002897205160000031
其中,
Figure FDA0002897205160000032
且:
Figure FDA0002897205160000033
4.根据权利要求1所述的电磁方向性互耦效应下相干干扰源测向方法,其特征在于:
所述MVP算法获得相干信源的到达角的步骤为:
代价函数写为:
V=tr(PAKKH),K=UsW1/2 (13)
代价函数的一阶导数为:
Figure FDA0002897205160000034
式中,矩阵D为A的一阶导数;
在解的过程中需要确定搜索的方向:
θk+1=θk-Δθ,Δθ=H-1V′ (15)
其中,θk=[θk,1,θk,2,…,θk,Q]T为第k次迭代的初始角度,k=1,2,…,K,θk,1表示第k次迭代第q个相干信源的来向,Δθ=[Δθ1,Δθ2,…,ΔθQ]为第k次迭代时Q个信源的信号搜索方向;
Figure FDA0002897205160000035
利用式(14)和(16)确定矩阵V′和H,从而由式(15)得到角度的搜索方向,进而得到更新的角度θk+1,当完成K次迭代后停止迭代,得到相干信源的到达角。
CN202110046094.6A 2021-01-14 2021-01-14 电磁方向性互耦效应下相干干扰源测向方法 Active CN112881971B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110046094.6A CN112881971B (zh) 2021-01-14 2021-01-14 电磁方向性互耦效应下相干干扰源测向方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110046094.6A CN112881971B (zh) 2021-01-14 2021-01-14 电磁方向性互耦效应下相干干扰源测向方法

Publications (2)

Publication Number Publication Date
CN112881971A CN112881971A (zh) 2021-06-01
CN112881971B true CN112881971B (zh) 2023-06-30

Family

ID=76045917

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110046094.6A Active CN112881971B (zh) 2021-01-14 2021-01-14 电磁方向性互耦效应下相干干扰源测向方法

Country Status (1)

Country Link
CN (1) CN112881971B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114679227B (zh) * 2022-03-25 2023-07-14 电子科技大学 一种测向误差的空间频域校正方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101375425B1 (ko) * 2013-02-28 2014-03-17 충북대학교 산학협력단 배열 안테나에서 전파의 도래각 추정 방법
CN108680891A (zh) * 2018-01-05 2018-10-19 大连大学 非均匀噪声条件下考虑互耦效应的doa估计方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101375425B1 (ko) * 2013-02-28 2014-03-17 충북대학교 산학협력단 배열 안테나에서 전파의 도래각 추정 방법
CN108680891A (zh) * 2018-01-05 2018-10-19 大连大学 非均匀噪声条件下考虑互耦效应的doa估计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于RARE-Cumulant的互耦校正和DOA估计;李新波;石要武;马彦;张涵;;吉林大学学报(工学版)(第04期);全文 *

Also Published As

Publication number Publication date
CN112881971A (zh) 2021-06-01

Similar Documents

Publication Publication Date Title
CN110113085B (zh) 一种基于协方差矩阵重构的波束形成方法及系统
CN108957387B (zh) 一种卫星信号二维到达角估计方法及系统
CN107870314B (zh) 基于极化敏感阵列的完备电磁分量加权融合测向优化方法
CN101149429A (zh) 适用于均匀圆阵的阵列互耦校正与信源测向方法
CN107576931B (zh) 一种基于协方差低维度迭代稀疏重构的相关/相干信号波达方向估计方法
CN111046591B (zh) 传感器幅相误差与目标到达角度的联合估计方法
CN112379327A (zh) 一种基于秩损估计的二维doa估计与互耦校正方法
CN104950297A (zh) 基于矩阵1范数拟合的阵元误差估计方法
CN109507635A (zh) 利用两个未知方位辅助源的阵列幅相误差估算方法
US11681006B2 (en) Method for jointly estimating gain-phase error and direction of arrival (DOA) based on unmanned aerial vehicle (UAV) array
CN106980104A (zh) 用于传感器阵列的信号波达方向自校正方法
CN112766304A (zh) 一种基于稀疏贝叶斯学习的机动阵列方位估计方法
CN112881971B (zh) 电磁方向性互耦效应下相干干扰源测向方法
CN109696651B (zh) 一种基于m估计的低快拍数下波达方向估计方法
CN112763972B (zh) 基于稀疏表示的双平行线阵二维doa估计方法及计算设备
CN113438738A (zh) 定位参数估计方法、装置、计算机设备和存储介质
Liu et al. Subspace techniques for multidimensional model order selection in colored noise
CN115826004B (zh) 一种基于二维角度及时差联合的三星协同直接定位方法
CN110286352B (zh) 一种基于秩损的非迭代混合信号源定位方法
CN108594165B (zh) 一种基于期望最大化算法的窄带信号波达方向估计方法
CN113821907B (zh) 针对大型平面天线阵列系统的幅度和相位自动校准方法
CN112327244B (zh) 一种基于l型阵列的二维非相干分布式目标参数估计方法
CN112051542B (zh) 基于多快照差分的sair rfi源地理定位方法和系统
CN112579972A (zh) 方向性电磁耦合效应下空域信息联合估计方法
CN112333629A (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