CN107315162B - 基于内插变换和波束形成的远场相干信号doa估计方法 - Google Patents

基于内插变换和波束形成的远场相干信号doa估计方法 Download PDF

Info

Publication number
CN107315162B
CN107315162B CN201710614333.7A CN201710614333A CN107315162B CN 107315162 B CN107315162 B CN 107315162B CN 201710614333 A CN201710614333 A CN 201710614333A CN 107315162 B CN107315162 B CN 107315162B
Authority
CN
China
Prior art keywords
array
covariance matrix
matrix
far
interpolation
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
CN201710614333.7A
Other languages
English (en)
Other versions
CN107315162A (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 CN201710614333.7A priority Critical patent/CN107315162B/zh
Publication of CN107315162A publication Critical patent/CN107315162A/zh
Application granted granted Critical
Publication of CN107315162B publication Critical patent/CN107315162B/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
    • G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及一种基于内插变换和波束形成的远场相干信号DOA估计方法，先使用内插矩阵
Figure DDA0002277178520000011
将非均匀线阵协方差矩阵
Figure DDA0002277178520000012
转化为虚拟阵列的协方差矩阵
Figure DDA0002277178520000013
将虚拟阵列的协方差矩阵
Figure DDA0002277178520000018
上的数据进行噪声预白化得到
Figure DDA0002277178520000014
再对
Figure DDA0002277178520000015
进行空间平滑处理来解相干，获得解相干后的数据协方差矩阵
Figure DDA0002277178520000016
构建代价函数对解相干后的数据协方差矩阵
Figure DDA0002277178520000017
进行处理，得到远场相干信号DOA的估计值。本发明在保证精度的前提下避免了计算复杂的特征分解和频谱搜索等操作，计算复杂度低，方法简单有效；同时，本发明将适用范围从均匀线阵推广到了任意线性阵列，从非相干信源推广到了相干信源，应用范围更广。

Description

基于内插变换和波束形成的远场相干信号DOA估计方法
技术领域
本发明属于阵列信号处理技术领域，具体涉及一种基于内插变换和波束形成的远场相干信号DOA估计方法。
背景技术
信源定位在雷达、声纳、无线通信等领域具有重要的应用，在实际应用中，由于信号的多径传播或军事上的智能干扰，往往会造成接收信号相干。同时，接收阵列的几何结构通常是由物理设施所决定，很难到达理想的均匀阵列。目前已经有很多用于解决相干信号源定位问题的方法，如基于子空间的估计方法和波束形成技术等。但是这些算法都无法直接用于非均匀阵列，并且还存在一些其他的缺陷和不足。以多重信号分类(MUSIC)为代表的基于子空间的算法使用了计算复杂的特征分解操作；而以Capon波束形成器为代表的波束形成技术估计精度有限且受信噪比(SNR)和快拍数影响很大，在低信噪比或小快拍数时算法性能急剧衰退。
发明内容
本发明的目的在于克服现有技术中存在的缺陷，提供一种可用于非均匀线阵的基于内插变换和波束形成的远场相干信号DOA估计方法。
为达到上述目的，本发明采用了以下技术方案。
包括以下步骤：
步骤一：使用内插矩阵
Figure GDA0002277178510000011
将非均匀线阵协方差矩阵
Figure GDA0002277178510000012
转化为虚拟阵列的协方差矩阵
Figure GDA0002277178510000013
将虚拟阵列的协方差矩阵
Figure GDA0002277178510000014
上的数据进行噪声预白化得到
Figure GDA0002277178510000015
步骤二：对
Figure GDA0002277178510000016
进行空间平滑处理来解相干，获得解相干后的数据协方差矩阵
Figure GDA0002277178510000017
步骤三：构建代价函数对解相干后的数据协方差矩阵
Figure GDA0002277178510000018
进行处理，得到远场相干信号DOA的估计值
Figure GDA0002277178510000021
其中，远场相干信号为入射到非均匀线阵上的p个相干窄带信号
Figure GDA0002277178510000022
非均匀线阵包含M个全向传感器阵元，M的取值范围为M≥1.5p，DOA的估计值
Figure GDA0002277178510000023
为远场相干信号的方位信息，其中
Figure GDA0002277178510000024
表示第i个远场相干信号的波达方向角，所述波达方向角为第i个远场相干信号相对于y轴的逆时针夹角。
进一步地，步骤一具体包括：
101、根据非均匀线阵接收的数据求得阵列协方差矩阵R的估计值
Figure GDA0002277178510000025
Figure GDA0002277178510000026
其中，N表示采样数，y(n)表示阵列输出信号，(·)H表示共轭转置；
102、使用内插变换得到虚拟阵列的协方差矩阵
Figure GDA0002277178510000027
Figure GDA0002277178510000028
其中，
Figure GDA0002277178510000029
表示内插矩阵；
103、将阵列协方差矩阵的估计值
Figure GDA00022771785100000210
划分为如下形式：
Figure GDA00022771785100000211
104、计算噪声方差估计值
Figure GDA00022771785100000212
Figure GDA00022771785100000213
其中
Figure GDA00022771785100000214
tr{·}表示求矩阵的迹，
Figure GDA00022771785100000215
表示广义逆；
105、噪声预白化后的虚拟阵列协方差矩阵
Figure GDA00022771785100000216
IM表示M×M的单位矩阵。
进一步地，内插矩阵
Figure GDA00022771785100000217
的计算方法为：
1021)将空间区域划分为K个子区间，在每个子区间上以Δθ为间隔定义一组角集：
Figure GDA0002277178510000031
1022)分别计算每组角集上真实阵列响应矩阵：
Figure GDA0002277178510000032
1023)设计虚拟阵列结构为均匀线性阵列，并计算其在每组角集上的响应矩阵：
Figure GDA0002277178510000033
1024)在每个子区间上，使用最小二乘法计算内插矩阵
Figure GDA0002277178510000034
使其满足：
Figure GDA0002277178510000035
由此得到每一个子区间上的内插矩阵
Figure GDA0002277178510000036
使用内插矩阵
Figure GDA0002277178510000037
进行统一表示。
进一步地，步骤二具体包括：
201)将虚拟阵列分别划分为L个部分重叠的前向子阵列和L个部分重叠的后向子阵列，每个子阵列有M0个阵元，计算第l个前向子阵列的协方差矩阵
Figure GDA0002277178510000038
和后向子阵列的协方差矩阵
Figure GDA0002277178510000039
Figure GDA00022771785100000310
Figure GDA00022771785100000311
其中(·)*表示共轭，
Figure GDA00022771785100000312
表示M0×M0的单位矩阵，A1表示虚拟阵列方向矩阵
Figure GDA00022771785100000313
的前M0行，D定义为对角矩阵：
Figure GDA00022771785100000314
利用噪声预白化后的虚拟阵列协方差矩阵
Figure GDA00022771785100000315
的元素求出上述前向子阵列协方差矩阵的估计值
Figure GDA00022771785100000316
和后项子阵列协方差矩阵的估计值
Figure GDA00022771785100000317
202)计算前向空间平滑协方差矩阵
Figure GDA00022771785100000318
和后向空间平滑协方差矩阵
Figure GDA00022771785100000319
Figure GDA00022771785100000320
Figure GDA0002277178510000041
203)计算解相干后的数据协方差矩阵：
Figure GDA0002277178510000042
进一步地，步骤三中构造代价函数具体包括以下步骤：
301、设定修正Capon波束形成器的代价函数fMCB(θ)为：
Figure GDA0002277178510000043
其中m表示修正系数；
302、根据式(13)构造多项式
Figure GDA0002277178510000044
其中
Figure GDA0002277178510000045
Figure GDA0002277178510000046
λ表示远场相干信号波长，j表示单位虚数，d表示虚拟阵列的阵元间隔，j2＝-1，通过求多项式
Figure GDA0002277178510000047
的p个零相位点来估计远场信号的波达方向角。
进一步地，m取值为大于1的整数。
与现有技术相比，本发明具有以下有益的技术效果：
本发明对任意线性阵列上远场相干信号方位信息进行估计，利用内插变换将任意线性阵列转化为虚拟均匀线阵，并且在虚拟阵列上采用空间平滑技术进行解相干，通过构建代价函数，估计得到信号的波达方向角。相比于已有的远场相干信号源定位算法，本发明在保证精度的前提下避免了计算复杂的特征分解和频谱搜索等操作，计算复杂度低，方法简单有效；同时，本发明将适用范围从均匀线阵推广到了任意线性阵列，从非相干信源推广到了相干信源，适用于更一般的阵列，应用范围更广。
进一步地，本发明采用了改进的Capon波束形成器，通过提升协方差矩阵的阶数来修正Capon波束形成器，估计得到信号的波达方向角，提高了估计精度，并且有效解决了传统Capon波束形成器在低信噪比或小快拍数时算法性能急剧衰退的问题。
附图说明
图1为阵列结构图。
图2为非均匀阵列的阵元设置图。
图3为波达方向角的估计均方根误差随信噪比变化曲线：快拍数设为500；带“△”实线：本发明方法(取修正系数m＝2)；带“*”实线：本发明方法(取修正系数m＝3)；带“□”实线：本发明方法(取修正系数m＝4)；带“○”虚线：经过内插变换和前后向平滑的传统Capon波束形成器；点划线：经过内插变换和前后向平滑后的MUSIC算法；虚线：克拉美罗下界(CRB)。
图4为波达方向角的估计均方根误差随快拍数变化曲线：信噪比设为5dB；带“△”实线：本发明方法(取修正系数m＝2)；带“*”实线：本发明方法(取修正系数m＝3)；带“□”实线：本发明方法(取修正系数m＝4)；带“○”虚线：经过内插变换和前后向平滑的传统Capon波束形成器；点划线：经过内插变换和前后向平滑后的MUSIC算法；虚线：克拉美罗下界(CRB)。
具体实施方式
下面结合附图和实施例对本发明做详细描述。
下文中，对于任意变量a，
Figure GDA0002277178510000051
表示该变量a的估计值。
DOA(Direction-Of-Arrival)：指波达方向。
本发明中将非均匀线阵内插为虚拟的均匀线阵，利用虚拟阵列上噪声预白化后的数据进行空间平滑处理来解相干，利用修正Capon波束形成器对解相干后的数据进行处理，从而估计出入射信号的波达方向角。入射信号为来自远场的p个相干窄带信号
Figure GDA0002277178510000052
即远场相干信号；
非均匀线阵包含M个全向传感器阵元，M的取值范围为M≥1.5p，角度
Figure GDA0002277178510000053
为入射信号的方位信息，其中θi表示第i个入射信号的波达方向角，波达方向角为第i个入射信号相对于y轴的逆时针夹角。
一种任意线性阵列上基于内插变换和波束形成的远场相干信号波达方向估计方法，具体实现步骤概括如下：
1)计算阵列协方差矩阵R的估计值
Figure GDA0002277178510000061
2)最小二乘法计算内插矩阵
Figure GDA0002277178510000062
3)由内插矩阵
Figure GDA0002277178510000063
和真实阵列的接收数据协方差矩阵估计值
Figure GDA0002277178510000064
计算虚拟均匀线阵上接收数据的协方差矩阵
Figure GDA0002277178510000065
4)对
Figure GDA0002277178510000066
进行噪声预白化得到
Figure GDA0002277178510000067
5)使用前后向平滑预处理
Figure GDA0002277178510000068
计算解相干后的数据协方差矩阵
Figure GDA0002277178510000069
Figure GDA00022771785100000610
6)由解相干后的数据协方差矩阵
Figure GDA00022771785100000611
构造代价函数，通过求解优化问题得到信号波达方向角的估计值
Figure GDA00022771785100000612
其中，代价函数：
Figure GDA00022771785100000613
其中(·)H表示共轭转置，m表示修正系数；
根据式(2)构造多项式
Figure GDA00022771785100000614
其中
Figure GDA00022771785100000615
Figure GDA00022771785100000620
λ表示入射信号波长，j表示单位虚数，d表上虚拟均匀阵列的阵元间隔，j2＝-1，通过求多项式
Figure GDA00022771785100000616
的p个零相位点来估计远场信号的波达方向角。
下面进行具体描述。
p个相干窄带信号
Figure GDA00022771785100000617
入射到非均匀线阵上，该非均匀线阵包含M个全向传感器阵元，M的取值范围为M≥1.5p，角度
Figure GDA00022771785100000618
为入射信号的方位信息，其中θi表示第i个入射信号相对于y轴的逆时针夹角(波达方向角)。参见图1。
令非均匀线阵的第一个阵元为参考阵元，阵列输出信号为
y(n)＝A(θ)s(n)+w(n) (1)
其中，A是阵列响应矩阵，
Figure GDA00022771785100000619
a是导向矢量，定义为
Figure GDA0002277178510000071
(·)T表示转置，λ是入射信号波长，di是第i个阵元与参考阵元的距离，j表示单位虚数，j2＝-1。假设入射信号全相干，则所有信号都可表示为一个复增益与第一个信号相乘的形式，即：si(n)＝βis1(n)，其中βi表示第i个信号的复衰减系数。
接收数据(即阵列输出信号)的阵列协方差矩阵为
Figure GDA0002277178510000072
其中Rs表示信号协方差，(·)H表示矩阵共轭转置，IM表示M×M的单位矩阵。
步骤1)计算阵列协方差矩阵估计值
Figure GDA0002277178510000073
的具体方法为：
根据真实阵列(非均匀线阵)接收数据求得阵列协方差矩阵R的估计值
Figure GDA0002277178510000074
Figure GDA0002277178510000075
其中，N表示采样数，y(n)表示阵列输出信号，(·)H表示共轭转置。
步骤2)计算内插矩阵
Figure GDA0002277178510000076
具体方法为：
a、将空间区域划分为个K子区间，在每个子区间上以Δθ为间隔定义一组角集：
Figure GDA0002277178510000077
b、分别计算每组角集上真实阵列响应矩阵：
Figure GDA0002277178510000078
c、设计虚拟阵列结构为均匀线性阵列，并计算其在每组角集上的响应矩阵：
Figure GDA0002277178510000079
d、在每个子区间上，使用最小二乘法计算内插矩阵
Figure GDA00022771785100000710
使其满足：
Figure GDA00022771785100000711
经过以上步骤，得到每一个子区间上的内插矩阵
Figure GDA00022771785100000712
使用内插矩阵
Figure GDA00022771785100000713
来代表这些矩阵，即此后的所有操作都是在每个子区间上使用对应的内插矩阵
Figure GDA00022771785100000714
分别进行的，使用
Figure GDA00022771785100000715
来统一表示。
步骤3)计算虚拟均匀线阵上接收数据的协方差矩阵
Figure GDA0002277178510000081
具体方法为：
Figure GDA0002277178510000082
其中
Figure GDA0002277178510000083
表示步骤1)计算的真实阵列输出信号协方差矩阵的估计值，
Figure GDA0002277178510000084
表示步骤2)计算的内插矩阵。
进一步展开(8)式：
Figure GDA0002277178510000085
说明本步骤得到的噪声预白化后的虚拟阵列协方差矩阵
Figure GDA0002277178510000086
中噪声不再是白噪声，因此需要后续步骤4)进行预白化。
步骤4)噪声预白化的具体方法为：
a、将真实阵列协方差矩阵估计值
Figure GDA0002277178510000087
划分为如下形式
Figure GDA0002277178510000088
b、计算噪声方差估计值
Figure GDA0002277178510000089
Figure GDA00022771785100000810
其中
Figure GDA00022771785100000811
tr{·}表示求矩阵的迹，
Figure GDA00022771785100000812
表示广义逆。
c、计算噪声预白化后的虚拟阵列协方差矩阵
Figure GDA00022771785100000813
步骤5)计算解相干后的数据协方差矩阵
Figure GDA00022771785100000814
的具体方法为：
a、将虚拟阵列分别划分为L个部分重叠的前向子阵列和L个部分重叠的后向子阵列，每个子阵列有M0个阵元，计算第l个前向子阵列的协方差矩阵
Figure GDA00022771785100000815
和后向子阵列的协方差矩阵
Figure GDA00022771785100000816
Figure GDA00022771785100000817
Figure GDA0002277178510000091
其中(·)*表示共轭，
Figure GDA0002277178510000092
表示M0×M0的单位矩阵，A1表示虚拟阵列方向矩阵
Figure GDA0002277178510000093
的前M0行，D定义为对角矩阵：
Figure GDA0002277178510000094
利用虚拟阵列的协方差矩阵
Figure GDA0002277178510000095
的元素可以求出上述前向和后项子阵列协方差矩阵的估计值
Figure GDA0002277178510000096
b、计算前向空间平滑协方差矩阵估计值
Figure GDA0002277178510000097
和后向空间平滑协方差矩阵估计值
Figure GDA0002277178510000098
Figure GDA0002277178510000099
Figure GDA00022771785100000910
c、综合前向平滑协方差矩阵和后向平滑协方差矩阵，得到前后向平滑协方差矩阵
Figure GDA00022771785100000911
Figure GDA00022771785100000912
步骤6)得到信号波达方向角的估计值
Figure GDA00022771785100000913
的具体方法为：
a、使用前后向平滑协方差矩阵
Figure GDA00022771785100000914
构造代价函数
对解相干后的数据协方差矩阵
Figure GDA00022771785100000915
进行特征值分解(EVD)可得：
Figure GDA00022771785100000916
其中，∑s＝[λ1,λ2,…,λp]为前p个大特征值组成的对角阵，
Figure GDA00022771785100000917
为剩余的小特征值组成的对角阵，并且
Figure GDA00022771785100000918
Us＝[u1,u2,…,up]称为信号子空间，
Figure GDA00022771785100000919
称为噪声子空间，λi，ui表示第i个特征值和其对应的特征向量；
以估计精度高而闻名的MUSIC算法的代价函数fMUSIC(θ)为：
Figure GDA00022771785100000920
其中，fi(θ)＝|aH(θ)ui|2。
本发明中提出的修正Capon波束形成器的代价函数fMCB(θ)为：
Figure GDA0002277178510000101
其中m表示修正系数；
将式(16)代入fMCB(θ)可得：
Figure GDA0002277178510000102
式(19)结合式(17)可得：
Figure GDA0002277178510000103
其中
Figure GDA0002277178510000104
由式(20)可见，本发明使用的代价函数由两部分构成，第一部分与信号子空间的信息相关，第二部分对应MUSIC算法的代价函数。当m增大时，第一部分信号子空间部分衰减，整个代价函数接近MUSIC算法的代价函数，因而估计精度上升。
选取m为大于1的整数，根据不同的m值可以构造不同的代价函数，m取值越大，估计精度越高，同时计算量变大，一般优选3或4。当m＝1时，该代价函数退化为标准Capon波束形成器的代价函数，即：
Figure GDA0002277178510000105
此时，也可以进行DOA的估算，但估算精度不高。
b、根据式(18)构造多项式
Figure GDA0002277178510000106
其中
Figure GDA0002277178510000107
Figure GDA0002277178510000108
λ表示入射信号波长，j表示单位虚数，d表示虚拟均匀阵列的阵元间隔，j2＝-1，通过求多项式
Figure GDA0002277178510000109
的p个零相位点来估计远场信号的波达方向角。
下面通过以下不同情形对上述方法的效果进行说明：
空间有两个波达方向未知的相干入射信号，其入射方向分别为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λ,]，参见图2。
仿真中加入经过同样的内插变换和前后向平滑预处理的传统Capon波束形成器和MUSIC算法作为对比，同时给出了克拉美罗下界(CRB)。每一个仿真结果都是经由1000次独立重复实验得到的。
由图3说明，本发明修正系数m取值分别为2、3和4时，随着修正参数m增大，本发明方法估计精度上升，当m＝4时本发明方法估计精度已经远高于经过内插变换和前后向平滑的传统Capon波束形成器(m＝1)，并且十分接近经过内插变换和前后向平滑后的高精度的MUSIC算法。同时本发明方法因为避免了MUSIC算法的特征分解操作，因而计算复杂度低于MUSIC算法。另外，随着信噪比的增大，各算法估计性能都有所提升。
由图4说明随着信噪比的增大，各算法估计性能都有所提升。同时和图3一致，本发明方法估计精度远高于经过内插变换和前后向平滑的传统Capon波束形成器，并且十分接近高精度的MUSIC算法。
本发明中，首先，将任意线性阵列进行内插变换得到虚拟的均匀线性阵列；然后，采用空间平滑技术在虚拟阵列上进行解相干；最后通过提升协方差矩阵的阶数来修正Capon波束形成器，估计得到信号的波达方向角。本发明从两方面降低了算法的计算复杂度，第一，构造出修正Capon波束形成器，在保证精度的前提下避免了传统子空间算法中计算复杂的特征分解操作；第二，将非均匀阵列内插为均匀阵列，进而使用求根算法估计波达方向角，避免了计算复杂的频谱搜素。此外，本发明将适用范围从均匀线阵推广到了任意线性阵列，从非相干信源推广到了相干信源。

Claims (5)

1.一种基于内插变换和波束形成的远场相干信号DOA估计方法，其特征在于，包括以下步骤：
步骤一：使用内插矩阵
Figure FDA0002277178500000011
将非均匀线阵协方差矩阵
Figure FDA0002277178500000012
转化为虚拟阵列的协方差矩阵
Figure FDA0002277178500000013
将虚拟阵列的协方差矩阵
Figure FDA0002277178500000014
上的数据进行噪声预白化得到
Figure FDA0002277178500000015
步骤二：对
Figure FDA0002277178500000016
进行空间平滑处理来解相干，获得解相干后的数据协方差矩阵
Figure FDA0002277178500000017
步骤三：构建代价函数对解相干后的数据协方差矩阵
Figure FDA0002277178500000018
进行处理，得到远场相干信号DOA的估计值
Figure FDA0002277178500000019
其中，远场相干信号为入射到非均匀线阵上的p个相干窄带信号
Figure FDA00022771785000000110
非均匀线阵包含M个全向传感器阵元，M的取值范围为M≥1.5p，DOA的估计值
Figure FDA00022771785000000111
为远场相干信号的方位信息，其中
Figure FDA00022771785000000112
表示第i个远场相干信号的波达方向角，所述波达方向角为第i个远场相干信号相对于y轴的逆时针夹角；步骤三中构造代价函数具体包括以下步骤：
301、设定修正Capon波束形成器的代价函数fMCB(θ)为：
Figure FDA00022771785000000113
其中m表示修正系数；
302、根据式(13)构造多项式
Figure FDA00022771785000000114
其中
Figure FDA00022771785000000115
Figure FDA00022771785000000116
λ表示远场相干信号波长，j表示单位虚数，d表示虚拟阵列的阵元间隔，j2＝-1，通过求多项式
Figure FDA00022771785000000117
的p个零相位点来估计远场信号的波达方向角。
2.根据权利要求1所述的一种基于内插变换和波束形成的远场相干信号DOA估计方法，其特征在于，步骤一具体包括：
101、根据非均匀线阵接收的数据求得阵列协方差矩阵R的估计值
Figure FDA00022771785000000118
Figure FDA00022771785000000119
其中，N表示采样数，y(n)表示阵列输出信号，(·)H表示共轭转置；
102、使用内插变换得到虚拟阵列的协方差矩阵
Figure FDA0002277178500000021
Figure FDA0002277178500000022
其中，
Figure FDA0002277178500000023
表示内插矩阵；
103、将阵列协方差矩阵的估计值
Figure FDA0002277178500000024
划分为如下形式：
Figure FDA0002277178500000025
104、计算噪声方差估计值
Figure FDA0002277178500000026
Figure FDA0002277178500000027
其中
Figure FDA0002277178500000028
tr{·}表示求矩阵的迹，
Figure FDA0002277178500000029
表示广义逆；
105、噪声预白化后的虚拟阵列协方差矩阵
Figure FDA00022771785000000210
IM表示M×M的单位矩阵。
3.根据权利要求2所述的一种基于内插变换和波束形成的远场相干信号DOA估计方法，其特征在于，最优内插矩阵
Figure FDA00022771785000000211
的计算方法为：
1021)将空间区域划分为K个子区间，在每个子区间上以Δθ为间隔定义一组角集：
Figure FDA00022771785000000212
1022)分别计算每组角集上真实阵列响应矩阵：
Figure FDA00022771785000000213
1023)设计虚拟阵列结构为均匀线性阵列，并计算其在每组角集上的响应矩阵：
Figure FDA00022771785000000214
1024)在每个子区间上，使用最小二乘法计算内插矩阵
Figure FDA00022771785000000215
使其满足：
Figure FDA00022771785000000216
由此得到每一个子区间上的内插矩阵
Figure FDA00022771785000000217
使用内插矩阵
Figure FDA00022771785000000218
进行统一表示。
4.根据权利要求1所述的一种基于内插变换和波束形成的远场相干信号DOA估计方法，其特征在于，步骤二具体包括：
201)将虚拟阵列分别划分为L个部分重叠的前向子阵列和L个部分重叠的后向子阵列，每个子阵列有M0个阵元，计算第l个前向子阵列的协方差矩阵
Figure FDA0002277178500000031
和后向子阵列的协方差矩阵
Figure FDA0002277178500000032
Figure FDA0002277178500000033
Figure FDA0002277178500000034
其中(·)*表示共轭，
Figure FDA0002277178500000035
表示M0×M0的单位矩阵，A1表示虚拟阵列方向矩阵
Figure FDA0002277178500000036
的前M0行，D定义为对角矩阵：
Figure FDA0002277178500000037
利用虚拟阵列协方差矩阵
Figure FDA0002277178500000038
的元素求出上述前向子阵列协方差矩阵的估计值
Figure FDA0002277178500000039
和后项子阵列协方差矩阵的估计值
Figure FDA00022771785000000310
202)计算前向空间平滑协方差矩阵
Figure FDA00022771785000000311
和后向空间平滑协方差矩阵
Figure FDA00022771785000000312
Figure FDA00022771785000000313
Figure FDA00022771785000000314
203)计算解相干后的数据协方差矩阵：
Figure FDA00022771785000000315
5.根据权利要求1所述的一种基于内插变换和波束形成的远场相干信号DOA估计方法，其特征在于，m取值为大于1的整数。
CN201710614333.7A 2017-07-25 2017-07-25 基于内插变换和波束形成的远场相干信号doa估计方法 Active CN107315162B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710614333.7A CN107315162B (zh) 2017-07-25 2017-07-25 基于内插变换和波束形成的远场相干信号doa估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710614333.7A CN107315162B (zh) 2017-07-25 2017-07-25 基于内插变换和波束形成的远场相干信号doa估计方法

Publications (2)

Publication Number Publication Date
CN107315162A CN107315162A (zh) 2017-11-03
CN107315162B true CN107315162B (zh) 2020-05-22

Family

ID=60178650

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710614333.7A Active CN107315162B (zh) 2017-07-25 2017-07-25 基于内插变换和波束形成的远场相干信号doa估计方法

Country Status (1)

Country Link
CN (1) CN107315162B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107870315B (zh) * 2017-11-06 2021-07-30 重庆邮电大学 一种利用迭代相位补偿技术估计任意阵列波达方向方法
CN108680891B (zh) * 2018-01-05 2022-02-22 大连大学 非均匀噪声条件下考虑互耦效应的doa估计方法
CN108459301B (zh) * 2018-01-18 2021-10-08 西安电子科技大学 一种基于异构阵的mimo雷达波形设计方法
CN110197112B (zh) * 2019-03-19 2023-07-18 上海大学 一种基于协方差修正的波束域Root-MUSIC方法
CN110031794B (zh) * 2019-04-16 2021-11-09 中国人民解放军国防科技大学 一种基于差分共性阵重构的相干信源doa估计方法
CN110531310B (zh) * 2019-07-25 2021-07-13 西安交通大学 基于子空间和内插变换的远场相干信号波达方向估计方法
CN111190136B (zh) * 2020-01-08 2023-03-24 华南理工大学 一种基于特定频率组合信号的一维doa估计方法
CN111337872B (zh) * 2020-02-21 2023-08-01 南京航空航天大学 一种用于相干信源测向的广义doa矩阵方法
CN111736111A (zh) * 2020-07-01 2020-10-02 哈尔滨理工大学 一种基于同心均匀分布双圆阵的相干源doa的估计方法
CN113777554A (zh) * 2021-08-26 2021-12-10 南京航空航天大学 一种基于求根Capon的二维DOA估计方法
CN114879131B (zh) * 2022-03-24 2024-07-16 西安电子科技大学 稀疏线阵结合内插虚拟变换技术的无网格doa估计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103353595A (zh) * 2013-06-18 2013-10-16 西安电子科技大学 基于阵列内插压缩感知的米波雷达测高方法
CN103383450A (zh) * 2013-06-25 2013-11-06 西安电子科技大学 共形阵列雷达幅相误差校正快速实现方法
CN103954931A (zh) * 2014-04-28 2014-07-30 西安交通大学 一种远场和近场混合信号源的定位方法
CN105022268A (zh) * 2015-07-09 2015-11-04 哈尔滨工程大学 一种线性约束虚拟天线波束形成方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2917180B1 (fr) * 2007-06-08 2010-05-14 Thales Sa Procede d'estimation des angles d'arrivees de sources coherentes par une technique de lissage spatial sur un reseau de capteurs quelconque

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103353595A (zh) * 2013-06-18 2013-10-16 西安电子科技大学 基于阵列内插压缩感知的米波雷达测高方法
CN103383450A (zh) * 2013-06-25 2013-11-06 西安电子科技大学 共形阵列雷达幅相误差校正快速实现方法
CN103954931A (zh) * 2014-04-28 2014-07-30 西安交通大学 一种远场和近场混合信号源的定位方法
CN105022268A (zh) * 2015-07-09 2015-11-04 哈尔滨工程大学 一种线性约束虚拟天线波束形成方法

Also Published As

Publication number Publication date
CN107315162A (zh) 2017-11-03

Similar Documents

Publication Publication Date Title
CN107315162B (zh) 基于内插变换和波束形成的远场相干信号doa估计方法
CN110927659B (zh) 互耦条件下任意阵列流形doa估计与互耦校准方法及系统
CN106054123B (zh) 一种稀疏l阵及其二维doa估计方法
CN110113085B (zh) 一种基于协方差矩阵重构的波束形成方法及系统
CN107340512B (zh) 一种基于子阵划分的近远场混合源被动定位方法
CN103837861B (zh) 基于特征子空间的子阵级线性约束自适应波束形成方法
CN107390197B (zh) 基于特征空间的雷达自适应和差波束测角方法
CN104991236B (zh) 一种单基地mimo雷达非圆信号相干源波达方向估计方法
CN112305495B (zh) 一种基于原子范数最小的互质阵列协方差矩阵重构方法
CN107576931B (zh) 一种基于协方差低维度迭代稀疏重构的相关/相干信号波达方向估计方法
CN104035069B (zh) 基于部分矫正对称均匀线阵的窄带近场信号源定位方法
CN104166136A (zh) 一种基于干扰子空间跟踪的高效自适应单脉冲测角方法
CN108872930B (zh) 扩展孔径二维联合对角化doa估计方法
CN107302391A (zh) 基于互质阵列的自适应波束成形方法
CN109709514B (zh) 一种基于阵列模型误差估计值校准的方法
CN105158741A (zh) 基于矩阵重构的自适应抗干扰多径多波束形成方法
CN115758087A (zh) 一种基于深度学习的非理想非均匀阵列单快拍测向方法
CN107102298B (zh) 基于迭代互耦校正的雷达协方差矩阵重构波束形成方法
CN106788655B (zh) 互耦条件下未知互耦信息的干扰相干稳健波束形成方法
CN113376568B (zh) 一种基于子空间正交补偿的圆阵doa估计方法
CN110535519A (zh) 一种基于空间平滑的稳健自适应波束形成方法
CN113625220B (zh) 一种多径信号波达方向和扩散角快速估计新方法
CN110531310A (zh) 基于子空间和内插变换的远场相干信号波达方向估计方法
CN115453487A (zh) 一种相控阵雷达鲁棒波束形成方法
CN114755628A (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