CN108008361B - 基于射频隐身的分布式mimo雷达干扰波形设计方法 - Google Patents
基于射频隐身的分布式mimo雷达干扰波形设计方法 Download PDFInfo
- Publication number
- CN108008361B CN108008361B CN201711082783.2A CN201711082783A CN108008361B CN 108008361 B CN108008361 B CN 108008361B CN 201711082783 A CN201711082783 A CN 201711082783A CN 108008361 B CN108008361 B CN 108008361B
- Authority
- CN
- China
- Prior art keywords
- interference
- radar
- distributed mimo
- mimo radar
- radio frequency
- 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 39
- 238000013461 design Methods 0.000 title claims abstract description 36
- 238000001228 spectrum Methods 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 32
- 238000013178 mathematical model Methods 0.000 claims description 7
- 238000005457 optimization Methods 0.000 claims description 7
- 230000005855 radiation Effects 0.000 claims description 7
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 2
- 230000005540 biological transmission Effects 0.000 abstract description 11
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 239000000243 solution Substances 0.000 description 11
- 238000001514 detection method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000035899 viability Effects 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000003471 anti-radiation Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
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
- 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/38—Jamming means, e.g. producing false echoes
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于射频隐身的分布式MIMO雷达干扰波形设计方法,首先在机载电子干扰***与分布式MIMO雷达进行电子对抗的情况下,根据先验知识,获取目标相对于分布式MIMO雷达的频率响应、雷达发射信号及接收机噪声功率谱;然后,以最小化机载电子干扰***总发射功率为目标,建立基于射频隐身的分布式MIMO雷达干扰波形设计模型,通过拉格朗日乘数法对模型进行求解;经二分法迭代计算,选取机载电子干扰***总发射功率最小的干扰波形作为最优解,即可得到符合约束条件的最小干扰总功率。本发明所达到的有益效果:本方法不仅保证了对分布式MIMO雷达***的干扰性能,而且有效提升了机载电子干扰***的射频隐身性能。
Description
技术领域
本发明涉及一种基于射频隐身的分布式MIMO雷达干扰波形设计方法,属于雷达干扰波形设计技术领域。
背景技术
分布式MIMO雷达***是一个由多部发射天线和多部接收天线组成并相互协作、紧密关联的有机整体,实现时域、频域和空域的协同工作,完成对目标的侦察探测、情报收集、识别跟踪、干扰压制、火力引导等功能。同一视场中的目标可被多部发射天线和接收天线探测。为了阻止分布式MIMO雷达对目标的探测跟踪,仅干扰***中的部分雷达接收机是不够的,必须对所有雷达接收机进行有效干扰才能掩护突防目标。而机载电子干扰***的干扰资源往往是有限的,为达到最佳干扰效果,需要对干扰发射波形进行合理设计。
常用的雷达干扰波形设计主要有两种准则:一是最小化信干噪比(Signal toInterference plus Noise Ratio,SINR),以降低雷达***的目标检测性能;二是最小化雷达接收回波与目标冲激响应之间的互信息(Mutual Information,MI),从而降低雷达***的目标参数估计性能。
然而,日益先进的无源探测***能够通过截获并处理电子干扰机发射的射频信号对其进行精确的探测、定位、分选和识别,对己方机载电子干扰***的战场生存能力构成严重威胁。射频隐身技术通过控制有源辐射***的辐射能量、波形参数等方法,可显著降低有源辐射***被敌方无源探测***截获、发现、分选、识别,以及被反辐射导弹攻击的概率,从而提高其自身及其搭载平台的战场生存力和突防能力。
目前,对射频隐身技术的研究主要集中在雷达与通信***领域,而针对分布式MIMO雷达的有源压制性干扰则未见公开报道。现有的关于有源压制性干扰的研究主要以干扰效能最大化为目标,但增大干扰***发射功率的同时也增加了敌方无源探测***对干扰信号的截获概率,从而大大降低了干扰***的射频隐身性能。因此,需要研究基于射频隐身的分布式MIMO雷达干扰波形设计问题。
发明内容
为解决现有技术的不足,本发明的目的在于提供一种基于射频隐身的分布式MIMO雷达干扰波形设计方法,能够通过干扰波形优化设计,最小化机载电子干扰***的总发射功率,提升其射频隐身性能。
为了实现上述目标,本发明采用如下的技术方案:
一种基于射频隐身的分布式MIMO雷达干扰波形设计方法,其特征是,包括如下步骤:
1)根据先验知识,确立目标频率响应、MIMO雷达发射信号和雷达接收机噪声功率谱,并确定干扰***的辐射参数与干扰性能MI门限参数;
2)以最小化机载电子干扰***总发射功率为目标,在满足一定干扰性能的条件下,建立基于射频隐身的分布式MIMO雷达干扰波形设计模型;
3)通过拉格朗日乘数法对模型进行求解,在满足干扰性能门限δMI的最小干扰总功率的情况下,确定拉格朗日乘子的最优解;
4)将步骤3)得到的最优解代入卡罗需-库恩-塔克条件的必要条件获取分布式MIMO雷达的最佳干扰发射波形,从而形成分布式MIMO雷达干扰波形设计方案。
前述的一种基于射频隐身的分布式MIMO雷达干扰波形设计方法,其特征是,所述步骤1)具体内容为:
获取分布式MIMO雷达发射信号矩阵S,目标相对于雷达***的频率响应h,雷达接收机处噪声n作为先验知识;其中,S=[S1,S2,…,SN],N为分布式MIMO雷达发射天线个数,***中第i部天线的发射信号Si为M×L维Toeplitz矩阵,M为雷达回波信号长度,L为目标频率响应长度,且有M=NL;目标相对于雷达的频率响应h为一个零均值复高斯随机矢量,且满足其中,Ξh为目标频率响应协方差矩阵;雷达接收机处噪声n为一个零均值复高斯随机矢量,且满足其中,Ξn为噪声协方差矩阵。
前述的一种基于射频隐身的分布式MIMO雷达干扰波形设计方法,其特征是,所述步骤2)具体内容为:
根据事先给定的干扰性能MI门限δMI,建立基于射频隐身的分布式MIMO雷达最优干扰波形设计数学模型:式中,上标(·)H表示矩阵的共轭转置,干扰信号J为一个零均值复高斯随机矢量,且满足其中,ΞJ为干扰协方差矩阵。
前述的一种基于射频隐身的分布式MIMO雷达干扰波形设计方法,其特征是,所述步骤2)处理后采用特征值分解,协方差矩阵Ξh、ΞJ、Ξn可分别分解如下:式中,Uh、UJ和Un分别为酉矩阵,对角矩阵Λh=diag[λ1,…,λNL],其中,λi、和分别为相应对角矩阵的特征值;经矩阵运算,将数学模型转化为:式中,为雷达发射信号矩阵S的特征值。
前述的一种基于射频隐身的分布式MIMO雷达干扰波形设计方法,其特征是,所述步骤3)具体包括如下步骤:
33)采用二分法迭代计算确定最优化的卡罗需-库恩-塔克条件下的拉格朗日乘子。
前述的一种基于射频隐身的分布式MIMO雷达干扰波形设计方法,其特征是,所述步骤32)包括如下内容:
前述的一种基于射频隐身的分布式MIMO雷达干扰波形设计方法,其特征是,所述步骤33)的具体内容为:
通过求解式(2),机载电子干扰***的最优干扰波形可表示为:ξ*是一个常数,其大小取决于给定的MI门限:经二分法迭代计算,将满足式(4)的ξ*值代入式(3)中,求得使机载电子干扰***总发射功率最小的一组干扰波形作为最优解,并最终确定***的干扰总功率。
本发明所达到的有益效果:
1)本方法所完成的主要任务是针对实际战场中机载电子干扰机与分布式MIMO雷达进行电子对抗的情况,在根据先验知识获取目标相对于分布式MIMO雷达的频率响应、雷达发射信号及噪声功率谱的基础上,以最小化***干扰总功率为目标,在满足一定干扰性能的条件下,对机载电子干扰***的发射波形进行自适应优化设计;
2)本方法在获取目标相对于分布式MIMO雷达的频率响应、雷达发射信号及噪声功率谱的基础上,以最小化机载电子干扰***的总发射功率为目标,在满足一定干扰性能的条件下,建立基于射频隐身的分布式MIMO雷达干扰波形设计模型,既保证了对分布式MIMO雷达***的干扰性能,还使机载电子干扰***具有最优的射频隐身性能;
3)与现有技术相比,本方法,不仅保证了对分布式MIMO雷达***的干扰性能,而且有效提升了机载电子干扰***的射频隐身性能。
附图说明
图1是分布式MIMO雷达干扰模型;
图2是分布式MIMO雷达干扰波形设计流程图;
图3是分布式MIMO雷达发射信号频谱、目标频率响应功率谱及噪声功率谱;
图4是分布式MIMO雷达干扰波形设计结果;
图5是不同干扰波形设计方法下干扰总功率对比图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
本发明首先在机载电子干扰***与分布式MIMO雷达进行电子对抗的情况下,根据先验知识,获取目标相对于分布式MIMO雷达的频率响应、雷达发射信号及接收机噪声功率谱;然后,以最小化机载电子干扰***总发射功率为目标,在满足一定干扰性能的条件下,建立基于射频隐身的分布式MIMO雷达干扰波形设计模型,并通过拉格朗日乘数法对模型进行求解。经二分法迭代计算,选取在满足一定干扰性能的条件下使得机载电子干扰***总发射功率最小的干扰波形作为最优解,将最优干扰发射波形代入卡罗需-库恩-塔克条件(KKT)的必要条件中,即可得到符合约束条件的最小干扰总功率。
具体的方法内容如下:
1)确定目标频率响应与MIMO雷达发射信号:
利用雷达发射信号、目标频率响应和接收机噪声等先验知识,假定分布式MIMO雷达干扰模型如图1所示。
2)以最小化机载电子干扰***总发射功率为目标,在满足一定干扰性能的条件下,建立基于射频隐身的分布式MIMO雷达干扰波形设计模型;
根据机载电子干扰***对分布式MIMO雷达干扰性能的要求,建立基于射频隐身的分布式MIMO雷达最优干扰波形设计数学模型,如下所示:
采用特征值分解,协方差矩阵Ξh、ΞJ、Ξn可分别分解如下:
其中,所有带“*”上标的变量分别表示各参数的最优解。
ξ*是一个常数,其大小取决于给定的MI门限:
4)将步骤3)得到的最优解代入卡罗需-库恩-塔克条件(KKT)的必要条件获取分布式MIMO雷达的最佳干扰发射波形,从而形成分布式MIMO雷达干扰波形设计方案。
这里给出一个实施例,本例子中,假设第2步中的参数如表1所示。
表1仿真参数设置
分布式MIMO雷达发射信号频谱、目标频率响应功率谱及噪声功率谱如图3所示。针对分布式MIMO雷达的干扰波形设计结果如图4所示。
基于射频隐身的分布式MIMO雷达干扰波形设计方法是根据MIMO雷达发射信号、目标相对于雷达的频率响应及雷达接收机处噪声功率谱计算所得的最优干扰波形。由图4可知,机载电子干扰***发射波形主要由MIMO雷达发射信号、目标相对于雷达的频率响应及雷达接收机噪声功率谱决定,干扰功率主要分配给雷达发射信号功率大、接收机噪声功率低的频段。为了在保证一定干扰性能的前提下最小化干扰***总发射功率,基于射频隐身的分布式MIMO雷达干扰波形设计方法根据注水原理进行干扰功率分配,即在雷达发射信号功率水平最大、接收机噪声功率最低所对应的频段处分配最多的功率。
图5给出了不同干扰波形设计方法下干扰总功率对比。由图5可知,随着对干扰性能要求的不断提高,机载电子干扰***发射功率逐渐增大。另外,最优干扰波形设计方法所得的干扰总功率明显小于均匀功率分配干扰波形设计方法,从而前者的射频隐身性能优于后者,这是由于均匀功率分配干扰波形是在没有任何关于雷达发射信号、目标频率响应和雷达接收机噪声功率谱等先验知识的情况下,将干扰功率均匀分配在整个频段上,因此,它具有更差的射频隐身性能。
由上述仿真结果可知,基于射频隐身的分布式MIMO雷达干扰波形设计方法,在根据先验知识获取目标相对于分布式MIMO雷达的频率响应、雷达发射信号及接收机噪声功率谱的基础上,以最小化机载电子干扰***总发射功率为目标,对干扰波形进行自适应优化设计,从而在保证一定干扰性能的条件下,有效提升机载电子干扰***的射频隐身性能。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (5)
1.一种基于射频隐身的分布式MIMO雷达干扰波形设计方法,其特征是,包括如下步骤:
1)根据先验知识,确立目标频率响应、MIMO雷达发射信号和雷达接收机噪声功率谱,并确定干扰***的辐射参数与干扰性能MI门限参数;
2)以最小化机载电子干扰***总发射功率为目标,在满足一定干扰性能的条件下,建立基于射频隐身的分布式MIMO雷达干扰波形设计模型;
3)通过拉格朗日乘数法对模型进行求解,在满足干扰性能门限δMI的最小干扰总功率的情况下,确定拉格朗日乘子的最优解;
4)将步骤3)得到的最优解代入卡罗需-库恩-塔克条件的必要条件获取分布式MIMO雷达的最佳干扰发射波形,从而形成分布式MIMO雷达干扰波形设计方案;
所述步骤1)具体内容为:
获取分布式MIMO雷达发射信号矩阵S,目标相对于雷达***的频率响应h,雷达接收机处噪声n作为先验知识;
其中,S=[S1,S2,…,SN],N为分布式MIMO雷达发射天线个数,***中第i部天线的发射信号Si为M×L维Toeplitz矩阵,M为雷达回波信号长度,L为目标频率响应长度,且有M=NL;
所述步骤2)具体内容为:
根据事先给定的干扰性能MI门限δMI,建立基于射频隐身的分布式MIMO雷达最优干扰波形设计数学模型:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711082783.2A CN108008361B (zh) | 2017-11-07 | 2017-11-07 | 基于射频隐身的分布式mimo雷达干扰波形设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711082783.2A CN108008361B (zh) | 2017-11-07 | 2017-11-07 | 基于射频隐身的分布式mimo雷达干扰波形设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108008361A CN108008361A (zh) | 2018-05-08 |
CN108008361B true CN108008361B (zh) | 2021-05-25 |
Family
ID=62051351
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711082783.2A Active CN108008361B (zh) | 2017-11-07 | 2017-11-07 | 基于射频隐身的分布式mimo雷达干扰波形设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108008361B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109164418B (zh) * | 2018-07-12 | 2021-01-05 | 南京航空航天大学 | 一种频谱共享环境下基于lpi的组网雷达稳健ofdm波形设计方法 |
CN109212494B (zh) * | 2018-07-12 | 2020-11-03 | 南京航空航天大学 | 一种针对组网雷达***的射频隐身干扰波形设计方法 |
CN109444832B (zh) * | 2018-10-25 | 2023-02-14 | 哈尔滨工程大学 | 基于多干扰效能值的群智能干扰决策方法 |
CN109669167B (zh) * | 2019-01-28 | 2023-04-07 | 南京航空航天大学 | 一种基于射频隐身的机载雷达发射波形选择方法 |
CN111142076B (zh) * | 2020-01-06 | 2023-09-05 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 提高雷达低截获性能的功率控制方法 |
CN111474526B (zh) * | 2020-04-24 | 2022-03-08 | 成都航空职业技术学院 | 一种机载stap杂波协方差矩阵的快速重建方法 |
CN111614432B (zh) * | 2020-05-26 | 2021-06-08 | 电子科技大学 | 一种设定阈值带惩罚项的干扰波形生成方法 |
CN111812599B (zh) * | 2020-06-23 | 2023-11-07 | 南京航空航天大学 | 博弈条件下基于低截获性能的组网雷达最优波形设计方法 |
CN114018101B (zh) * | 2021-12-08 | 2022-03-04 | 成都市克莱微波科技有限公司 | 一种弹载干扰装置、***及方法 |
CN115184926B (zh) * | 2022-09-13 | 2022-11-15 | 中国电子科技集团公司信息科学研究院 | 利用相参mimo雷达的分布式协同探测***及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006258449A (ja) * | 2005-03-15 | 2006-09-28 | Yokohama Rubber Co Ltd:The | レドーム |
CN106526546A (zh) * | 2016-12-09 | 2017-03-22 | 南京航空航天大学 | 一种针对雷达与通信联合***的雷达干扰功率分配方法 |
CN106685441A (zh) * | 2016-12-09 | 2017-05-17 | 南京航空航天大学 | 一种基于合作博弈的射频隐身雷达组网功率分配方法 |
CN106680780A (zh) * | 2016-12-09 | 2017-05-17 | 南京航空航天大学 | 频谱共享环境下基于射频隐身的雷达最优波形设计方法 |
-
2017
- 2017-11-07 CN CN201711082783.2A patent/CN108008361B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006258449A (ja) * | 2005-03-15 | 2006-09-28 | Yokohama Rubber Co Ltd:The | レドーム |
CN106526546A (zh) * | 2016-12-09 | 2017-03-22 | 南京航空航天大学 | 一种针对雷达与通信联合***的雷达干扰功率分配方法 |
CN106685441A (zh) * | 2016-12-09 | 2017-05-17 | 南京航空航天大学 | 一种基于合作博弈的射频隐身雷达组网功率分配方法 |
CN106680780A (zh) * | 2016-12-09 | 2017-05-17 | 南京航空航天大学 | 频谱共享环境下基于射频隐身的雷达最优波形设计方法 |
Non-Patent Citations (2)
Title |
---|
Research on MIMO Radar RF Stealth Algorithm in Searching Mode;Shuangye Zhao 等;《2014 IEEE International Conference on Signal Processing,Communications and Computing》;20140808;全文 * |
基于协同噪声干扰的机载雷达射频隐身性能优化算法;时晨光 等;《数据采集与处理》;20170331;第32卷(第2期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108008361A (zh) | 2018-05-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108008361B (zh) | 基于射频隐身的分布式mimo雷达干扰波形设计方法 | |
CN105158756B (zh) | 集中式mimo雷达射频隐身时多目标跟踪波束指向方法 | |
CN108732543B (zh) | 一种基于射频隐身的机载组网雷达辐射参数联合优化方法 | |
CN106896351B (zh) | 一种基于非合作博弈的雷达组网功率控制方法 | |
CN106909779B (zh) | 基于分布式处理的mimo雷达克拉美罗界计算方法 | |
CN107907863B (zh) | 基于雷达-通信频谱共享的组网雷达波形设计方法 | |
CN106291481A (zh) | 基于射频隐身的分布式mimo雷达资源联合优化方法 | |
CN106526546B (zh) | 一种针对雷达与通信联合***的雷达干扰功率分配方法 | |
CN111812599B (zh) | 博弈条件下基于低截获性能的组网雷达最优波形设计方法 | |
CN107944597B (zh) | 一种面对先进无源探测***的编队雷达资源管理方法 | |
CN106501778B (zh) | 一种基于射频隐身性能优化的dmrs稳健波形设计方法 | |
CN108614261B (zh) | 一种组网雷达***多目标跟踪下的辐射参数控制方法 | |
CN109459727B (zh) | 基于射频隐身的雷达-通信联合***最优波形设计方法 | |
CN107918116B (zh) | 一种基于雷达与通信***共存的多目标雷达波形设计方法 | |
Ma et al. | A novel two-dimensional sparse-weight NLMS filtering scheme for passive bistatic radar | |
CN113203993B (zh) | 针对组网雷达的自适应多干扰机协同压制干扰方法 | |
CN116299408B (zh) | 一种多雷达自主协同探测***及探测方法 | |
Karimi et al. | Adaptive OFDM waveform design for cognitive radar in signal-dependent clutter | |
CN105891799B (zh) | 适用于机械扫描雷达的有源干扰侦察方法 | |
CN103852749A (zh) | 提高mimo-stap检测性能的稳健波形优化方法 | |
CN104020459A (zh) | 一种提高mimo-stap检测性能的波形优化方法 | |
CN108037487B (zh) | 一种基于射频隐身的分布式mimo雷达发射信号优化设计方法 | |
CN113567934B (zh) | 一种基于协同的雷达探测方法 | |
CN110927692A (zh) | 海杂波场景下搜索雷达射频隐身工作方式的求解方法 | |
CN110149134A (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 |