基于模糊数学和MATLAB图形用户界面的雷达干扰效果评估方法与流程

基于模糊数学和matlab图形用户界面的雷达干扰效果评估方法
技术领域
1.本发明属于雷达干扰效果评估领域，涉及一种基于模糊数学和matlab图形用户界面的雷达干扰效果评估方法。


背景技术：

2.雷达对抗在现代战争中发挥着重要作用。典型的雷达对抗系统主要包含电子对抗侦察、干扰、指挥控制3个子系统。其中，雷达对抗指挥控制子系统需要根据侦察信息以及干扰效果调整干扰策略，以达到最佳干扰效果。因此，对于干扰效果评估技术进行研究，准确地判断战场环境和精确地评估雷达干扰效果，对战争的结果具有重要影响，具有重要的军事应用价值。
3.雷达干扰效果评估是一个多指标综合评估问题，它需要从多个角度和层次反映雷达干扰效能的水平。现代的雷达对抗技术是以体系对抗为代表的，对抗的双方都是众多装备的集合，特别是信息化的引入，进一步提高了双方综合作战的性能，使它们的相互作用不仅仅是数量的结果叠加，而可能引发许多相关性的新效果，如遮盖/欺骗复合干扰等，在雷达的时、频、空等多维目标检测域形成非常复杂的检测局面。也正是因此，至今没有一个全面的统一的实用准则和方法。因此，对雷达干扰效果评估的研究具有重要的意义。
4.早期的评估因子法，其基本思想是将一些对干扰效果有影响的主要因素通过数学处理组合成评估因子，以此来评估干扰效果。评估因子法利用了对干扰效果有影响的一些重要因素，能够简单地对干扰效果进行评估。但在一般情况下，该方法考虑的因素数量较少，不足以全面反映干扰效果。为了解决干扰效果评估的模糊性问题，研究人员将模糊数学的工具引入到干扰效果评估中，形成了模糊综合评判法。模糊综合评判法可以较好地适应干扰效果评估中影响因素众多的特点，能够将分析得出的可能影响干扰效果的所有因素全部引入到因素集中，使干扰效果的评估更加全面。
5.传统的干扰效果评估方法不具有易操作性，将matlab图形用户界面与雷达干扰效果评估方法结合，提升了使用matlab对雷达干扰效果做出评估的数据可视化效果，能够反映出雷达干扰方案的干扰效果。


技术实现要素：

6.为了解决能够在干扰因素众多的情况下，利用matlab图形用户界面的易操作性，完成雷达干扰效果评估的问题,本发明提出如下技术方案:一种基于模糊数学和matlab图形用户界面的雷达干扰效果评估方法，包括以下步骤：
7.步骤1：建立雷达干扰因素集，其中的干扰因素至少包括干扰时间、干扰空间、干扰频率、干扰功率和干扰样式；
8.步骤2：建立评价集，将一实数集划分为若区间，每一区间对应表示一个雷达干扰效果的评估等级和该等级对应的评估效果；
9.步骤3：创建用于获取用户输入的各项参数的matlab图形用户界面，输入的参数包括干扰因素及对应的干扰效益值，干扰因素的权重值；
10.步骤4：通过各干扰因素的干扰效益值和隶属度函数，对各输入的干扰因素进行单因素评判，并将其写成矩阵形式，而得到模糊关系矩阵；
11.步骤5：通过权重和模糊关系矩阵，计算得到模糊综合评估矩阵；
12.步骤6：对模糊综合评估矩阵去模糊处理，得到雷达干扰效果的评估结果，将评估结果与评价集对比，最终获取雷达干扰效果的评估结果等级并显示在matlab图形用户界面上。
13.进一步的，步骤2中，采用实数集v＝{v|0≤v≤1}来描述评价集，将实数集v划分为若干个区间来表示不同的等级。
14.进一步的，将实数集v划分为六个区间来表示不同的等级，对应的等级依次为“很差、差、一般、较好、好、很好”，v越大，评判等级越高，评判结果越好。
15.进一步的，步骤2中建立评价集，将[0,1]之间的数划分为若干区间，每一区间对应表示一个等级。
[0016]
进一步的，步骤3中，创建的matlab图形用户界面包括三个模块：
[0017]
干扰因素数据输入模块，用于获取目标雷达和己方雷达的因素集中的参数数据，所述雷达干扰因素集至少包括干扰时间、干扰空间、干扰频率、干扰功率和干扰样式，并获取所述雷达干扰因素集中各干扰因素的干扰效益值；
[0018]
权重输入模块，用于获取用户输入的所述干扰因素的权重值，表明所述干扰因素对干扰效果的影响程度；
[0019]
评估结果显示模块，用于根据所述干扰因素的干扰效益值和隶属度函数，计算得到模糊关系矩阵，并结合用户输入的权重获得干扰因素的模糊综合矩阵，接下来进行去模糊处理得到雷达干扰效果的评估结果，将评估结果与评价集作对比，获取雷达干扰效果的评估等级，并将最终的评估效果予以显示在matlab图形用户界面上。
[0020]
进一步的，步骤4中，采用梯形隶属度函数对雷达干扰因素集中各个干扰因素进行单因素评判。
[0021]
进一步的，步骤6中，采用最大隶属度原则对数据进行去模糊处理。
[0022]
有益效果：本发明创建了matlab图形用户界面，增强了各功能模块间的联系，从干扰方出发，首先通过分析雷达干扰要素，建立合理的雷达干扰因素集；并结合每种因素的属性通过预先设定的不同计算方法获得所述雷达干扰因素集中各干扰因素的干扰效益值；确定隶属度函数，对各个因素进行单因素评价得到模糊关系矩阵；获取用户输入的所述干扰因素的权重值；根据所述模糊关系矩阵和权重获得干扰因素的模糊综合矩阵，并根据最大隶属度原则获得最终的雷达干扰效果评估结果。本发明的评估方法可以对不同干扰条件下的干扰效果进行综合评价，反映出雷达干扰方案的干扰效果，提升了使用matlab对雷达干扰效果做出评估的数据可视化能力，使雷达干扰效果评估过程更加具有可操作性。
附图说明
[0023]
图1为本发明的流程图；
[0024]
图2是根据本发明实施例的matlab图形用户界面静态文本参数设置；
[0025]
图3是根据本发明实施例的matlab图形用户界面外观图；
[0026]
图4是根据本发明实施例的matlab图形用户界面测试图。
具体实施方式
[0027]
如图1所示，本发明基于模糊数学和matlab图形用户界面的雷达干扰效果评估方法，包括以下步骤：
[0028]
步骤1：建立因素集，分析影响干扰效果的因素，获取目标雷达和雷达干扰方案的行为参数建立雷达干扰因素集u＝{u1,u2,
···
,u
n
}，表示n个因素组成的有序集。本发明所述雷达干扰因素集设定n＝5，因素集中包括的因素为干扰时间、干扰空间、干扰频率、干扰功率和干扰样式共五个因素。确定matlab图形用户界面的输入和输出；创建matlab图形用户界面。
[0029]
在本发明的实施例中，根据不同干扰因素的干扰原理不同，通过预先设定的不同计算方法获得所述雷达干扰因素集中各干扰因素的干扰效益值，五个干扰因素的效益值计算方法如下：
[0030]
干扰时间效益值：对于干扰时间因素，本发明采用干扰时间效益函数e
t
来表示干扰机对雷达的干扰压制效果。e
t
为在雷达的威胁时间内，干扰机对雷达能有效干扰的时间段对整体干扰效果的影响程度。在本发明的实施例中，设雷达的信号脉宽为[t
r1
,t
r2
]，干扰机对雷达实施有效干扰的时长为[t
j1
,t
j2
]，通过公式(1)计算干扰时间效益值。
[0031][0032]
干扰空间效益值：对于干扰空间因素，本发明采用空域影响效益函数e
h
来计算干扰空间效益值，获取目标雷达的探测覆盖范围ω
r
和干扰机干扰范围ω
j
，用δω
j
表示目标雷达的探测覆盖范围和干扰机干扰范围的交叠区域，即空域有效干扰区域，δω
j
的取值范围是[0：ω
r
]，通过公式(2)计算干扰空间效益值。
[0033][0034]
干扰频率效益值：对于干扰频率因素，通过干扰机的干扰频率与目标雷达的工作频率两者的对准程度来评估干扰机对目标雷达的干扰压制效果，采用频率对准效益函数e
f
来计算干扰频率效益值，设侦察到的雷达工作频率范围为f
r1
‑
f
r2
，干扰机的频率覆盖范围为f
j1
‑
f
j2
，通过公式(3)计算干扰空间效益值。
[0035]
[0036]
干扰功率效益值：对于干扰功率因素，采用功率压制效益函数e
p
来表示干扰机对雷达的干扰压制效果。在雷达接收机端，只有当干信比达到要求的压制系数时，才能对回波信号实施干扰，设获得的目标雷达的信号功率p
rs
，干扰机的干扰信号功率为p
rj
，由雷达方程可知：
[0037][0038]
其中，p
s
和p
j
分别表示目标雷达和干扰机的发射功率，g
t
是雷达主瓣增益，g
a
是干扰机的增益，σ是雷达的散射面积，a是雷达天线的有效接收面积，r是雷达与干扰机间的距离，k
j
为正常工作的干信比，通过公式(5)计算干扰功率效益值。
[0039][0040]
雷达干扰样式效益值：分析干扰样式因素时，用干扰样式与雷达工作体制的匹配程度来表示干扰效果，在干扰样式选择上，需要从实际的效果出发，虽然实际环境中敌方抗干扰方式不能够准确得知，但存在着依靠侦查设备分析出雷达信号频率等先验信息后，采用一定方法对敌雷达的信号样式进行可能性的判断，一般通过大量的实验和专家的经验给出一个介于0和1之间的数来描述所选择的干扰样式和侦察到的雷达工作体制之间的匹配程度，应当理解的是，在本实施例中干扰样式的干扰效益值通过用户输入获得。
[0041]
步骤2：建立评价集，将[0,1]之间的的数划分为若干段，每一段对应地表示一个等级，用实数集描述为v＝{v|0≤v≤1}，根据v的值，将其划分为以下6个评价等级，见表1。
[0042]
表1评价级别的划分
[0043][0044]
步骤3：如图2、3所示，创建matlab图形用户界面，用于获取用户输入的各项参数，并计算得到雷达干扰效果评估结果，创建的matlab图形用户界面包括3个模块：
[0045]
干扰因素数据输入模块，用于获取目标雷达和己方雷达的因素集中的参数数据，所述雷达干扰因素集至少包括干扰时间、干扰空间、干扰频率、干扰功率和干扰样式；并分别通过预先设定的不同的计算方法获得所述雷达干扰因素集中各干扰因素的干扰效益值；
[0046]
权重输入模块，用于获取用户输入的所述干扰因素的权重值，表明所述因素对干扰效果的影响大小；
[0047]
评估结果显示模块，用于根据所述干扰因素的干扰效益值和隶属度函数，计算得到模糊关系矩阵，并结合用户输入的权重获得干扰因素的模糊综合矩阵，采用最大隶属度原则获得雷达干扰效果评估结果，并显示最终的干扰效果评估结果。
[0048]
步骤4：确定隶属度函数，隶属函数的确定有多种方法，在评估中可根据因素对干
扰效果产生影响的原理，依据其数学表达形式，并结合实际情况，选取适当的函数。本发明选用梯形隶属度函数对干扰因素进行评估。隶属度函数表达式如公式(6)：
[0049][0050]
其中f
ij
(a
i
)是第i个因素(用数值a
i
表征)对第j个评价等级的隶属函数；x1，x2，x3，x4是函数中的参数，有0≤x1＜x2＜x3＜x4≤1，通过隶属度函数对各个因素进行单因素评判，将其写成矩阵形式，得到模糊关系矩阵如公式(7)。
[0051]
r＝(r
ij
)
n
×
m
ꢀꢀ
(7)
[0052]
其中0≤i≤n，0≤j≤m。
[0053]
步骤5：确定权重集并计算综合评估结果，通过权重和模糊关系矩阵计算模糊综合评估结果,如公式(8)。
[0054][0055]
采用最大隶属度原则对数据结果进行去模糊处理,得到雷达干扰效果的评估结果。
[0056]
下面结合具体实施例，进一步说明本发明的有益效果。
[0057]
假设被干扰雷达是单频脉冲多普勒雷达，它的工作时间为2μs
‑
4μs，工作频率范围是7.5khz
‑
12.5khz，功率范围是300w
‑
1200w，天线的主瓣增益范围为5
‑
10db，干扰机距雷达的距离10km，取雷达正常工作所需的干信比界限k
j
＝0.5。干扰机的工作时间为1μs
‑
3μs，工作频率范围是(f
j1
,f
j2
)＝5khz
‑
15khz，脉冲功率为1000w，连续波功率为100w，具有拖距、拖速、连续杂波、选择距离等多种干扰样式，散射面积σ＝10m2,干扰机天线增益是10db。
[0058]
首先按照之前说明的步骤，计算各个干扰因素的效益值，并对各个因素进行单因素评判，在本实施例中，以干扰时间因素为例，按照公式(1)计算干扰时间效益值。其中，目标雷达的工作时间范围是2μs
‑
4μs，干扰机的工作时间范围是1μs
‑
3μs，计算得e
t
＝0.5。
[0059]
按照公式(6)的隶属度函数，对各个干扰因素的效益值做单因素评价，在本实施例中，同样以干扰时间因素效益值为例，得到r
1,1≤j≤6
＝[0,0,0.5,1,0.5,0]，将各因素的单因素评价结果写成矩阵形式，得到模糊关系矩阵如公式(9)所示。
[0060][0061]
在本实施例中，输入权重a＝[0.1,0.1,0.2,0.2,0.4]，利用公式(8)计算，得到模糊综合矩阵如式(10)。
[0062][0063]
采用最大隶属度原则对数据结果进行去模糊处理，得到雷达干扰效果的评估结
果，在v
j
＝[0.4 0.075 0.15 0.325 0.45 0.2]中，最大值0.45对应评价集v中的5级(差)。
[0064]
接下来用本发明创建的matlab用户图形界面验证评估结果是否一致。
[0065]
如图4所示，在干扰因素数据输入模块的干扰时间编辑框，输入目标雷达的工作时间和干扰机的干扰时间，在本实施例中，分别输入目标雷达的工作时间2μs
‑
4μs，干扰机的干扰时间1μs
‑
3μs。
[0066]
在干扰因素数据输入模块的干扰空间编辑框，输入目标雷达的探测覆盖范围和干扰机干扰范围，在本实施例中，分别输入目标雷达的探测覆盖范围0
°‑5°
，干扰机干扰范围0
°‑4°
。
[0067]
在干扰因素数据输入模块的干扰频率编辑框，输入侦察到的雷达工作频率范围以及干扰机的频率覆盖范围，在本实施例中，分别输入雷达工作频率范围7.5khz
‑
12.5khz，干扰机的频率覆盖范围5khz
‑
15khz。
[0068]
干扰因素数据输入模块的干扰功率编辑框，输入目标雷达的信号功率和干扰机的干扰信号功率，在本实施例中，输入雷达信号功率和干扰机干扰信号功率的比值1。
[0069]
在干扰因素数据输入模块的干扰样式编辑框，输入用户得到的专家经验性干扰样式效益值，在本实施例中，输入经验性干扰样式效益值0.3。
[0070]
在权重编辑框，输入用户既定权重，在本实施例中，设权重集为a＝[0.1,0.1,0.2,0.2,0.4]。
[0071]
输入以上数据后，点击“计算”按钮，评估结果显示为0.45，属于评价集中的5级(差)，和之前的计算结果一致，表明matlab图形用户界面能够实现雷达效果评估的可操作性，对雷达干扰的干扰效果做出准确评估。
[0072]
以上所述，仅为本发明创造较佳的具体实施方式，但本发明创造的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明创造披露的技术范围内，根据本发明创造的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明创造的保护范围之内。