CN114781191B - 一种复杂电磁环境的雷达探测能力仿真方法及装置 - Google Patents

一种复杂电磁环境的雷达探测能力仿真方法及装置 Download PDF

Info

Publication number
CN114781191B
CN114781191B CN202210677568.1A CN202210677568A CN114781191B CN 114781191 B CN114781191 B CN 114781191B CN 202210677568 A CN202210677568 A CN 202210677568A CN 114781191 B CN114781191 B CN 114781191B
Authority
CN
China
Prior art keywords
radar
interference
power
target
calculating
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
Application number
CN202210677568.1A
Other languages
English (en)
Other versions
CN114781191A (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.)
Aerospace Hongtu Information Technology Co Ltd
Original Assignee
Aerospace Hongtu Information Technology Co Ltd
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 Aerospace Hongtu Information Technology Co Ltd filed Critical Aerospace Hongtu Information Technology Co Ltd
Priority to CN202210677568.1A priority Critical patent/CN114781191B/zh
Publication of CN114781191A publication Critical patent/CN114781191A/zh
Application granted granted Critical
Publication of CN114781191B publication Critical patent/CN114781191B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本申请提供了一种复杂电磁环境的雷达探测能力仿真方法及装置，涉及雷达探测技术领域，具体为：对预设高度层的平面进行网格化，获取多个网格点；将雷达能够探测到的所有网格点作为目标；基于地空传播损耗模型计算目标回波功率；计算所有干扰雷达的邻频干扰功率的总和，作为总干扰功率；将目标回波功率减去雷达接收机噪声和总干扰功率，得到所述目标的雷达信干噪比；若所述目标的雷达信干噪比大于雷达接收灵敏度，则确定雷达在干扰环境下能够探测到所述目标；对预设高度层的所有目标的雷达信干噪比进行绘制，得到预设高度层的雷达威力覆盖图。本申请在雷达探测中考虑了无意干扰即其它雷达的干扰影响，提高了复杂电磁环境的雷达探测能力仿真的真实性。

Description

一种复杂电磁环境的雷达探测能力仿真方法及装置
技术领域
本申请涉及雷达探测技术领域，尤其是涉及一种复杂电磁环境的雷达探测能力仿真方法及装置。
背景技术
当多雷达同时开机时，会产生相互干扰，影响雷达探测能力。这些干扰主要包括邻频干扰、互调干扰等，会使雷达接收机输出信噪比变差，严重影响对小目标的检测。
目前未有考虑无意干扰对雷达探测能力影响的技术方案。
发明内容
有鉴于此，本申请提供了一种复杂电磁环境的雷达探测能力仿真方法及装置，以解决上述技术问题。
第一方面，本申请提供了一种复杂电磁环境的雷达探测能力仿真方法，包括：
对预设高度层的平面进行网格化，获取多个网格点；将雷达能够探测到的所有网格点作为目标；
基于地空传播损耗模型计算目标回波功率；
计算所有干扰雷达的邻频干扰功率的总和，作为总干扰功率；
将目标回波功率减去雷达接收机噪声和总干扰功率，得到所述目标的雷达信干噪比；若所述目标的雷达信干噪比大于雷达接收灵敏度，则确定雷达在干扰环境下能够探测到所述目标；
对预设高度层的所有目标的雷达信干噪比进行绘制，得到预设高度层的雷达威力覆盖图。
进一步，将雷达能够探测到的所有网格点作为目标；包括：
步骤S1：计算雷达和一个网格点的通视距离
Figure DEST_PATH_IMAGE002
：
Figure DEST_PATH_IMAGE004
其中，
Figure DEST_PATH_IMAGE006
为等效地球半径，
Figure DEST_PATH_IMAGE008
为雷达天线高度；
Figure DEST_PATH_IMAGE010
为网格点高度；
步骤S2：判断通视距离
Figure 161277DEST_PATH_IMAGE002
是否小于雷达的最大探测距离，若为是，进入步骤S3，否 则，确定雷达不能够探测该网格点处的目标；
步骤S3：将雷达和网格点之间的水平距离等间隔分段，对每个分段的端点计算遮 蔽角
Figure DEST_PATH_IMAGE012
：
Figure DEST_PATH_IMAGE014
其中，
Figure DEST_PATH_IMAGE016
为端点的障碍物海拔高度，
Figure DEST_PATH_IMAGE018
为雷达天线到障碍物顶端斜距；
步骤S4：将遮蔽角
Figure DEST_PATH_IMAGE019
加上一个修正量
Figure DEST_PATH_IMAGE021
，得到修正后的遮蔽角
Figure DEST_PATH_IMAGE023
：
Figure DEST_PATH_IMAGE025
其中，
Figure DEST_PATH_IMAGE027
，
Figure DEST_PATH_IMAGE029
为雷达波长；
进一步，基于地空传播损耗模型计算目标回波功率，包括：
目标回波功率
Figure DEST_PATH_IMAGE031
为：
Figure DEST_PATH_IMAGE033
其中，
Figure DEST_PATH_IMAGE035
为初始目标回波功率：
Figure DEST_PATH_IMAGE037
其中，
Figure DEST_PATH_IMAGE039
是大气传输损耗，根据地空传播损耗模型计算得出；
上述公式中，
Figure DEST_PATH_IMAGE041
和
Figure DEST_PATH_IMAGE043
分别为雷达发射功率和天线增益，
Figure DEST_PATH_IMAGE045
为目标雷达截面积，
Figure DEST_PATH_IMAGE047
为雷达与目标的距离。
进一步，所有干扰雷达的邻频干扰功率的计算步骤包括：
对于一个干扰雷达，首先判断下式是否成立：
Figure DEST_PATH_IMAGE049
或者
Figure DEST_PATH_IMAGE051
其中，
Figure DEST_PATH_IMAGE053
和
Figure DEST_PATH_IMAGE055
是干扰雷达发射频率和带宽，
Figure DEST_PATH_IMAGE057
和
Figure DEST_PATH_IMAGE059
是雷 达的发射频率和带宽；
若成立，则所述干扰雷达的邻频干扰对雷达无效，否则，所述干扰雷达的邻频干扰对雷达有影响，计算重叠带宽：
若
Figure DEST_PATH_IMAGE061
且
Figure DEST_PATH_IMAGE063
则重叠带宽
Figure DEST_PATH_IMAGE065
为：
Figure DEST_PATH_IMAGE067
若
Figure DEST_PATH_IMAGE068
且
Figure DEST_PATH_IMAGE070
则重叠带宽
Figure 503044DEST_PATH_IMAGE065
为：
Figure DEST_PATH_IMAGE072
若
Figure DEST_PATH_IMAGE074
且
Figure 319691DEST_PATH_IMAGE063
则重叠带宽
Figure 56702DEST_PATH_IMAGE065
为：
Figure DEST_PATH_IMAGE076
若
Figure 642404DEST_PATH_IMAGE074
且
Figure DEST_PATH_IMAGE077
则重叠带宽
Figure 636031DEST_PATH_IMAGE065
为：
Figure DEST_PATH_IMAGE079
基于重叠带宽
Figure 2290DEST_PATH_IMAGE065
，计算所述干扰雷达对雷达的邻频干扰功率
Figure DEST_PATH_IMAGE081
：
Figure DEST_PATH_IMAGE083
其中，
Figure DEST_PATH_IMAGE085
为干扰雷达的输出功率，
Figure DEST_PATH_IMAGE087
为所述干扰雷达天线增益；
Figure DEST_PATH_IMAGE089
为所述干 扰雷达与雷达之间的直线距离；
Figure DEST_PATH_IMAGE091
为所述干扰雷达辐射的干扰信号带宽；
Figure DEST_PATH_IMAGE093
为所述干 扰雷达信号在距离
Figure DEST_PATH_IMAGE095
上的单程传输损耗；
Figure DEST_PATH_IMAGE097
为雷达干扰发射损耗；
Figure DEST_PATH_IMAGE099
为雷达天线在 干扰方向的增益值，其中，
Figure DEST_PATH_IMAGE101
;
Figure DEST_PATH_IMAGE103
雷达和目标间的方位角与雷达和干扰机 间的方位角的差；
Figure DEST_PATH_IMAGE105
为雷达和目标间的仰角与雷达和干扰机间的仰角的差；
将所有对雷达有邻频干扰的干扰雷达的邻频干扰功率相加，得到所有干扰雷达的邻频干扰功率。
进一步，雷达天线在干扰方向的增益值
Figure 586068DEST_PATH_IMAGE099
的计算步骤包括：
主瓣内的方位方向图函数
Figure DEST_PATH_IMAGE107
为：
Figure DEST_PATH_IMAGE109
其中，参数
Figure DEST_PATH_IMAGE111
，
Figure DEST_PATH_IMAGE113
为单程半功率波束宽度；
主瓣外的方位方向图函数
Figure DEST_PATH_IMAGE115
为：
Figure DEST_PATH_IMAGE117
其中，参数
Figure DEST_PATH_IMAGE119
；
俯仰方向图函数
Figure DEST_PATH_IMAGE121
如下：
Figure DEST_PATH_IMAGE123
其中，
Figure DEST_PATH_IMAGE125
为低仰角的单程半功率波束宽度，
Figure DEST_PATH_IMAGE127
为最大点的仰角；
则有：
Figure DEST_PATH_IMAGE129
。
进一步，干扰雷达信号在距离
Figure DEST_PATH_IMAGE131
上的单程传输损耗
Figure DEST_PATH_IMAGE133
为：
Figure DEST_PATH_IMAGE135
其中，
Figure DEST_PATH_IMAGE137
为干扰雷达的天线高度。
进一步，干扰雷达信号在距离
Figure DEST_PATH_IMAGE138
上的单程传输损耗
Figure 728074DEST_PATH_IMAGE133
为：
Figure DEST_PATH_IMAGE139
其中，
Figure 63241DEST_PATH_IMAGE137
为干扰雷达的天线高度。
当干扰雷达的个数为2时，所述方法还包括：计算两个干扰雷达的互调干扰功率；所述总干扰功率为所有干扰雷达的邻频干扰功率的总和与互调干扰功率的和。
进一步，计算两个干扰雷达的互调干扰功率；包括：
获取干扰雷达A和干扰雷达B的频点
Figure DEST_PATH_IMAGE141
和
Figure DEST_PATH_IMAGE143
，其中干扰雷达A与雷达的距离小 于干扰雷达B与雷达的距离；
计算两个三阶互调干扰频点为
Figure DEST_PATH_IMAGE145
和
Figure DEST_PATH_IMAGE147
；
分别计算两个三阶互调干扰频点与雷达频点的差的绝对值，只要有一个差的绝对 值小于设定阈值，则判断雷达受到干扰雷达的互调干扰，互调干扰功率
Figure DEST_PATH_IMAGE149
为：
Figure DEST_PATH_IMAGE151
其中，
Figure DEST_PATH_IMAGE153
和
Figure DEST_PATH_IMAGE155
分别为干扰雷达A和干扰雷达B的发射功率，
Figure DEST_PATH_IMAGE157
为两信号三阶互 调常数，
Figure DEST_PATH_IMAGE159
是两个干扰雷达偏离雷达频点的平均值。
进一步，雷达接收机噪声功率的计算步骤包括：
计算雷达接收机噪声平均值
Figure DEST_PATH_IMAGE161
：
Figure DEST_PATH_IMAGE163
其中，
Figure DEST_PATH_IMAGE165
为波尔茨曼常数，基准温度
Figure DEST_PATH_IMAGE167
，
Figure DEST_PATH_IMAGE169
为雷达工作瞬时带宽，
Figure DEST_PATH_IMAGE171
为噪声系数；
则雷达接收机瞬时噪声功率
Figure DEST_PATH_IMAGE173
为：
Figure DEST_PATH_IMAGE175
其中，
Figure DEST_PATH_IMAGE177
和
Figure DEST_PATH_IMAGE179
均为服从标准正态分布的随机数。
第二方面，本申请实施例提供了一种复杂电磁环境的雷达探测能力仿真装置，包括：
目标确定单元，用于对预设高度层的平面进行网格化，获取多个网格点；将雷达能够探测到的所有网格点作为目标；
第一计算单元，用于基于地空传播损耗模型计算目标回波功率；
第二计算单元，用于计算所有干扰雷达的邻频干扰功率的总和，作为总干扰功率；
第三计算单元，用于将目标回波功率减去雷达接收机噪声和总干扰功率，得到所述目标的雷达信干噪比；若所述目标的雷达信干噪比大于雷达接收灵敏度，则确定雷达在干扰环境下能够探测到所述目标；
雷达威力覆盖图绘制单元，用于对预设高度层的所有目标的雷达信干噪比进行绘制，得到预设高度层的雷达威力覆盖图。
本申请在雷达探测中考虑了无意干扰即其它雷达的干扰影响，提高了复杂电磁环境的雷达探测能力仿真的真实性。
附图说明
为了更清楚地说明本申请具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本申请的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本申请实施例提供的复杂电磁环境的雷达探测能力仿真方法的流程图；
图2为本申请实施例提供的复杂电磁环境的雷达探测能力仿真装置的功能结构图；
图3为本申请实施例提供的电子设备的结构图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本申请一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本申请实施例的组件可以以各种不同的配置来布置和设计。
因此，以下对在附图中提供的本申请的实施例的详细描述并非旨在限制要求保护的本申请的范围，而是仅仅表示本申请的选定实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
首先对本申请实施例的设计思想进行简单介绍。
当多雷达同时开机时，会产生相互干扰，影响雷达探测能力。这些干扰主要包括邻频干扰、互调干扰等，会使雷达接收机输出信噪比变差，严重影响对小目标的检测。而现有的雷达功能仿真对这方面考虑较少，所以有必要对这些干扰进行量化分析。
为解决上述技术问题，本申请提供了一种复杂电磁环境的雷达探测能力仿真方法，其应用场景为：地面和空中部署多部雷达，雷达的性能参数可配置，如位置、技术参数等。多雷达工作时互相干扰，分析干扰对雷达探测性能的影响，同时考虑通视和传播损耗的影响；可以得出雷达在无意干扰的情况下，是否能探测到目标，同时计算出回波信干噪比，最终绘出雷达在不同高度层的威力覆盖图，用不同颜色来区分信干噪比大小。
本申请在无意干扰功率计算时，考虑邻频干扰和互调干扰；邻频干扰考虑主雷达和干扰雷达带宽的重合度；互调干扰主要考虑三阶互调干扰，可通过经验公式得到干扰功率。两种干扰都用到了雷达干扰方程和损耗模型。雷达干扰方程中，主要是计算雷达天线在干扰雷达方向上的增益。干扰功率损耗值可根据应用场景，调用不同损耗模型计算。计算出各自干扰功率后，直接相加即可得到总干扰功率。
此外，本申请在计算雷达回波功率时，针对各种损耗，大多数仿真都是直接给出近似值。而本申请建立了精确的损耗数学模型，依据不同情况给出不同损耗值，计算精度有了很大提高。
在介绍了本申请实施例的应用场景和设计思想之后，下面对本申请实施例提供的技术方案进行说明。
如图1所示，本申请实施例提供了一种复杂电磁环境的雷达探测能力仿真方法，包括如下步骤：
步骤101：对预设高度层的平面进行网格化，获取多个网格点；将雷达能够探测到的所有网格点作为目标；
其中，预设高度层的平面垂直于雷达所在的地平面，与雷达的距离为预设的高度H。为了研究雷达在预设高度层平面的探测能力，首先要选取多个平面上的点作为目标，首先对预设高度层的平面进行网格化，获取多个网格点，但是并不是每个网格点所在的目标都能被雷达探测到，因此对于每个网格点计算视线通视和电磁通视分布情况；地球曲率对电磁波的传播有遮蔽作用，导致在低仰角情况下，雷达对目标的探测受视距影响。当雷达和目标之间是通视的，还需要考虑地形地物对雷达探测范围的影响，其对目标的遮挡通常用遮蔽角来衡量。遮蔽角度定义为雷达天线和障碍物顶端的连线和水平线之间的夹角，它是某一方位上雷达的最小仰角。
由于电磁波在障碍物附近传播会产生绕射现象，当目标和雷达之间仰角与遮蔽角接近时，绕射会产生很大的衰减，故有必要对遮蔽角加上一个修正量。
在本实施例中，将雷达能够探测到的所有网格点作为目标；包括：
步骤S1：计算雷达和一个网格点的通视距离
Figure 854479DEST_PATH_IMAGE002
：
Figure DEST_PATH_IMAGE180
其中，
Figure 933294DEST_PATH_IMAGE006
为等效地球半径，
Figure 493588DEST_PATH_IMAGE008
为雷达天线高度；
Figure 265235DEST_PATH_IMAGE010
为网格点高度；
步骤S2：判断通视距离
Figure 278190DEST_PATH_IMAGE002
是否小于雷达的最大探测距离，若为是，进入步骤S3，否 则，确定雷达不能够探测该网格点处的目标；
步骤S3：将雷达和网格点之间的水平距离等间隔分段，对每个分段的端点计算遮 蔽角
Figure 160696DEST_PATH_IMAGE012
：
Figure DEST_PATH_IMAGE181
其中，
Figure 575496DEST_PATH_IMAGE016
为端点的障碍物海拔高度，
Figure 518045DEST_PATH_IMAGE018
为雷达天线到障碍物顶端斜距；
步骤S4：将遮蔽角
Figure 955979DEST_PATH_IMAGE019
加上一个修正量
Figure 937448DEST_PATH_IMAGE021
，得到修正后的遮蔽角
Figure DEST_PATH_IMAGE182
：
Figure DEST_PATH_IMAGE183
其中，
Figure 472335DEST_PATH_IMAGE027
，
Figure 585784DEST_PATH_IMAGE029
为雷达波长；
步骤S5：判断雷达与目标间的仰角是否大于每个修正后的遮蔽角
Figure 573332DEST_PATH_IMAGE182
，若为是，则 确定雷达能够探测该网格点处的目标，否则，确定雷达不能够探测该网格点处的目标。
步骤102：基于地空传播损耗模型计算目标回波功率；
雷达电波在空间传播，会有各种能量损耗，包括大气分子和水汽及由水汽凝聚成的降落颗粒引起的损耗，电波绕过地形表面或障碍物的绕射引起的损耗。实际雷达工作时的传播衰减与雷达的作用距离以及目标高度有关。实际仿真中，把各个典型频段传播损耗预先利用地空传播损耗模型计算出结果，使用时直接调用即可。
目标回波功率
Figure 797640DEST_PATH_IMAGE031
为：
Figure DEST_PATH_IMAGE184
其中，
Figure 921454DEST_PATH_IMAGE035
为初始目标回波功率：
Figure DEST_PATH_IMAGE185
其中，
Figure DEST_PATH_IMAGE186
是大气传输损耗，根据地空传播损耗模型计算得出；
上述公式中，
Figure DEST_PATH_IMAGE187
和
Figure DEST_PATH_IMAGE188
分别为雷达发射功率和天线增益，
Figure 894220DEST_PATH_IMAGE045
为目标雷达截面积，
Figure 369064DEST_PATH_IMAGE047
为雷达与目标的距离。
步骤103：计算所有干扰雷达的邻频干扰功率的总和，作为总干扰功率；
邻频干扰指多个雷达辐射源和己方雷达频段接近，空间相隔较近，发射电波能进入己方雷达接收机内部，从而对有用信号造成较大干扰。所以有必要给出定量化的分析。
在本实施例中，所有干扰雷达的邻频干扰功率的计算步骤包括：
对于一个干扰雷达，首先判断下式是否成立：
Figure DEST_PATH_IMAGE189
或者
Figure 957915DEST_PATH_IMAGE051
其中，
Figure DEST_PATH_IMAGE190
和
Figure DEST_PATH_IMAGE191
是干扰雷达发射频率和带宽，
Figure 201814DEST_PATH_IMAGE057
和
Figure 391487DEST_PATH_IMAGE059
是雷 达的发射频率和带宽；
若成立，则所述干扰雷达的邻频干扰对雷达无效，否则，所述干扰雷达的邻频干扰对雷达有影响，计算重叠带宽：
若
Figure 353627DEST_PATH_IMAGE068
且
Figure 919737DEST_PATH_IMAGE063
则重叠带宽
Figure DEST_PATH_IMAGE192
为：
Figure 283723DEST_PATH_IMAGE067
若
Figure 644297DEST_PATH_IMAGE068
且
Figure 595197DEST_PATH_IMAGE077
则重叠带宽
Figure DEST_PATH_IMAGE193
为：
Figure 964999DEST_PATH_IMAGE072
若
Figure 917911DEST_PATH_IMAGE074
且
Figure 714966DEST_PATH_IMAGE063
则重叠带宽
Figure 323802DEST_PATH_IMAGE193
为：
Figure 559611DEST_PATH_IMAGE076
若
Figure 632610DEST_PATH_IMAGE074
且
Figure 600566DEST_PATH_IMAGE077
则重叠带宽
Figure 257550DEST_PATH_IMAGE193
为：
Figure 969154DEST_PATH_IMAGE079
基于重叠带宽
Figure 834341DEST_PATH_IMAGE193
，计算所述干扰雷达对雷达的邻频干扰功率
Figure 35516DEST_PATH_IMAGE081
：
Figure DEST_PATH_IMAGE194
其中：
Figure 681261DEST_PATH_IMAGE081
的单位为W；
Figure 930976DEST_PATH_IMAGE085
为干扰雷达的输出功率，
Figure 542349DEST_PATH_IMAGE087
为所述干扰雷达天线增 益；
Figure 852107DEST_PATH_IMAGE089
为所述干扰雷达与雷达之间的直线距离；
Figure 985148DEST_PATH_IMAGE091
为所述干扰雷达辐射的干扰信号带 宽，单位为Hz；
Figure 38555DEST_PATH_IMAGE093
为所述干扰雷达信号在距离
Figure 940652DEST_PATH_IMAGE095
上的单程传输损耗；
Figure 483629DEST_PATH_IMAGE097
为雷达干扰发 射损耗；
Figure 776070DEST_PATH_IMAGE099
为雷达天线在干扰方向的增益值，其中，
Figure 217457DEST_PATH_IMAGE101
;
Figure 911743DEST_PATH_IMAGE103
雷达和目 标间的方位角与雷达和干扰机间的方位角的差；
Figure 625621DEST_PATH_IMAGE105
为雷达和目标间的仰角与雷达和干扰 机间的仰角的差；
将所有对雷达有邻频干扰的干扰雷达的邻频干扰功率相加，得到所有干扰雷达的邻频干扰功率。
其中，
Figure 670938DEST_PATH_IMAGE133
采用地地传播损耗模型计算。地地传播主要考虑地形的影响，如下式 所示：
Figure 66147DEST_PATH_IMAGE135
其中，
Figure 677257DEST_PATH_IMAGE137
为干扰雷达的天线高度。可见，路径损耗与传播距离的四次方成正比，随 传播距离的增加损耗迅速增大。
将所有对雷达有邻频干扰的干扰雷达的邻频干扰功率相加，得到所有干扰雷达的邻频干扰功率。
在邻频干扰基础上，若相邻雷达的频点和主雷达频点比较接近，很容易对主雷达产生互调干扰。它产生的原因是因为多个频率分量进入接收机后，由于雷达接收机非线性器件相互作用，产生的接近主雷达频点的其它频率分量，而这种频率分量接收机很难抑制，从而对有用信号形成比较大的干扰。这里主要考虑三阶互调干扰，因为这种干扰产生的能量最大，能对有用信号产生很大危害。
在本实施例中，只考虑两个干扰雷达的情况下产生的互调干扰，具体的，计算两个干扰雷达的互调干扰功率；包括：
获取干扰雷达A和干扰雷达B的频点
Figure 765299DEST_PATH_IMAGE141
和
Figure 94649DEST_PATH_IMAGE143
，其中干扰雷达A与雷达的距离小 于干扰雷达B与雷达的距离；
计算两个三阶互调干扰频点为
Figure 293549DEST_PATH_IMAGE145
和
Figure 696849DEST_PATH_IMAGE147
；
分别计算两个三阶互调干扰频点与雷达频点的差的绝对值，只要有一个差的绝对 值小于设定阈值，则判断雷达受到干扰雷达的互调干扰，互调干扰功率
Figure 519573DEST_PATH_IMAGE149
为：
Figure 273903DEST_PATH_IMAGE151
其中，
Figure 135548DEST_PATH_IMAGE153
和
Figure 658934DEST_PATH_IMAGE155
分别为干扰雷达A和干扰雷达B的发射功率，
Figure 88778DEST_PATH_IMAGE157
为两信号三阶互 调常数，
Figure 127141DEST_PATH_IMAGE159
是两个干扰雷达偏离雷达频点的平均值。
此时，所述总干扰功率为所有干扰雷达的邻频干扰功率的总和与互调干扰功率的和。
步骤104：将目标回波功率减去雷达接收机噪声和总干扰功率，得到所述目标的雷达信干噪比；若所述目标的雷达信干噪比大于雷达接收灵敏度，则确定雷达在干扰环境下能够探测到所述目标；
噪声包括接收机产生的噪声和大气噪声参数，具体的，雷达接收机噪声功率的计算步骤包括：
计算雷达接收机噪声平均值
Figure 667844DEST_PATH_IMAGE161
：
Figure 872167DEST_PATH_IMAGE163
其中，
Figure 207333DEST_PATH_IMAGE165
为波尔茨曼常数，基准温度
Figure 998572DEST_PATH_IMAGE167
，
Figure 77386DEST_PATH_IMAGE169
为雷达工作瞬时带宽，
Figure 575364DEST_PATH_IMAGE171
为噪声系数；
则雷达接收机瞬时噪声功率
Figure 143748DEST_PATH_IMAGE173
为：
Figure DEST_PATH_IMAGE195
其中，
Figure 422283DEST_PATH_IMAGE177
和
Figure 304788DEST_PATH_IMAGE179
均为服从标准正态分布的随机数。
将目标回波功率，雷达接收机噪声和总干扰功率均取dB值后，计算目标回波功率减去雷达接收机噪声和总干扰功率后的值，作为所述目标的雷达信干噪比（单位为dB）。
步骤105：对预设高度层的所有目标的雷达信干噪比进行绘制，得到预设高度层的雷达威力覆盖图。
此外，可以用颜色来区分不同位置的信干噪比大小。取几个典型高度层，重复同样步骤，就能够得到多个雷达整体威力覆盖图。
基于上述实施例，本申请实施例提供了一种复杂电磁环境的雷达探测能力仿真装置，参阅图2所示，本申请实施例提供的复杂电磁环境的雷达探测能力仿真装置200至少包括：
目标确定单元201，用于对预设高度层的平面进行网格化，获取多个网格点；将雷达能够探测到的所有网格点作为目标；
第一计算单元202，用于基于地空传播损耗模型计算目标回波功率；
第二计算单元203，用于计算所有干扰雷达的邻频干扰功率的总和，作为总干扰功率；
第三计算单元204，用于将目标回波功率减去雷达接收机噪声和总干扰功率，得到所述目标的雷达信干噪比；若所述目标的雷达信干噪比大于雷达接收灵敏度，则确定雷达在干扰环境下能够探测到所述目标；
雷达威力覆盖图绘制单元205，用于对预设高度层的所有目标的雷达信干噪比进行绘制，得到预设高度层的雷达威力覆盖图。
需要说明的是，本申请实施例提供的复杂电磁环境的雷达探测能力仿真装置200解决技术问题的原理与本申请实施例提供的复杂电磁环境的雷达探测能力仿真方法相似，因此，本申请实施例提供的复杂电磁环境的雷达探测能力仿真装置200的实施可以参见本申请实施例提供的复杂电磁环境的雷达探测能力仿真方法的实施，重复之处不再赘述。
如图3所示，本申请实施例提供的电子设备300至少包括：处理器301、存储器302和存储在存储器302上并可在处理器301上运行的计算机程序，处理器301执行计算机程序时实现本申请实施例提供的复杂电磁环境的雷达探测能力仿真方法。
本申请实施例提供的电子设备300还可以包括连接不同组件（包括处理器301和存储器302）的总线303。其中，总线303表示几类总线结构中的一种或多种，包括存储器总线、外围总线、局域总线等。
存储器302可以包括易失性存储器形式的可读介质，例如随机存储器（RandomAccess Memory，RAM）3021和/或高速缓存存储器3022，还可以进一步包括只读存储器（ReadOnly Memory，ROM）3023。
存储器302还可以包括具有一组（至少一个）程序模块3025的程序工具3024，程序模块3025包括但不限于：操作子系统、一个或者多个应用程序、其它程序模块以及程序数据，这些示例中的每一个或某种组合中可能包括网络环境的实现。
电子设备300也可以与一个或多个外部设备303（例如键盘、遥控器等）通信，还可以与一个或者多个使得用户能与电子设备300交互的设备通信（例如手机、电脑等），和/或，与使得电子设备300与一个或多个其它电子设备300进行通信的任何设备（例如路由器、调制解调器等）通信。这种通信可以通过输入/输出（Input /Output，I/O）接口305进行。并且，电子设备300还可以通过网络适配器306与一个或者多个网络（例如局域网（Local AreaNetwork，LAN），广域网（Wide Area Network，WAN）和/或公共网络，例如因特网）通信。如图3所示，网络适配器306通过总线303与电子设备300的其它模块通信。应当理解，尽管图3中未示出，可以结合电子设备300使用其它硬件和/或软件模块，包括但不限于：微代码、设备驱动器、冗余处理器、外部磁盘驱动阵列、磁盘阵列（Redundant Arrays of IndependentDisks，RAID）子系统、磁带驱动器以及数据备份存储子系统等。
需要说明的是，图3所示的电子设备300仅仅是一个示例，不应对本申请实施例的功能和使用范围带来任何限制。
本申请实施例还提供了一种计算机可读存储介质，该计算机可读存储介质存储有计算机指令，该计算机指令被处理器执行时实现本申请实施例提供的复杂电磁环境的雷达探测能力仿真方法。
此外，尽管在附图中以特定顺序描述了本申请方法的操作，但是，这并非要求或者暗示必须按照该特定顺序来执行这些操作，或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地，可以省略某些步骤，将多个步骤合并为一个步骤执行，和/或将一个步骤分解为多个步骤执行。
尽管已描述了本申请的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本申请范围的所有变更和修改。
最后应说明的是：以上各实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述各实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的范围。

Claims (6)

1.一种复杂电磁环境的雷达探测能力仿真方法，其特征在于，包括：
对预设高度层的平面进行网格化，获取多个网格点；将雷达能够探测到的所有网格点作为目标；
计算目标回波功率；
计算所有干扰雷达的邻频干扰功率的总和，作为总干扰功率；
将目标回波功率减去雷达接收机噪声和总干扰功率，得到所述目标的雷达信干噪比；若所述目标的雷达信干噪比大于雷达接收灵敏度，则确定雷达在干扰环境下能够探测到所述目标；
对预设高度层的所有目标的雷达信干噪比进行绘制，得到预设高度层的雷达威力覆盖图；
计算目标回波功率，包括：
目标回波功率
Figure 694779DEST_PATH_IMAGE001
为：
Figure 917949DEST_PATH_IMAGE002
其中，
Figure 717278DEST_PATH_IMAGE003
为初始目标回波功率：
Figure 899385DEST_PATH_IMAGE004
其中，
Figure 798071DEST_PATH_IMAGE005
是大气传输损耗，根据地空传播损耗模型计算得出；
其中，
Figure 570855DEST_PATH_IMAGE006
和
Figure 783662DEST_PATH_IMAGE007
分别为雷达发射功率和天线增益，
Figure 895974DEST_PATH_IMAGE008
为目标雷达截面积，
Figure 293458DEST_PATH_IMAGE009
为雷达与目标的距离；
所有干扰雷达的邻频干扰功率的计算步骤包括：
对于一个干扰雷达，首先判断下式是否成立：
Figure 412592DEST_PATH_IMAGE010
或者
Figure 491407DEST_PATH_IMAGE011
其中，
Figure 786122DEST_PATH_IMAGE012
和
Figure 26610DEST_PATH_IMAGE013
是干扰雷达发射频率和带宽，
Figure 649353DEST_PATH_IMAGE014
和
Figure 328596DEST_PATH_IMAGE015
是雷达的发射频率和带宽；
若成立，则所述干扰雷达的邻频干扰对雷达无效，否则，所述干扰雷达的邻频干扰对雷达有影响，计算重叠带宽：
若
Figure 415500DEST_PATH_IMAGE016
且
Figure 13841DEST_PATH_IMAGE017
则重叠带宽
Figure 982934DEST_PATH_IMAGE018
为：
Figure 403551DEST_PATH_IMAGE019
若
Figure 282645DEST_PATH_IMAGE016
且
Figure 130515DEST_PATH_IMAGE020
则重叠带宽
Figure 852484DEST_PATH_IMAGE018
为：
Figure 732584DEST_PATH_IMAGE021
若
Figure 794081DEST_PATH_IMAGE022
且
Figure 875169DEST_PATH_IMAGE017
则重叠带宽
Figure 756538DEST_PATH_IMAGE018
为：
Figure 191061DEST_PATH_IMAGE023
若
Figure 169382DEST_PATH_IMAGE022
且
Figure 359054DEST_PATH_IMAGE024
则重叠带宽
Figure 194897DEST_PATH_IMAGE018
为：
Figure 292166DEST_PATH_IMAGE025
基于重叠带宽
Figure 62676DEST_PATH_IMAGE018
，计算所述干扰雷达对雷达的邻频干扰功率
Figure 95354DEST_PATH_IMAGE026
：
Figure 13632DEST_PATH_IMAGE027
其中，
Figure 117854DEST_PATH_IMAGE028
为干扰雷达的输出功率，
Figure 195400DEST_PATH_IMAGE029
为所述干扰雷达天线增益；
Figure 726876DEST_PATH_IMAGE030
为所述干扰雷达与雷达之间的直线距离；
Figure 132449DEST_PATH_IMAGE031
为所述干扰雷达辐射的干扰信号带宽；
Figure 446887DEST_PATH_IMAGE032
为所述干扰雷达信号在距离
Figure 191989DEST_PATH_IMAGE033
上的单程传输损耗；
Figure 956683DEST_PATH_IMAGE034
为雷达干扰发射损耗；
Figure 708607DEST_PATH_IMAGE035
为雷达天线在干扰方向的增益值，其中，
Figure 154632DEST_PATH_IMAGE036
;
Figure 816557DEST_PATH_IMAGE037
雷达和目标间的方位角与雷达和干扰机间的方位角的差；
Figure 627519DEST_PATH_IMAGE038
为雷达和目标间的仰角与雷达和干扰机间的仰角的差；
将所有对雷达有邻频干扰的干扰雷达的邻频干扰功率相加，得到所有干扰雷达的邻频干扰功率；
当干扰雷达的个数为2时，所述方法还包括：计算两个干扰雷达的互调干扰功率；所述总干扰功率为所有干扰雷达的邻频干扰功率的总和与互调干扰功率的和；
计算两个干扰雷达的互调干扰功率；包括：
获取干扰雷达A和干扰雷达B的频点
Figure 679788DEST_PATH_IMAGE039
和
Figure 991821DEST_PATH_IMAGE040
，其中干扰雷达A与雷达的距离小于干扰雷达B与雷达的距离；
计算两个三阶互调干扰频点为
Figure 632887DEST_PATH_IMAGE041
和
Figure 677066DEST_PATH_IMAGE042
；
分别计算两个三阶互调干扰频点与雷达频点的差的绝对值，只要有一个差的绝对值小于设定阈值，则判断雷达受到干扰雷达的互调干扰，互调干扰功率
Figure 278949DEST_PATH_IMAGE043
为：
Figure 66776DEST_PATH_IMAGE044
其中，
Figure 578660DEST_PATH_IMAGE045
和
Figure 590478DEST_PATH_IMAGE046
分别为干扰雷达A和干扰雷达B的发射功率，
Figure 148499DEST_PATH_IMAGE047
为两信号三阶互调常数，
Figure 398739DEST_PATH_IMAGE048
是两个干扰雷达偏离雷达频点的平均值。
2.根据权利要求1所述的复杂电磁环境的雷达探测能力仿真方法，其特征在于，将雷达能够探测到的所有网格点作为目标；包括：
步骤S1：计算雷达和一个网格点的通视距离
Figure 889763DEST_PATH_IMAGE049
：
Figure 275745DEST_PATH_IMAGE050
其中，
Figure 727586DEST_PATH_IMAGE051
为等效地球半径，
Figure 653954DEST_PATH_IMAGE052
为雷达天线高度；
Figure 202747DEST_PATH_IMAGE053
为网格点高度；
步骤S2：判断通视距离
Figure 946581DEST_PATH_IMAGE054
是否小于雷达的最大探测距离，若为是，进入步骤S3，否则，确定雷达不能够探测该网格点处的目标；
步骤S3：将雷达和网格点之间的水平距离等间隔分段，对每个分段的端点计算遮蔽角
Figure 948035DEST_PATH_IMAGE055
：
Figure 943673DEST_PATH_IMAGE056
其中，
Figure 19076DEST_PATH_IMAGE057
为端点的障碍物海拔高度，
Figure 746861DEST_PATH_IMAGE058
为雷达天线到障碍物顶端斜距；
步骤S4：将遮蔽角
Figure 32348DEST_PATH_IMAGE055
加上一个修正量
Figure 956311DEST_PATH_IMAGE059
，得到修正后的遮蔽角
Figure 10855DEST_PATH_IMAGE060
：
Figure 909541DEST_PATH_IMAGE061
其中，
Figure 823270DEST_PATH_IMAGE062
，
Figure 98393DEST_PATH_IMAGE063
为雷达波长；
步骤S5：判断雷达与目标间的仰角是否大于每个修正后的遮蔽角
Figure 7444DEST_PATH_IMAGE060
，若为是，则确定雷达能够探测该网格点处的目标，否则，确定雷达不能够探测该网格点处的目标。
3.根据权利要求1所述的复杂电磁环境的雷达探测能力仿真方法，其特征在于，雷达天线在干扰方向的增益值
Figure 529561DEST_PATH_IMAGE035
的计算步骤包括：
主瓣内的方位方向图函数
Figure 727324DEST_PATH_IMAGE064
为：
Figure 602876DEST_PATH_IMAGE065
其中，参数
Figure 507378DEST_PATH_IMAGE066
，
Figure 13446DEST_PATH_IMAGE067
为单程半功率波束宽度；
主瓣外的方位方向图函数
Figure 760822DEST_PATH_IMAGE068
为：
Figure 36470DEST_PATH_IMAGE069
其中，参数
Figure 185692DEST_PATH_IMAGE070
；
俯仰方向图函数
Figure 862661DEST_PATH_IMAGE071
如下：
Figure 972699DEST_PATH_IMAGE072
其中，
Figure 190054DEST_PATH_IMAGE073
为低仰角的单程半功率波束宽度，
Figure 131465DEST_PATH_IMAGE074
为最大点的仰角；
则有：
Figure 166286DEST_PATH_IMAGE075
。
4.根据权利要求3所述的复杂电磁环境的雷达探测能力仿真方法，其特征在于，干扰雷达信号在距离
Figure 825937DEST_PATH_IMAGE076
上的单程传输损耗
Figure 581404DEST_PATH_IMAGE077
为：
Figure 580584DEST_PATH_IMAGE078
其中，
Figure 864935DEST_PATH_IMAGE079
为干扰雷达的天线高度。
5.根据权利要求1所述的复杂电磁环境的雷达探测能力仿真方法，其特征在于，雷达接收机噪声功率的计算步骤包括：
计算雷达接收机噪声平均值
Figure 198833DEST_PATH_IMAGE080
：
Figure 961252DEST_PATH_IMAGE081
其中，
Figure 939573DEST_PATH_IMAGE082
为波尔茨曼常数，基准温度
Figure 129246DEST_PATH_IMAGE083
，
Figure 170014DEST_PATH_IMAGE084
为雷达工作瞬时带宽，
Figure 532862DEST_PATH_IMAGE085
为噪声系数；
则雷达接收机瞬时噪声功率
Figure 568951DEST_PATH_IMAGE086
为：
Figure 116476DEST_PATH_IMAGE087
其中，
Figure 769174DEST_PATH_IMAGE088
和
Figure 873397DEST_PATH_IMAGE089
均为服从标准正态分布的随机数。
6.一种复杂电磁环境的雷达探测能力仿真装置，其特征在于，包括：
目标确定单元，用于对预设高度层的平面进行网格化，获取多个网格点；将雷达能够探测到的所有网格点作为目标；
第一计算单元，用于计算目标回波功率；
第二计算单元，用于计算所有干扰雷达的邻频干扰功率的总和，作为总干扰功率；
第三计算单元，用于将目标回波功率减去雷达接收机噪声和总干扰功率，得到所述目标的雷达信干噪比；若所述目标的雷达信干噪比大于雷达接收灵敏度，则确定雷达在干扰环境下能够探测到所述目标；
雷达威力覆盖图绘制单元，用于对预设高度层的所有目标的雷达信干噪比进行绘制，得到预设高度层的雷达威力覆盖图；
计算目标回波功率，包括：
目标回波功率
Figure 436096DEST_PATH_IMAGE001
为：
Figure 29889DEST_PATH_IMAGE002
其中，
Figure 373145DEST_PATH_IMAGE003
为初始目标回波功率：
Figure 939780DEST_PATH_IMAGE004
其中，
Figure 747199DEST_PATH_IMAGE005
是大气传输损耗，根据地空传播损耗模型计算得出；
其中，
Figure 449576DEST_PATH_IMAGE006
和
Figure 217812DEST_PATH_IMAGE007
分别为雷达发射功率和天线增益，
Figure 460574DEST_PATH_IMAGE008
为目标雷达截面积，
Figure 60183DEST_PATH_IMAGE009
为雷达与目标的距离；
所有干扰雷达的邻频干扰功率的计算步骤包括：
对于一个干扰雷达，首先判断下式是否成立：
Figure 120412DEST_PATH_IMAGE090
或者
Figure 438261DEST_PATH_IMAGE011
其中，
Figure 750293DEST_PATH_IMAGE012
和
Figure 876512DEST_PATH_IMAGE013
是干扰雷达发射频率和带宽，
Figure 920692DEST_PATH_IMAGE014
和
Figure 788153DEST_PATH_IMAGE015
是雷达的发射频率和带宽；
若成立，则所述干扰雷达的邻频干扰对雷达无效，否则，所述干扰雷达的邻频干扰对雷达有影响，计算重叠带宽：
若
Figure 575981DEST_PATH_IMAGE016
且
Figure 71553DEST_PATH_IMAGE017
则重叠带宽
Figure 348951DEST_PATH_IMAGE018
为：
Figure 641392DEST_PATH_IMAGE019
若
Figure 905014DEST_PATH_IMAGE016
且
Figure 130459DEST_PATH_IMAGE020
则重叠带宽
Figure 782020DEST_PATH_IMAGE018
为：
Figure 483129DEST_PATH_IMAGE021
若
Figure 409497DEST_PATH_IMAGE022
且
Figure 958290DEST_PATH_IMAGE091
则重叠带宽
Figure 452856DEST_PATH_IMAGE018
为：
Figure 454310DEST_PATH_IMAGE092
若
Figure 449948DEST_PATH_IMAGE022
且
Figure 765830DEST_PATH_IMAGE020
则重叠带宽
Figure 493614DEST_PATH_IMAGE018
为：
Figure 44681DEST_PATH_IMAGE025
基于重叠带宽
Figure 516114DEST_PATH_IMAGE018
，计算所述干扰雷达对雷达的邻频干扰功率
Figure 446024DEST_PATH_IMAGE026
：
Figure 672606DEST_PATH_IMAGE027
其中，
Figure 648652DEST_PATH_IMAGE028
为干扰雷达的输出功率，
Figure 845147DEST_PATH_IMAGE029
为所述干扰雷达天线增益；
Figure 223039DEST_PATH_IMAGE030
为所述干扰雷达与雷达之间的直线距离；
Figure 354943DEST_PATH_IMAGE031
为所述干扰雷达辐射的干扰信号带宽；
Figure 490389DEST_PATH_IMAGE032
为所述干扰雷达信号在距离
Figure 100362DEST_PATH_IMAGE093
上的单程传输损耗；
Figure 254132DEST_PATH_IMAGE034
为雷达干扰发射损耗；
Figure 760199DEST_PATH_IMAGE035
为雷达天线在干扰方向的增益值，其中，
Figure 507576DEST_PATH_IMAGE036
;
Figure 124502DEST_PATH_IMAGE037
雷达和目标间的方位角与雷达和干扰机间的方位角的差；
Figure 149090DEST_PATH_IMAGE038
为雷达和目标间的仰角与雷达和干扰机间的仰角的差；
将所有对雷达有邻频干扰的干扰雷达的邻频干扰功率相加，得到所有干扰雷达的邻频干扰功率；
当干扰雷达的个数为2时，所述装置还用于：计算两个干扰雷达的互调干扰功率；所述总干扰功率为所有干扰雷达的邻频干扰功率的总和与互调干扰功率的和；
计算两个干扰雷达的互调干扰功率；包括：
获取干扰雷达A和干扰雷达B的频点
Figure 622796DEST_PATH_IMAGE039
和
Figure 795152DEST_PATH_IMAGE094
，其中干扰雷达A与雷达的距离小于干扰雷达B与雷达的距离；
计算两个三阶互调干扰频点为
Figure 137140DEST_PATH_IMAGE041
和
Figure 140868DEST_PATH_IMAGE042
；
分别计算两个三阶互调干扰频点与雷达频点的差的绝对值，只要有一个差的绝对值小于设定阈值，则判断雷达受到干扰雷达的互调干扰，互调干扰功率
Figure 988738DEST_PATH_IMAGE043
为：
Figure 586073DEST_PATH_IMAGE044
其中，
Figure 544802DEST_PATH_IMAGE045
和
Figure 403036DEST_PATH_IMAGE046
分别为干扰雷达A和干扰雷达B的发射功率，
Figure 346109DEST_PATH_IMAGE047
为两信号三阶互调常数，
Figure 493057DEST_PATH_IMAGE048
是两个干扰雷达偏离雷达频点的平均值。
CN202210677568.1A 2022-06-16 2022-06-16 一种复杂电磁环境的雷达探测能力仿真方法及装置 Active CN114781191B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210677568.1A CN114781191B (zh) 2022-06-16 2022-06-16 一种复杂电磁环境的雷达探测能力仿真方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210677568.1A CN114781191B (zh) 2022-06-16 2022-06-16 一种复杂电磁环境的雷达探测能力仿真方法及装置

Publications (2)

Publication Number Publication Date
CN114781191A CN114781191A (zh) 2022-07-22
CN114781191B true CN114781191B (zh) 2022-09-09

Family

ID=82422178

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210677568.1A Active CN114781191B (zh) 2022-06-16 2022-06-16 一种复杂电磁环境的雷达探测能力仿真方法及装置

Country Status (1)

Country Link
CN (1) CN114781191B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116224260A (zh) * 2023-05-06 2023-06-06 成都众享天地网络科技有限公司 一种基于干扰的雷达三维威力范围计算方法
CN117315182B (zh) * 2023-11-30 2024-03-22 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) 一种用于确定在地形影响下的雷达探测范围的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102298140A (zh) * 2011-05-16 2011-12-28 中国人民解放军海军航空工程学院 基于探测概率几何分布曲线的雷达网抗干扰效能评估方法
CN113406579A (zh) * 2021-06-07 2021-09-17 哈尔滨工程大学 一种基于深度强化学习的伪装干扰波形生成方法
CN114265066A (zh) * 2021-12-27 2022-04-01 中国人民解放军63891部队 运动船舶的sar成像模拟方法、装置及电子设备
CN114527436A (zh) * 2022-02-16 2022-05-24 西安电子科技大学 基于最优干扰规划模型的集群协同干扰规划方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101551811B1 (ko) * 2014-05-19 2015-09-10 최수호 레이더 장치 및 그의 주파수 간섭 제거방법

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102298140A (zh) * 2011-05-16 2011-12-28 中国人民解放军海军航空工程学院 基于探测概率几何分布曲线的雷达网抗干扰效能评估方法
CN113406579A (zh) * 2021-06-07 2021-09-17 哈尔滨工程大学 一种基于深度强化学习的伪装干扰波形生成方法
CN114265066A (zh) * 2021-12-27 2022-04-01 中国人民解放军63891部队 运动船舶的sar成像模拟方法、装置及电子设备
CN114527436A (zh) * 2022-02-16 2022-05-24 西安电子科技大学 基于最优干扰规划模型的集群协同干扰规划方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
《ITU-R P.2108-1 建议书(09/2021)地物损耗的预测》;国际电信联盟;《https://www.itu.int/dms_pubrec/itu-r/rec/p/R-REC-P.2108-1-202109-I!!PDF-C.pdf》;20220315;全文 *
《一种基于局部特性的红外弱小目标检测方法》;王宇翔 等;《空军雷达学院学报》;20120415;第26卷(第2期);全文 *
《雷达信号传播建模》;Monster0112;《https://blog.csdn.net/weixin_49905873/article/details/116274099》;20210527;全文 *
地面防空预警雷达探测概率仿真模型研究;毛滔等;《空军预警学院学报》;20131215(第06期);全文 *
复杂电磁环境下雷达探测能力分析及仿真;刘己斌等;《系统工程理论与实践》;20080515(第05期);全文 *
复杂电磁环境下雷达探测能力的定量描述;胡进;《航天电子对抗》;20170228(第01期);全文 *

Also Published As

Publication number Publication date
CN114781191A (zh) 2022-07-22

Similar Documents

Publication Publication Date Title
CN114781191B (zh) 一种复杂电磁环境的雷达探测能力仿真方法及装置
CN114781190B (zh) 一种雷达探测能力仿真方法及装置
US7372774B1 (en) System for detecting, tracking, and reconstructing signals in spectrally competitive environments
CN107121670B (zh) 一种基于合成孔径雷达的反无人机防御方法
CN106872831B (zh) 微波载荷卫星内高灵敏接收设备天线带内干扰源确定方法
CN113253213B (zh) 高频地波雷达电离层回波空间分布特性测量方法及系统
CN106886021A (zh) 高分辨率星载sar成像质量提升方法
CN107422312B (zh) 欺骗干扰信号发射功率估计方法
CN107390186A (zh) 压制干扰信号发射功率估计方法
Thomas et al. Response of high frequency radar to meteor backscatter
CN115032605A (zh) 一种基于目标能力画像的雷达目标体系分析方法
CN112036011B (zh) 一种用于水下航行器的甚低频波通信传输分析方法及系统
Lee et al. A Survey and analysis on a troposcatter propagation model based on ITU-R recommendations
CN114374979B (zh) 5GHz频段RLAN台站与雷达系统隔离距离确定方法及装置
CN114244457A (zh) 一种多平台电磁干扰问题的快速预测方法
Penzin et al. Real-time techniques for interpretation of ionospheric backscatter sounding data
Cacciola et al. Impact of transmit interference on receive sensitivity in a bi-static active array system
Lange Performance prediction and improvement of a Bistatic Passive coherent location Radar.
Breton et al. Terrain-Scatter Augmented Vertical Plane Model for Radio Path Loss Estimation in Complex Terrain
JP7300654B2 (ja) 電波遮蔽損失算出方法、電波遮蔽損失算出装置、及び電波遮蔽損失算出プログラム
Riddolls Practical two-dimensional antenna arrays for high-latitude over-the-horizon radar systems
Hampton The impact of evaporative ducting on covert communications
Shuang-xun et al. Study of NLMS algorithm used in the cancellation of navigation interference in the H I observation
Saavedra et al. A new approach to OTH main parameters determination
Xun et al. Laddered Hierarchy Analysis Method for Electromagnetic Compatibility of UAV Based on Relevant Frequency Spectrum Curve

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