CN101839975B - 警戒雷达抗同型干扰方法 - Google Patents
警戒雷达抗同型干扰方法 Download PDFInfo
- Publication number
- CN101839975B CN101839975B CN2010101577164A CN201010157716A CN101839975B CN 101839975 B CN101839975 B CN 101839975B CN 2010101577164 A CN2010101577164 A CN 2010101577164A CN 201010157716 A CN201010157716 A CN 201010157716A CN 101839975 B CN101839975 B CN 101839975B
- Authority
- CN
- China
- Prior art keywords
- antibody
- radar
- homotype
- sigma
- fitness
- 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.)
- Expired - Fee Related
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公布了一种基于克隆选择的警戒雷达抗同型干扰方法。该方法基于克隆选择框架构造启发式搜索,求解一个有约束的整数二次规划问题,从任意的初始解出发,通过克隆、变异、选择逐代进化,直至达到最大进化代数或最优解无改善,得到各部雷达选用各套参数配置的最优频度,设计最优工作序列表,实现抗同型干扰同时不降低每部雷达性能的目标。具有良好的普适性、稳健性和较高的效率,是一种实际可行、可靠、高效的方法,并且程序简单,容易在嵌入式数字处理平台上实现,具有较高的实用价值。
Description
技术领域
本发明属于雷达与电子对抗技术领域,尤其涉及一种基于克隆选择的警戒雷达抗同型干扰方法。
背景技术
干扰与抗干扰问题自雷达诞生之日起就是雷达与电子对抗技术领域最为重要的课题之一,相关的技术一直伴随着雷达自身的发展而发展。现代空间电磁环境错综复杂、瞬息万变,雷达面对的干扰源和干扰类型也日趋多样,其中同型雷达相互间干扰(简称为同型干扰)是较为严重的一种。同型干扰是多部同类型雷达在近距离内同时工作带来的相互间干扰。雷达是通过发射一定形式的电磁波、接收目标回波来工作的,当多部雷达在近距离内同时工作、发射波形参数相同或者近似的电磁波时,一部雷达发射的电磁波会被周围雷达接收,同时也会接收到周围雷达发射的电磁波,这些电磁波的强度可能大于甚至远远大于目标回波的强度,掩盖甚至完全淹没目标,从而形成严重的同型干扰。同型干扰分为同步干扰和异步干扰,在雷达显示画面上,前者表现为同心圆,后者表现为螺旋线,这些干扰可使雷达发现和跟踪目标的能力大大降低乃至完全失效。特别是在舰艇编队中,由于同型舰艇往往配备同型雷达,并且舰艇之间通常相距较近,同型干扰尤为常见,一般来说编队规模越大,同型雷达数量越多,同型干扰越严重。
为使多部同型雷达能够在近距离内同时正常工作,必须设法尽可能地减小同型干扰。几十年来,各种抗同型干扰方法在国内外雷达与电子对抗技术领域的专著中不断被提出和改进。国内南京船舶雷达研究所、海军大连舰艇学院等单位开展了相关研究,较有代表性的如陈正禄等于2006年、2008年、2009年在《雷达与对抗》期刊上发表的文章。目前,较为成熟、应用较多的抑制或消除同型干扰的办法大都可以概括为以下两类:
1)不同雷达采用不同波形参数,在接收机中配合采用相应的信号处理技术,实现抑制干扰、增强目标。这些技术通常关注的波形参数有雷达工作频率(载频)、脉冲重复频率(重频)、调频斜率(对于线性调频信号)等。例如,若不同雷达采用不同载频,则接收机采用相应的带通滤波;若不同雷达采用不同重频,则接收机采用相应的相关积累;若不同雷达采用不同调频斜率,则接收机采用相应的匹配滤波。
2)各部雷达采用时分或空分方式协同工作,从源头上减少干扰的产生。这类办法通常着眼于同时、近距、天线相互指向等造成强干扰的关键因素,有针对性地予以去除。例如,采用时分工作方式,即各部雷达根据一定的顺序规则依次工作,保证任一时刻只有一部雷达处于工作状态;采用空分工作方式,即增大各部雷达之间的相互距离、减小干扰强度,或者使用低副瓣天线(或应用抑制天线副瓣的信号处理技术)、通过空间分集避免天线主瓣对指。
上述措施作为雷达抗同型干扰方法、特别是对于舰艇编队雷达抗同型干扰来说或多或少都存在问题。错开载频结合带通滤波可以消除或减少进入接收机的干扰能量,但是很大程度上限制了雷达性能的发挥,同时需要占用较多的频谱资源,并且由于实际滤波器非理想、系统非线性等因素,干扰很难彻底滤除;错开重频结合相关积累可以先将同步干扰变为异步干扰,再进一步通过方位向反异步处理消除异步干扰,但是这种处理影响到系统相参积累,将导致系统信噪比和雷达探测概率的下降;不同雷达采用不同的调制编码方式或参数,如线性调频下采用不同的调频斜率等,一般都涉及比较复杂的技术,需要重新设计装备或对现役装备进行较大的技术改造;错开各部雷达的工作时间或拉开各部雷达的相互距离均与雷达协同、组网工作的战术要求相悖,大大降低系统效能,这些措施实际上是在回避问题而没有解决问题;使用低副瓣天线(或应用抑制天线副瓣的信号处理技术),通过空间分集避免天线主瓣对指,需要各部雷达接受统一调度,即失去工作的独立性,同时也需要重新设计天线,不适合现役雷达。
现役雷达通常都具备多套可选的参数配置,如多个工作频点,上述抗同型干扰方法一般都要限制特定雷达只能在其中的特定范围内选择,这使得每部雷达的性能都得不到充分发挥,特别是参数捷变体制的雷达(如频率捷变雷达)工作时需要在各套参数配置间捷变,限制越严则其性能损失越大。为充分发挥每部雷达的性能,应允许每部雷达独立工作,并且能够充分利用各套可选的参数配置。只要各部雷达在选用各套参数配置时尽量避开形成强同型干扰的条件,即可实现抗同型干扰同时不降低每部雷达性能的目标。因此,需要为各部雷达设计最优的工作序列,即选用各套参数配置的次序表,使得平均同型干扰最小。文献检索表明,目前国内外关于这个问题还没有专门的论述,特别是当系统中雷达数量和可选参数配置数目较多时,上述问题的解空间很大,目前还没有提出有效的计算方法,如果直接进行全局的搜索,计算效率非常低。
人工免疫系统是模拟生物免疫系统构建的一种新型人工智能工具。1959年,F.M.Burnet在《The Clonal Selection Theory of Acquired Immunity》一书中提出了著名的克隆选择学说。克隆选择理论认为,抗体的大量复制、变异和抗原对抗体的选择是生物免疫系统进化的基本过程,在一代代克隆选择的过程中,优秀的抗体不断产生并被保留下来。利用基于这个框架构造的启发式搜索求解非线性优化问题,具有稳健快速地收敛到全局最优解的性能。上述设计雷达最优工作序列表的问题可以建模为一个有约束的整数二次规划问题,因此可以基于克隆选择框架构造启发式搜索进行求解,获得使平均同型干扰最小的雷达最优工作序列表,实现抗同型干扰同时不降低每部雷达性能的目标。
经过文献检索,没有发现克隆选择方法在雷达抗同型干扰方面应用的论文。
经过对美国专利商标局USPTO、欧洲专利局EPO和日本专利局JPO的检索,没有找到类似的专利。
发明内容
本发明给出一种基于克隆选择的警戒雷达抗同型干扰方法。该方法基于克隆选择框架构造启发式搜索,求解一个有约束的整数二次规划问题,获得使平均同型干扰最小的雷达最优工作序列表(即选用各套参数配置的最优次序表),实现抗同型干扰同时不降低每部雷达性能的目标。
设系统中有P部雷达同时工作,其中第p(p=1,2,...,P)部雷达有Qp个可变参数,第q(q=1,2,...,Qp)个参数有个可能值,则理论上第p部雷达最多有套参数配置(实际上由于雷达各个参数之间的配合需要,可用的参数配置可能少于这个数目)。用Mp表示第p部雷达实际可选参数配置的数目,理论分析表明,在各部雷达相互独立工作的情况下,平均同型干扰与每部雷达选用每套参数配置的频度有关,而与具体次序无关,归一化的平均同型干扰可以表示为
或
其中,表示雷达p1,p2分别采用各自的第i,j套参数配置时雷达p1对p2的同型干扰(归一化至[0,1]区间),式(1)中的分别表示雷达p1,p2选用各自的第i,j套参数配置的频度,式(2)中的N,为整数,分别为的离散化。因此设计使平均同型干扰最小的雷达最优工作序列表,关键在于确定各部雷达选用各套参数配置的最优频度,即使式(1)或式(2)最小的或(p=1,2,...,P,i=1,2,...,Mp)。基于克隆选择框架构造启发式搜索,求解以最小化式(2)为目标的有约束整数二次规划问题,具体步骤如下:
1.根据实际情况确定优化目标式(2)中的各个参量,包括系统中雷达数量P、各部雷达可选参数配置数目Mp(p=1,2,...,P)、频度量化数N以及归一化同型干扰系数(与两部雷达本身的特性和各自选用的参数配置等因素有关,应经过实际测试确定并归一化至[0,1]区间)。
此外还应包括其他由实际情况决定的约束条件,例如选用每套参数配置的频度均不低于某个给定值等。
3.任意给定一组满足约束条件的初始解,该组解称为第1代抗体群,记为常数NA为抗体群规模,即其中抗体的数目,应根据实际情况选定,一般来说抗体群规模越大,越有利于搜索到全局最优解,但是每一代计算时间也越长。每一代抗体群规模都是相同的。
并记录其最大值。抗体适应度越大,即平均归一化同型干扰越小。
并构造克隆抗体群适应度越大的抗体被克隆得越多,意味着在较好的解周围进行较多搜索,有利于提高计算效率。常数NC为克隆抗体群规模,也应根据实际情况选定,一般来说克隆抗体群规模越大,越有利于搜索到局部最优解,但是每一代计算时间也越长。每一代克隆抗体群规模都是相同的。由于式(5)中使用了向上取整操作克隆抗体总可能略大于NC。
6.对克隆抗体群中的每个抗体
随机选取p0∈{1,2,...,P}并随机选构造变异抗体
满足
其中,常数α>0控制选择的随机性,一般来说其值越小越有利于搜索到局部最优解,其值越大越有利于跳出局部最优解。
9.记录第k+1代抗体群的适应度最大值。
10.k:=k+1。若k=Kmax,或连续Kstop代抗体群适应度最大值没有改善,则终止计算,当前代抗体群中适应度最大的抗体即为求得的最优解;否则返回步骤5继续进化。
确定了各部雷达选用各套参数配置的最优频度之后,可以任意安排具体选用次序、设计最优工作序列表,平均同型干扰不变。例如,若从反辐射源信号侦察的角度考虑,雷达参数应具有较强的随机性,因此最优工作序列表可以按照完全随机排列的方式生成,这相当于各部雷达按照与最优频度相一致的概率分布随机选用各套参数配置。
与现有技术相比,本发明的优点如下:
1)采用本发明设计得到最优工作序列表,能够保证各部雷达在满足一定频度要求的条件下、充分利用系统参数配置独立工作,充分发挥每部雷达的性能,同时最大限度地减小同型干扰。
2)本发明建立的模型具有一般性,对雷达数量、类型、参数种类均没有限制,且可与现有技术兼容,普适性较强。
3)本发明提出基于克隆选择求解雷达抗同型干扰最优工作序列表的方法,充分发挥启发式搜索算法稳健快速地收敛到全局最优解的性能,计算效率高。
具体实施方式
以下以某型号舰载远程警戒雷达为例,说明本发明方法的有效性,具体如下:
某编队装备了2部某型号舰载远程警戒雷达,相互独立工作,采用相同的工作序列表。雷达载频可以在8个频点(由低到高依次编号为频点1,2,...,8)上捷变,每经过固定的时间间隔跳频一次,跳频表长度为100个时间间隔(跳频100次)。当两部雷达工作在任意相同频点上时,相互之间的归一化同型干扰近似为1;当两部雷达工作在任意相邻频点上时,相互之间的归一化同型干扰近似为0.5;其他情况下两部雷达相互之间无同型干扰。需要设计最优工作序列表(即最优跳频表),使得平均归一化同型干扰最小。设计步骤如下:
1.确定系统参量:雷达数量P=2,频度量化数即跳频表长度N=100,可选参数配置即可选工作频点数目M1=M2=8,由于两部雷达完全相同,相互同型干扰也相同,即归一化同型干扰系数满足和因此可简记为一个对称矩阵
优化目标式(2)可简化为
3.初始化第1代抗体群:
n=1,2,...,50
即所有频点大致平均选用,共50个这样的抗体。
4.用式(12)及式(4)计算第1代抗体群中每个抗体对应的平均归一化同型干扰及其适应度(所有抗体均相同)。
6.变异:对于克隆抗体群中的每个抗体[n1,...,n8],随机选取i,j∈{1,2,...,8},在和不变的条件下随机重新分配ni和nj。
10.k:=k+1。若k=50,或连续5代抗体群适应度最大值没有改善,则终止计算,当前代抗体群中适应度最大的抗体即为求得的最优解;否则返回步骤5继续进化。
采用上述方法求得的最优解为[20 5 15 10 10 15 5 20],即在长度为100的跳频表中频点1,2,...,8出现的次数依次为20,5,15,10,10,15,5,20次时平均同型干扰最小,这与全局搜索得到的结果一致。按照完全随机排列的方式生成的最优跳频表见下表。
表1某型号舰载远程警戒雷达最优跳频表
为验证上述方法的效率和稳健性,进行了20次重复试验,正确(求得上述最优解)率为100%,进化代数介于15至20代之间,每一代大约搜索1000个解,总计搜索不到20000个解,而如果直接进行全局的搜索,根据有重复组合数计算公式,需要搜索个解。
在采用本发明之前,编队一般采用时分和频分的方法控制多部同型警戒雷达的工作。如果采用时分方法,整个编队同一时刻只能有一部雷达工作,不能充分发挥全部雷达的功能;如果采用频分方法,每部雷达只能在个别频点上工作,不能充分发挥每部雷达的性能。采用本发明给出的最优跳频表,各部雷达可以同时独立工作,大大提高了编队探测防区空域的数据更新率,每部雷达的性能也得到了充分发挥。
由上面的实例可知,本发明给出的方法具有较高的效率和良好的稳健性,是一种实际可行、可靠、高效的方法,并且该方法程序简单,容易在嵌入式数字处理平台上实现,具有较高的实用价值。
Claims (3)
1.一种基于克隆选择的警戒雷达抗同型干扰方法,其特征在于,包括以下步骤:
步骤1:将归一化的平均同型干扰式(1)作为优化目标式,根据实际情况确定各个参量,
其中,P为系统中雷达数量,各部雷达可选参数配置数目为Mp(p=1,2,...,P),N为频度量化数,为归一化同型干扰系数,该系数需经过实际测试确定并归一化至[0,1]区间,和为整数,分别为的离散化,分别表示雷达p1,p2选用各自的第i,j套参数配置的频度;
步骤4:用式(1)计算所述第1代抗体群中每个抗体对应的平均归一化同型干扰计算抗体适应度
并记录其最大值;
并构造克隆抗体群 其中,常数NC为克隆抗体群规模;
步骤6:对克隆抗体群中的每个抗体
满足
在和不变的条件下两个分量的重新分配也是随机的;
其中,常数α>0控制选择的随机性;
步骤10:k:=k+1,若k=Kmax,或连续Kstop代抗体群适应度最大值没有改善,则终止计算,当前代抗体群中适应度最大的抗体即为求得的最优解;否则返回步骤5继续进化。
2.如权利要求1所述的方法,其特征在于,若所述步骤2中还列出了解的其他约束条件,则所述步骤6中构造变异抗体时应同时予以考虑,保证变异抗体对应可行解。
3.如权利要求1所述的方法,其特征在于,确定各部雷达选用各套参数配置的最优频度之后,按照完全随机排列的方式生成最优工作序列表,平均同型干扰不变。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101577164A CN101839975B (zh) | 2010-04-21 | 2010-04-21 | 警戒雷达抗同型干扰方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101577164A CN101839975B (zh) | 2010-04-21 | 2010-04-21 | 警戒雷达抗同型干扰方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101839975A CN101839975A (zh) | 2010-09-22 |
CN101839975B true CN101839975B (zh) | 2012-07-25 |
Family
ID=42743475
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101577164A Expired - Fee Related CN101839975B (zh) | 2010-04-21 | 2010-04-21 | 警戒雷达抗同型干扰方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101839975B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106714196A (zh) * | 2016-12-26 | 2017-05-24 | 全球能源互联网研究院 | 一种无线通信设备部署方法及装置 |
CN112204417A (zh) * | 2019-10-09 | 2021-01-08 | 深圳市大疆创新科技有限公司 | 雷达的抗干扰方法、设备、系统及存储介质 |
CN112363121B (zh) * | 2020-11-03 | 2023-03-03 | 中国人民解放军海军航空大学 | 基于回波预处理和相参积累的c&i干扰抑制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4068233A (en) * | 1976-08-13 | 1978-01-10 | Raytheon Company | Radar system having interference rejection |
CN2118313U (zh) * | 1991-10-28 | 1992-10-07 | 四川大学 | 同频异步干扰脉冲抑制器 |
EP1983353A1 (en) * | 2007-04-20 | 2008-10-22 | IDS Ingegneria Dei Sistemi S.p.A. | Radar method and device with verification of presence of active services in the band of frequency |
-
2010
- 2010-04-21 CN CN2010101577164A patent/CN101839975B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4068233A (en) * | 1976-08-13 | 1978-01-10 | Raytheon Company | Radar system having interference rejection |
CN2118313U (zh) * | 1991-10-28 | 1992-10-07 | 四川大学 | 同频异步干扰脉冲抑制器 |
EP1983353A1 (en) * | 2007-04-20 | 2008-10-22 | IDS Ingegneria Dei Sistemi S.p.A. | Radar method and device with verification of presence of active services in the band of frequency |
Non-Patent Citations (4)
Title |
---|
刘冬利等.时域多脉冲相关法抗雷达同频干扰.《现代雷达》.2009,第31卷(第06期),12-14. * |
薛春祥等.舰船雷达同频干扰来源分析及抗同频干扰的方法.《雷达与对抗》.2008,(第01期),1-4,9. * |
邱卫新.一种有效的反同型雷达相互间干扰的方法.《雷达与对抗》.2002,(第03期),5-7,30. * |
陈正禄等.舰载脉冲压缩体制雷达的抗同频干扰技术研究.《雷达与对抗》.2006,(第02期),19-22,36. * |
Also Published As
Publication number | Publication date |
---|---|
CN101839975A (zh) | 2010-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107451393B (zh) | 基于随机变邻域搜索算法的护士排班方法和装置 | |
CN101839975B (zh) | 警戒雷达抗同型干扰方法 | |
CN1174259C (zh) | 传送压缩地震探测数据的方法和系统 | |
CN109194508A (zh) | 基于区块链的数据处理方法和装置 | |
CN106656308B (zh) | 空间信息网络中任务的规划方法和装置 | |
CN105911532B (zh) | 基于深度协同的合成孔径雷达回波并行模拟方法 | |
US20220308249A1 (en) | Method and system for generating logarithmic non-uniform pseudo-random electromagnetic exploration signal | |
CN104063749A (zh) | 一种基于滚动时域控制的成像卫星自主任务规划算法 | |
CN105046412A (zh) | 一种被动相控阵雷达多站联合资源调度与分配方法 | |
CN106951991A (zh) | 一种核电站检修路径智能规划系统及其应用方法 | |
CN103036251B (zh) | 一种提高电网风电接纳能力的源网协调调度方法 | |
CN105372723A (zh) | 基于卷积神经网络模型的太阳耀斑预报方法 | |
CN109904878A (zh) | 一种多风电场发电时序模拟场景构建方法 | |
Ou et al. | An interval chance-constrained fuzzy modeling approach for supporting land-use planning and eco-environment planning at a watershed level | |
CN106789347A (zh) | 一种基于告警数据实现告警关联和网络故障诊断的方法 | |
CN109343059A (zh) | 一种正交非线性调频信号生成方法及装置 | |
CN106022652A (zh) | 一种森林碳汇经营方案的处理方法和装置 | |
Nozu | A super asperity model for the 2011 off the Pacific coast of Tohoku earthquake | |
CN104376363A (zh) | 一种基于改进的免疫遗传算法的多相正交码生成方法 | |
CN109038631A (zh) | 一种分布式发电和储能的容量配置方法 | |
Prokopenko et al. | Optimization of use of resource in multi-position radar systems | |
CN109752758B (zh) | 一种地震数据分解方法、系统及存储介质与终端 | |
Tran et al. | Structures of association rule set | |
CN105117470B (zh) | 数据存储方法和装置 | |
CN102253366B (zh) | 基于信道脉冲响应自相关序列的定位方法 |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120725 Termination date: 20140421 |