一种基于全频带谱分析的雷达压制式干扰检测方法技术

本发明专利技术属于雷达信号处理技术领域，公开了一种基于全频带谱分析的雷达压制式干扰检测方法，包括：得到雷达一个扫描周期内L个脉冲分别对应的频谱累加值；根据雷达一个扫描周期内L个脉冲分别对应的频谱累加值计算每个脉冲的频谱累加值的中间值、频谱累加值的平均值，得到该脉冲对应的加权中值滤波后的值，组成雷达一个扫描周期内L个脉冲的加权中值滤波结果；确定雷达一个扫描周期内L个脉冲组成的干扰区域集合，干扰区域集合包含J个干扰区域；获取J个干扰区域内包含的脉冲分别对应的方位码，确定J个干扰区域内包含的脉冲分别对应的干扰方位、干扰强度、干扰类型以及干扰频点，能够实现雷达工作频带的频谱功率分析，识别率高，易于工程实现，并提高了雷达的抗干扰能力。

【技术实现步骤摘要】

本专利技术涉及雷达信号处理
，尤其涉及一种基于全频带谱分析的雷达压制式干扰检测方法。
技术介绍
现代战争中电磁环境日益复杂，雷达抗干扰面临越来越多的挑战，压制式干扰是对雷达威胁较大的一类干扰；压制式干扰是通过释放强噪声干扰，其从天线旁瓣进入雷达主通道，与目标回波叠加在一起，使雷达接收机过载、饱和，严重影响雷达的正常检测性能。雷达抗干扰能力已成为雷达系统设计的一项核心战术指标，干扰检测性能也已成为雷达抗干扰能力的重要组成部分。良好的干扰检测性能使得雷达能够根据干扰检测结果及时自动选取最优的抗干扰措施，进而使得雷达的整体抗干扰性能得到充分发挥。雷达干扰检测功能实现的方法很多，但多数方法较复杂，不利于工程实现，具有很大的应用局限性；对于传统机械扫描雷达来说，受到自身体制局限性，通常采用在主通道外增加一路辅助通道进行干扰检测，用此方法检测得到的信息往往存在虚警，不够准确；此外，由于主通道和辅助通道工作在当前雷达工作频点下，所以只有干扰处在雷达当前工作频点时才能够被有效检测到，大大影响了干扰检测功能的完备性。
技术实现思路
针对以上现有技术存在的不足，本专利技术的目的在于提出一种基于全频带谱分析的雷达压制式干扰检测方法，该方法能够实现雷达所有工作频带的频谱功率分析，方法简单，实用性强，识别率高，易于工程实现，并提高了雷达的抗干扰能力。为达到上述目的，本专利技术的实施例采用如下技术方案予以实现。一种基于全频带谱分析的雷达压制式干扰检测方法，所述方法包括如下步骤：步骤1，获取雷达发射第l个脉冲后的回波射频信号，对所述回波射频信号依次进行模数转换和正交采样，得到第l个脉冲覆盖雷达工作频带的基带数据；从所述基带数据中均匀截取M段子基带数据，对每段子基带数据依次进行时域加窗和N点快速傅里叶变换，得到该段子基带数据的N个频点数据，从而得到第l个脉冲M段子基带数据的M×N个频点数据；×表示乘号，其中，l初值为1，l＝1，...，L，L为雷达一个扫描周期内包含的脉冲个数，且为自然数，M≥1，N≥1；步骤2，分别计算第l个脉冲M段子基带数据的M×N个频点数据的模值平方，将M段子基带数据中每段子基带数据对应的第n个频点的模值平方进行累加，得到第n个频点对应的频谱值，令n依次取1，…，N，从而得到第l个脉冲N个频点分别对应的频谱值；将所述N个频点分别对应的频谱值进行累加，得到第l个脉冲的频谱累加值；M和N为自然数步骤3，令l加1，并依次重复步骤1至步骤2，从而得到雷达一个扫描周期内L个脉冲分别对应的频谱累加值；步骤4，获取加权中值滤波线性窗，并根据雷达一个扫描周期内L个脉冲分别对应的频谱累加值计算每个脉冲的频谱累加值的中间值、频谱累加值的平均值，然后根据每个脉冲的频谱累加值的中间值、频谱累加值的平均值得到该脉冲对应的加权中值滤波后的值，从而得到L个脉冲对应的加权中值滤波后的值，并组成雷达一个扫描周期内L个脉冲的加权中值滤波结果；步骤5，根据所述雷达一个扫描周期内L个脉冲的加权中值滤波结果，确定雷达一个扫描周期内的干扰区域集合，所述干扰区域集合包含J个干扰区域，J为自然数；步骤6，获取J个干扰区域内包含的脉冲分别对应的方位码，根据J个干扰区域内包含的脉冲分别对应的方位码、J个干扰区域内包含的脉冲分别对应的加权中值滤波后的值，确定J个干扰区域内包含的脉冲分别对应的干扰方位、干扰强度、干扰类型以及干扰频点。本专利技术的有益效果：第一，本专利技术方法结构清晰，易于硬件实现，实用性高；第二，本专利技术方法无需通过切换雷达工作频点，就能够实现雷达工作频带的全频带干扰检测，应用范围大大增强；第三，中值滤波方法的使用，有利于减少脉冲噪声的带来的检测虚警，自适应门限的应用，提高应用的灵活性，提高系统适应环境变化的能力；第四，本专利技术参数可变，能够根据使用过程中的外界环境变化和应用需求进行灵活调节，有效提高方法使用的灵活性、复用性和可移植性。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术一种基于全频带谱分析的雷达压制式干扰检测方法流程示意图；图2a为切比雪夫窗时域波形示意图；图2b为切比雪夫窗频域波形示意图；图3为加切比雪夫窗与不加切比雪夫窗分别进行快速傅里叶变换后得到的频谱示意图；图4为加权中值滤波算法的流程示意图；图5为利用自适应门限划分干扰区域流程示意图；图6a为中值滤波前频谱累加和的曲线示意图；图6b为中值滤波后频谱累加和的曲线示意图；图7a为方位融合后的干扰标志位示意图；图7b为更新后的自适应门限示意图；图8为某脉冲32频点功率分布示意图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。本专利技术实施例提供一种基于全频带谱分析的雷达压制式干扰检测方法，参照图1，所述方法包括如下步骤：步骤1，获取雷达发射第l个脉冲后的回波射频信号，对所述回波射频信号依次进行模数转换和正交采样，得到第l个脉冲覆盖雷达工作频带的基带数据；从所述基带数据中均匀截取M段子基带数据，对每段子基带数据依次进行时域加窗和N点快速傅里叶变换，得到该段子基带数据的N个频点数据，从而得到第l个脉冲M段子基带数据的M×N个频点数据；×表示乘号，其中，l初值为1，l＝1，...，L，L为雷达一个扫描周期内包含的脉冲个数，且为自然数，M≥1，N≥1。对所述第l个脉冲的基带数据中截取的M段数据依次进行时域加窗和快速傅里叶变换，依次得到加窗后第l个脉冲的M段数据yl1，…，ylM和第l个脉冲的M段N点频域结果Yll，…，YlM；其中，M≥1，N≥1，每段数据包含N个频点，则第l个脉冲的基带数据中截取的M段数据包含M×N个频点。步骤1的子步骤为：(1a)i∈{1，…M本文档来自技高网...

【技术保护点】
一种基于全频带谱分析的雷达压制式干扰检测方法，其特征在于，所述方法包括如下步骤：步骤1，获取雷达发射第l个脉冲后的回波射频信号，对所述回波射频信号依次进行模数转换和正交采样，得到第l个脉冲覆盖雷达工作频带的基带数据；从所述基带数据中均匀截取M段子基带数据，对每段子基带数据依次进行时域加窗和N点快速傅里叶变换，得到该段子基带数据的N个频点数据，从而得到第l个脉冲M段子基带数据的M×N个频点数据；×表示乘号，其中，l初值为1，l＝1，...，L，L为雷达一个扫描周期内包含的脉冲个数，且为自然数，M≥1，N≥1；步骤2，分别计算第l个脉冲M段子基带数据的M×N个频点数据的模值平方，将M段子基带数据中每段子基带数据对应的第n个频点的模值平方进行累加，得到第n个频点对应的频谱值，令n依次取1，…，N，从而得到第l个脉冲N个频点分别对应的频谱值；将所述N个频点分别对应的频谱值进行累加，得到第l个脉冲的频谱累加值；步骤3，令l加1，并依次重复步骤1至步骤2，从而得到雷达一个扫描周期内L个脉冲分别对应的频谱累加值；步骤4，获取加权中值滤波线性窗，并根据雷达一个扫描周期内L个脉冲分别对应的频谱累加值计算每个脉冲的频谱累加值的中间值、频谱累加值的平均值，然后根据每个脉冲的频谱累加值的中间值、频谱累加值的平均值得到该脉冲对应的加权中值滤波后的值，从而得到L个脉冲对应的加权中值滤波后的值，并组成雷达一个扫描周期内L个脉冲的加权中值滤波结果；步骤5，根据所述雷达一个扫描周期内L个脉冲的加权中值滤波结果，确定雷达一个扫描周期内的干扰区域集合，所述干扰区域集合包含J个干扰区域，J为自然数；步骤6，获取J个干扰区域内包含的脉冲分别对应的方位码，根据J个干扰区域内包含的脉冲分别对应的方位码、J个干扰区域内包含的脉冲分别对应的加权中值滤波后的值，确定J个干扰区域内包含的脉冲分别对应的干扰方位、干扰强度、干扰类型以及干扰频点。...

【技术特征摘要】
1.一种基于全频带谱分析的雷达压制式干扰检测方法，其特征在于，所述方法包括如下步骤：步骤1，获取雷达发射第l个脉冲后的回波射频信号，对所述回波射频信号依次进行模数转换和正交采样，得到第l个脉冲覆盖雷达工作频带的基带数据；从所述基带数据中均匀截取M段子基带数据，对每段子基带数据依次进行时域加窗和N点快速傅里叶变换，得到该段子基带数据的N个频点数据，从而得到第l个脉冲M段子基带数据的M×N个频点数据；×表示乘号，其中，l初值为1，l＝1，...，L，L为雷达一个扫描周期内包含的脉冲个数，且为自然数，M≥1，N≥1；步骤2，分别计算第l个脉冲M段子基带数据的M×N个频点数据的模值平方，将M段子基带数据中每段子基带数据对应的第n个频点的模值平方进行累加，得到第n个频点对应的频谱值，令n依次取1，…，N，从而得到第l个脉冲N个频点分别对应的频谱值；将所述N个频点分别对应的频谱值进行累加，得到第l个脉冲的频谱累加值；步骤3，令l加1，并依次重复步骤1至步骤2，从而得到雷达一个扫描周期内L个脉冲分别对应的频谱累加值；步骤4，获取加权中值滤波线性窗，并根据...

【专利技术属性】
技术研发人员：苏涛，于祥龙，邱瑾，徐丹，许士杰，贾树雄，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61