CN103353587B - 基于改进的自适应数字波束形成的干涉测角方法 - Google Patents
基于改进的自适应数字波束形成的干涉测角方法 Download PDFInfo
- Publication number
- CN103353587B CN103353587B CN201310234572.1A CN201310234572A CN103353587B CN 103353587 B CN103353587 B CN 103353587B CN 201310234572 A CN201310234572 A CN 201310234572A CN 103353587 B CN103353587 B CN 103353587B
- Authority
- CN
- China
- Prior art keywords
- submatrix
- matrix
- data
- delta
- beam forming
- 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
Links
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于改进的自适应数字波束形成的干涉测角方法,主要解决现有技术自适应数字波束形成运算量大,难以满足测角***实时性的问题。其实现步骤是:(1)将中心对称阵列划分为上下左右四个子阵;(2)计算任意一个子阵的波束形成权值矢量;(3)计算其它三个子阵的波束形成权值矢量;(4)利用求得的各子阵的波束形成权值矢量对各子阵的采样数据进行加权,并进行波束增强;(5)根据波束增强后输出数据的最大值所在的距离门,选出各子阵的加权后数据计算俯仰角与方位角。本发明在不影响测角精度的前提下,降低了自适应数字波束形成的运算量,提高了测角***的实时性,可用于目标信息提取。
Description
技术领域
本发明属于信号处理技术领域,特别涉及一种雷达信号处理中的干涉测角方法,可用于目标信息提取。
背景技术
信号处理是要把记录在某种媒体上的信号进行处理,以便抽取出有用信息的过程,针对不同的用途和方法,信号处理现已经发展为很多分支,如智能信号处理、语音信号处理、阵列信号处理等。
阵列信号处理作为信号处理领域的一个重要分支,它的应用涉及雷达、声纳、通信等多种应用领域。阵列信号处理的目的是通过对阵列接收的信号进行处理,增强所需要的有用信号,抑制无用的干扰和噪声,并提取有用的信号特征。阵列信号处理最主要的研究内容包括波达方向DOA估计和数字波束形成DBF。
在已有的波达方向估计方法中,干涉法是最简单、最容易实现的一类方法,以其精度高、速度快的特点在电子战中得到广泛应用。干涉法测角实质就是利用辐射信号在接收天线上形成的相位差来确定辐射源的方向。
数字波束形成是在原来模拟波束形成原理的基础上,引入数字信号处理方法后建立的一门新技术,它将阵列所有阵元上的信号以一定方式进行加权合并,以便检测从一个特定角度到达的信号,可以认为是形成了一个波束。自适应数字波束形成是根据环境的变化以自适应的方式对接收数据进行加权。
基于传统的自适应数字波束形成的干涉测角方法是将一个阵列分为几个子阵,每个子阵分别进行自适应数字波束形成,然后根据每个子阵的等效相位中心的相位差利用干涉法进行测角。由于每个子阵要分别进行自适应数字波束形成,要多次实现矩阵求逆或矩阵特征值分解等运算复杂度很高的算法,因而在不增加硬件成本的前提下,难以满足测角***实时性的要求。
发明内容
本发明的目的在于针对基于现有的自适应数字波束形成的干涉测角方法的不足,提出一种基于改进的自适应数字波束形成的干涉测角方法,以降低自适应数字波束形成的运算量,提高测角***的实时性。
为实现上述目的,本发明基于改进的自适应数字波束形成的干涉测角方法,包括如下步骤:
(1)将中心对称的阵列划分为上下左右四个子阵U、D、L、R,每个子阵所含阵元数相同;
(2)计算任意一个子阵的波束形成权值矢量;
(3)计算其它三个子阵的波束形成权值矢量:
(3a)根据子阵间的几何关系,分别计算其它三个子阵与步骤(2)所选的子阵之间的权值变换矩阵;
(3b)通过这三个子阵中的每个子阵对应的权值变换矩阵与步骤(2)所选的子阵的波束形成权值矢量相乘,得到这三个子阵中的每个子阵的波束形成权值矢量;
(4)分别对上下左右各子阵的采样数据进行加权,并将各子阵加权后的输出数据按相同距离门直接相加,得到波束增强后的一路数据;
(5)在波束增强后的一路单个脉冲数据中找到最大值,记录最大值对应的距离门,然后分别找到此距离门对应的波束增强之前的上下左右各子阵的加权后输出值yu、yd、yl和yr,代入如下干涉法测角公式得到测角结果,即俯仰角与方位角
其中,arcsin(·)表示反正弦函数,phase(·)表示取相角函数,表示俯仰向法线指向斜率,表示方位向法线指向斜率,Ly与Lx分别表示上子阵U与右子阵R的等效相位中心到坐标原点的距离,θB与分别表示波束指向的俯仰角和方位角。
本发明与现有技术相比,具有如下优点:
本发明由于只需要利用采样数据对中心对称的阵列的四个子阵中的一个子阵的波束形成权值矢量进行计算,然后根据子阵间的几何关系,分别得到其它三个子阵与已求子阵之间的权值变换矩阵,再由每个权值变换矩阵跟已求子阵的波束形成权值矢量直接相乘,分别得到这三个子阵的波束形成权值矢量,因而避免了矩阵求逆或者特征值分解等运算复杂度高的算法的多次实现,降低了自适应数字波束形成的运算量,提高了测角***的实时性,利于硬件实现。
附图说明
图1是本发明的实现总流程图;
图2是本发明中使用的阵列排布、阵元编号及坐标系;
图3是本发明中计算一个子阵波束形成权值矢量的子流程图。
具体实施方式
参照图1,本发明给出如下四种实施例:
每个实施例均使用如图2所示的阵列和坐标系,阵元间距为载波波长的一半;设置一个目标,俯仰角为3°、方位角θ为2°;再设置一个欺骗式干扰,俯仰角为40°、方位角为30°;波束指向的俯仰角为3°、方位角为2°;信噪比SNR=20dB,干噪比INR=30dB。
实施例1
步骤1.划分子阵。
将图2所示的中心对称阵列划分为上子阵U、下子阵D、左子阵L、右子阵R,这四个子阵共用中心位置的阵元,每个子阵均包括6个阵元,由图2所示阵元编号可知上子阵包括阵元{1,2,3,4,5,6},下子阵包括阵元{6,7,8,9,10,11},左子阵包括阵元{12,13,14,15,16,6},右子阵包括阵元{6,17,18,19,20,21}。
步骤2.计算上子阵U的波束形成权值矢量。
参照图3,本步骤的具体实现如下:
(2.1)对上子阵U进行512次快拍采样,构成采样数据矩阵Xu为:
其中,是第i次快拍采样数据,1≤i≤512,再计算采样协方差矩阵Su为:
其中,H表示共轭转置;
(2.2)对采样协方差矩阵Su进行特征值分解,得到
其中,是由采样协方差矩阵Su的6个特征矢量所构成的矩阵,Λu是由采样协方差矩阵Su的6个特征值构成的对角矩阵;
(2.3)根据干扰是欺骗式干扰,且个数为1,干噪比INR=30dB,信噪比SNR=20dB这些先验信息,在采样协方差矩阵Su的6个特征矢量中,选出构成干扰子空间的特征矢量为采样协方差矩阵Su的6个特征值中最大值对应的特征矢量,假设为
(2.4)根据选出的构成干扰子空间的特征矢量计算干扰子空间矩阵Ju为:
(2.5)根据干扰子空间矩阵Ju,计算干扰子空间的正交补空间矩阵为:
其中I为6维的单位矩阵;
(2.6)设上子阵U的阵元{1,2,3,4,5,6}的坐标分别为p1,p2,p3,p4,p5,p6,计算上子阵U的导向矢量au为:
其中,T表示转置,j为虚数单位,k为波束指向对应的波数矢量;
(2.7)根据干扰子空间的正交补空间矩阵和导向矢量au,计算上子阵U的波束形成权值矢量wu为:
步骤3.分别计算下子阵D、左子阵L、右子阵R的波束形成权值矢量。
(3.1)设下子阵D的阵元{6,7,8,9,10,11}的坐标分别为p6,p7,p8,p9,p10,p11,左子阵L的阵元{12,13,14,15,16,6}的坐标分别为p12,p13,p14,p15,p16,p6,右子阵R的阵元{6,17,18,19,20,21}的坐标分别为p6,p17,p18,p19,p20,p21;
(3.2)计算下子阵D的波束形成权值矢量wd:
(3.2.1)计算下子阵D的每个阵元与上子阵U的对应阵元之间的坐标差:
(3.2.2)计算下子阵D与上子阵U之间的权值变换矩阵Cdu为:
其中,T表示转置,j为虚数单位,k为波束指向对应的波数矢量;
(3.2.3)根据下子阵D与上子阵U之间的权值变换矩阵Cdu和上子阵U的波束形成权值矢量wu,计算下子阵D的波束形成权值矢量wd为:
wd=Cdu·wu; [10]
(3.3)计算左子阵L的波束形成权值矢量wl:
(3.3.1)计算左子阵L的每个阵元与上子阵U的对应阵元之间的坐标差:
(3.3.2)计算左子阵L与上子阵U之间的权值变换矩阵Clu为:
(3.3.3)根据左子阵L与上子阵U之间的权值变换矩阵Clu和上子阵U的波束形成权值矢量wu,计算左子阵L的波束形成权值矢量wl为:
wl=Clu·wu; [13]
(3.4)计算右子阵R的波束形成权值矢量wr:
(3.4.1)计算右子阵R的每个阵元与上子阵U的对应阵元之间的坐标差:
(3.4.2)计算右子阵R与上子阵U之间的权值变换矩阵Cru为:
(3.4.3)根据右子阵R与上子阵U之间的权值变换矩阵Cru和上子阵U的波束形成权值矢量wu,计算右子阵R的波束形成权值矢量wr为:
wr=Cru·wu。 [16]
步骤4.加权与波束增强。
(4.1)根据求得的各子阵的波束形成权值矢量wu,wd,wl,wr,分别对上子阵采样数据矩阵Xu、下子阵采样数据矩阵Xd、左子阵采样数据矩阵Xl、右子阵采样数据矩阵Xr进行加权,得到各子阵加权后的输出数据矢量Yu,Yd,Yl,Yr分别为:
其中,H表示共轭转置,Yu为上子阵加权后的输出数据矢量,Yd为下子阵加权后的输出数据矢量,Yl为左子阵加权后的输出数据矢量,Yr为右子阵加权后的输出数据矢量;
(4.2)将各子阵加权后的输出数据矢量Yu,Yd,Yl,Yr相加合成一路数据,得到增强后数据矢量Yen为:
Yen=Yu+Yd+Yl+Yr; [18]
假设各子阵加权后的输出数据矢量Yu,Yd,Yl,Yr分别为:
则增强后数据矢量Yen为:
步骤5.干涉法测角。
(5.1)找到增强后数据矢量Yen中所有元素的最大值ymax,假设ymax是增强后数据矢量Yen的第t个元素,1≤t≤512;
(5.2)在各子阵的加权后输出数据矢量Yu,Yd,Yl,Yr中,分别找到上子阵加权后输出数据矢量Yu的第t个元素yu、下子阵加权后输出数据矢量Yd的第t个元素yd、左子阵加权后输出数据矢量Yl的第t个元素yl和右子阵加权后输出数据矢量Yr的第t个元素yr;
(5.3)根据(5.2)得到的各元素值yu、yd、yl、yr,求得俯仰角和方位角
其中,arcsin(·)表示反正弦函数,phase(·)表示取相角函数,俯仰向法线指向斜率ky0=2.5π,方位向法线指向斜率kx0=2.5π,θB与分别表示波束指向的俯仰角和方位角。
实施例2
步骤A.与实施例1的步骤1相同。
步骤B.计算下子阵D的波束形成权值矢量。
参照图3,本步骤的具体实现如下:
(B1)对下子阵D进行512次快拍采样,构成采样数据矩阵Xd为:
其中,是第i次快拍采样数据,1≤i≤512,再计算采样协方差矩阵Sd为:
其中,H表示共轭转置;
(B2)对采样协方差矩阵Sd进行特征值分解,得到
其中,是由采样协方差矩阵Sd的6个特征矢量所构成的矩阵,Λd是由采样协方差矩阵Sd的6个特征值构成的对角矩阵;
(B3)根据干扰是欺骗式干扰,且个数为1,干噪比INR=30dB,信噪比SNR=20dB这些先验信息,在采样协方差矩阵Sd的6个特征矢量中,选出构成干扰子空间的特征矢量为采样协方差矩阵Sd的6个特征值中最大值对应的特征向量,假设为
(B4)根据选出的构成干扰子空间的特征矢量计算干扰子空间矩阵Jd为:
(B5)根据干扰子空间矩阵Jd,计算干扰子空间的正交补空间矩阵为:
其中I为6维的单位矩阵;
(B6)设下子阵D的阵元{6,7,8,9,10,11}的坐标分别为p12,p13,p14,p15,p16,p6,计算下子阵D的导向矢量ad为:,
其中,T表示转置,j为虚数单位,k为波束指向对应的波数矢量;
(B7)根据干扰子空间的正交补空间矩阵和导向矢量ad,计算下子阵D的波束形成权值矢量wd为:
步骤C.分别计算上子阵U、左子阵L、右子阵R的波束形成权值矢量。
(C1)设上子阵U的阵元{1,2,3,4,5,6}的坐标分别为p1,p2,p3,p4,p5,p6,左子阵L的阵元{12,13,14,15,16,6}的坐标分别为p12,p13,p14,p15,p16,p6,右子阵R的阵元{6,17,18,19,20,21}的坐标分别为p6,p17,p18,p19,p20,p21;
(C2)计算上子阵U的波束形成权值矢量wu:
(C2-1)计算上子阵U的每个阵元与下子阵D的对应阵元之间的坐标差:
(C2-2)计算上子阵U与下子阵D之间的权值变换矩阵Cud为:
其中,T表示转置,j为虚数单位,k为波束指向对应的波数矢量;
(C2-3)根据上子阵U与下子阵D之间的权值变换矩阵Cud和下子阵D的波束形成权值矢量wd,计算上子阵U的波束形成权值矢量wu为:
wu=Cud·wd; [31]
(C3)计算左子阵L的波束形成权值矢量wl:
(C3-1)计算左子阵L的每个阵元与下子阵D的对应阵元之间的坐标差:
(C3-2)计算左子阵L与下子阵D之间的权值变换矩阵Cld为:
(C3-3)根据左子阵L与下子阵D之间的权值变换矩阵Cld和下子阵D的波束形成权值矢量wd,计算左子阵L的波束形成权值矢量wl为:
wl=Cld·wd; [34]
(C4)计算右子阵R的波束形成权值矢量wr:
(C4-1)计算右子阵R的每个阵元与下子阵D的对应阵元之间的坐标差:
(C4-2)计算右子阵R与下子阵D之间的权值变换矩阵Crd为:
(C4-3)根据右子阵R与下子阵D之间的权值变换矩阵Crd和下子阵D的波束形成权值矢量wd,计算右子阵R的波束形成权值矢量wr为:
wr=Crd·wd。 [37]
步骤D.与实施例1的步骤4相同。
步骤E.与实施例1的步骤5相同。
实施例3
步骤一.与实施例1的步骤1相同。
步骤二.计算左子阵L的波束形成权值矢量。
参照图3,本步骤的具体实现如下:
2.1)对左子阵L进行512次快拍采样,构成采样数据矩阵Xl为:
其中,是第i次快拍采样数据,1≤i≤512,再计算采样协方差矩阵Sl为:
其中,H表示共轭转置;
2.2)对采样协方差矩阵Sl进行特征值分解,得到
其中,是由采样协方差矩阵Sl的6个特征矢量所构成的矩阵,Λl是由采样协方差矩阵Sl的6个特征值构成的对角矩阵;
2.3)根据干扰是欺骗式干扰,且个数为1,干噪比INR=30dB,信噪比SNR=20dB这些先验信息,在采样协方差矩阵Sl的6个特征矢量中,选出构成干扰子空间的特征矢量为采样协方差矩阵Sl的6个特征值中最大值对应的特征向量,假设为
2.4)根据选出的构成干扰子空间的特征矢量计算干扰子空间矩阵Jl为:
2.5)根据干扰子空间矩阵Jl,计算干扰子空间的正交补空间矩阵为:
其中I为6维的单位矩阵;
2.6)设左子阵L的阵元{12,13,14,15,16,6}的坐标分别为p12,p13,p14,p15,p16,p6,计算左子阵L的导向矢量al为:
其中,T表示转置,j为虚数单位,k为波束指向对应的波数矢量;
2.7)根据干扰子空间的正交补空间矩阵和导向矢量al,计算左子阵L的波束形成权值矢量wl为:
步骤三.分别计算上子阵U、下子阵D、右子阵R的波束形成权值矢量。
3.1)设上子阵U的阵元{1,2,3,4,5,6}的坐标分别为p1,p2,p3,p4,p5,p6,下子阵D的阵元{6,7,8,9,10,11}的坐标分别为p6,p7,p8,p9,p10,p11,右子阵R的阵元{6,17,18,19,20,21}的坐标分别为p6,p17,p18,p19,p20,p21;
3.2)计算上子阵U的波束形成权值矢量wu:
3.2.1)计算上子阵U的每个阵元与左子阵L的对应阵元之间的坐标差:
3.2.2)计算上子阵U与左子阵L之间的权值变换矩阵Cul为:
其中,T表示转置,j为虚数单位,k为波束指向对应的波数矢量;
3.2.3)根据上子阵U与左子阵L之间的权值变换矩阵Cul和左子阵L的波束形成权值矢量wl,计算上子阵U的波束形成权值矢量wu为:
wu=Cul·wl; [47]
3.3)计算下子阵D的波束形成权值矢量wd:
3.3.1)计算下子阵D的每个阵元与左子阵L的对应阵元之间的坐标差:
3.3.2)计算下子阵D与左子阵L之间的权值变换矩阵Cdl为
3.3.3)根据下子阵D与左子阵L之间的权值变换矩阵Cdl和左子阵L的波束形成权值矢量wl,计算下子阵D的波束形成权值矢量wd为:
wd=Cdl·wl; [50]
3.4)计算右子阵R的波束形成权值矢量wr:
3.4.1)计算右子阵R的每个阵元与左子阵L的对应阵元之间的坐标差:
3.4.2)计算右子阵R与左子阵L之间的权值变换矩阵Crl为:
3.4.3)根据右子阵R与左子阵L之间的权值变换矩阵Crl和左子阵L的波束形成权值矢量wl,计算右子阵R的波束形成权值矢量wr为:
wr=Crl·wl。 [53]
步骤四.与实施例1的步骤4相同。
步骤五.与实施例1的步骤5相同。
实施例4
第一步.与实施例1的步骤1相同。
第二步.计算右子阵R的波束形成权值矢量。
参照图3,本步骤的具体实现如下:
(2a)对右子阵R进行512次快拍采样,构成采样数据矩阵Xr为:
其中,是第i次快拍采样数据,1≤i≤512,再计算采样协方差矩阵Sr为:
其中,H表示共轭转置;
(2b)对采样协方差矩阵Sr进行特征值分解,得到
其中,是由采样协方差矩阵Sr的6个特征矢量所构成的矩阵,Λr是由采样协方差矩阵Sr的6个特征值构成的对角矩阵;
(2c)根据干扰是欺骗式干扰,且个数为1,干噪比INR=30dB,信噪比SNR=20dB这些先验信息,在采样协方差矩阵Sr的6个特征矢量中,选出构成干扰子空间的特征矢量为采样协方差矩阵Sr的6个特征值中最大值对应的特征向量,假设为
(2d)根据选出的构成干扰子空间的特征矢量计算干扰子空间矩阵Jr为:
(2e)根据干扰子空间矩阵Jr,计算干扰子空间的正交补空间矩阵为:
其中I为6维的单位矩阵;
(2f)设右子阵R的阵元{6,17,18,19,20,21}的坐标分别为p6,p17,p18,p19,p20,p21,计算右子阵R的导向矢量ar为:
其中,T表示转置,j为虚数单位,k为波束指向对应的波数矢量;
(2g)根据干扰子空间的正交补空间矩阵和导向矢量ar,计算右子阵R波束形成权值矢量wr为:
第三步.分别计算上子阵U、下子阵D、左子阵L的波束形成权值矢量。
(3a)设上子阵U的阵元{1,2,3,4,5,6}的坐标分别为p1,p2,p3,p4,p5,p6,下子阵D的阵元{6,7,8,9,10,11}的坐标分别为p6,p7,p8,p9,p10,p11,左子阵L的阵元{12,13,14,15,16,6}的坐标分别为p12,p13,p14,p15,p16,p6;
(3b)计算上子阵U的波束形成权值矢量wu:
(3b1)计算上子阵U的每个阵元与右子阵R的对应阵元之间的坐标差:
(3b2)计算上子阵U与右子阵R之间的权值变换矩阵Cur为:
其中,T表示转置,j为虚数单位,k为波束指向对应的波数矢量;
(3b3)根据上子阵U与右子阵R之间的权值变换矩阵Cur和右子阵R波束形成权值矢量wr,计算上子阵U的波束形成权值矢量wu为:
wu=Cur·wr; [63]
(3c)计算下子阵D的波束形成权值矢量wd:
(3c1)计算下子阵D的每个阵元与右子阵R的对应阵元之间的坐标差:
(3c2)计算下子阵D与右子阵R之间的权值变换矩阵Cdr为:
(3c3)根据下子阵D与右子阵R之间的权值变换矩阵Cdr和右子阵R波束形成权值矢量wr,计算下子阵D的波束形成权值矢量wd为:
wd=Cdr·wr; [66]
(3d)计算左子阵L的波束形成权值矢量wl:
(3d1)计算左子阵L的每个阵元与右子阵R的对应阵元之间的坐标差:
(3d2)计算左子阵L与右子阵R之间的权值变换矩阵Clr为:
(3d3)根据左子阵L与右子阵R之间的权值变换矩阵Clr和右子阵R波束形成权值矢量wr,计算左子阵L的波束形成权值矢量wl为:
wl=Clr·wr。 [69]
第四步.与实施例1的步骤4相同。
第五步.与实施例1的步骤5相同。
本发明的效果可通过以下测角结果与理论分析进一步说明:
1.测角结果
将四种实施例的测角结果与角度理论值及传统方法的测角结果进行比较,比较结果如下表所示:
表1四种实施例的测角结果与角度理论值及传统方法的测角结果对比
由表1所示的比较结果可知,本发明方法是行之有效的。
2.合理性证明
假设有A,B两个子阵,每个子阵均包括M个阵元,A子阵的各阵元坐标分别为pA1,pA2,…,pAM,B子阵的各阵元坐标分别为pB1,pB2,…,pBM,则A子阵的导向矢量vA为:
其中,j为虚数单位,T表示转置,k为波数矢量,B子阵的导向矢量vB为:
从而可得
其中阵元坐标差△pAi,Bi=pBi-pAi,i=1,2,…,M,记
则有
vB=CAB·vA, [74]
其中,H表示共轭转置,I为M维的单位矩阵;
设{rA1,rA2,…,rAm}为构成A子阵的干扰子空间的一组导向矢量,并设矩阵RA=[rA1,rA2,…,rAm];对A子阵的采样协方差矩阵进行特征值分解,设A子阵的采样协方差矩阵的一组特征矢量{sA1,sA2,…,sAm}构成A子阵的干扰子空间,并设矩阵SA=[sA1,sA2,…,sAm],则有
设{rB1,rB2,…,rBm}为构成B子阵干扰子空间的一组导向矢量,并设矩阵RB=[rB1,rB2,…,rBm];对B子阵的采样协方差矩阵进行特征值分解,设B子阵的采样协方差矩阵的一组特征矢量{sB1,sB2,…,sBm}构成B子阵的干扰子空间,并设矩阵SB=[sB1,sB2,…,sBm],则有
并且
RB=CAB·RA; [78]
于是得到B子阵的波束形成权值矢量wB为:
由上式可知,如果已经求得矩阵CAB,B子阵的波束形成权值矢量可由矩阵CAB与A子阵的波束形成权值矢量直接相乘得到,CAB即为B子阵与A子阵之间的权值变换矩阵,CAB的具体形式如式[73]所示。
由上述证明过程可知,如果已经求得中心对称阵列的一个子阵的波束形成权值矢量,计算中心对称阵列的其它三个子阵的波束形成权值矢量,可以根据子阵间的几何关系,分别得到其它三个子阵与已求子阵之间的权值变换矩阵,再由每个权值变换矩阵跟已求子阵的波束形成权值矢量直接相乘,得到这三个子阵的波束形成权值矢量,即本发明方法在理论上是合理的。
3.运算量分析
本发明的运算量可由包括的乘法次数和加法次数来说明,假设本发明所使用的阵列的四个子阵的每个子阵均包括M个阵元,其运算量分析如下:
(3.1)本发明方法比传统方法增加了分别计算三个子阵对应的权值变换矩阵、三个子阵对应的权值变换矩阵分别与已求子阵的波束形成权值矢量相乘,下面分别分析增加的这两个运算过程的运算量:
(3.1a)计算一个子阵对应的权值变换矩阵包括3·M次乘法和2·M·(M-1)次加法,因此分别计算三个子阵对应的权值变换矩阵共包括3·3·M=9·M次乘法和3·2·M·(M-1)=6·M2-6·M次加法;
(3.1b)一个子阵对应的权值变换矩阵与已求子阵的波束形成权值矢量相乘包括4·M次乘法和2·M次加法,因此三个子阵对应的权值变换矩阵分别与已求子阵的波束形成权值矢量相乘共包括3·4·M=12·M次乘法和3·2·M=6·M次加法。
根据(3.1)的分析,本发明方法比传统方法增加的运算量总共是:21·M次乘法和6·M2次加法。
(3.2)本发明方法比传统方法减少了分别计算三个子阵的采样协方差矩阵、分别对三个子阵的采样协方差矩阵进行特征值分解、分别计算三个子阵的干扰子空间矩阵、分别计算三个子阵的干扰子空间的正交补空间矩阵、三个子阵的每个子阵的干扰子空间的正交补空间矩阵与这个子阵的导向矢量相乘,下面分别分析减少的这五个运算过程的运算量:
(3.2a)计算一个子阵的采样协方差矩阵包括4·N·M2+1次乘法和2·N·M2+N-1次加法,因此分别计算三个子阵的采样协方差矩阵共包括3·(4·N·M2+1)次乘法和3·(2·N·M2+N-1)次加法;
(3.2b)由于矩阵特征值分解的具体算法的不同,运算量也不完全相同,对一个子阵的采样协方差矩阵进行特征值分解的运算量数量级为О(M3),假设包括M3次乘法和M3次加法,因此分别对三个子阵的采样协方差矩阵进行特征值分解共包括3·M3次乘法和3·M3次加法;
(3.2c)计算一个子阵的干扰子空间矩阵包括4·M2次乘法和2·M2次加法,因此分别计算三个子阵的干扰子空间矩阵共包括3·4·M2=12·M2次乘法和3·2·M2=6·M2次加法;
(3.2d)计算一个子阵的干扰子空间的正交补空间矩阵包括M次加法运算,因此分别计算三个子阵的干扰子空间的正交补空间矩阵共包括3·M次加法运算;
(3.2e)一个子阵的干扰子空间的正交补空间矩阵与这个子阵的导向矢量相乘包括4·M2次乘法和2·M2+2·M·(M-1)次加法,因此三个子阵的每个子阵的干扰子空间的正交补空间矩阵与这个子阵的导向矢量相乘共包括3·4·M2=12·M2次乘法和3·(2·M2+2·M·(M-1))=12·M2-6·M次加法。
根据(3.2)的分析,本发明方法比传统方法减少的运算量总共是:3·M3+12·(N+2)·M2+3次乘法和3·M3+6·(N+3)·M2+3·(N-M-1)次加法。
综上(3.1)和(3.2)对运算量的分析,本发明方法比传统方法增加的运算量相对于减少的运算量可以忽略不计,又因为传统方法需要利用采样数据分别计算四个子阵的波束形成权值矢量,而本发明方法只需要利用采样数据计算一个子阵的波束形成权值矢量,所以本发明方法的自适应数字波束形成的运算量约为传统方法的自适应数字波束形成的运算量的1/4。由此可知,本发明方法显著降低了自适应数字波束形成的运算量。
Claims (3)
1.一种基于改进的自适应数字波束形成的干涉测角方法,包括如下步骤:
(1)将中心对称的阵列划分为上下左右四个子阵U、D、L、R,每个子阵所含阵元数相同;
(2)计算任意一个子阵的波束形成权值矢量;
(3)计算其它三个子阵的波束形成权值矢量:
(3a)根据子阵间的几何关系,分别计算其它三个子阵与步骤(2)所选的子阵之间的权值变换矩阵:
(3a1)记步骤(2)所选的子阵为子阵1,其它三个子阵分别为子阵2、3、4,并设四个子阵均包含M个阵元,子阵1的阵元坐标分别为p11,p12,…,p1M,子阵2的阵元坐标分别为p21,p22,…,p2M,子阵3的阵元坐标分别为p31,p32,…,p3M,子阵4的阵元坐标分别为p41,p42,…,p4M;
(3a2)计算子阵2的每个阵元与子阵1的对应阵元之间的坐标差:
Δp1i,2i=p2i-p1i;
其中i=1,2,…,M,再计算子阵2与子阵1之间的权值变换矩阵:
其中k为波束指向的波数矢量,T表示转置,j为虚数单位;
(3a3)计算子阵3的每个阵元与子阵1的对应阵元之间的坐标差:
Δp1i,3i=p3i-p1i;
再计算子阵3与子阵1之间的权值变换矩阵:
(3a4)计算子阵4的每个阵元与子阵1的对应阵元之间的坐标差:
Δp1i,4i=p4i-p1i;
再计算子阵4与子阵1之间的权值变换矩阵:
(3b)通过这三个子阵中的每个子阵对应的权值变换矩阵与步骤(2)所选的子阵的波束形成权值矢量相乘,得到这三个子阵中的每个子阵的波束形成权值矢量;
(4)分别对上下左右各子阵的采样数据进行加权,并将各子阵加权后的输出数据按相同距离门直接相加,得到波束增强后的一路数据;
(5)在波束增强后的一路单个脉冲数据中找到最大值,记录最大值对应的距离门,然后分别找到此距离门对应的波束增强之前的上下左右各子阵的加权后输出值yu、yd、yl和yr,代入如下干涉法测角公式得到测角结果,即俯仰角与方位角
其中,arcsin(·)表示反正弦函数,phase(·)表示取相角函数,表示俯仰向法线指向斜率,表示方位向法线指向斜率,Ly与Lx分别表示上子阵U与右子阵R的等效相位中心到坐标原点的距离,θB与分别表示波束指向的俯仰角和方位角。
2.按权利要求1所述的方法,其特征在于步骤(2)所述的计算任意一个子阵的波束形成权值矢量,按如下步骤进行:
(2a)对所选子阵进行N次快拍采样,构成采样数据矩阵X,根据该矩阵X计算采样协方差矩阵S:
S=X·XH/N,
其中,H表示共轭转置;
(2b)对采样协方差矩阵S进行特征值分解,得到
S=QΛQH,
其中,Q=[v1,v2,…,vm]是由采样协方差矩阵S的m个特征向量vi所构成的矩阵,Λ是由采样协方差矩阵S的m个特征值λi构成的对角矩阵,m是采样协方差矩阵S的维数,i=1,2,…,m;
(2c)根据干扰的类型、强度和个数这些先验信息和采样协方差矩阵S的m个特征值{λ1,λ2,…,λm}之间的大小关系,在采样协方差矩阵S的m个特征向量{v1,v2,…,vm}中选出构成干扰子空间的r个特征矢量{v1,v2,…,vr},1≤r≤m;
(2d)根据选出的构成干扰子空间的r个特征矢量{v1,v2,…,vr},计算干扰子空间矩阵J为:
(2e)根据求得的干扰子空间矩阵J,计算干扰子空间的正交补空间矩阵J⊥为:
J⊥=I-J,
其中,I是与干扰子空间矩阵J相同维数的单位矩阵;
(2f)根据所选子阵的阵元坐标,设为p1,p2,…,pm,和波束指向对应的波数矢量k,计算所选子阵的导向矢量a1为:
其中m为阵元个数,与采样协方差矩阵S的维数相同;
(2g)根据求得的干扰子空间的正交补空间矩阵J⊥和导向矢量a1,计算所选子阵的波束形成权值矢量w1为:
w1=J⊥·a1。
3.按权利要求1所述的方法,其特征在于步骤(4)所述的分别对上下左右各子阵的采样数据进行加权,是分别将各子阵求得的波束形成权值矢量取共轭转置后再与各子阵的采样数据相乘,即
其中Dataup_in,Datadown_in,Dataleft_in,Dataright_in分别是上下左右各子阵的采样数据,wup,wdown,wleft,wright分别是上下左右各子阵的波束形成权值矢量,Dataup_out,Datadown_out,Dataleft_out,Dataright_out分别是上下左右各子阵加权后的输出数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310234572.1A CN103353587B (zh) | 2013-06-13 | 2013-06-13 | 基于改进的自适应数字波束形成的干涉测角方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310234572.1A CN103353587B (zh) | 2013-06-13 | 2013-06-13 | 基于改进的自适应数字波束形成的干涉测角方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103353587A CN103353587A (zh) | 2013-10-16 |
CN103353587B true CN103353587B (zh) | 2015-01-28 |
Family
ID=49309974
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310234572.1A Active CN103353587B (zh) | 2013-06-13 | 2013-06-13 | 基于改进的自适应数字波束形成的干涉测角方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103353587B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106855622B (zh) * | 2015-12-08 | 2019-06-14 | 中国航空工业集团公司雷华电子技术研究所 | 一种子阵级相控阵雷达的测角方法 |
CN106054122B (zh) * | 2016-05-30 | 2018-09-25 | 西安电子科技大学 | 基于数字信号处理机的时域宽带信号频域闭环测向方法 |
CN107064884B (zh) * | 2017-01-05 | 2020-01-31 | 西安电子科技大学 | 基于规则重叠子阵的自适应波束形成方法 |
CN106950527B (zh) * | 2017-04-14 | 2020-01-31 | 中国电子科技集团公司第二十九研究所 | 一种多基线干涉仪测向体制下的脉冲信号分选方法 |
CN107340495A (zh) * | 2017-06-28 | 2017-11-10 | 西安电子科技大学 | 一种基于阵列雷达的目标波达方向快速估计方法 |
CN112448748A (zh) | 2019-08-30 | 2021-03-05 | 中兴通讯股份有限公司 | 一种实现波束对准的方法和装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1503487A (zh) * | 2002-11-25 | 2004-06-09 | 深圳市中兴通讯股份有限公司 | 码滤波自适应波束形成方法 |
CN1855760A (zh) * | 2005-04-28 | 2006-11-01 | 上海原动力通信科技有限公司 | 一种波束赋形方法 |
EP2096457A2 (en) * | 2008-02-29 | 2009-09-02 | Universität Karlsruhe, (Technische Hochschule) | Digital beam forming using frequency-modulated signals |
CN101644773A (zh) * | 2009-03-20 | 2010-02-10 | 中国科学院声学研究所 | 一种实时频域超分辨方位估计方法及装置 |
CN102142609A (zh) * | 2010-12-16 | 2011-08-03 | 哈尔滨工业大学 | 具有低旁瓣特性的子阵级自适应数字波束形成器 |
CN102608588A (zh) * | 2012-03-14 | 2012-07-25 | 西安电子科技大学 | 基于子带分解的宽带子阵自适应波束形成方法 |
CN102707258A (zh) * | 2012-06-05 | 2012-10-03 | 西安交通大学苏州研究院 | 基于l型传感器阵列的信号方位角和仰角联合估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6697009B2 (en) * | 2001-06-15 | 2004-02-24 | Lockheed Martin Corporation | Adaptive digital beamforming architecture for target detection and angle estimation in multiple mainlobe and sidelobe jamming |
-
2013
- 2013-06-13 CN CN201310234572.1A patent/CN103353587B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1503487A (zh) * | 2002-11-25 | 2004-06-09 | 深圳市中兴通讯股份有限公司 | 码滤波自适应波束形成方法 |
CN1855760A (zh) * | 2005-04-28 | 2006-11-01 | 上海原动力通信科技有限公司 | 一种波束赋形方法 |
EP2096457A2 (en) * | 2008-02-29 | 2009-09-02 | Universität Karlsruhe, (Technische Hochschule) | Digital beam forming using frequency-modulated signals |
CN101644773A (zh) * | 2009-03-20 | 2010-02-10 | 中国科学院声学研究所 | 一种实时频域超分辨方位估计方法及装置 |
CN102142609A (zh) * | 2010-12-16 | 2011-08-03 | 哈尔滨工业大学 | 具有低旁瓣特性的子阵级自适应数字波束形成器 |
CN102608588A (zh) * | 2012-03-14 | 2012-07-25 | 西安电子科技大学 | 基于子带分解的宽带子阵自适应波束形成方法 |
CN102707258A (zh) * | 2012-06-05 | 2012-10-03 | 西安交通大学苏州研究院 | 基于l型传感器阵列的信号方位角和仰角联合估计方法 |
Non-Patent Citations (1)
Title |
---|
一种稳健的自适应波束形成器;郭庆华等;《电子与信息学报》;20040131;第26卷(第1期);146-150 * |
Also Published As
Publication number | Publication date |
---|---|
CN103353587A (zh) | 2013-10-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103353587B (zh) | 基于改进的自适应数字波束形成的干涉测角方法 | |
CN110927659B (zh) | 互耦条件下任意阵列流形doa估计与互耦校准方法及*** | |
Liu et al. | An eigenstructure method for estimating DOA and sensor gain-phase errors | |
Jiang et al. | Mixed near-field and far-field sources localization using the uniform linear sensor array | |
CN108375763B (zh) | 一种应用于多声源环境的分频定位方法 | |
CN110045323B (zh) | 一种基于矩阵填充的互质阵稳健自适应波束形成算法 | |
CN103344940B (zh) | 低复杂度的doa估计方法及*** | |
Zhang et al. | Two-dimensional direction of arrival estimation for coprime planar arrays via polynomial root finding technique | |
Zhang et al. | A rank-reduction based 2-D DOA estimation algorithm for three parallel uniform linear arrays | |
CN109254272B (zh) | 一种共点式极化mimo雷达的两维角度估计方法 | |
CN109471063B (zh) | 基于延迟快拍的均匀线列阵高分辨波达方向估计方法 | |
CN110954861A (zh) | 一种基于增强型嵌套阵列的doa估计方法 | |
CN106405487A (zh) | 一种通用基于扩展esprit技术的空间谱估计方法 | |
CN105005038A (zh) | 一种改进的声矢量阵相干源doa估计算法 | |
CN108872930B (zh) | 扩展孔径二维联合对角化doa估计方法 | |
He et al. | Polarization difference smoothing for direction finding of coherent signals | |
Ge et al. | A low complexity algorithm for direction of arrival estimation with direction-dependent mutual coupling | |
He et al. | Two-dimensional direction finding of acoustic sources by a vector sensor array using the propagator method | |
Huang et al. | Unitary transformations for spherical harmonics MUSIC | |
CN110196417B (zh) | 基于发射能量集中的双基地mimo雷达角度估计方法 | |
CN104931937A (zh) | 基于协方差矩阵归一化的子阵级正交投影波束形成方法 | |
Zhang et al. | An efficient 2-D DOA estimation method for unfolded coprime L-shaped array | |
CN112733333A (zh) | 互质面阵中一种基于多项式求根的二维测向估计方法 | |
CN116299150B (zh) | 一种均匀面阵中降维传播算子的二维doa估计方法 | |
CN112016037A (zh) | 一种互质面阵中基于降维Capon求根的二维测向估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |