申请/专利权人：西安电子科技大学

申请日：2017-08-04

公开（公告）日：2020-10-16

公开（公告）号：CN107390197B

主分类号：G01S7/41(20060101)

分类号：G01S7/41(20060101)

优先权：

专利状态码：有效-授权

法律状态：2020.10.16#授权;2017.12.22#实质审查的生效;2017.11.24#公开

摘要：本发明公开了一种基于特征空间的雷达自适应和差波束测角方法，其主要思路为：确定均匀线阵，该均匀线阵包括M个阵元，均匀线阵的检测范围内存在J+1个信号源，J+1个信号源向均匀线阵发射J+1个入射信号，所述J+1个入射信号中包含目标信号；进而确定t时刻均匀线阵接收的M维信号；将M个阵元的均匀线阵划分为L个子阵，并根据t时刻均匀线阵接收的M维信号，得到t时刻均匀线阵的L维降维信号输出；计算得到L个子阵接收J+1个入射信号的和波束权值；计算得到L个子阵接收J+1个入射信号的差波束权值Wdiff；根据L个子阵接收J+1个入射信号的和波束权值Wsum、L个子阵接收J+1个入射信号的差波束权值Wdiff，以及t时刻均匀线阵的L维降维信号输出，计算得到目标信号的实际来波方向。

主权项：1.一种基于特征空间的雷达自适应和差波束测角方法，其特征在于，包括以下步骤：步骤1，确定均匀线阵，该均匀线阵包括M个阵元，均匀线阵的检测范围内存在J+1个信号源，J+1个信号源向均匀线阵发射J+1个入射信号，所述J+1个入射信号中包含目标信号；进而确定t时刻均匀线阵接收的M维信号；其中，t表示时间变量，M表示均匀线阵包含的阵元个数，M、J分别为大于0的正整数；所述t时刻均匀线阵接收的M维信号，其确定过程为：J+1个入射信号包括一个目标信号和J个干扰信号，且J+1个入射信号分别为窄带信号，第i个入射信号的来波方向为θi，i∈{1,2,…,J+1}，且i＝1时对应目标信号，i≠1时分别对应干扰信号；J+1个入射信号的来波方向不同，分别为θ1,θ2,…,θi,…,θJ+1；J+1个信号源向均匀线阵发射J+1个入射信号，每个阵元分别接收J+1个不相关的窄带信号，其中将第m个阵元接收的J+1个不相关的窄带信号记为smt，smt＝{sm,1t,sm,2t,…,sm,it,…,sm,J+1t}，m∈{1,2,…,M}，sm,it表示第m个阵元接收的第i个窄带信号；则t时刻均匀线阵接收的M维信号为xt，其表达式为： nt＝[n1t,n2t,…,nmt,…,nMt]T其中，xmt表示t时刻均匀线阵中第m个阵元接收的信号，所述t时刻均匀线阵中第m个阵元接收的信号为第m个阵元接收的J+1个不相关的窄带信号；sit为t时刻第i个入射信号的复包络；aθi为第i个入射信号的导向矢量，nt为t时刻均匀线阵中M个阵元的高斯白噪声，nmt表示t时刻均匀线阵中第m个阵元的高斯白噪声，上标T表示转置操作，d表示均匀线阵的阵元间距，θi表示第i个入射信号的来波方向，e表示指数函数，上标j表示虚数单位，sin为正弦函数；λi表示第i个入射信号的波长，且每一个入射信号的波长取值分别相等；t表示时间变量，M表示均匀线阵包含的阵元个数，M为大于0的正整数；步骤2，将M个阵元的均匀线阵划分为L个子阵，并根据t时刻均匀线阵接收的M维信号，得到t时刻均匀线阵的L维降维信号输出；其中，L表示将M个阵元的均匀线阵划分后包含的子阵个数，L为大于0的正整数；所述t时刻均匀线阵的L维降维信号输出，其得到过程为：将M个阵元的均匀线阵划分为L个子阵，每个子阵包含的阵元个数均为g，g为正整数，g≥1，表示向上取整；进而计算得到划分L个子阵后第i个入射信号的导向矢量为 其中，g表示每个子阵包含的阵元个数，d表示均匀线阵的阵元间距，θi表示第i个入射信号的来波方向，上标T表示转置操作，e表示指数函数，上标j表示虚数单位，sin为正弦函数；λi表示第i个入射信号的波长，且每一个入射信号的波长取值分别相等；进而得到同一来波方向的条件下第i个入射信号的导向矢量其表达式为： 进而计算得到L个子阵的L×M维降维矩阵T： 其中，θ1表示目标信号的来波方向，目标信号的来波方向θ1为目标信号的期望来波方向；e表示指数函数，上标j表示虚数单位，sin为正弦函数，λ表示每一个入射信号的波长，g表示每个子阵包含的阵元个数，g为正整数，g≥1；然后根据L个子阵的L×M维降维矩阵T，将t时刻均匀线阵接收的M维信号xt降到L维，得到t时刻均匀线阵的L维降维信号输出mt：mt＝Txt＝[m1t,m2t,…,ml't,…,mLt]T其中，ml't表示t时刻均匀线阵第l'个子阵的降维信号输出，l'＝1,…L，L表示将M个阵元的均匀线阵划分后包含的子阵个数，L为大于0的正整数；步骤3，计算得到L个子阵接收J+1个入射信号的和波束权值Wsum；其得到过程为：3a计算得到L个子阵接收J+1个入射信号的协方差矩阵R，并对其作特征值分解，得到： 其中，L个子阵接收J+1个入射信号的协方差矩阵R包含L个特征值，且λ1≥λ2≥…≥λJ+1λJ+2＝…＝λL＝σn2，表示每个子阵接收J+1个入射信号的高斯白噪声功率，且将λ1,λ2,…,λi',…,λJ+1记为J+1个大特征值，将λJ+2,λJ+3,…,λi”,…,λL记为L-J-1个小特征值，vi'表示第i'个特征值对应的特征矢量，i'＝1,…,J+1，vi”表示第i”个特征值对应的特征矢量，i”＝J+2,…,L，λi'表示第i'个特征值，i”＝J+2,J+3,…,L，λi”表示第i”个特征值，上标H表示共轭转置操作；分别令由J+1个大特征值构成的对角矩阵为Ds，令由L-J-1个小特征值构成的对角矩阵为Dn，令J+1个大特征值对应的特征向量构成的矩阵为Vs，令L-J-1个小特征值对应的特征向量构成的矩阵为Vn，其表达式分别为：Ds＝diagλ1,λ2,…,λi',…,λJ+1，Dn＝diagλJ+2,λJ+3,…,λi”,…,λLVs＝[v1,v2,…,vi',…,vJ+1]，Vn＝[vJ+2,vJ+3,…,vi”,…,vL]其中，λi”表示第i”个特征值，Ds表示由J+1个大特征值构成的对角矩阵，Dn表示由L-J-1个小特征值构成的对角矩阵，diag表示对角矩阵，Vs表示J+1个大特征值对应的特征向量构成的矩阵，记为信号子空间；Vn表示L-J-1个小特征值对应的特征向量构成的矩阵，记为噪声子空间；3b计算L个子阵接收J+1个入射信号的自适应权矢量为W0：W0＝μR-1aθ1其中，μ为设定系数，μ≠0；R表示L个子阵接收J+1个入射信号的协方差矩阵，上标-1为求逆操作，aθ1表示目标信号的导向矢量；从而计算得到L个子阵接收J+1个入射信号的和波束权值Wsum，其表达式为：Wsum＝VsVsHW0＝μVsVsHR-1aθ1；步骤4，计算得到L个子阵接收J+1个入射信号的差波束权值Wdiff；其得到过程为：将L个子阵接收J+1个入射信号的差波束权值记为Wdiff，并根据L个子阵接收J+1个入射信号的和波束权值Wsum，确定单脉冲比曲线，其表达式为： 其中，Γ是单脉冲比曲线的斜率，aθ1表示目标信号的导向矢量，θ1表示目标信号的来波方向，θ表示目标信号的实际来波方向，目标信号的来波方向θ1为目标信号的期望来波方向；Wdiff表示L个子阵接收J+1个入射信号的差波束权值，Wdiff∈Cm×1，Wsum表示L个子阵接收J+1个入射信号的和波束权值，Wsum∈Cm×1，Cm×1表示m×1维的复数矩阵，Sk表示L个子阵在k时刻接收J+1个入射信号的复包络，nk表示L个子阵在k时刻接收J+1个入射信号的高斯白噪声，k表示时间变量；Γθ1-θ表示单脉冲比， 表示L个子阵接收J+1个入射信号的和波束， 表示L个子阵接收J+1个入射信号的差波束，Re表示取实部操作；4a以单脉冲比曲线表达式作为求L个子阵接收J+1个入射信号的差波束的约束条件，以L个子阵接收J+1个入射信号的差波束输出信号功率最小为目标函数，构造以下形式来推导出L个子阵接收J+1个入射信号的差波束权值： 其中，ReWdiffHC＝g'为约束条件函数，C表示单脉冲比曲线表达式中 的归一化形式对应的复数矩阵，g'表示单脉冲比曲线表达式中Γθ1-θ的形式对应的实数矩阵，g'＝[Γ△θ0-Γ△θ],Γ是单脉冲比曲线的斜率，θ表示目标信号的实际来波方向，aθ+△θ表示L个子阵接收的J+1个入射信号中偏离目标实际方向△θ的导向矢量，aθ表示目标信号的实际来波方向θ的导向矢量，aθ-△θ表示L个子阵接收的J+1个入射信号中偏离目标实际方向-△θ的导向矢量，△θ表示目标信号偏离L个子阵接收J+1个入射信号波束中心的角度，aθ1表示目标信号的导向矢量，θ1表示目标信号的来波方向，目标信号的来波方向θ1为目标信号的期望来波方向；θ表示目标信号的实际来波方向，Wdiff表示L个子阵接收J+1个入射信号的差波束权值，Wsum表示L个子阵接收J+1个入射信号的和波束权值，Sk表示L个子阵在k时刻接收J+1个入射信号的复包络，nk表示L个子阵在k时刻接收J+1个入射信号的高斯白噪声，k表示时间变量；Γθ1-θ表示单脉冲比，表示L个子阵接收J+1个入射信号的和波束，表示L个子阵接收J+1个入射信号的差波束，Re表示取实部操作，min表示取最小值操作，上标H表示共轭转置操作；4b运用拉格朗日乘子法构造函数： 其中，f表示根据约束条件函数ReWdiffHC＝g'与得到的函数，ζ表示拉格朗日乘子，R表示L个子阵接收J+1个入射信号的协方差矩阵，Wdiff表示L个子阵接收J+1个入射信号的差波束权值；计算L个子阵接收J+1个入射信号的差波束权值Wdiff的复梯度其表达式为： 其中，表示求偏导，Re表示取实部操作，Im表示取虚部操作；从而计算得到L个子阵接收J+1个入射信号的差波束权值Wdiff，其计算表达式为：Wdiff＝R-1C[ReCHR-1C]-Tg'T；其中，上标-1为求逆操作，Re表示取实部操作；步骤5，根据L个子阵接收J+1个入射信号的和波束权值Wsum、L个子阵接收J+1个入射信号的差波束权值Wdiff，以及t时刻均匀线阵的L维降维信号输出，计算得到目标信号的实际来波方向。

全文数据：基于特征空间的雷达自适应和差波束测角方法技术领域[0001]本发明属于雷达技术领域，特别涉及一种基于特征空间的雷达自适应和差波束测角方法，适用于有干扰的条件下雷达检测目标的实际位置。背景技术[0002]阵列信号处理算法因其宽广的应用前景，近些年来发展相当快速，自适应波束形成也受到了更多的关注;在理想条件下，自适应波束形成能够有效保留期望信号，抑制干扰信号和杂波，使得阵列输出信干噪比最大。然而实际系统却存在各种各样的误差，包括自适应样本有限引起的协方差估计误差、约束导向矢量的指向误差和各种系统误差，如阵元位置误差、阵元的幅相误差、阵元之间的互耦、通道频率特性失配等，这时自适应波束性能会大大下降，甚至完全失效，尤其是当协方差矩阵中含有期望信号时，这种现象更为明显。[0003]传统的自适应波束形成方法绝大多数是基于阵元级的天线阵列来进行研究的，形成阵元级波束形成方法，但是随着大型雷达阵列的推广使用，阵元级波束形成方法的弊端也逐渐凸显出来:运算量相当大，系统实现软硬件成本也相当高;和、差波束方向图因为在实现角度估计和跟踪方面发挥着独特的作用，在相控阵雷达以及合成孔径声纳方面应用广泛，使得这个技术一直以来得到了很多学者的重视并且有很多成果问世。但是存在干扰时，自适应信号处理在抑制干扰的同时影响了和、差波束的权矢量，使和差方向图在主瓣附近产生畸变，致使自适应鉴角曲线与普通和差鉴角曲线有差异，给测角带来了偏差，从而导致和差测角性能下降甚至失效。发明内容[0004]针对上述现有技术存在的不足，本发明的目的在于提出一种基于特征空间的雷达自适应和差波束测角方法，该种基于特征空间的雷达自适应和差波束测角方法能够运用子阵来减少阵列的复杂度，同时在干扰存在的情况下依然能够维持单脉冲测角，保持良好的测角性能；而且运用一种新的方法求得差波束权值，易得出目标的实际方向，易于工程实现。[0005]本发明的主要思路:对均匀线阵的接收信号进行降维处理，得到子阵接收信号，对子阵接收信号的协方差矩阵进行特征值分解，得到信号子空间和噪声子空间；在输出期望信号不变的情况下使均匀线阵输出的方差最小，求出降维后子阵的最优权值，再将此最优权值向信号子空间作投影，得到基于特征空间的自适应和波束最优权值;接着考虑存在干扰的情况下，为了保持单脉冲测角性能，提出新的差波束权值的线性关系作为求差波束权值的约束条件，以差波束输出信号功率最小为目标函数，求得差波束权值和单脉冲比曲线，从而易得出目标的实际方向。[0006]为达到上述技术目的，本发明采用如下技术方案予以实现。[0007]—种基于特征空间的雷达自适应和差波束测角方法，包括以下步骤：[0008]步骤1，确定均匀线阵，该均匀线阵包括M个阵元，均匀线阵的检测范围内存在J+1个信号源，J+1个信号源向均匀线阵发射J+1个入射信号，所述J+1个入射信号中包含目标信号;进而确定t时刻均匀线阵接收的M维信号；[0009]其中，t表示时间变量，M表示均匀线阵包含的阵元个数，M、J分别为大于0的正整数；[0010]步骤2,将M个阵元的均匀线阵划分为L个子阵，并根据t时刻均匀线阵接收的M维信号，得到t时刻均匀线阵的L维降维信号输出；[0011]其中，L表示将M个阵元的均匀线阵划分后包含的子阵个数，L为大于0的正整数；[0012]步骤3，计算得到L个子阵接收J+1个入射信号的和波束权值Wsum;[0013]步骤4，计算得到L个子阵接收J+1个入射信号的差波束权值Wdiff;[00M]步骤5,根据L个子阵接收J+1个入射信号的和波束权值W_、L个子阵接收J+1个入射信号的差波束权值Wdlff，以及t时刻均匀线阵的L维降维信号输出，计算得到目标信号的实际来波方向。[0015]本发明与现有技术相比具有以下优点：[0016]现有和差波束测角方法在有干扰的条件下，自适应信号处理在抑制干扰的同时影响了和、差波束的权矢量，使和差方向图发生畸变，致使自适应鉴角曲线与普通和差鉴角曲线有差异，给测角带来了较大的偏差，从而导致和差测角性能下降甚至失效;而本发明运用了子阵，从而减小了阵列的复杂度，不仅运算量小、收敛速度快，而且大大地减少了硬件与软件的成本;关键是在存在干扰的条件下，仍能维持单脉冲测角方法，稳健性较好，保持正确的测角性能，能够有效精确地在干扰环境下测得目标的实际位置。附图说明[0017]下面结合附图和具体实施方式对本发明作进一步详细说明。[0018]图1为本发明的一种基于特征空间的雷达自适应和差波束测角方法流程图；[0019]图2为本发明方法的和差波束方向图；[0020]图3为本发明方法得到的单脉冲比曲线图；[0021]图4为在幅相误差条件下本发明方法的和差波束方向图；[0022]图5为在幅相误差条件下本发明方法得到的单脉冲比曲线图。具体实施方式[0023]参照图1，为本发明的一种基于特征空间的雷达自适应和差波束测角方法流程图；其中所述基于特征空间的雷达自适应和差波束测角方法，包括以下步骤：[0024]步骤1，确定均匀线阵，该均匀线阵包括M个阵元，阵元间距为d，M个阵元的方向图为各向同性，均匀线阵的检测范围内存在J+1个信号源，J+1个信号源向均匀线阵发射J+1个入射信号，J+1个入射信号包括一个目标信号和J个干扰信号，且J+1个入射信号分别为窄带信号，且波长都为λ，第i个入射信号的来波方向为0i，ie{1，2,…，J+1}，J+1个入射信号的来波方向不同，分别为θ1，θ2,…，θί，···，θ』+1。[0025]J+1个信号源向均匀线阵发射J+1个入射信号，每个阵元分别接收J+1个不相关的窄带信号，其中将第m个阵元接收的J+1个不相关的窄带信号记为sm⑴，[0026]Smt={Sm,lt，Sm,2⑴，…，Sm,i⑴，…，SmJ+lt}，me{1，2，…，M}，Sm,i⑴表示第m个阵元接收的第i个窄带信号;则t时刻均匀线阵接收的M维信号为Xt，其表达式为：[0030]其中，ie{l，2，〜，J+l}，且i=l时对应目标信号，i辛1时分别对应干扰信号；Xmt表示t时刻均勾线阵中第m个阵元接收的信号，所述t时刻均勾线阵中第m个阵元接收的信号为第m个阵元接收的J+1个不相关的窄带信号；Sl⑴为t时刻第i个入射信号的复包络；aQ1为第i个入射信号的导向矢量，η⑴为t时刻均匀线阵中M个阵元的高斯白噪声，nmt表示t时刻均匀线阵中第m个阵元的高斯白噪声，上标T表示转置操作，d表示均匀线阵的阵元间距，Q1表示第i个入射信号的来波方向，e表示指数函数，上标j表示虚数单位，sin为正弦函数;M表不第i个入射信号的波长，且每一个入射信号的波长取值分别相等;t表不时间变量，M表示均匀线阵包含的阵元个数，M为大于0的正整数。[0031]步骤2,将M个阵元的均匀线阵划分为L个子阵，每个子阵包含的阵元个数均为g，g为正整数：表示向上取整;进而计算得到划分L个子阵后第i个入射信号的导向矢量为[0033]其中，划分L个子阵后第i个入射信号的导向矢量钱⑩包含L个元素，以第1个子阵为参考单元，则L个元素中第1个元素为1，且L个元素之间的间距为gd;g表示每个子阵包含的阵元个数，g为正整数，g多l，d表示均匀线阵的阵元间距，Q1表示第i个入射信号的来波方向，上标T表示转置操作，e表示指数函数，上标j表示虚数单位，sin为正弦函数;Ai表示第i个入射信号的波长，且每一个入射信号的波长取值分别相等。[0034]每个子阵的导向矢量均以对应子阵中第一个阵元为参考，因此得到同一来波方向的条件下第i个入射信号的导向矢量:，其表达式为：[0036]由于每一个入射信号的波长取值分别相等，则令λ表不每一个入射信号的波长;进而计算得到L个子阵的LXM维降维矩阵T:[0038]其中1表不目标信号的来波方向，目标信号的来波方向01为目标·?信号的期望来波方向；e表示指数函数，上标j表示虚数单位，sin为正弦函数，λ表示每一个入射信号的波长，g表示每个子阵包含的阵元个数，g为正整数，g多1;然后根据L个子阵的LXM维降维矩阵T，将t时刻均匀线阵接收的M维信号X⑴降到L维，得到t时刻均匀线阵的L维降维信号输出mt:[0039]m⑴=Tx⑴=[mit，m2t，…，mrt，…，mLt]τ[0040]其中，mr⑴表示t时刻均匀线阵第Γ个子阵的降维信号输出，Γ=I，…L。[0041]步骤3,现对子阵进行处理，利用基于特征空间的自适应波束形成方法求出和波束权值。[0042]3a—个目标信号，J个干扰信号，在有限次快拍下，计算得至IjL个子阵接收J+1个入射信号的协方差矩阵R，并对其作特征值分解，得到：[0043][0044]其中，L个子阵接收J+1个入射信号的协方差矩阵R包含L个特征值，且彡…彡表示每个子阵接收J+1个入射信号的高斯白噪声功率，且将A1,入1，…，λί，，···，λ:+!记为J+1个大特征值，将λ:+2，λ:+3,…，λί”，···人记为L-J-I个小特征值，Vi1表示第i’个特征值对应的特征矢量，i’=1，…，J+1，vr表示第i”个特征值对应的特征矢量，i”=J+2,…，L，Ai’表示第i’个特征值，i”=J+2,J+3,…，L，Ar表示第i”个特征值，上标H表示共辄转置操作。[0045]分别令由J+1个大特征值构成的对角矩阵为Ds，令由L-J-1个小特征值构成的对角矩阵为Dn，令J+1个大特征值对应的特征向量构成的矩阵为Vs，令L-J-I个小特征值对应的特征向量构成的矩阵为Vn，其表达式分别为：[0048]其中，V’表示第i”个特征值，Ds表示由J+1个大特征值构成的对角矩阵，Dn表示由L-J-I个小特征值构成的对角矩阵，diag表示对角矩阵，Vs表示J+1个大特征值对应的特征向量构成的矩阵，记为信号子空间；Vn表示L-J-1个小特征值对应的特征向量构成的矩阵，记为噪声子空间。[0049]3b在线性约束最小方差准则中，L个子阵接收J+1个入射信号的自适应权矢量为Wo:[0050]W0=Ur1aθι[0051]其中，μ为设定系数，取值为任意非零常数;R表示L个子阵接收J+1个入射信号的协方差矩阵，上标_1为求逆操作，aθι表不目标彳目号的导向矢量;此自适应权矢量Wo是由彳目号子空间和噪声子空间构成。[0052]在理想情况下，目标信号位于信号子空间内，因此有，因此，L个子阵接收J+I个入射信号的自适应权矢量Wo仅为信号子空间的分量，噪声子空间的分量为0;从而计算得到L个子阵接收J+1个入射信号的和波束权值Wsum，其表达式为：[0054]步骤4，将L个子阵接收J+1个入射信号的差波束权值记为Wdiff，并根据L个子阵接收J+1个入射信号的和波束权值Wsum，确定单脉冲比曲线，其表达式为：[0056]其中，Γ是单脉冲比曲线的斜率，a^表示目标信号的导向矢量，Q1表示目标信号的来波方向，Θ表示目标信号的实际来波方向，目标信号的来波方向Θ1为目标信号的期望来波方向;Wdiff表示L个子阵接收J+1个入射信号的差波束权值，表示L个子阵接收J+I个入射信号的和波束权值：^表不mX1维的复数矩阵，Sk表不L个子阵在k时刻接收J+1个入射信号的复包络，nk表示L个子阵在k时刻接收J+1个入射信号的高斯白噪声，k表示时间变量；Γθ^θ表示单脉冲比，[0057]表不L个子阵接收J+1个入射信号的和波束，[0058]ί表示L个子阵接收J+1个入射信号的差波束，Re表示取实部操作。[0059]L个子阵接收J+1个入射信号的和波束为L个子阵接收的J+1个入射信号在目标信号方向上形成的波束，L个子阵接收J+1个入射信号的差波束为L个子阵接收的J+1个入射信号在目标信号方向形成零陷后的波束。[0060]在存在干扰的情况下，保持单脉冲测角性能，要求L个子阵接收J+1个入射信号的差波束最优权值满足单脉冲比曲线表达式等式。[0061]4a以单脉冲比曲线表达式作为求L个子阵接收J+1个入射信号的差波束的约束条件，以L个子阵接收J+1个入射信号的差波束输出信号功率最小为目标函数，构造以下形式来推导出L个子阵接收J+1个入射信号的差波束权值：[0064]其中，ReWdiffflC=g’为约束条件函数，C表示单脉冲比曲线表达式中[0065]的归一化形式对应的复数矩阵，表示mX3维的复数矩阵，g’表示单脉冲比曲线表达式中rQ1-Q的形式对应的实数矩阵，θ表示目标信号的实际来波方向，aΘ+ΛΘ表不L个子阵接收的J+1个入射信号中偏离目标实际方向ΛΘ的导向矢量，aΘ表不目标信号的实际来波方向Θ的导向矢量，aΘ-ΛΘ表不L个子阵接收的J+1个入射信号中偏离目标实际方向-ΛΘ的导向矢量，ΛΘ表不目标信号偏离L个子阵接收J+I个入射信号波束中心的角度，aθι表不目标信号的导向矢量，0:表不目标信号的来波方向，目标信号的来波方向θι为目标信号的期望来波方向；Θ表示目标信号的实际来波方向，Wdiff表不L个子阵接收J+1个入射信号的差波束权值，WdiffeCmxl，Wsum表不L个子阵接收J+1个入射信号的和波束权值，WsumeCmxl，Cmxl表示mX1维的复数矩阵，Sk表示L个子阵在k时刻接收J+1个入射信号的复包络，nk表示L个子阵在k时刻接收J+1个入射信号的高斯白噪声，k表示时间变量；Γθ^θ表示单脉冲比，表示L个子阵接收J+1个入射信号的和波束[0066]表不L个子阵接收J+1个入射信号的差波束，Re表不取实部操作，min表示取最小值操作，上标H表示共辄转置操作。[0067]而单脉冲比曲线表达式中Γθ^θ的形式对应的实数矩阵g’为：[0068][0069]其中，Rlx3表示IX3维的实数矩阵。[0070]4b求得L个子阵接收J+1个入射信号的差波束权值。[0071]运用拉格朗日乘子法构造函数：[0073]其中，f表示根据约束条件函数ReWdlffflC=g’与得到的函数，ζ表示拉格朗日乘子，GeR3xl,R3xl表示3Χ1维的实数矩阵，R表示L个子阵接收J+1个入射信号的协方差矩阵，Wdiff表示L个子阵接收J+1个入射信号的差波束权值。[0074]为了解决上述问题，引进L个子阵接收J+1个入射信号的差波束权值Wdiff的复梯度，其表达式为：[0075][0076]其中，3表不求偏导，Re表不取实部操作，Im表不取虚部操作。[0077]从而计算得到L个子阵接收J+1个入射信号的差波束权值Wdlff，其计算表达式为：[0078][0079]其中，上标-1为求逆操作，Re表示取实部操作。[0080]步骤5,根据L个子阵接收J+1个入射信号的和波束权值Wsum和L个子阵接收J+1个入射信号的差波束权值Wdlff，得到L个子阵接收J+1个入射信号的差波束信号输出5^与1^个子阵接收J+1个入射信号的和波束信号输出7Σ的比值，并取所述比值的实部，得到L个子阵的自适应单脉冲比，进而得到目标信号偏离L个子阵接收J+1个入射信号波束中心的角度ΛΘ：[0083]其中，上标H表示共辄转置操作，m⑴表示t时刻L个子阵接收的J+1个入射信号，Γ表示单脉冲比曲线的斜率，t表示时间变量。[0084]根据目标信号的来波方向9:和目标信号偏离L个子阵接收J+1个入射信号波束中心的角度ΛΘ，计算得到目标信号的实际来波方向θ,θιθ^ΛΘ;得到所述目标信号的实际来波方向的同时也抑制了J个干扰信号的来波方向。[0085]通过以下仿真试验对本发明效果作进一步验证说明。[0086]一仿真条件:仿真采用阵元间距为0.6波长的40元均匀线阵，阵元天线方向图各向同性;均匀线阵接收数据中包含加性高斯白噪声，目标信号的方向为〇°，干扰信号与目标信号不相干，每个子阵的维数相同均为1X4。[0087]均匀线阵发射的信号为线性调频信号，目标信号带宽为50ΜΗΖ;均匀线阵发射的信号基带频率为10GHZ，采样频率为100ΜΗΖ，信噪比为OdB，5个互不相关且与目标信号不相关的远场干扰信号分别从-8°、-10°、_19°、19°、36°入射，干噪比为31.6^，采样快拍数为4096〇[0088]二仿真内容：[0089]采用本发明方法得到的和差波束方向图，仿真结果如图2所示，图2为本发明方法的和差波束方向图；由本发明方法得到的单脉冲比曲线图，仿真结果如图3所示，图3为本发明方法得到的单脉冲比曲线图；在有幅相误差条件下，采用本发明方法得到的和差波束方向图，仿真结果如图4所示，图4为在幅相误差条件下本发明方法的和差波束方向图；在有幅相误差条件下，由本发明方法得到的单脉冲比曲线图，仿真结果如图5所示，图5为在幅相误差条件下本发明方法得到的单脉冲比曲线图。[0090]三仿真结果分析：[0091]从图2和图4可以看出，采用本发明在干扰方向形成很深的零陷，能够有效的抑制干扰，且在幅相误差条件下依然在干扰方向能形成很深的零陷，说明稳健性好。[0092]从图3、图5可以看出，在无幅相误差和有幅相误差的条件下，得出的单脉冲比曲线方向图和普通和差方向图单脉冲比曲线相一致，说明运用这种新的方法求得差波束权值，测出目标偏离波束中心的角度，从而得出目标的实际方向的性能很好。[0093]综上所述，仿真实验验证了本发明的正确性，有效性和可靠性。[0094]显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围；这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。

权利要求：1.一种基于特征空间的雷达自适应和差波束测角方法，其特征在于，包括以下步骤：步骤1，确定均匀线阵，该均匀线阵包括M个阵元，均匀线阵的检测范围内存在J+1个信号源，J+1个信号源向均匀线阵发射J+1个入射信号，所述J+1个入射信号中包含目标信号；进而确定t时刻均匀线阵接收的M维信号；其中，t表示时间变量，M表示均匀线阵包含的阵元个数，M、J分别为大于O的正整数；步骤2,将M个阵元的均匀线阵划分为L个子阵，并根据t时刻均匀线阵接收的M维信号，得到t时刻均匀线阵的L维降维信号输出；其中，L表示将M个阵元的均匀线阵划分后包含的子阵个数，L为大于O的正整数；步骤3,计算得到L个子阵接收J+1个入射信号的和波束权值Wsum;步骤4，计算得到L个子阵接收J+1个入射信号的差波束权值Wdiff;步骤5,根据L个子阵接收J+1个入射信号的和波束权值Wsum、L个子阵接收J+1个入射信号的差波束权值Wdlff，以及t时刻均匀线阵的L维降维信号输出，计算得到目标信号的实际来波方向。2.如权利要求1所述的一种基于特征空间的雷达自适应和差波束测角方法，其特征在于，在步骤1中，所述t时刻均匀线阵接收的M维信号，其确定过程为：J+1个入射信号包括一个目标信号和J个干扰信号，且J+1个入射信号分别为窄带信号，第i个入射信号的来波方向为9i，ie{1，2，〜，J+1}，且i=l时对应目标信号，i辛1时分别对应干扰信号;J+1个入射信号的来波方向不同，分别为θ1，θ2,…，θί,…，0j+1;J+1个信号源向均匀线阵发射J+1个入射信号，每个阵元分别接收J+1个不相关的窄带信号，其中将第m个阵元接收的J+1个不相关的窄带信号记为Sm⑴，Sm⑴={Sm,l⑴，Sm,2⑴，…，Sm,i⑴，…，Sm,J+lt}，me{1，2,…，M}，Sm,i⑴表示第Π1个阵元接收的第i个窄带信号;则t时刻均匀线阵接收的M维信号为Xt，其表达式为：其中，xmt表示t时刻均匀线阵中第m个阵元接收的信号，所述t时刻均匀线阵中第m个阵元接收的信号为第m个阵元接收的J+1个不相关的窄带信号；Slt为t时刻第i个入射信号的复包络;aQ1为第i个入射信号的导向矢量，nt为t时刻均匀线阵中M个阵元的高斯白噪声，nmt表示t时刻均匀线阵中第m个阵元的高斯白噪声，上标T表示转置操作，d表示均匀线阵的阵元间距，S1表示第i个入射信号的来波方向，e表示指数函数，上标j表示虚数单位，sin为正弦函数;M表示第i个入射信号的波长，且每一个入射信号的波长取值分别相等;t表示时间变量，M表示均匀线阵包含的阵元个数，M为大于O的正整数。3.如权利要求2所述的一种基于特征空间的雷达自适应和差波束测角方法，其特征在于，在步骤2中，所述t时刻均匀线阵的L维降维信号输出，其得到过程为：将M个阵元的均匀线阵划分为L个子阵，每个子阵包含的阵元个数均为g，g为正整数，g，「Ί表示向上取整;进而计算得到划分L个子阵后第i个入射信号的导向矢量为m:其中，g表示每个子阵包含的阵元个数，d表示均匀线阵的阵元间距，Q1表示第i个入射信号的来波方向，上标T表示转置操作，e表示指数函数，上标j表示虚数单位，sin为正弦函数；M表示第i个入射信号的波长，且每一个入射信号的波长取值分别相等；进而得到同一来波方向的条件下第i个入射信号的导向矢量其表达式为：进而计算得到L个子阵的LXM维降维矩阵T:其中，Q1表不目标信号的来波方向，目标信号的来波方向Q1为目标信号的期望来波方向；e表示指数函数，上标j表示虚数单位，sin为正弦函数，λ表示每一个入射信号的波长，g表示每个子阵包含的阵元个数，g为正整数，g$l;然后根据L个子阵的LXM维降维矩阵T，将t时刻均匀线阵接收的M维信号X⑴降到L维，得到t时刻均匀线阵的L维降维信号输出mt:mt=Txt=[mit，m2t，…，mrt，…，mLt]τ其中，mrt表示t时刻均匀线阵第Γ个子阵的降维信号输出，Γ=1，一UL表示将M个阵元的均匀线阵划分后包含的子阵个数，L为大于O的正整数。4.如权利要求3所述的一种基于特征空间的雷达自适应和差波束测角方法，其特征在于，在步骤3中，所述L个子阵接收J+1个入射信号的和波束权值Wsum，其得到过程为：3a计算得到L个子阵接收J+1个入射信号的协方差矩阵R，并对其作特征值分解，得到：其中，L个子阵接收J+1个入射信号的协方差矩阵R包含L个特征值，且表示每个子阵接收J+1个入射信号的高斯白噪声功率，且将λι,A1,…，λί，，…，λ:+!记为J+1个大特征值，将λ:+2，λ:+3，…，λί”，…，\[^记为L-J-I个小特征值，vr表不第i’个特征值对应的特征矢量，i’=1，…，J+I，vi”表示第i”个特征值对应的特征矢量，i”=J+2,…,L，Ai，表示第i’个特征值，i”=J+2，J+3，…，L，A1，，表示第i”个特征值，上标H表示共辄转置操作；分别令由J+1个大特征值构成的对角矩阵为Ds，令由L-J-I个小特征值构成的对角矩阵为〇„，令J+1个大特征值对应的特征向量构成的矩阵为Vs，令L-J-1个小特征值对应的特征向量构成的矩阵为Vn，其表达式分别为：Ds=diagλι，入2，…，λί’，…，Aj+i，Dn=diagAj+2，Aj+3，…，λί”，…，λ〇Vs=[VI，V2,···，Vi，，…，VJ+l]，Vn=[VJ+2，VJ+3,···，Vi”，…，VL]其中，V表示第i”个特征值，Ds表示由J+1个大特征值构成的对角矩阵，Dn表示由L-J-I个小特征值构成的对角矩阵，diag表示对角矩阵，Vs表示J+1个大特征值对应的特征向量构成的矩阵，记为信号子空间；Vn表示L-J-I个小特征值对应的特征向量构成的矩阵，记为噪声子空间；3b计算L个子阵接收J+1个入射信号的自适应权矢量为W0:Wo=yR_1aθι其中，μ为设定系数，μ辛〇;R表示L个子阵接收J+1个入射信号的协方差矩阵，上标-1为求逆操作，aQ1表不目标信号的导向矢量；从而计算得到L个子阵接收J+1个入射信号的和波束权值Wsum，其表达式为：5.如权利要求4所述的一种基于特征空间的雷达自适应和差波束测角方法，其特征在于，在步骤4中，所述L个子阵接收J+1个入射信号的差波束权值Wdlff，其得到过程为：将L个子阵接收J+1个入射信号的差波束权值记为Wdiff，并根据L个子阵接收J+1个入射信号的和波束权值Ws™，确定单脉冲比曲线，其表达式为：其中，Γ是单脉冲比曲线的斜率，aQ1表示目标信号的导向矢量，Q1表示目标信号的来波方向，Θ表示目标信号的实际来波方向，目标信号的来波方向Θ1为目标信号的期望来波方向;Wdiff表示L个子阵接收J+1个入射信号的差波束权值表示L个子阵接收J+1个入射信号的和波束权值表不mX1维的复数矩阵，Sk表不L个子阵在k时刻接收J+1个入射信号的复包络，nk表示L个子阵在k时刻接收J+1个入射信号的高斯白噪声，k表示时间变量；Γθ^θ表示单脉冲比，表不L个子阵接收J+1个入射信号的和波束，表示L个子阵接收J+1个入射信号的差波束，Re表示取实部操作；4a以单脉冲比曲线表达式作为求L个子阵接收J+1个入射信号的差波束的约束条件，以L个子阵接收J+1个入射信号的差波束输出信号功率最小为目标函数，构造以下形式来推导出L个子阵接收J+1个入射信号的差波束权值：其中，ReWdiff11C=g’为约束条件函数，C表示单脉冲比曲线表达式中的归一化形式对应的复数矩阵，g’表示单脉冲比曲线表达式中ΓQ1-Q的形式对应的实数矩阵，g’=[rA0O-ΓΛΘ]，Γ是单脉冲比曲线的斜率，Θ表示目标信号的实际来波方向，aΘ+ΛΘ表示L个子阵接收的J+1个入射信号中偏离目标实际方向ΛΘ的导向矢量，aΘ表不目标信号的实际来波方向Θ的导向矢量，aΘ-ΛΘ表不L个子阵接收的J+1个入射信号中偏离目标实际方向-ΛΘ的导向矢量，ΛΘ表不目标信号偏离L个子阵接收J+I个入射信号波束中心的角度，aθι表不目标信号的导向矢量，0:表不目标信号的来波方向，目标信号的来波方向θι为目标信号的期望来波方向；Θ表示目标信号的实际来波方向，Wdiff表示L个子阵接收J+1个入射信号的差波束权值，Wsum表示L个子阵接收J+1个入射信号的和波束权值，Sk表不L个子阵在k时刻接收J+1个入射信号的复包络，nk表不L个子阵在k时刻接收J+1个入射信号的高斯白噪声，k表示时间变量；Γθ^θ表示单脉冲比，表不L个子阵接收J+1个入射信号的和波束，表示L个子阵接收J+1个入射信号的差波束，Re表示取实部操作，min表示取最小值操作，上标H表示共辄转置操作；4b运用拉格朗日乘子法构造函数：其中，f表示根据约束条件函数得到的函数，ζ表示拉格朗日乘子，R表示L个子阵接收J+1个入射信号的协方差矩阵，Wdlff表示L个子阵接收J+1个入射信号的差波束权值；计算L个子阵接收J+1个入射信号的差波束权值Wdlff的复梯度，其表达式为：其中，δ表示求偏导，Re表示取实部操作，Im表示取虚部操作；从而计算得到L个子阵接收J+1个入射信号的差波束权值Wdlff,其计算表达式为：其中，上标-1为求逆操作，Re表示取实部操作。6.如权利要求5所述的一种基于特征空间的雷达自适应和差波束测角方法，其特征在于，在步骤5中，所述目标信号的实际来波方向，其得到过程为：根据L个子阵接收J+1个入射信号的和波束权值Wsum和L个子阵接收J+1个入射信号的差波束权值Wdiff，得到L个子阵接收J+1个入射信号的差波束信号输出yA与L个子阵接收J+1个入射信号的和波束信号输出7Σ的比值，并取所述比值的实部，得到L个子阵的自适应单脉冲比，进而得到目标信号偏离L个子阵接收J+1个入射信号波束中心的角度ΛΘ:其中，上标H表不共辄转置操作，mt表不t时刻L个子阵接收的J+1个入射信号，Γ表不单脉冲比曲线的斜率，t表示时间变量；根据目标信号的来波方向Θ#Ρ目标信号偏离L个子阵接收J+1个入射信号波束中心的角度ΛΘ，计算得到目标信号的实际来波方向θ，Θ=Θ^ΛΘ。

百度查询： 西安电子科技大学 基于特征空间的雷达自适应和差波束测角方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD