CN112612019B - 一种基于频控阵相位中心的主动抗干扰的方法 - Google Patents
一种基于频控阵相位中心的主动抗干扰的方法 Download PDFInfo
- Publication number
- CN112612019B CN112612019B CN202011380267.XA CN202011380267A CN112612019B CN 112612019 B CN112612019 B CN 112612019B CN 202011380267 A CN202011380267 A CN 202011380267A CN 112612019 B CN112612019 B CN 112612019B
- Authority
- CN
- China
- Prior art keywords
- fda
- phase center
- phase
- array
- radar
- 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
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000005457 optimization Methods 0.000 claims abstract description 19
- 230000005855 radiation Effects 0.000 claims abstract description 18
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 16
- 239000002245 particle Substances 0.000 claims abstract description 13
- 238000001514 detection method Methods 0.000 claims abstract description 11
- 239000013598 vector Substances 0.000 claims description 23
- 239000002131 composite material Substances 0.000 claims description 12
- 239000008186 active pharmaceutical agent Substances 0.000 claims description 8
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 2
- 230000005684 electric field Effects 0.000 claims 1
- 230000007123 defense Effects 0.000 abstract 1
- 230000006870 function Effects 0.000 description 8
- 230000008859 change Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000010485 coping Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000030279 gene silencing Effects 0.000 description 1
- 210000000987 immune system Anatomy 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000010587 phase diagram Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S2013/0236—Special technical features
- G01S2013/0245—Radar with phased array antenna
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于频控阵相位中心的主动抗干扰的方法,具体按照以下步骤实施:步骤1、构造频控阵雷达系统;步骤2、推导FDA辐射信号的相位中心;步骤3、通过调节频率增量,控制相位中心偏差;步骤4、通过一种免疫粒子群优化算法寻找最优频偏,使相位中心偏差最大。本发明的方法使相位中心偏差最大,参数可控,在提高自身探测精度的同时,实现对敌方干扰机的角度欺骗,从而实现基于频控阵相位中心的主动抗干扰,为电子战中的雷达对抗提供了新思路,与传统的被动抗干扰方法相比,具有实现主导对抗的积极防御,参数可控,避免出现干扰发生变化穷于应对等情况的优点。
Description
技术领域
本发明涉及雷达信号处理技术,具体涉及一种基于频控阵相位中心的主动抗干扰的方法。
背景技术
随着电子对抗系统的发展,特别是能够产生多维、灵活调制的欺骗性干扰信号的有源干扰机的出现和发展,对雷达系统构成了严重威胁。因此,迫切需要发展新的雷达抗干扰技术。由于在远场条件下波束的距离角度相关性,FDA被用于在距离范围内识别假目标。但是,现有技术中大多都是针对特定场景,对产生的假目标进行被动抗干扰,被动抗干扰容易出现受干扰变化而穷于应对的被动情况。
发明内容
本发明提供了一种基于频控阵相位中心的主动抗干扰的方法,解决了现有技术对假目标进行被动抗干扰,易受干扰变化而穷于应对的问题。
本发明采用以下技术方案:一种基于频控阵相位中心的主动抗干扰的方法,其特征在于,具体按照以下步骤实施:
步骤1、构造频控阵雷达系统;
步骤2、推导FDA辐射信号的相位中心;
步骤3、通过调节频率增量,控制相位中心偏差;
步骤4、通过一种免疫粒子群优化算法寻找最优频偏,使相位中心偏差最大。
本发明的特点还在于,
步骤1的具体步骤为:
假设一部由N个阵元组成的均匀线阵频控阵雷达,如图2所示,其中频控阵阵元数为Q,C为光速,相邻两个阵元之间的频率增量为Δf,初始载频f0,第q个阵元处的频率为fq=f0+Δfq,阵元间距为d,目标位置为(θ’,R’),目标方位角为θ,目标与频控阵之间的距离为R;第q个阵元发射的信号表示为:
Stq(t)=exp[j2πfq(t-Rq/c)] (1)
Aq表示第q个阵元辐射信号的幅度,在远场目标位置处的合成场强为:
在远场处,各阵元发射的信号几乎平行,即Rq=R-(q-1)d sinθ,对式(2)进行展开可得:
每个单元的信号辐射幅值是相等的,假设Aq=A=1,由于Δfq<<f0,Δfq(q-1)dsinθ/c可忽略,
为简化计算,假设:
at(t),aR(R),aθ(θ)分别包含每个元素的唯一包含t,R和θ的项,a0(t,R)包含每个元素的公共项;式(4)写作⊙表示阿达玛积;
通过推导FDA的全向合成信号模型,当需要探测位于(θ0,R0)的目标时,合成信号指向(θ0,R0),引入导向矢量:wR(R0)=a0(0,-R0)⊙aR(-R0),wθ(θ0)=aθ(-θ0),
wR(R0)和wθ(θ0)分别表示距离角度向量,则无方向性合成信号表示为:
EFDA(R,θ,R0,θ0,t)=(a0(t,R)⊙w)T·(at(t)⊙aθ(θ)⊙aR(R)) (6)
将场强用幅度和相位两部分来表示:
其相位信息:
ΨFDA(R,θ,R0,θ0,t)=angle(EFDA(R,θ,R0,θ0,t)) (8)
angle为相位角求解函数,AF=|EFDA(R,θ,R0,θ0,t)|为幅值,AF≥0。
步骤2假设相位中心M,坐标原点O,dM为从M到O的偏差向量,目标位于(θ’,R’),第q个阵元到场源的距离为:
Rq′=R′-(q-1)d sinθ+dMsinθ (9)
结合式(9)和式(4),合成信号的场强:
将其用导向矢量表示为:
E′FDA(R′,θ′,t)=a′0 T(t,R′,θ′)·(a′t(t)⊙a′θ(θ′)⊙a′R(R′)) (11)
为了探测远场(θ0,R0)处的目标,导向矢量为:
无方向性的合成信号为:
E′FDA(R′,θ′,R0,θ0,t)=(a′0(t,R′,θ′)⊙w′)T·(a′t(t)⊙a′θ(θ′)⊙a′R(R′)) (13)
AF′=|E′FDA(R′,θ′,R0,θ0,t)|表示以M为参考点的FDA信号幅度信息,上式写为:
假设FDA辐射的为球面波,以相位中心为参考点,在相同距离R′0处的相位为常数值:
C=Ψ′FDA(R′0,θ′,R0,θ0,t)=angle(E′FDA(R′0,θ′,R0,θ0,t)) (15)
结合式(5)和式(12),式(13)写为:
结合式(15)和式(16),式(15)写成:
式(17)是在假设线性FDA辐射球面波的基础上引入的,其实线性FDA辐射的波是准球面波;但是,方位角θ′在一个小范围内变化,而范围R′是固定的,所以式(17)成立,基于上述分析,当目标处于(θ′0,R′0),
当方位角变化δθ→0,
推出干扰机(θ′0,R′0)处的相位中心偏差为:
式(19)给出了相位中心相对于阵列原点的偏移距离,看出不同位置的相中心可能是不同的,由于相位方向图随不同雷达参数的变化而变化,所以雷达参数的变化也会导致相位中心的变化。
从步骤2中的式(19)看出干扰机(θ′0,R′0)处的相位中心偏差与ΨFDA有关,而ΨFDA与频偏有关,通过调节FDA雷达的频率增量序列Δf=[Δf1 Δf2 … Δfq … ΔfQ],(θ′0,R′0)处的相位分布也会改变,进而改变相位中心dM(R′0,θ′0);
在不考虑FDA雷达的波束多样性的情况下,假设导向矢量为:
无方向合成信号表示为:
EFDA(R′,θ′,t)=a0 T(t,R)·(at(t)⊙aθ(θ′)⊙aR(R′)) (21)
FDA的相位分布写为:
ΨFDA(R′,θ′,t)=angle(EFDA(R′,θ′,t)) (22)
对于固定时刻,t=0,式(22)写为:
通过10000次蒙特卡洛实验,证明选择合适的频偏使相位中心偏差较大。
步骤2的作用为推导出相位中心的偏差与阵列原点,目标位置因素有关,并且由于ΨFDA随雷达参数的改变而改变,通过控制雷达参数来调节相位中心,其中假设了FDA辐射的波为球面波,实际上FDA辐射的为准球面波,但如果当距离R′固定,方位角θ′在很小的方位角内变化,式(17)成立。
FDA雷达的功能包括探测目标和自我保护两个方面,因此,要求FDA雷达在保证探测能力的同时,对干扰源实施角度欺骗,在不影响正常工作的情况下实现主动抗干扰。
步骤3假设一个场景,在固定时刻,要求雷达对干扰源实施主动抗干扰,同时保持对目标的探测,从而在正常工作中达到自我保护的目的;如式(16)所示,导向矢量通过控制指向(θ0,R0),所以雷达探测固定在(θ0,R0)位置出的目标,同时对C处的干扰源实施主动抗干扰来防止自身被准确定位,因此,要求FDA雷达发射的电磁波到达C点处的相位中心与自身雷达阵列天线偏离的距离(假设用db表示)满足下式:
d′M(θ′0,R′0)≥db (24)
从式(19)得出,当时刻固定,t=0,假设目标固定在(θ0,R0),通过调节Δf使C处的相位中心满足式(20);
d′M(Δf)≥db (26)
主动抗干扰问题转化为了最优化问题,
Δfmax和Δfmin分别为阵列的最大和最小频偏,最优化问题转化为最小值问题;
步骤3中的蒙特卡洛实验中设置的FDA载频为1GHZ,阵元间距d为0.15m,阵元数目为10,时间固定在零时刻,干扰机位置为(30°,250km)。
步骤4中采用免疫粒子群优化算法寻找最优频偏,其算法复杂度为O[DS·M·Q+(MaxT-DS)M(Q+N)]。
与现有技术相比,本发明具有以下优点:(1)本发明提出了一种基于FDA相位中心的主动抗干扰方法,该方法所完成的主要任务是推导出了FDA波束的相位中心,通过比较相位辐射函数与阵列原点和相位中心之间的关系,可以推导出相位中心与远场相位的关系。在此基础上,通过寻找最优频偏进行适当的相位中心调节既可以实现准确的探测,也有利于FDA雷达的角度欺骗,可以用于主动抗干扰。同时,为了解决优化表达式中的周期性问题,提出了一种免疫粒子群优化(PSO-IMMU)算法,以提高求解速度和精度。
(2)与现有技术相比,本发明提出基于FDA相位中心的主动抗干扰方法的,即可提高雷达探测精度,也能提高FDA雷达对敌方干扰机的角度欺骗效果,实现了参数可控,大大提高了雷达主动抗干扰能力,可避免出现当干扰变化而穷于应对的被动情况,提高了武器系统主导对抗的积极防御的能力。
附图说明:
图1是本发明一种基于频控阵相位中心的主动抗干扰的方法的流程图;
图2是本发明频控阵雷达发射阵列模型图;
图3是本发明相位中心偏移示意图;
图4是本发明10000次蒙特卡洛实验的结果图;
图5是本发明雷达、目标和干扰源的位置图;
图6(a)是本发明优化FDA的波束图;
图6(b)是本发明优化FDA的相位辐射图
图7(a)是本发明基本FDA的波束图;
图7(b)是本发明基本FDA的相位辐射图。
具体实施方式:
下面结合附图和具体实施方式对本发明进行详细说明。
本发明提出了一种基于频控阵相位中心的主动抗干扰的方法,从实现雷达主动抗干扰,提高战场雷达生存能力的实际军事需求出发,如图1所示,具体按照以下步骤实施:
步骤1、构造频控阵雷达系统;步骤1的作用为推导出FDA的信号模型包含波束方向图,相位方向图两部分内容,为后面相位中心的推导作理论准备。
步骤1的具体步骤为:
假设一部由N个阵元组成的均匀线阵频控阵雷达,如图2所示,其中频控阵阵元数为Q,C为光速,相邻两个阵元之间的频率增量为Δf,初始载频f0,第q个阵元处的频率为fq=f0+Δfq,阵元间距为d,目标位置为(θ’,R’),目标方位角为θ,目标与频控阵之间的距离为R。
第q个阵元发射的信号可以表示为:
Stq(t)=exp[j2πfq(t-Rq/c)] (1)
Aq表示第q个阵元辐射信号的幅度,在远场目标位置处的合成场强为:
在远场处,各阵元发射的信号几乎平行,即Rq=R-(q-1)d sinθ,对式(2)进行展开可得:
每个单元的信号辐射幅值是相等的,假设Aq=A=1,由于Δfq<<f0,Δfq(q-1)dsinθ/c可忽略,
为简化计算,假设:
at(t),aR(R),aθ(θ)分别包含每个元素的唯一包含t,R和θ的项,a0(t,R)包含每个元素的公共项。式(4)可以写作⊙表示阿达玛积。
通过推导FDA的全向合成信号模型,当需要探测位于(θ0,R0)的目标时,合成信号指向(θ0,R0),引入导向矢量:wR(R0)=a0(0,-R0)⊙aR(-R0),wθ(θ0)=aθ(-θ0),
wR(R0)和wθ(θ0)分别表示距离角度向量,则无方向性合成信号表示为:
EFDA(R,θ,R0,θ0,t)=(a0(t,R)⊙w)T·(at(t)⊙aθ(θ)⊙aR(R)) (6)
将场强用幅度和相位两部分来表示:
其相位信息:
ΨFDA(R,θ,R0,θ0,t)=angle(EFDA(R,θ,R0,θ0,t)) (8)
angle为相位角求解函数,AF=|EFDA(R,θ,R0,θ0,t)|为幅值,AF≥0。
步骤1根据基本FDA阵列结构推导出场强表达式,场强由幅度相位两部分组成,而敌方干扰机一般根据接收到的我方雷达辐射信号的信息来判断我方雷达的位置进而实施干扰。为了保护我方雷达,采用无线电静默的方式,但是当雷达不发射信号工作时,又无法探测目标,所以角度欺骗是更灵活的抗干扰手段。
目前,敌方干扰机探测我方雷达方向主要有两种方法:振幅法和相位法。虽然这两种方法的过程不同,但是这两种方法的原理是相似的,FDA雷达的信号是以球面波的形式辐射的,而在远场中,FDA信号的方向可以计算为同相波前的法向。实际上,法线方向和雷达所在的水平面会有一个交点,这个交点就是虚辐射源。对于球面波,虚拟辐射源也可以表示为相位中心。
步骤2、推导FDA辐射信号的相位中心;
如图3所示,假设相位中心M,坐标原点O,dM为从M到O的偏差向量,目标位于(θ’,R’),第q个阵元到场源的距离为:
Rq′=R′-(q-1)d sinθ+dM sinθ (9)
结合式(9)和式(4),合成信号的场强:
将其用导向矢量表示为:
E′FDA(R′,θ′,t)=a′0 T(t,R′,θ′)·(a′t(t)⊙a′θ(θ′)⊙a′R(R′)) (11)
为了探测远场(θ0,R0)处的目标,导向矢量为:
无方向性的合成信号为:
E′FDA(R′,θ′,R0,θ0,t)=(a′0(t,R′,θ′)⊙w′)T·(a′t(t)⊙a′θ(θ′)⊙a′R(R′)) (13)
AF′=|E′FDA(R′,θ′,R0,θ0,t)|表示以M为参考点的FDA信号幅度信息,上式写为:
假设FDA辐射的为球面波,以相位中心为参考点,在相同距离R′0处的相位为常数值:
C=Ψ′FDA(R′0,θ′,R0,θ0,t)=angle(E′FDA(R′0,θ′,R0,θ0,t)) (15)
结合式(5)和式(12),式(13)可以写为:
结合式(15)和式(16),式(15)可以写成:
式(17)是在假设线性FDA辐射球面波的基础上引入的,其实线性FDA辐射的波是准球面波。但是,方位角θ′在一个小范围内变化,而范围R′是固定的,所以式(17)成立。基于上述分析,当目标处于(θ′0,R′0),
当方位角变化δθ→0,
可以推出干扰机(θ′0,R′0)处的相位中心偏差为:
式(19)给出了相位中心相对于阵列原点的偏移距离,可以看出不同位置的相中心可能是不同的,由于相位方向图随不同雷达参数的变化而变化,所以雷达参数的变化也会导致相位中心的变化。
步骤2的作用为推导出相位中心的偏差与阵列原点,目标位置因素有关,并且由于ΨFDA随雷达参数的改变而改变,通过控制雷达参数来调节相位中心。
步骤3、通过调节频率增量,控制相位中心偏差;
从步骤2中的式(19)看出干扰机(θ′0,R′0)处的相位中心偏差与ΨFDA有关,而ΨFDA与频偏有关,通过调节FDA雷达的频率增量序列Δf=[Δf1 Δf2 … Δfq … ΔfQ],(θ′0,R′0)处的相位分布也会改变,进而改变相位中心dM(R′0,θ′0)。
在不考虑FDA雷达的波束多样性的情况下,假设导向矢量为:
无方向合成信号表示为:
EFDA(R′,θ′,t)=a0 T(t,R)·(at(t)⊙aθ(θ′)⊙aR(R′)) (21)
FDA的相位分布写为:
ΨFDA(R′,θ′,t)=angle(EFDA(R′,θ′,t)) (22)
对于固定时刻,t=0,式(22)写为:
通过10000次蒙特卡洛实验,证明选择合适的频偏可以使相位中心偏差较大。
步骤3中FDA雷达的功能包括探测目标和自我保护两个方面,因此,要求FDA雷达在保证探测能力的同时,对干扰源实施角度欺骗,在不影响正常工作的情况下实现主动抗干扰。
假设一个场景如图5所示,在固定时刻,要求雷达对干扰源实施主动抗干扰(角度欺骗),同时保持对目标的探测,从而在正常工作中达到自我保护的目的。如式(16)所示,导向矢量可以通过控制指向(θ0,R0),所以雷达可以探测固定在(θ0,R0)位置出的目标,同时对C处的干扰源实施主动抗干扰来防止自身被准确定位,因此,要求FDA雷达发射的电磁波到达C点处的相位中心与自身雷达阵列天线偏离的距离(假设用db表示)满足下式:
d′M(θ′0,R′0)≥db (24)
从式(19)得出,当时刻固定,t=0,假设目标固定在(θ0,R0),通过调节Δf使C处的相位中心满足式(20)。
d′M(Δf)≥db (26)
主动抗干扰问题转化为了最优化问题,
Δfmax和Δfmin分别为阵列的最大和最小频偏,最优化问题转化为最小值问题;
步骤3中的蒙特卡洛实验中设置的FDA载频为1GHZ,阵元间距d为0.15m,阵元数目为10,时间固定在零时刻,干扰机位置为(30°,250km)。结果如图4所示,可以看出当频偏改变时,大多数的相位中心偏差较小,但一旦找到了最优频偏,相位中心偏差将很大。因此证明了可以通过控制频偏来调节相位中心,从而使干扰机定位不准,进而实现主动抗干扰。
步骤4、通过一种免疫粒子群优化算法寻找最优频偏,使相位中心偏差最大。
步骤4中采用免疫粒子群优化算法,这种算法结合粒子群优化算法具有的全局寻优能力和免疫系统的免疫信息处理机制,实现简单,改善了粒子群优化算法摆脱局部极值点的能力,提高了算法进化过程中的收敛速度和精度。步骤4算法复杂度为O[DS·M·Q+(MaxT-DS)M(Q+N)],而步骤3中的复杂度为O{[(Δfmax-Δfmin)/δ]Q},δ为寻找次数,复杂度大大降低,步骤4解决了优化表达式中的周期性问题,提高了求解速度和精度,降低了运算复杂度。
实施例1
仿真实例:设时间固定,载频f0为3GHz,最大频偏Δfmax为10KHz,最小频偏Δfmin为0,阵元数目Q为10,阵元间距d为0.05m,光速c为3×108m/s,相位中心的最小偏差距离db为100m,惩罚因子r为109,学习因子c1,c2为2,粒子数M为50,惯性权重ω为0.8,最大迭代次数MaxDT为500,检验间隔DS为10,目标与干扰机的位置分别为(45°,200km),(30°,250km),Δf1=0,在进行第297次迭代时,达到最大的相位偏差119.567m。频率增量序列Δf=(06.74494.39005.32902.83485.09957.11355.08524.20073.5184)KHz。均在阵列的最大频率增量和最小频率增量之间,满足优化要求。
根据改进的免疫粒子群优化算法得到频率增量序列,FDA的波束图和相位辐射图如图6所示,称其为优化FDA。如图7所示,给出了频率递增序列下FDA的波束图和相位辐射图,Δf=(0123456789)KHz,称其为基本FDA,图6(a)和图7(a)分别为优化FDA和基本FDA的波束图,目标和干扰源的位置如图所示,对比图6(a)和图7(a)可以看出,优化FDA和基本FDA都能使固定目标处于更高的波束能量下,而只有优化FDA能使固定干扰源在更低的波束能量下。这说明优化FDA和基本FDA都具有指向固定目标的能力,但由于接收到更多的波束能量,基本FDA比优化FDA更容易被固定干扰源检测到。此外,基本FDA的波束在距离约500公里处也是主波束,这意味着存在距离模糊,图6(b)和图7(b)分别为优化FDA和基本FDA的相位辐射图,对比图6(b)和图7(b)可以看出,优化FDA和基本FDA都能使固定目标处于一个规律的相位分布,而只有优化FDA才能使固定干扰源处于畸变的相位分布。这说明优化FDA和基本FDA由于相位分布规律,都能检测到固定的目标,但由于相位分布畸变,基本FDA可以在干扰源处实现相位中心偏移,从而实现主动抗干扰。综合分析图6和图7,发现优化FDA在保证目标探测能力的同时,通过降低波束能量和畸变相位分布,实现对干扰源的主动抗干扰。
Claims (3)
1.一种基于频控阵相位中心的主动抗干扰的方法,其特征在于,具体按照以下步骤实施:
步骤1、构造频控阵雷达系统;
步骤2、推导FDA辐射信号的相位中心;
步骤3、通过调节频率增量,控制相位中心偏差;
步骤4、通过一种免疫粒子群优化算法寻找最优频偏,使相位中心偏差最大;
所述步骤1的具体步骤为:
假设一部均匀线阵频控阵雷达,其中频控阵阵元数为Q,C为光速,相邻两个阵元之间的频率增量为Δf,初始载频f0,第q个阵元处的频率为fq=f0+Δfq,阵元间距为d,目标位置为(θ’,R’),目标方位角为θ,目标与频控阵之间的距离为R;第q个阵元发射的信号表示为:
Stq(t)=exp[j2πfq(t-Rq/c)] (1)
Aq表示第q个阵元辐射信号的幅度,在远场目标位置处的合成场强为:
在远场处,各阵元发射的信号几乎平行,即Rq=R-(q-1)dsinθ,对式(2)进行展开可得:
每个单元的信号辐射幅值是相等的,假设Aq=A=1,由于Δfq<<f0,Δfq(q-1)d sinθ/c可忽略,
为简化计算,假设:
at(t),aR(R),aθ(θ)分别包含每个元素的唯一包含t,R和θ的项,a0(t,R)包含每个元素的公共项;式(4)写作⊙表示阿达玛积;
通过推导FDA的全向合成信号模型,当需要探测位于(θ0,R0)的目标时,合成信号指向(θ0,R0),引入导向矢量:wR(R0)=a0(0,-R0)⊙aR(-R0),wθ(θ0)=aθ(-θ0),
wR(R0)和wθ(θ0)分别表示距离角度向量,则无方向性合成信号表示为:
EFDA(R,θ,R0,θ0,t)=(a0(t,R)⊙w)T·(at(t)⊙aθ(θ)⊙aR(R)) (6)
W为距离向量wR(R0)与角度向量wθ(θ0)的阿达码积;
将场强用幅度和相位两部分来表示:
其相位信息:
ΨFDA(R,θ,R0,θ0,t)=angle(EFDA(R,θ,R0,θ0,t)) (8)
angle为相位角求解函数,AF=|EFDA(R,θ,R0,θ0,t)|为幅值,AF≥0;
所述步骤2假设相位中心M,坐标原点O,dM为从M到O的偏差向量,目标位于(θ’,R’),第q个阵元到场源的距离为:
Rq′=R′-(q-1)d sinθ+dMsinθ (9)
结合式(9)和式(4),合成信号的场强:
类比式(6),目标(θ’,R’)处的电场强度将其用导向矢量表示为:
为了探测远场(θ0,R0)处的目标,导向矢量为:
无方向性的合成信号为:
E′FDA(R′,θ′,R0,θ0,t)=(a0′(t,R′,θ′)⊙w′)T·(a′t(t)⊙a′θ(θ′)⊙a′R(R′)) (13)
AF′=E′FDA(R′,θ′,R0,θ0,t)表示以M为参考点的FDA信号幅度信息,上式写为:
假设FDA辐射的为球面波,以相位中心为参考点,在相同距离R0′处的相位为常数值:
C=Ψ′FDA(R′0,θ′,R0,θ0,t)=angle(E′FDA(R′0,θ′,R0,θ0,t)) (15)
结合式(5)和式(12),式(13)写为:
结合式(15)和式(16),式(15)写成:
式(17)是在假设线性FDA辐射球面波的基础上引入的,其实线性FDA辐射的波是准球面波;但是,方位角θ′在一个小范围内变化,而范围R′是固定的,所以式(17)成立,基于上述分析,当目标处于(θ′0,R′0),
当方位角变化δθ→0,
推出干扰机(θ′0,R′0)处的相位中心偏差为:
式(19)给出了相位中心相对于阵列原点的偏移距离,看出不同位置的相位中心可能是不同的,由于相位方向图随不同雷达参数的变化而变化,所以雷达参数的变化也会导致相位中心的变化;
所述从步骤2中的式(19)看出干扰机(θ′0,R′0)处的相位中心偏差与ΨFDA有关,而ΨFDA与频偏有关,通过调节FDA雷达的频率增量序列Δf=[Δf1 Δf2…Δfq…ΔfQ],(θ′0,R′0)处的相位分布也会改变,进而改变相位中心dM(R′0,θ′0);
在不考虑FDA雷达的波束多样性的情况下,假设导向矢量为:
无方向合成信号表示为:
EFDA(R′,θ′,t)=a0 T(t,R)·(at(t)⊙aθ(θ′)⊙aR(R′)) (21)
FDA的相位分布写为:
ΨFDA(R′,θ′,t)=angle(EFDA(R′,θ′,t)) (22)
对于固定时刻,t=0,式(22)写为:
通过10000次蒙特卡洛实验,证明选择合适的频偏使相位中心偏差较大;
所述步骤3假设一个场景,在固定时刻,要求雷达对干扰源实施主动抗干扰,同时保持对目标的探测,从而在正常工作中达到自我保护的目的;如式(16)所示,导向矢量通过控制指向(θ0,R0),所以雷达探测固定在(θ0,R0)位置出的目标,同时对C处的干扰源实施主动抗干扰来防止自身被准确定位,因此,要求FDA雷达发射的电磁波到达C点处的相位中心与自身雷达阵列天线偏离的距离db满足下式:
d′M(θ′0,R′0)≥db (24)
从式(19)得出,当时刻固定,t=0,假设目标固定在(θ0,R0),通过调节Δf使C处的相位中心满足式(20);
d′M(Δf)≥db (26)
主动抗干扰问题转化为了最优化问题,
Δfmax和Δfmin分别为阵列的最大和最小频偏,最优化问题转化为最小值问题;
2.根据权利要求1所述的一种基于频控阵相位中心的主动抗干扰的方法,其特征在于,所述步骤3中的蒙特卡洛实验中设置的FDA载频为1GHZ,阵元间距d为0.15m,阵元数目为10,时间固定在零时刻,干扰机位置为(30°,250km)。
3.根据权利要求1所述的一种基于频控阵相位中心的主动抗干扰的方法,其特征在于,所述步骤4中采用免疫粒子群优化算法,步骤4算法复杂度为O[DS·M·Q+(MaxDT-DS)·M·Q];其中,DS为检验间隔,M为粒子数,Q为阵元数目,MaxDT为最大迭代次数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011380267.XA CN112612019B (zh) | 2020-11-30 | 2020-11-30 | 一种基于频控阵相位中心的主动抗干扰的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011380267.XA CN112612019B (zh) | 2020-11-30 | 2020-11-30 | 一种基于频控阵相位中心的主动抗干扰的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112612019A CN112612019A (zh) | 2021-04-06 |
CN112612019B true CN112612019B (zh) | 2024-02-02 |
Family
ID=75229785
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011380267.XA Active CN112612019B (zh) | 2020-11-30 | 2020-11-30 | 一种基于频控阵相位中心的主动抗干扰的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112612019B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113791392B (zh) * | 2021-08-05 | 2023-11-07 | 中国人民解放军空军工程大学 | 一种基于频控阵对全向比幅单脉冲测向系统的欺骗方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105589070A (zh) * | 2015-12-11 | 2016-05-18 | 桂林电子科技大学 | 基于频率分集阵列的雷达目标成像方法 |
CN107037410A (zh) * | 2017-04-17 | 2017-08-11 | 电子科技大学 | 一种对雷达进行干扰的方法、装置及频控阵干扰机 |
CN108761414A (zh) * | 2018-05-25 | 2018-11-06 | 电子科技大学 | 一种基于频控阵的s形干扰波束测试方法 |
CN109143275A (zh) * | 2018-07-30 | 2019-01-04 | 中国电子科技集团公司第五十四研究所 | 一种基于粒子群的小型化阵列天线抗干扰实现方法 |
CN109768584A (zh) * | 2018-12-24 | 2019-05-17 | 深圳供电局有限公司 | 基于免疫粒子群算法的微网自治调频调压方法和装置 |
CN110109069A (zh) * | 2019-05-14 | 2019-08-09 | 桂林电子科技大学 | 一种去时间依赖性频控阵点状干扰波束形成方法 |
CN110471040A (zh) * | 2019-08-22 | 2019-11-19 | 电子科技大学 | 一种基于fda天线的逆合成孔径雷达干扰方法 |
CN111830482A (zh) * | 2020-07-10 | 2020-10-27 | 西安电子科技大学 | 基于捷变ofdm的fda雷达目标定位方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB201803239D0 (en) * | 2018-02-28 | 2018-04-11 | Secr Defence | A radio or sonic wave detector, transmitter, reciver and method thereof |
US10955543B2 (en) * | 2019-12-20 | 2021-03-23 | Intel Corporation | Station detection using frequency-diverse array signals |
-
2020
- 2020-11-30 CN CN202011380267.XA patent/CN112612019B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105589070A (zh) * | 2015-12-11 | 2016-05-18 | 桂林电子科技大学 | 基于频率分集阵列的雷达目标成像方法 |
CN107037410A (zh) * | 2017-04-17 | 2017-08-11 | 电子科技大学 | 一种对雷达进行干扰的方法、装置及频控阵干扰机 |
CN108761414A (zh) * | 2018-05-25 | 2018-11-06 | 电子科技大学 | 一种基于频控阵的s形干扰波束测试方法 |
CN109143275A (zh) * | 2018-07-30 | 2019-01-04 | 中国电子科技集团公司第五十四研究所 | 一种基于粒子群的小型化阵列天线抗干扰实现方法 |
CN109768584A (zh) * | 2018-12-24 | 2019-05-17 | 深圳供电局有限公司 | 基于免疫粒子群算法的微网自治调频调压方法和装置 |
CN110109069A (zh) * | 2019-05-14 | 2019-08-09 | 桂林电子科技大学 | 一种去时间依赖性频控阵点状干扰波束形成方法 |
CN110471040A (zh) * | 2019-08-22 | 2019-11-19 | 电子科技大学 | 一种基于fda天线的逆合成孔径雷达干扰方法 |
CN111830482A (zh) * | 2020-07-10 | 2020-10-27 | 西安电子科技大学 | 基于捷变ofdm的fda雷达目标定位方法 |
Non-Patent Citations (3)
Title |
---|
Guangming Li 等.An Anti-jamming Method of ISAR Imaging with FDA-MIMO Radar .AND FDA.2018,全文. * |
王博 ; 谢军伟 ; 葛佳昂 ; 张晶 ; .FDA发射干扰机对无源雷达干涉仪测向系统的欺骗机理.华南理工大学学报(自然科学版).2020,(第01期),全文. * |
王博 ; 谢军伟 ; 葛佳昂 ; 张晶 ; .FDA对测向时差组合定位的欺骗研究.空军工程大学学报(自然科学版).2019,(第06期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN112612019A (zh) | 2021-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7504982B2 (en) | Anti-Missile system and method | |
US7961133B2 (en) | System and method for diverting a guided missile | |
Zhou et al. | Performance evaluation of radar and decoy system counteracting antiradiation missile | |
CN108614261B (zh) | 一种组网雷达系统多目标跟踪下的辐射参数控制方法 | |
CN108196231B (zh) | 一种基于新型频控阵技术的s形干扰波束实现方法 | |
Ge et al. | A cognitive active anti-jamming method based on frequency diverse array radar phase center | |
CN112949100A (zh) | 一种机载雷达抗主瓣干扰的方法 | |
CN112612019B (zh) | 一种基于频控阵相位中心的主动抗干扰的方法 | |
CN111352077B (zh) | 基于多比例分式规划的低截获频控阵mimo雷达系统的设计方法 | |
CN110133604A (zh) | 一种基于多元合成技术的机载自卫欺骗式干扰方法 | |
CN111175731A (zh) | 一种多雷达目标单机无源定位方法 | |
CN109324507B (zh) | 针对隐身目标动态rcs的雷达发射功率自适应控制方法 | |
CN115063544A (zh) | 一种面向随队支援干扰的雷达探测威力范围计算方法 | |
Zhouyan et al. | Research on track deception technology based on jamming of flexible lines by jammers | |
CN111427015A (zh) | 基于飞行载体的四源非均匀线阵反向交叉眼干扰方法 | |
CN116125400A (zh) | 一种基于对消隐身的同步慢闪烁角度欺骗干扰方法 | |
CN110109069A (zh) | 一种去时间依赖性频控阵点状干扰波束形成方法 | |
CN108896970A (zh) | 一种基于多元矢量合成技术的雷达干扰方法 | |
CN111090079B (zh) | 基于无源传感器协同的雷达组网辐射间隔优化控制方法 | |
Zhang et al. | Performance evaluation for UAV-based distributed jamming system: an illustrative example | |
US11662428B2 (en) | Radar system for jammed environments | |
Zhao et al. | Frequency diversity array MIMO track-before-detect in coherent repeated interference | |
Song et al. | Research on anti-deceptive jamming of multi station fusion of netted radar | |
CN117491954A (zh) | 三机编队闪烁干扰抑制地面单脉冲雷达跟踪性能的方法 | |
RU190395U1 (ru) | Радиолокационная станция кругового обзора |
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 |