CN107944597A - 一种面对先进无源探测系统的编队雷达资源管理方法 - Google Patents

一种面对先进无源探测系统的编队雷达资源管理方法 Download PDF

Info

Publication number
CN107944597A
CN107944597A CN201710996643.XA CN201710996643A CN107944597A CN 107944597 A CN107944597 A CN 107944597A CN 201710996643 A CN201710996643 A CN 201710996643A CN 107944597 A CN107944597 A CN 107944597A
Authority
CN
China
Prior art keywords
mrow
msub
radar
station
mtd
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.)
Granted
Application number
CN201710996643.XA
Other languages
English (en)
Other versions
CN107944597B (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201710996643.XA priority Critical patent/CN107944597B/zh
Publication of CN107944597A publication Critical patent/CN107944597A/zh
Application granted granted Critical
Publication of CN107944597B publication Critical patent/CN107944597B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Educational Administration (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种面对敌方先进无源探测系统的编队雷达资源管理方法,首先根据典型敌方先进无源探测系统和编队雷达系统性能,确定典型敌方系统的截获概率建模方法;确定我方每部雷达发射功率和驻留时间的上、下界;以最优化我方编队雷达系统在功率域与时间域的被截获性能为目标,构建最优化被截获性能的编队雷达驻留时间与功率资源联合优化模型,经数值计算,再得到各雷达驻留时间和发射功率的最优解,得到利于射频隐身的编队雷达功率与时间资源管理结果。本发明面对敌方先进无源探测系统,平衡了编队雷达系统在功率域与时间域的被截获概率,提升了编队雷达系统面对敌方先进无源探测系统时的射频隐身性能。

Description

一种面对先进无源探测系统的编队雷达资源管理方法
技术领域
本发明涉及编队雷达资源管理方法,具体涉及一种面对先进无源探测系统的 编队雷达资源管理方法。
背景技术
随着现代战场中电子对抗的日益激烈,雷达的生存环境受到了严重的威胁。 通过对雷达资源进行有效管理能显著降低雷达被探测、发现、识别和攻击的概率, 是提高雷达及其运载平台的战场生存能力和作战效能的重要保证。相比目标外形 特征减缩和红外特征减缩,雷达资源管理并非无限制的减缩雷达在功率域与时间 域的辐射特征,而是在满足设备功能、性能要求的基础上对辐射功率与时间进行 有效控制,提高其对抗敌方无源探测系统的性能。
基于雷达资源管理理论,目前可采用的对抗敌方无源探测系统的优化策略主 要有两大类:最小辐射能量策略和最大信号不确定性策略。最小辐射能量策略要 求在任何时间都应以系统所需的最小能量向外辐射,该策略通过主动辐射源的辐 射功率管理、辐射时间优化和低旁瓣天线设计,降低系统的辐射能量和旁瓣功率。 目前关于单个机载雷达的辐射能量控制策略已相对较成熟。
随着计算机技术、通信技术和微波集成电路的快速发展,以及现代战争复杂 性的日益提高,越来越多的传感器被纳入一体化网络参与协同作战。同时,面对 日益复杂的战场电磁环境,综合利用多传感器的信息在空间域进行多传感器信息 融合不仅可以提高系统的可靠性和生存能力,而且可以尽可能全面、准确地获取 信息。编队雷达系统是未来网络化斗争发展的必然趋势,它是利用多部雷达同步 地发射正交的波形,同时使用多部雷达接收回波信号,并集中处理的一种新型雷 达体制。它主要利用目标雷达截面积的空间分集增益来提高检测性能。在提高雷 达探测威力、抑制干扰和对抗无源探测系统等方面具有巨大的潜力。
编队雷达作为研究的新领域,许多文献主要着眼于编队雷达的探测性能,而 以对抗敌方先进无源探测系统作为优化目标的研究相对较少。为提高编队雷达面 对敌方无源探测系统时的被截获性能,有文献通过最小辐射能量控制策略,优化 每一时刻的雷达发射功率,达到降低雷达截获因子的目的。但是编队雷达作为一 种组网雷达系统,可控参数很多,单一优化一个参量,对雷达性能的改变并不明 显,也会造成其他资源的不必要浪费。
发明内容
发明目的:为解决现有技术存在的问题,本发明的目的是提供一种在满足编 队雷达跟踪性能要求的前提下,通过动态地优化各部雷达的发射功率和驻留时间, 达到最优化编队雷达系统被截获性能的目的,以提升系统对抗无源探测系统的能 力的编队雷达资源管理方法。
技术方案:一种面对先进无源探测系统的编队雷达资源管理方法,包括如下 步骤:
步骤一:根据战场环境中典型敌方先进无源探测系统和编队雷达系统性能, 确定典型敌方先进无源探测系统的截获概率建模方法;
步骤二:确定我方每部雷达发射功率和驻留时间的上、下界;
步骤三:以最优化我方编队雷达系统在功率域与时间域的被截获性能为目标, 在满足目标跟踪过程中检测性能的条件下,面对先进无源探测系统,构建最优化 被截获性能的编队雷达驻留时间与功率资源联合优化模型;
步骤四:经数值计算,得到在满足目标跟踪过程中检测性能条件下,使得编 队雷达系统在功率域与时间域的被截获性能最优的各雷达的驻留时间和发射功 率作为最优解,进而可得到利于射频隐身的编队雷达功率与时间资源管理结果。
所述步骤二具体包括:
步骤2.1:确定编队雷达系统中每部雷达的最小和大发射功率Pt min和Pt max,最 短和最长驻留时间并将作为优化模型的资源优 化区间;i代表第i部雷达,相应的Pti和τei就代表第i部雷达某一时刻的发射功率 与驻留时间;
步骤2.2:根据编队雷达系统目标跟踪过程中检测性能的要求,确定目标回 波的二元检测门限γth
步骤2.3:根据给定的目标反射系数的方差Rg、传播损耗因子pij、雷达接收 机的噪声方差Rθ以及发射信号脉冲重复频率fr参量,计算每一时刻编队雷达系 统目标跟踪过程中检测性能的二元检测门限γth,并计算对应于γth的回波检测信 噪比门限作为优化模型的约束条件;
步骤2.4:根据先进无源探测系统在功率域的虚警概率Pfa、检测信噪比SNRi和时间域的短时平均搜索时间以及雷达的性能参数,计算编队雷达系统截获 性能pai=(wdPd-wfPf)2,其中,Pd和Pf分别代表敌方先进无源探测系统的功率 域与时间域截获概率,wd和wf分别对应于Pd和Pf的经验加权系数。并将 作为优化模型的目标函数;
步骤2.5:根据步骤1确定的第i部雷达的发射功率和驻留时间区间、步骤 3确定的约束条件以及步骤4确定的目标函数,构建面对敌方先进无源探测系统 的编队雷达资源联合优化模型;
步骤2.6:对步骤2.5建立的优化模型进行求解,得到当前时刻使得编队雷 达系统被截获性能pai最优的发射功率Pti *和驻留时间解,并循环求解出目标 跟踪过程中满足检测性能要求的所有时刻的发射功率和驻留时间的解集。
进一步的,步骤2.3中编队雷达系统目标跟踪过程中回波二元检测门限γth的 数学表达式为:
其中,njk(t0)~N(0,Rθ/fr),Rθ为雷达接收机的噪声方差,gjk为目标反射系数的 方差,pjk为传播损耗因子,fr为雷达发射信号的脉冲重复频率,Nt和Nt分别为 雷达系统的发射机和接收机个数,xk是雷达发射的信号,τkj是雷达信号从第k部 雷达信号经目标反射到第i部雷达的时间延迟。
步骤2.4中编队雷达系统面对先进无源探测系统的功率域被截获概率的数 学表达式为:
Pd=maxPdi
编队雷达系统面对先进无源探测系统的时间域被截获概率的数学表达式为:
Pf=maxPfi
其中,Pfa和SNRi分别是先进无源探测系统的虚警概率和检测信噪比,是 短时平均搜索时间,I0是零阶修正贝赛尔函数。
步骤2.5中构建的面对敌方先进无源探测系统的编队雷达被截获性能的优 化模型为:
为优化目标,为约束条件,采用最小二乘算法进行 计算,求得使目标函数pai最优的一组解Pti *即为当前时刻发射功率Pti和 驻留时间τei的一组最优解。
有益效果
与现有技术相比,本发明具有以下效果:1、通过将实际战场中雷达的发射 功率和驻留时间建模为上、下界已知的不确定集合,以先进无源探测系统的功率 域与时间域的一般性能作为先验知识,以最优化系统的截获性能为目标,在满足 一定目标跟踪性能的条件下建立对抗先进无源探测系统的驻留时间与功率资源 联合优化模型;因此,本发明能够在满足编队雷达跟踪过程中的检测性能要求的 前提下,通过动态地优化各部雷达的发射功率和驻留时间,以实现最优化编队雷 达对抗敌方先进无源探测系统的性能;不仅保证了系统在目标跟踪过程中的检测 性能,而且使系统具备对抗敌方先进无源探测系统的最优性能。2、本发明不仅 仅考虑了系统在目标跟踪过程中的功率域与时间域资源管理问题,同时也是实现 了编队雷达系统资源的有效利用。
附图说明
图1为编队雷达驻留时间与功率资源联合优化方法流程图;
图2为目标跟踪场景;
图3为编队雷达系统与目标之间的距离关系;
图4(a)-4(b)为跟踪过程中编队雷达系统的跟踪误差和最大化采样间隔;
图5(a)-5(b)为跟踪过程中编队雷达系统得最优功率和驻留时间分配结 果。
具体实施方法
下面结合具体实施方式和附图对本发明的技术方案作进一步说明。
本实施例的一种面对先进无源探测系统的编队雷达资源管理方法,首先根据战场环境中典型敌方先进无源探测系统和编队雷达系统性能,确定典型敌方先进无源 探测系统的截获概率建模方法;确定我方每部雷达发射功率和驻留时间的上、下 界;然后以最优化我方编队雷达系统在功率域与时间域的被截获性能为目标,在 满足目标跟踪过程中检测性能的条件下,面对先进无源探测系统,构建最优化被 截获性能的编队雷达驻留时间与功率资源联合优化模型;经数值计算,得到在满 足目标跟踪过程中检测性能条件下使得编队雷达系统截获性能pai最优的各雷达 的驻留时间和发射功率作为最优解,进而可得到编队雷达系统当前时刻 对抗敌方先进无源探测系统的最优性能。
如图1所示,具体包括以下步骤:
1、确定发射功率和驻留时间的优化区间
雷达的发射功率和驻留时间的上、下界不仅和雷达系统的性能参数有关,还 和当前战场环境中目标距雷达的距离有关。首先,根据编队雷达系统的性能参数, 再根据预测的目标距雷达的距离,确定编队雷达系统中每部雷达的最小和大发射 功率Pt min和Pt max,最短和最长驻留时间并将作为 优化模型的资源优化区间;;
2、建立约束条件
根据给定的目标反射系数的方差Rg、传播损耗因子pij、雷达接收机的噪声 方差Rθ以及发射信号脉冲重复频率fr参量,计算每一时刻编队雷达系统目标跟 踪过程中检测性能的二元检测门限γth,并计算对应于γth的回波检测信噪比门限 作为优化模型的约束条件;
编队雷达系统目标跟踪过程中回波二元检测门限γth的数学表达式为:
其中,njk(t0)~N(0,Rθ/fr),Rθ为雷达接收机的噪声方差,gjk为目标反射系数的 方差,pjk为传播损耗因子,fr为雷达发射信号的脉冲重复频率,Nt和Nt分别为 雷达系统的发射机和接收机个数,xk是雷达发射的信号,τkj是雷达信号从第k部 雷达信号经目标反射到第i部雷达的时间延迟。
3、建立优化模型的目标函数
根据先进无源探测系统在功率域的虚警概率Pfa、检测信噪比SNRi和时间域 的短时平均搜索时间以及雷达的性能参数,计算编队雷达系统截获性能 pai=(wdPd-wfPf)2,其中,Pd和Pf分别代表敌方先进无源探测系统的功率域与 时间域截获概率,wd和wf分别对应于Pd和Pf的经验加权系数。并将作 为优化模型的目标函数。编队雷达系统面对先进无源探测系统的功率域被截获概 率的数学表达式为:
Pd=max Pdi (2)
编队雷达系统面对先进无源探测系统的时间域被截获概率的数学表达式为:
Pf=maxPfi (3)
其中,Pfa和SNRi分别是先进无源探测系统的虚警概率和检测信噪比,是 短时平均搜索时间,I0是零阶修正贝赛尔函数。
4、建立驻留时间与功率资源联合优化模型
根据步骤1确定的第i部雷达的发射功率和驻留时间区间、步骤2确定的约 束条件以及步骤3确定的目标函数,构建编队雷达系统的驻留时间与功率资源联 合优化模型:
5、获得发射功率和驻留时间的最优解
为优化目标,为约束条件,采用最小二乘算法进行 计算,求得使目标函数pai最优的一组解Pti *即为当前时刻发射功率Pti和 驻留时间τei的一组最优解。循环求解出目标跟踪过程中满足检测性能要求的所 有时刻的发射功率和驻留时间的解集。
6、仿真结果
本实施例针对一个二维平面中运动的目标场景进行了仿真;仿真中,假设 Nt=4;初始时刻雷达位置分布如表1所示。
编队雷达系统中每部雷达最大发射功率为最小发射功率为每部雷达最大驻留时间为最小驻留时间为
表1
雷达 位置
Radar1 [0,0]km
Radar2 [40,0]km
Radar3 [0,30]km
Radar4 [40,30]km
采用卡尔曼滤波算法对单目标进行跟踪。目标跟踪场景如图2所示,图2 中显示了四部雷达、单个目标的位置与航迹。图3中显示了四部雷达跟踪单个目 标的距离关系。通过卡尔曼滤波算法,图4(a)、图4(b)显示了计算出的跟踪 过程中编队雷达系统的跟踪误差和最大化采样间隔。经数值计算,图5(a)、图 5(b)显示了对抗敌先进无源探测系统时的各个时刻每部雷达的发射功率和驻留 时间分配情况。
由上述仿真结果可知,本发明在保证目标跟踪性能的前提下,通过动态地调 整各部雷达的发射功率和驻留时间,能有效地优化编队雷达系统的截获性能。并 且,在整个目标跟踪过程中,各部雷达并非时刻采用最大发射功率和驻留时间进 行工作,而是通过合理的分配雷达的资源,实现了雷达资源的有效利用。
上面结合附图所描述的本发明优选具体实施例仅用于说明本发明的实施方 式,而不是作为对前述发明目的和所附权利要求内容和范围的限制,凡是依据本 发明的技术实质对以上实施例所做的任何简单修改、等同变化与修饰,均仍属本 发明技术和权利保护范畴。

Claims (5)

1.一种面对先进无源探测系统的编队雷达资源管理方法,其特征在于,包括如下步骤:
步骤一:根据战场环境中典型敌方先进无源探测系统和编队雷达系统性能,确定典型敌方先进无源探测系统的截获概率建模方法;
步骤二:确定我方每部雷达发射功率和驻留时间的上、下界;
步骤三:以最优化我方编队雷达系统在功率域与时间域的被截获性能为目标,在满足目标跟踪过程中检测性能的条件下,面对先进无源探测系统,构建最优化被截获性能的编队雷达驻留时间与功率资源联合优化模型;
步骤四:经数值计算,得到在满足目标跟踪过程中检测性能条件下,使得编队雷达系统在功率域与时间域的被截获性能最优的各雷达的驻留时间和发射功率作为最优解,进而得到利于射频隐身的编队雷达功率与时间资源管理结果。
2.根据权利要求1所述的编队雷达资源管理方法,其特征在于:所述步骤二具体包括:
步骤2.1:确定编队雷达系统中每部雷达的最小和大发射功率Pt min和Pt max,最短和最长驻留时间并将作为优化模型的资源优化区间,i代表第i部雷达,相应的Pti和τei就代表第i部雷达某一时刻的发射功率与驻留时间;
步骤2.2:根据编队雷达系统目标跟踪过程中检测性能的要求,确定目标回波的二元检测门限γth
步骤2.3:根据给定的目标反射系数的方差Rg、传播损耗因子pij、雷达接收机的噪声方差Rθ以及发射信号脉冲重复频率fr参量,计算每一时刻编队雷达系统目标跟踪过程中检测性能的二元检测门限γth,并计算对应于γth的回波检测信噪比门限作为优化模型的约束条件;
步骤2.4:根据先进无源探测系统在功率域的虚警概率Pfa、检测信噪比SNRi和时间域的短时平均搜索时间以及雷达的性能参数,计算编队雷达系统截获性能pai=(wdPd-wfPf)2,并将作为优化模型的目标函数;其中,Pd和Pf分别代表敌方先进无源探测系统的功率域与时间域截获概率,wd和wf分别对应于Pd和Pf的经验加权系数。
步骤2.5,根据步骤2.1确定的第i部雷达的发射功率和驻留时间区间、步骤2.3确定的约束条件以及步骤2.4确定的目标函数,构建面对敌方先进无源探测系统的编队雷达资源联合优化模型;
步骤2.6,对步骤2.5建立的优化模型进行求解,得到当前时刻使得编队雷达系统被截获性能pai最小的最优发射功率Pti *和驻留时间解,并循环求解出目标跟踪过程中满足检测性能要求的所有时刻的发射功率和驻留时间的解集。
3.根据权利要求2所述的编队雷达资源管理方法,其特征在于:所述步骤2.3中编队雷达系统目标跟踪过程中回波二元检测门限γth的数学表达式为:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>H</mi> <mn>0</mn> </msub> <mo>:</mo> <mi>&amp;gamma;</mi> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>t</mi> </msub> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>t</mi> </msub> </munderover> <mo>|</mo> <msub> <mi>n</mi> <mrow> <mi>j</mi> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>H</mi> <mn>1</mn> </msub> <mo>:</mo> <mi>&amp;gamma;</mi> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>t</mi> </msub> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>t</mi> </msub> </munderover> <mo>|</mo> <msub> <mi>y</mi> <mrow> <mi>j</mi> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced>
<mrow> <msub> <mi>y</mi> <mrow> <mi>j</mi> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>m</mi> <mi>a</mi> <mi>x</mi> <mo>{</mo> <mo>&amp;Integral;</mo> <msub> <mi>y</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msubsup> <mi>x</mi> <mi>k</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <mi>&amp;tau;</mi> <mo>-</mo> <mi>t</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>&amp;tau;</mi> <mo>=</mo> <msub> <mi>p</mi> <mrow> <mi>j</mi> <mi>k</mi> </mrow> </msub> <msub> <mi>g</mi> <mrow> <mi>j</mi> <mi>k</mi> </mrow> </msub> <msqrt> <msub> <mi>P</mi> <mrow> <mi>t</mi> <mi>k</mi> </mrow> </msub> </msqrt> <mo>&amp;Integral;</mo> <msub> <mi>x</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>&amp;tau;</mi> <mo>-</mo> <msub> <mi>&amp;tau;</mi> <mrow> <mi>k</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msubsup> <mi>x</mi> <mi>k</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <mi>&amp;tau;</mi> <mo>-</mo> <mi>t</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>&amp;tau;</mi> <mo>+</mo> <msub> <mi>n</mi> <mrow> <mi>j</mi> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>}</mo> </mrow>
其中,njk(t0)~N(0,Rθ/fr),Rθ为雷达接收机的噪声方差,gjk为目标反射系数的方差,pjk为传播损耗因子,fr为雷达发射信号的脉冲重复频率,Nt和Nt分别为雷达系统的发射机和接收机个数,xk是雷达发射的信号,τkj是雷达信号从第k部雷达信号经目标反射到第i部雷达的时间延迟。
4.根据权利要求2所述的面对敌方先进无源探测系统的编队雷达资源管理方法,其特征在于:所述步骤2.4中编队雷达系统面对先进无源探测系统的功率域被截获概率的数学表达式为:
Pd=maxPdi
<mrow> <msub> <mi>P</mi> <mrow> <mi>d</mi> <mi>i</mi> </mrow> </msub> <mo>=</mo> <mn>2</mn> <msubsup> <mo>&amp;Integral;</mo> <msqrt> <mrow> <mo>-</mo> <mi>ln</mi> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mrow> <mi>f</mi> <mi>a</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </msqrt> <mi>&amp;infin;</mi> </msubsup> <mi>y</mi> <mi> </mi> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mo>(</mo> <mrow> <msup> <mi>y</mi> <mn>2</mn> </msup> <mo>+</mo> <msub> <mi>SNR</mi> <mi>i</mi> </msub> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <msub> <mi>I</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mn>2</mn> <mi>y</mi> <msqrt> <mrow> <msub> <mi>SNR</mi> <mi>i</mi> </msub> </mrow> </msqrt> <mo>)</mo> </mrow> <mi>d</mi> <mi>y</mi> </mrow>
编队雷达系统面对先进无源探测系统的时间域被截获概率的数学表达式为:
Pf=maxPfi
<mrow> <msub> <mi>P</mi> <mi>f</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mfrac> <mn>1</mn> <msub> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mi>c</mi> </msub> </mfrac> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>x</mi> <mo>/</mo> <msub> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mi>c</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>x</mi> <mo>&gt;</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>x</mi> <mo>&lt;</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>
其中,Pfa和SNRi分别是先进无源探测系统的虚警概率和检测信噪比,是短时平均搜索时间,I0是零阶修正贝赛尔函数。
5.根据权利要求2所述的面对敌方先进无源探测系统的编队雷达资源管理方法,其特征在于:所述步骤2.5中构建的面对敌方先进无源探测系统的编队雷达被截获性能的优化模型为:
<mrow> <mfenced open = "" close = "}"> <mtable> <mtr> <mtd> <mrow> <munder> <mi>min</mi> <msub> <mi>N</mi> <mi>t</mi> </msub> </munder> <msub> <mi>p</mi> <mrow> <mi>a</mi> <mi>i</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>s</mi> <mo>.</mo> <mi>t</mi> <mo>.</mo> <mo>:</mo> <mn>0</mn> <mo>&lt;</mo> <msub> <mi>P</mi> <mi>f</mi> </msub> <mo>&amp;le;</mo> <msub> <mi>P</mi> <mrow> <mi>f</mi> <mi>t</mi> <mi>h</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>0</mn> <mo>&lt;</mo> <msub> <mi>P</mi> <mi>d</mi> </msub> <mo>&amp;le;</mo> <msub> <mi>P</mi> <mrow> <mi>d</mi> <mi>t</mi> <mi>h</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>P</mi> <mi>t</mi> <mi>min</mi> </msubsup> <mo>&lt;</mo> <msub> <mi>P</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mi>t</mi> <mi>max</mi> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>&amp;tau;</mi> <mrow> <mi>e</mi> <mi>i</mi> </mrow> <mi>min</mi> </msubsup> <mo>&amp;le;</mo> <msub> <mi>&amp;tau;</mi> <mrow> <mi>e</mi> <mi>i</mi> </mrow> </msub> <mo>&amp;le;</mo> <msubsup> <mi>&amp;tau;</mi> <mrow> <mi>e</mi> <mi>i</mi> </mrow> <mi>max</mi> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>SNR</mi> <mrow> <mi>n</mi> <mi>e</mi> <mi>t</mi> </mrow> <mrow> <mi>t</mi> <mi>h</mi> </mrow> </msubsup> <mo>&amp;le;</mo> <msub> <mi>SNR</mi> <mrow> <mi>n</mi> <mi>e</mi> <mi>t</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>.</mo> </mrow>
为优化目标,为约束条件,采用最小二乘算法进行计算,求得使目标函数pai最优的一组解Pti *即为当前时刻发射功率Pti和驻留时间τei的一组最优解。
CN201710996643.XA 2017-10-19 2017-10-19 一种面对先进无源探测系统的编队雷达资源管理方法 Active CN107944597B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710996643.XA CN107944597B (zh) 2017-10-19 2017-10-19 一种面对先进无源探测系统的编队雷达资源管理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710996643.XA CN107944597B (zh) 2017-10-19 2017-10-19 一种面对先进无源探测系统的编队雷达资源管理方法

Publications (2)

Publication Number Publication Date
CN107944597A true CN107944597A (zh) 2018-04-20
CN107944597B CN107944597B (zh) 2021-07-27

Family

ID=61935596

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710996643.XA Active CN107944597B (zh) 2017-10-19 2017-10-19 一种面对先进无源探测系统的编队雷达资源管理方法

Country Status (1)

Country Link
CN (1) CN107944597B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109308344A (zh) * 2018-08-30 2019-02-05 西北工业大学 基于空情航迹共享的机载雷达信号被截获概率建模方法
CN109829251A (zh) * 2019-03-12 2019-05-31 中国人民解放军空军预警学院 基于算子的无源雷达目标截获概率柔性仿真建模方法
CN110412515A (zh) * 2019-08-19 2019-11-05 南京航空航天大学 基于射频隐身的组网雷达多目标跟踪发射功率分配方法
CN110794371A (zh) * 2019-10-30 2020-02-14 南京航空航天大学 基于改进递推法的雷达采样间隔控制方法
CN111323772A (zh) * 2020-02-19 2020-06-23 南京航空航天大学 基于射频隐身的有/无人机混合编队协同任务分配方法
CN112766775A (zh) * 2021-01-27 2021-05-07 中国人民解放军海军工程大学 微波武器在舰艇防空反导体系贡献率评估方法
CN113917399A (zh) * 2021-09-30 2022-01-11 中国船舶重工集团公司第七二四研究所 一种机扫被动探测雷达自动化扫频策略设定方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104346537A (zh) * 2014-11-17 2015-02-11 西安电子科技大学 机载雷达射频隐身性能评测方法
CN106291481A (zh) * 2016-07-27 2017-01-04 南京航空航天大学 基于射频隐身的分布式mimo雷达资源联合优化方法
CN107064882A (zh) * 2017-01-22 2017-08-18 南京航空航天大学 无源协同下基于射频隐身的雷达组网资源控制方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104346537A (zh) * 2014-11-17 2015-02-11 西安电子科技大学 机载雷达射频隐身性能评测方法
CN106291481A (zh) * 2016-07-27 2017-01-04 南京航空航天大学 基于射频隐身的分布式mimo雷达资源联合优化方法
CN107064882A (zh) * 2017-01-22 2017-08-18 南京航空航天大学 无源协同下基于射频隐身的雷达组网资源控制方法

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109308344A (zh) * 2018-08-30 2019-02-05 西北工业大学 基于空情航迹共享的机载雷达信号被截获概率建模方法
CN109829251A (zh) * 2019-03-12 2019-05-31 中国人民解放军空军预警学院 基于算子的无源雷达目标截获概率柔性仿真建模方法
CN110412515A (zh) * 2019-08-19 2019-11-05 南京航空航天大学 基于射频隐身的组网雷达多目标跟踪发射功率分配方法
CN110412515B (zh) * 2019-08-19 2023-02-03 南京航空航天大学 基于射频隐身的组网雷达多目标跟踪发射功率分配方法
CN110794371A (zh) * 2019-10-30 2020-02-14 南京航空航天大学 基于改进递推法的雷达采样间隔控制方法
CN111323772A (zh) * 2020-02-19 2020-06-23 南京航空航天大学 基于射频隐身的有/无人机混合编队协同任务分配方法
CN111323772B (zh) * 2020-02-19 2022-04-19 南京航空航天大学 基于射频隐身的有/无人机混合编队协同任务分配方法
CN112766775A (zh) * 2021-01-27 2021-05-07 中国人民解放军海军工程大学 微波武器在舰艇防空反导体系贡献率评估方法
CN112766775B (zh) * 2021-01-27 2022-05-03 中国人民解放军海军工程大学 微波武器在舰艇防空反导体系贡献率评估方法
CN113917399A (zh) * 2021-09-30 2022-01-11 中国船舶重工集团公司第七二四研究所 一种机扫被动探测雷达自动化扫频策略设定方法

Also Published As

Publication number Publication date
CN107944597B (zh) 2021-07-27

Similar Documents

Publication Publication Date Title
CN107944597A (zh) 一种面对先进无源探测系统的编队雷达资源管理方法
CN106291481B (zh) 基于射频隐身的分布式mimo雷达资源联合优化方法
CN106896351B (zh) 一种基于非合作博弈的雷达组网功率控制方法
CN108008361B (zh) 基于射频隐身的分布式mimo雷达干扰波形设计方法
CN108732543B (zh) 一种基于射频隐身的机载组网雷达辐射参数联合优化方法
CN107064882B (zh) 无源协同下基于射频隐身的雷达组网资源控制方法
CN111812599B (zh) 博弈条件下基于低截获性能的组网雷达最优波形设计方法
CN106501778B (zh) 一种基于射频隐身性能优化的dmrs稳健波形设计方法
CN111090078B (zh) 一种基于射频隐身的组网雷达驻留时间优化控制方法
CN105425225B (zh) 一种被动雷达低空目标检测方法
CN107907863B (zh) 基于雷达-通信频谱共享的组网雷达波形设计方法
CN108614261B (zh) 一种组网雷达系统多目标跟踪下的辐射参数控制方法
CN110412515A (zh) 基于射频隐身的组网雷达多目标跟踪发射功率分配方法
CN111812628B (zh) 基于非合作博弈的雷达通信一体化组网系统功率控制方法
CN112118621B (zh) 一种机载雷达通信一体化系统辐射功率优化设计方法
CN110007278A (zh) 一种基于斯塔科尔伯格博弈的多基雷达射频隐身优化方法
CN105891799B (zh) 适用于机械扫描雷达的有源干扰侦察方法
CN104020459A (zh) 一种提高mimo-stap检测性能的波形优化方法
CN103852749A (zh) 提高mimo-stap检测性能的稳健波形优化方法
Shi et al. Low probability of intercept optimization for radar network based on mutual information
CN110927692A (zh) 海杂波场景下搜索雷达射频隐身工作方式的求解方法
CN109212494B (zh) 一种针对组网雷达系统的射频隐身干扰波形设计方法
CN107589409B (zh) 一种分置天线mimo雷达分布式低通信量检测融合方法
CN113238198A (zh) 基于时分复用的雷达通信一体化系统辐射时间控制方法
CN103678906B (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