CN110531310B - 基于子空间和内插变换的远场相干信号波达方向估计方法 - Google Patents

基于子空间和内插变换的远场相干信号波达方向估计方法 Download PDF

Info

Publication number
CN110531310B
CN110531310B CN201910677931.8A CN201910677931A CN110531310B CN 110531310 B CN110531310 B CN 110531310B CN 201910677931 A CN201910677931 A CN 201910677931A CN 110531310 B CN110531310 B CN 110531310B
Authority
CN
China
Prior art keywords
array
correlation
virtual
real
far
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
CN201910677931.8A
Other languages
English (en)
Other versions
CN110531310A (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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong 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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201910677931.8A priority Critical patent/CN110531310B/zh
Publication of CN110531310A publication Critical patent/CN110531310A/zh
Application granted granted Critical
Publication of CN110531310B publication Critical patent/CN110531310B/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

本发明公开基于子空间和内插变换的远场相干信号波达方向估计方法，具体的计算每个子区间上的相关矢量；计算最优内插矩阵；内插矩阵以及真实阵列的相关矢量，求解虚拟矩阵的相关矢量；然后构建真是阵列和虚拟阵列的相关矩阵；再构建线性算子和正交投影算子；最后计算多项式根，求解得到远场相干信号波达方向估计，构造线性阵列子空间，求解相关矢量，在保证精度的前提下避免了传统子空间算法中计算复杂的特征分解操作；将非均匀阵列内插为均匀阵列，进而使用求根算法估计波达方向角，避免了计算复杂的频谱搜素；此外，本发明将适用范围从均匀线阵推广到了任意线性阵列，并且当传感器数量较大或需要实时估计时仍然有效。

Description

基于子空间和内插变换的远场相干信号波达方向估计方法
技术领域
本发明属于阵列信号处理技术领域，具体涉及一种基于子空间和内插变换的远场相干信号波达方向估计方法。
背景技术
信源定位在雷达、声纳、无线通信等领域具有重要的应用，在实际应用中，由于信号的多径传播或军事上的智能干扰，往往会造成接收信号相干。近几十年，以多重信号分类法(Multiple Signal Classification，MUSIC)为代表的基于子空间DOA估计算法由于其良好的估计精度而备受青睐。然而，这些传统的子空间算法都基于协方差矩阵的特征分解，计算复杂，难以满足大规模阵列或实时估计的要求。为避免特征分解这一计算密集型操作，研究者们提出了很多新型子空间算法，其中，无特征分解的子空间算法(Subspace-BasedMethods Without Eigendecomposition，SUMWE)算法在减小计算量的同时，保证了传统子空间算法的估计精度，总体性能较优。然而，SUMWE算法在应用过程中受限于阵列的几何结构，均匀线性阵列是该算法性能优秀的必备条件。但是，实际应用中均匀线阵这一理想条件却很难满足。
发明内容
本发明的目的在于克服现有技术中存在的缺陷，提供一种可用于非均匀线阵的基于子空间和内插变换的远场相干信号波达方向估计方法，对每个子阵列利用内插变换，将任意线性阵列转化为虚拟均匀线阵，从而实现SUWWE算法的推广，降低了算法的计算复杂度。
为达到上述目的，本发明采用了以下技术方案：一种基于子空间和内插变换的远场相干信号波达方向估计方法，包括以下步骤：
S100，构建真实阵列的子阵列α和子阵列β，用真实阵列前M-1个阵元构建子阵列α，后M-1个阵元构建子阵列β，其中，所述真实阵列为非均匀线阵，所述非均匀线阵包括M 个全向传感器阵元；计算子阵列α接收数据yα(n)和最后一个阵元的接受数据yM(n)的相关矢量
Figure BDA0002143889350000021
计算子阵列β的接收数据yβ(n)和第一个阵元的接受数据y1(n)的相关矢量
Figure BDA0002143889350000022
构建虚拟阵列的子阵列
Figure BDA0002143889350000023
和子阵列
Figure BDA0002143889350000024
用虚拟阵列的前M-1个阵元构建子阵列
Figure BDA0002143889350000025
和虚拟阵列的后M-1个阵元构建子阵列
Figure BDA0002143889350000026
所述虚拟阵列为为首位阵元与真实阵列首位阵元位置相同的均匀线性阵列，虚拟阵列的阵元数与真实阵列的阵元数相等；
S200，计算S100所构建的子阵列α和子阵列β的最优内插矩阵Bα和Bβ；
S300，利用S200所得最优内插矩阵Bα和Bβ以及S100所得真实阵列的相关矢量
Figure BDA0002143889350000027
和
Figure BDA0002143889350000028
求解虚拟阵列的相关矢量
Figure BDA0002143889350000029
和
Figure BDA00021438893500000210
得到：
Figure BDA00021438893500000211
Figure BDA00021438893500000212
S400，根据S100所得相关矢量
Figure BDA00021438893500000213
和
Figure BDA00021438893500000214
以及S300所得虚拟阵列的相关矢量
Figure BDA00021438893500000215
和
Figure BDA00021438893500000216
分别构建真实阵列相关矩阵Фα和Фβ以及虚拟阵列的相关矩阵
Figure BDA00021438893500000217
和
Figure BDA00021438893500000218
S500，基于S400所得真实阵列相关矩阵Фα和Фβ虚拟阵列的相关矩阵
Figure BDA00021438893500000219
和
Figure BDA00021438893500000220
构建线性算子T和正交投影算子ΠQ；
S600，利用S500所得线性算子T和正交投影算子ΠQ构建多项式
Figure BDA00021438893500000221
并求解所述多项式的根，其中，
Figure BDA00021438893500000222
为一组关于z的矢量；由所述多项式的根求解出波达方向角的估计值
Figure BDA00021438893500000223
即远场相干信号波达方向估计。
远场相干信号为入射到非均匀线阵上的p个相干窄带信号
Figure BDA00021438893500000224
M的取值范围为M＞p，角度
Figure BDA00021438893500000225
为入射相干窄带信号的方位信息，其中θi表示第i个入射相干窄带信号的波达方向角，所述波达方向角为第i个入射相干窄带信号相对于y轴的逆时针夹角。
S100中，真实阵列上接收数据y(n)为：
y(n)＝A(θ)s(n)+ω(n) (3)
其中n＝1，2，...N，N表示采样数，A(θ)＝[a(θ1)，a(θ2)，...，a(θp)]为方向矩阵；ω(n)为噪声矢量，s(n)为远场相干信号；
子阵列α接收数据yα(n)和最后一个阵元的接受数据yM(n)的相关矢量
Figure BDA0002143889350000031
为：
Figure BDA0002143889350000032
子阵列β的接收数据yβ(n)和第一个阵元的接受数据y1(n)的相关矢量
Figure BDA0002143889350000033
为：
Figure BDA0002143889350000034
其中(.)H表示共轭转置。
S200中，内插矩阵为Bα和Bβ的计算方法为：
S201，将真实阵列空间区域划分为Ls个子空间，在每个子空间上以Δθ为间隔定义一组角集，得到Ls组角集；
Figure BDA0002143889350000035
S203，分别计算S2011所述每组角集上真实阵列的子阵列α和β的响应矩阵：
Figure BDA0002143889350000036
Figure BDA0002143889350000037
其中m＝1，2，...，Ls；
S203，设计虚拟阵列结构为首位阵元与真实阵列首位阵元位置相同的均匀线性阵列，并计算其在S2011所述每组角集上的响应矩阵：
Figure BDA0002143889350000038
Figure BDA0002143889350000039
其中m＝1，2，...，Ls；
S204，在每个子阵列上，使用最小二乘法计算内插矩阵Bα和Bβ，使其满足：
Figure BDA0002143889350000041
Figure BDA0002143889350000042
其中
Figure BDA0002143889350000043
Figure BDA0002143889350000044
虚拟阵列的相关矢量
Figure BDA0002143889350000045
和
Figure BDA0002143889350000046
Figure BDA0002143889350000047
额
Figure BDA0002143889350000048
均表示虚拟阵列相关系数，Mv为虚拟阵列阵元数，真实阵列的相关矢量
Figure BDA0002143889350000049
和
Figure BDA00021438893500000410
r1，2，r1，3，...，r1，M和 rM，1，rM，2，...，rM，M-1均表示真实阵列相关系数，M为真实阵列阵元数；
S400中利用S300所述相关矢量中的元素构造相关矩阵Фα、
Figure BDA00021438893500000411
Фβ和
Figure BDA00021438893500000412
Figure BDA00021438893500000413
Figure BDA00021438893500000414
其中，
Figure BDA00021438893500000415
Figure BDA00021438893500000416
l＝1，2，...，L-1，L＝Mv-p-1。
S500具体包括以下步骤：
S501，对S400所得相关矩阵Фα、
Figure BDA00021438893500000417
Фβ和
Figure BDA00021438893500000418
进行分块：
Figure BDA00021438893500000419
Figure BDA00021438893500000420
其中，L＝Mv-p-1；
S502，用S401所得相关矩阵构造线性算子T：
Figure BDA0002143889350000051
其中
Figure BDA0002143889350000052
S503，构造正交投影算子ΠQ：
ΠQ＝Q(QHQ)-1QH (14)
其中Q＝[TT，-IL-p-1]T，(.)T表示转置。
根据S500所得求解多项式
Figure BDA0002143889350000053
的根，其中，
Figure BDA0002143889350000054
选出 p个距离单位圆最近的根，利用计算式
Figure BDA0002143889350000055
得到各远场相干信号的波达方向。
与现有技术相比，本发明具有以下有益的技术效果：本发明对任意线性阵列上远场相干信号方位信息进行估计，对每个子阵列利用内插变换将任意线性阵列转化为虚拟均匀线阵，从而实现SUWWE算法的推广，估计得到信号的波达方向角；相比于已有的远场相干信号源定位算法，本发明估计结果接近基于特征分解的精度很高的FBSS-based Root-MUSIC 算法，同时，计算复杂度和其他避免特征分解的算法相当，显著低于需要特征分解的子空间算法；在保证估计精度的同时避免了特征分解这一计算密集型操作，将传感器阵列分为两个子阵列求解相关矢量，降低了计算复杂度，总的来说，SUMWE算法是在保证估计精度的情况下，大大减小了计算量，更加符合实际需求；
同时，本发明将适用范围从均匀线阵推广到了任意线性阵列，从非相干信源推广到了相干信源，应用范围更广。
附图说明
图1为阵列结构图。
图2为阵列子空间划分图。
图3为虚拟阵列的阵元设置图。
图4为波达方向角的估计均方根误差随信噪比变化曲线：快拍数设为600；带“○”实线：本发明方法；带“□”实线：经过内插变换和前后向平滑的传统Capon波束形成器；带“△”实线：经过内插变换、前后向平滑的改进Capon波束形成器(m＝2)；带“*”实线：经过内插变换和前后向平滑的root-MUSIC算法；虚线：克拉美罗下界(CRB)。
图5为波达方向角的估计均方根误差随快拍数变化曲线：信噪比设为5dB；带“○”实线：本发明方法；带“□”实线：经过内插变换和前后向平滑的传统Capon波束形成器；带“△”实线：经过内插变换、前后向平滑的改进Capon波束形成器(m＝2)；带“*”实线：经过内插变换和前后向平滑的root-MUSIC算法；虚线：克拉美罗下界(CRB)。
具体实施方式
下面结合附图和实施例对本发明做详细描述。
对真实阵列划分子空间，采用内插变换技术，将任意阵列通过内插矩阵转化为虚拟的均匀线阵，从而实现SUWWE算法的推广；所述远场相干信号为入射到非均匀线阵上的p个相干窄带信号
Figure BDA0002143889350000061
所述非均匀线阵包含M个全向传感器阵元，M的取值范围为M＞p，角度
Figure BDA0002143889350000062
为入射信号的方位信息，其中θi表示第i个入射信号的波达方向角，所述波达方向角为第i个入射信号相对于y轴的逆时针夹角，如图1所示。
基于子空间和内插变换的远场相干信号波达方向估计方法，包括以下步骤：
1)根据真实阵列的几何结构，划分子阵列，求出每个子阵列上的相关矢量
Figure BDA0002143889350000063
和
Figure BDA0002143889350000064
Figure BDA0002143889350000065
Figure BDA0002143889350000066
2)根据真实阵列的相关矢量
Figure BDA0002143889350000067
和
Figure BDA0002143889350000068
利用内插矩阵Bα和Bβ求解获得虚拟阵列的相关矢量
Figure BDA0002143889350000069
和
Figure BDA00021438893500000610
并基于真实阵列和虚拟阵列的相关矢量中的元素分别构建相关矩阵Фα、
Figure BDA00021438893500000611
Фβ和
Figure BDA00021438893500000612
Figure BDA0002143889350000071
Figure BDA0002143889350000072
3)构建线性算子T和正交投影算子ΠQ，求解多项式
Figure BDA0002143889350000073
最接近单位圆的p的根，解得各远场相干信号波达方向
Figure BDA0002143889350000074
Figure BDA0002143889350000075
其中(·)H表示共轭转置；
一种基于子空间和内插变换的远场相干信号波达方向估计方法，包括以下步骤：
S100，构建真实阵列的子阵列α和子阵列β，计算子阵列α接收数据yα(n)和最后一个阵元的接受数据yM(n)的相关矢量
Figure BDA0002143889350000076
计算子阵列β的接收数据yβ(n)和第一个阵元的接受数据y1(n)的相关矢量
Figure BDA0002143889350000077
构建虚拟阵列的子阵列
Figure BDA0002143889350000078
和子阵列
Figure BDA0002143889350000079
S200，计算S100所构建的子阵列α和子阵列β的最优内插矩阵Bα和Bβ；
S300，利用S200所得最优内插矩阵Bα和Bβ以及S100所得真实阵列的相关矢量
Figure BDA00021438893500000710
和
Figure BDA00021438893500000711
求解虚拟阵列的相关矢量
Figure BDA00021438893500000712
和
Figure BDA00021438893500000713
S400，根据S100所得相关矢量
Figure BDA00021438893500000714
和
Figure BDA00021438893500000715
以及S300所得虚拟阵列的相关矢量
Figure BDA00021438893500000716
和
Figure BDA00021438893500000717
分别构建真实阵列相关矩阵Фα和Фβ以及虚拟阵列的相关矩阵
Figure BDA00021438893500000718
和
Figure BDA00021438893500000719
S500，基于S400所得真实阵列相关矩阵Фα和Фβ虚拟阵列的相关矩阵
Figure BDA00021438893500000720
和
Figure BDA00021438893500000721
构建线性算子T和正交投影算子ΠQ；
S600，利用S500所得线性算子T和正交投影算子ΠQ构建多项式
Figure BDA00021438893500000722
并求解所述多项式的根，其中，
Figure BDA00021438893500000723
为一组关于z的矢量；由所述多项式的根求解出波达方向角的估计值
Figure BDA00021438893500000724
即远场相干信号波达方向估计。
下面进行具体描述。
p个相干窄带信号
Figure BDA0002143889350000081
入射到非均匀线阵上，该非均匀线阵包含M个全向传感器阵元，M的取值范围为M＞p，角度
Figure BDA0002143889350000082
为入射信号的方位信息，其中θi表示第i个入射信号相对于y轴的逆时针夹角(波达方向角)，参见图1。
将所述非均匀阵列的第一个阵元作为参考阵元，阵列输出信号为
y(n)＝A(θ)s(n)+ω(n) (1)
其中，n＝1，2，...，N，N表示采样数，A(θ)＝[a(θ1)，a(θ2)，...，a(θp)]为方向矩阵，a为导向矢量，其定义为
Figure BDA0002143889350000083
(.)T表示转置，λ是入射信号波长，di是第i个阵元与参考阵元的距离，j表示单位虚数，j2＝-1；ω(n)为噪声矢量，s(n)为远场相干信号。
S100计算每个子阵列上的相关矢量
Figure BDA0002143889350000084
和
Figure BDA0002143889350000085
的具体方法为：
S101，如图2所示，利用真实阵列前M-1个阵元构建子阵列α，后M-1个阵元构成子阵列β，子阵列α接收数据yα(n)和子阵列β的接收数据yβ(n)分别如下：
yα(n)＝[y1(n)，y2(n)，...，yM-1(n)]
yβ(n)＝[y2(n)，y3(n)，...，yM(n)] (2)
S102，计算子阵列α接收数据yα(n)和最后一个阵元的接受数据yM(n)的相关矢量
Figure BDA0002143889350000086
计算子阵列β的接收数据yβ(n)和第一个阵元的接受数据y1(n)的相关矢量
Figure BDA0002143889350000087
Figure BDA0002143889350000088
Figure BDA0002143889350000089
其中(.)H表示共轭转置；
S200最小二乘法计算最优内插矩阵Bα和Bβ具体方法为：
S201，将阵列空间区域划分为Ls个子空间，在每个子空间上以Δθ为间隔定义一组角集：
Figure BDA0002143889350000091
S202，分别计算每组角集上真实阵列的子阵列α和β的响应矩阵：
Figure BDA0002143889350000092
Figure BDA0002143889350000093
其中m＝1，2，...，Ls；
S203，如图3所示，设计虚拟阵列结构为首位阵元与真实阵列首位阵元位置相同均匀线性阵列，并计算其在每组角集上的响应矩阵：
Figure BDA0002143889350000094
Figure BDA0002143889350000095
其中m＝1，2，...，Ls；
S204，在每个子空间上，使用最小二乘法计算内插矩阵Bα和Bβ，使其满足：
Figure BDA0002143889350000096
Figure BDA0002143889350000097
其中
Figure BDA0002143889350000098
Figure BDA0002143889350000099
S300，利用S200所得内插矩阵Bα和Bβ以及S100所得真实阵列的相关矢量
Figure BDA00021438893500000910
和
Figure BDA00021438893500000911
求解虚拟阵列的相关矢量
Figure BDA00021438893500000912
和
Figure BDA00021438893500000913
Figure BDA00021438893500000914
Figure BDA00021438893500000915
得到为均匀分布的虚拟阵列的相关矢量。
S400，构建相关矩阵Фα、
Figure BDA00021438893500000916
Фβ和
Figure BDA00021438893500000917
具体过程为：
S401，真实阵列相关矢量
Figure BDA00021438893500000918
虚拟阵列相关矢量
Figure BDA0002143889350000101
和
Figure BDA0002143889350000102
Mv为虚拟阵列阵元数，分别利用虚拟阵列相关矢量和真实阵列相关矢量中元素构造相关矩阵Фα、
Figure BDA0002143889350000103
Фβ和
Figure BDA0002143889350000104
Figure BDA0002143889350000105
Figure BDA0002143889350000106
其中
Figure BDA0002143889350000107
Figure BDA0002143889350000108
l＝1，2，...，L-1，L＝Mv-p-1；
S402，对相关矩阵Фα、
Figure BDA0002143889350000109
Фβ和
Figure BDA00021438893500001010
进行分块：
Figure BDA00021438893500001011
Figure BDA00021438893500001012
S500，构造线性算子T和正交投影算子ΠQ具体过程为：
Figure BDA00021438893500001013
ΠQ＝Q(QHQ)-1QH (15)
其中
Figure BDA00021438893500001014
Q＝[TT，-IL-p-1]T，(.)T表示转置；
S600，计算多项式
Figure BDA00021438893500001015
的根，所述多项式的p个峰值对应p个远场信号的波达方向，求解波达方向角的估计值
Figure BDA00021438893500001016
具体过程为：
Figure BDA00021438893500001017
选出p个最接近单位圆的根z1，z2，，...，zp，得到各远场相干信号波达方向
Figure BDA00021438893500001018
下面通过以下不同情形对本发明所述方法的效果进行解释说明：
空间有两个波达方向未知的相干入射信号，其入射方向分别为8°和20°。非均匀线阵含有10个阵元，其构造方式为：增加一个水平偏移量Δdi至一个10阵元均匀线阵(ULA)。均匀线阵的阵元间隔为d＝λ/2，水平偏移向量为：Δd＝[0，0，0.1λ，-0.1λ，0.1λ，-0.1λ，0.1λ，-0.1λ，0.1λ，-0.1λ]T，参见图3。
仿真中加入Interpolated-FBSS-root-SCB、Interpolated-FBSS-root-MCB以及Interpolated-FBSS-root-MUSIC等方法作为对比，同时给出了CRB界；每一个仿真结果都是经由1000次独立重复实验得到的。
由图4可以看出，所有的方法都相对准确地估计了远场相干信号的DOA，得益于插值变换和阵列去相关，RMSE随着信噪比的增加而显著降低，由图可以看出本发明提出的方法在低信噪比下完全优于Interpolated-FBSS-root-SCB方法和Interpolated-FBSS-root-MCB方法，随着信噪比的增加，该方法的分辨率接近于以其高精度著称的方法，但是计算量远低于该方法。
由图5说明本发明提出的方法在高快拍数下估计效果优于Interpolated-FBSS-root-MCB 方法，其RMSE随快拍数的增加单调下降，且逐渐逼近Interpolated-FBSS-root-MUSIC方法。
综上，本发明首先将线性阵列分为两个子阵列；然后，对每个子阵列利用内插变换，将任意线性阵列转化为虚拟均匀线阵，从而实现SUWWE算法的推广；最后，通过多项式求根估计得到信号的波达方向角；本发明从两方面降低了算法的计算复杂度，第一，构造线性阵列子空间，求解相关矢量，在保证精度的前提下避免了传统子空间算法中计算复杂的特征分解操作；第二，将非均匀阵列内插为均匀阵列，进而使用求根算法估计波达方向角，避免了计算复杂的频谱搜素，此外，本发明将适用范围从均匀线阵推广到了任意线性阵列，并且当传感器数量较大或需要实时估计时仍然有效。

Claims (8)

1.基于子空间和内插变换的远场相干信号波达方向估计方法，其特征在于，包括以下步骤：
S100，构建真实阵列的子阵列α和子阵列β，用真实阵列前M-1个阵元构建子阵列α，后M-1个阵元构建子阵列β，其中，所述真实阵列为非均匀线阵，所述非均匀线阵包括M个全向传感器阵元；计算子阵列α接收数据yα(n)和最后一个阵元的接收数据yM(n)的相关矢量
Figure FDA0002944346150000011
计算子阵列β的接收数据yβ(n)和第一个阵元的接收数据y1(n)的相关矢量
Figure FDA0002944346150000012
构建虚拟阵列的子阵列
Figure FDA0002944346150000013
和子阵列
Figure FDA0002944346150000014
用虚拟阵列的前M-1个阵元构建子阵列
Figure FDA0002944346150000015
和虚拟阵列的后M-1个阵元构建子阵列
Figure FDA0002944346150000016
所述虚拟阵列为首位阵元与真实阵列首位阵元位置相同的均匀线性阵列，虚拟阵列的阵元数与真实阵列的阵元数相等；
S200，计算S100所构建的子阵列α和子阵列β的最优内插矩阵Bα和Bβ；
S300，利用S200所得最优内插矩阵Bα和Bβ以及S100所得真实阵列的相关矢量
Figure FDA0002944346150000017
和
Figure FDA0002944346150000018
求解虚拟阵列的相关矢量
Figure FDA0002944346150000019
和
Figure FDA00029443461500000110
得到：
Figure FDA00029443461500000111
Figure FDA00029443461500000112
S400，根据S100所得相关矢量
Figure FDA00029443461500000113
和
Figure FDA00029443461500000114
以及S300所得虚拟阵列的相关矢量
Figure FDA00029443461500000115
和
Figure FDA00029443461500000116
分别构建真实阵列相关矩阵Φα和Φβ以及虚拟阵列的相关矩阵
Figure FDA00029443461500000117
和
Figure FDA00029443461500000118
S500，基于S400所得真实阵列相关矩阵Φα和Φβ虚拟阵列的相关矩阵
Figure FDA00029443461500000119
和
Figure FDA00029443461500000120
构建线性算子T和正交投影算子ΠQ；
S600，利用S500所得线性算子T和正交投影算子ΠQ构建多项式
Figure FDA00029443461500000121
并求解所述多项式的根，其中，
Figure FDA00029443461500000122
为一组关于z的矢量；由所述多项式的根求解出波达方向角的估计值
Figure FDA00029443461500000123
即远场相干信号波达方向估计。
2.根据权利要求1所述的基于子空间和内插变换的远场相干信号波达方向估计方法，其特征在于，远场相干信号为入射到非均匀线阵上的p个相干窄带信号
Figure FDA0002944346150000021
3.根据权利要求2所述的基于子空间和内插变换的远场相干信号波达方向估计方法，其特征在于，M的取值范围为M＞p，角度
Figure FDA0002944346150000022
为入射相干窄带信号的方位信息，其中θi表示第i个入射相干窄带信号的波达方向角，所述波达方向角为第i个入射相干窄带信号相对于y轴的逆时针夹角。
4.根据权利要求1所述的基于子空间和内插变换的远场相干信号波达方向估计方法，其特征在于，S100中，真实阵列上接收数据y(n)为：
y(n)＝A(θ)s(n)+ω(n) (3)
其中n＝1，2，...N，N表示采样数，A(θ)＝[a(θ1)，a(θ2)，...，a(θp)]为方向矩阵；ω(n)为噪声矢量，s(n)为远场相干信号；
子阵列α接收数据yα(n)和最后一个阵元的接收数据yM(n)的相关矢量
Figure FDA0002944346150000023
为：
Figure FDA0002944346150000024
子阵列β的接收数据yβ(n)和第一个阵元的接收数据y1(n)的相关矢量
Figure FDA0002944346150000029
为：
Figure FDA0002944346150000026
其中(.)H表示共轭转置。
5.根据权利要求4所述的基于子空间和内插变换的远场相干信号波达方向估计方法，其特征在于，S200中，内插矩阵为Bα和Bβ的计算方法为：
S201，将真实阵列空间区域划分为Ls个子空间，在每个子空间上以Δθ为间隔定义一组角集，得到Ls组角集；
Figure FDA0002944346150000027
S203，分别计算S2011所述每组角集上真实阵列的子阵列α和β的响应矩阵：
Figure FDA0002944346150000028
Figure FDA0002944346150000031
其中m＝1，2，...，Ls；
S203，设计虚拟阵列结构为首位阵元与真实阵列首位阵元位置相同的均匀线性阵列，并计算其在S2011所述每组角集上的响应矩阵：
Figure FDA0002944346150000032
Figure FDA0002944346150000033
其中m＝1，2，...，Ls；
S204，在每个子阵列上，使用最小二乘法计算内插矩阵Bα和Bβ，使其满足：
Figure FDA0002944346150000034
Figure FDA0002944346150000035
其中
Figure FDA0002944346150000036
Figure FDA0002944346150000037
6.根据权利要求5所述的基于子空间和内插变换的远场相干信号波达方向估计方法，其特征在于，虚拟阵列的相关矢量
Figure FDA0002944346150000038
和
Figure FDA0002944346150000039
Figure FDA00029443461500000310
和
Figure FDA00029443461500000311
均表示虚拟阵列相关系数，Mv为虚拟阵列阵元数，真实阵列的相关矢量
Figure FDA00029443461500000312
和
Figure FDA00029443461500000313
r1，2，r1，3，…，r1，M和rM，1，rM，2，…，rM，M-1均表示真实阵列相关系数，M为真实阵列阵元数；
S400中利用S300所述相关矢量中的元素构造相关矩阵Φα、
Figure FDA00029443461500000314
Φβ和
Figure FDA00029443461500000315
Figure FDA00029443461500000316
Figure FDA00029443461500000317
其中，
Figure FDA00029443461500000318
Figure FDA00029443461500000319
l＝1，2，…，L-1，L＝Mv-p-1。
7.根据权利要求6所述的基于子空间和内插变换的远场相干信号波达方向估计方法，其特征在于，S500具体包括以下步骤：
S501，对S400所得相关矩阵Φα、
Figure FDA0002944346150000041
Φβ和
Figure FDA0002944346150000042
进行分块：
Figure FDA0002944346150000043
Figure FDA0002944346150000044
其中，L＝Mv-p-1；
S502，用S401所得相关矩阵构造线性算子T：
Figure FDA0002944346150000045
其中
Figure FDA0002944346150000046
S503，构造正交投影算子∏Q：
ΠQ＝Q(QHQ)-1QH (14)
其中Q＝[TT，-IL-p-1]T，(.)T表示转置。
8.根据权利要求1所述的基于子空间和内插变换的远场相干信号波达方向估计方法，其特征在于，根据S500所得求解多项式
Figure FDA0002944346150000048
的根，其中，
Figure FDA0002944346150000049
选出p个距离单位圆最近的根，利用计算式
Figure FDA0002944346150000047
得到各远场相干信号的波达方向。
CN201910677931.8A 2019-07-25 2019-07-25 基于子空间和内插变换的远场相干信号波达方向估计方法 Active CN110531310B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910677931.8A CN110531310B (zh) 2019-07-25 2019-07-25 基于子空间和内插变换的远场相干信号波达方向估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910677931.8A CN110531310B (zh) 2019-07-25 2019-07-25 基于子空间和内插变换的远场相干信号波达方向估计方法

Publications (2)

Publication Number Publication Date
CN110531310A CN110531310A (zh) 2019-12-03
CN110531310B true CN110531310B (zh) 2021-07-13

Family

ID=68660537

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910677931.8A Active CN110531310B (zh) 2019-07-25 2019-07-25 基于子空间和内插变换的远场相干信号波达方向估计方法

Country Status (1)

Country Link
CN (1) CN110531310B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111669697B (zh) * 2020-05-25 2021-05-18 中国科学院声学研究所 一种多通道信号的相干声与环境声提取方法及系统
CN111711918B (zh) * 2020-05-25 2021-05-18 中国科学院声学研究所 一种多通道信号的相干声与环境声提取方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008105748A1 (en) * 2007-02-26 2008-09-04 Temel Engin Tuncer Method and apparatus for the joint detection of the number of signal sources and their direction of arrivals
CN108710101A (zh) * 2018-04-10 2018-10-26 贵州理工学院 一种利用正交性获取相干信号子空间的doa估计算法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10641861B2 (en) * 2000-06-02 2020-05-05 Dennis J. Dupray Services and applications for a communications network
CN102175990B (zh) * 2011-01-27 2013-03-13 西安交通大学 基于龙伯格观测器与子空间更新的波达方向跟踪方法及装置
RU2607977C1 (ru) * 2015-06-30 2017-01-11 Александр Игоревич Колотыгин Способ создания модели объекта
CN107315162B (zh) * 2017-07-25 2020-05-22 西安交通大学 基于内插变换和波束形成的远场相干信号doa估计方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008105748A1 (en) * 2007-02-26 2008-09-04 Temel Engin Tuncer Method and apparatus for the joint detection of the number of signal sources and their direction of arrivals
CN108710101A (zh) * 2018-04-10 2018-10-26 贵州理工学院 一种利用正交性获取相干信号子空间的doa估计算法

Also Published As

Publication number Publication date
CN110531310A (zh) 2019-12-03

Similar Documents

Publication Publication Date Title
CN107315162B (zh) 基于内插变换和波束形成的远场相干信号doa估计方法
CN108508423B (zh) 基于异型阵的子阵数字和差单脉冲测角方法
CN109471086B (zh) 基于多采样快拍和集阵列信号离散傅里叶变换的互质mimo雷达波达方向估计方法
CN109765521B (zh) 一种基于子阵划分的波束域成像方法
Uttam et al. Superresolution of coherent sources in real-beam data
CN107870315B (zh) 一种利用迭代相位补偿技术估计任意阵列波达方向方法
CN110531310B (zh) 基于子空间和内插变换的远场相干信号波达方向估计方法
CN107703478B (zh) 基于互相关矩阵的扩展孔径二维doa估计方法
CN113126087B (zh) 一种星载干涉成像高度计天线
CN110687538A (zh) 一种基于近场聚焦的超波束形成方法
CN109270483B (zh) 三维阵虚拟扩展相干源二维波达方向估计方法
CN114779236A (zh) 一种基于空间平滑music改进的米波雷达低仰角测高方法
CN109541526A (zh) 一种利用矩阵变换的圆环阵方位估计方法
CN109541573A (zh) 一种弯曲水听器阵列的阵元位置校准方法
CN116299150B (zh) 一种均匀面阵中降维传播算子的二维doa估计方法
CN114563760B (zh) 一种基于sca阵型的二阶超波束形成方法、设备及介质
CN114487985A (zh) 一种基于差-和信号的波束锐化方法及系统
CN114814830A (zh) 一种基于鲁棒主成分分析降噪的米波雷达低仰角测高方法
CN114755628A (zh) 非均匀噪声下声矢量传感器阵列波达方向估计方法
CN111830458B (zh) 一种平行线阵单快拍二维测向方法
CN112698263A (zh) 一种基于正交传播算子的单基地互质mimo阵列doa估计算法
CN114152913A (zh) 一种互耦均匀线阵doa估计方法
CN114167346A (zh) 基于协方差矩阵拟合阵元扩展的doa估计方法及系统
CN111537974A (zh) 一种去耦合的空时自适应信号处理方法
CN114142900B (zh) 一种基于lcmv准则的多通道自适应波束形成方法

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