一种电磁频谱信号能量异常检测方法 |
|||||||
申请号 | CN202311110359.X | 申请日 | 2023-08-30 | 公开(公告)号 | CN117147966B | 公开(公告)日 | 2024-05-07 |
申请人 | 中国人民解放军军事科学院系统工程研究院; | 发明人 | 汪生; 方旖; 刘传文; 张琳; | ||||
摘要 | 本 发明 公开了一种电磁 频谱 信号 能量 异常检测 方法,该方法包括:在执行重要监测任务前,获取电磁频谱信号;对电磁频谱信号进行处理,得到k个典型 频率 的频谱,对k个典型频率的频谱进行处理,得到信号能量 波动 门 限;在执行重要监测任务时,实时计算电磁频谱信号的能量波动方差;将能量波动方差与信号能量波动门限进行频谱能量波动变化对比检测,得到信号能量异常检测结果;在信号能量异常时,对能量波动方差进行处理,得到能量异常程度信息;根据能量异常程度信息进行告警处理。本发明实现了对电磁频谱信号数据的客观分析,以及信号能量异常的自动检测,提高了导致信号变化的重要事件发现效率,增加了目标意图分析的维度。 | ||||||
权利要求 | 1.一种电磁频谱信号能量异常检测方法,其特征在于,所述方法包括: |
||||||
说明书全文 | 一种电磁频谱信号能量异常检测方法技术领域[0001] 本发明涉及电磁频谱空间异常检测技术领域,尤其涉及一种电磁频谱信号能量异常检测方法。 背景技术[0002] 掌握电磁频谱信号活动情况是电磁频谱空间领域的一项重要活动,其途径和方式的优势决定了分析的效率和准确性。目前,电磁频谱信号活动情况的掌握,针对电子目标主要依赖于电子侦察,针对小范围辐射源主要依赖于频谱监测等设备,但随着频谱监测设备能力的提升,也能发现更远的信号,包括电子目标的信号。频谱监测设备的优势是能从时、空、频、能、谱等多维度对电磁频谱信号进行分析和处理,并且能域是频谱信号的重要要素,能够弥补传统电子侦察方式在能域分析处理上的薄弱之处。通过频谱监测设备观测到的电磁频谱信号剧烈变化,最直接的表现就是信号能量大幅度波动,比如:信号出现、消失,信号功率增强、降低等,对这类信号能量大幅度波动统一定义为信号能量异常,这些异常往往预示着重要事件的发生。现有技术对电磁频谱信号能量异常的发现,往往依赖于人的经验,具体体现在以下两方面: [0003] 人工分析电磁频谱信号能量幅值变化剧烈程度判断异常。该类方法主要是利用频谱监测设备测量的频率、带宽、调制方式、方向、信号电平等参数,与历史数据进行比对,并观察信号变化幅度和持续时间,继而人工判断出信号是否发生异常。这种方式判断信号能量异常的正确性取决于人的经验,并需要人工长时值守才能发现频谱能量异常。 [0004] 设置电磁频谱信号能量门限判断超标异常。该类方法是基于频谱监测设备获取的信号频谱,依据经验设置能量幅度门限,或者依据统计信号历史平均变化情况设置门限,超过该门限则表示发生异常。这种方式较前面的人工分析方式有所进步,但仍然需要根据人的经验或者数据平均统计情况设置信号门限,无法保证门限的准确性;并且电磁环境发生变化,或者监测设备发生搬动,均需要重新设置门限,无法自动适应性调整门限。 发明内容[0005] 本发明所要解决的技术问题在于,提供一种电磁频谱信号能量异常检测方法,应用于对电磁频谱信号变化的检测处理,能够基于频谱监测设备获取的电磁频谱信号谱变化的特点,进行电磁频谱信号能量异常的客观自动检测与发现。 [0006] 为了解决上述技术问题,本发明实施例公开了电磁频谱信号能量异常检测方法,所述方法包括: [0007] S1,在执行重要监测任务前,获取电磁频谱信号; [0008] S2,对所述电磁频谱信号进行处理,得到k个典型频率的频谱,k为整数; [0009] S3,对所述k个典型频率的频谱进行处理,得到信号能量波动门限; [0010] S4,在执行重要监测任务时,实时计算电磁频谱信号的能量波动方差; [0011] S5,将所述能量波动方差与所述信号能量波动门限进行频谱能量波动变化对比检测,得到信号能量异常检测结果; [0012] S6,在信号能量异常时,对所述能量波动方差进行处理,得到能量异常程度信息; [0013] S7,根据所述能量异常程度信息进行告警处理。 [0014] 作为一种可选的实施方式,本发明实施例中,所述对所述k个典型频率的频谱进行处理,得到信号能量波动门限,包括: [0015] S31,对所述k个典型频率的频谱进行处理,得到频谱能量方差矩阵; [0016] S32,对所述频谱能量方差矩阵进行映射,得到频谱能量方差图像信息; [0017] S33,对所述频谱能量方差图像信息进行处理,得到图像总的模糊熵; [0018] S34,对所述图像总的模糊熵进行处理,得到信号能量波动门限。 [0019] 作为一种可选的实施方式,本发明实施例中,所述对所述k个典型频率的频谱进行处理,得到频谱能量方差矩阵,包括: [0020] 利用方差计算模型,对所述k个典型频率的频谱进行处理,得到频谱能量方差矩阵; [0021] 所述方差计算模型为: [0022] [0024] 作为一种可选的实施方式,本发明实施例中,所述对所述频谱能量方差矩阵进行映射,得到频谱能量方差图像信息,包括: [0025] 对所述频谱能量方差矩阵进行处理,得到图像灰度值; [0026] 所述图像灰度值为: [0027] [0028] 其中,mi=1,2,…,M,M为采集的电磁频谱信号的帧数, 为第mi帧频谱能量方差矩阵,L=256,*表示相乘,ximi为第i个典型频率第mi帧的图像灰度值; [0029] k个典型频率图像灰度值构成频谱能量方差图像信息; [0030] 所述频谱能量方差图像信息为: [0031] {xi1,xi2,…,xiM} [0032] 其中i=1,2,…,k。 [0033] 作为一种可选的实施方式,本发明实施例中,所述对所述频谱能量方差图像信息进行处理,得到图像总的模糊熵,包括: [0034] S331,对所述频谱能量方差图像信息进行处理,得到模糊集矩阵; [0035] S332,对所述频谱能量方差图像信息进行处理,得到图像直方图; [0036] S333,预设阈值,根据所述阈值对图像直方图进行划分,得到目标图像和背景图像; [0038] S335,对所述像素属于目标图像的条件概率和像素属于背景图像的条件概率进行处理,得到图像总的模糊熵。 [0039] 作为一种可选的实施方式,本发明实施例中,所述对所述图像总的模糊熵进行处理,得到信号能量波动门限,包括: [0040] S341,对所述图像总的模糊熵进行处理,得到最大模糊熵和平均模糊熵; [0041] S342,对所述最大模糊熵和平均模糊熵进行处理,得到启发函数; [0042] S343,对所述启发函数进行处理,得到最优模糊熵; [0043] S344,对所述最优模糊熵进行处理,得到图像分割阈值; [0044] S345,对所述图像分割阈值进行处理,得到信号能量波动门限。 [0045] 作为一种可选的实施方式,本发明实施例中,所述实时计算电磁频谱信号的能量波动方差,包括: [0046] 利用实时滑窗模型,计算电磁频谱信号的能量波动方差; [0047] 所述实时滑窗模型为: [0048] [0049] 其中,ni=t0,t0+1,t0+2,…,N+t0‑1,窗的宽度为N,步长表示为1, 表示第i个典型频率t0窗内电磁频谱信号的能量波动方差,Sni(E)为第i个典型频率t0窗内的频谱,为k个典型频率的频谱均值。 [0050] 作为一种可选的实施方式,本发明实施例中,所述将所述能量波动方差与所述信号能量波动门限进行频谱能量波动变化对比检测,得到信号能量异常检测信息,包括: [0051] S51,记录当前窗口内的能量波动方差和当前窗口1个步长前窗口内的能量波动方差,得到第一能量波动方差和第二能量波动方差; [0052] S52,将所述第一能量波动方差和第二能量波动方差与所述信号能量波动门限进行频谱能量波动变化对比检测,得到信号能量异常检测信息。 [0053] 作为一种可选的实施方式,本发明实施例中,所述将所述第一能量波动方差和第二能量波动方差与所述信号能量波动门限进行频谱能量波动变化对比检测,得到信号能量异常检测结果,包括: [0054] S521,如果所述信号能量波动门限大于等于第一能量波动方差,且小于第二能量波动方差,则当前时间窗口信号能量大幅度增加,信号能量异常检测结果为信号能量异常; [0055] S522,如果所述信号能量波动门限大于等于第二能量波动方差,且小于第一能量波动方差,则当前时间窗口信号能量大幅度下降,信号能量异常检测结果为信号能量异常。 [0056] 作为一种可选的实施方式,本发明实施例中,所述在信号能量异常时,对所述能量波动方差进行处理,得到能量异常程度信息,包括: [0057] 对所述能量波动方差进行处理,得到第一能量波动方差和第二能量波动方差; [0058] 利用能量异常程度计算模型,对所述第一能量波动方差和第二能量波动方差进行处理,得到能量异常程度信息; [0059] 所述能量异常程度计算模型为: [0060] [0061] 其中,E1为能量异常程度信息, 为第二能量波动方差, 为第一能量波动方差,i=1,2,…,k,ni=t0,t0+1,t0+2,…,N+t0‑1,窗的宽度为N,t0为窗的位置,k为典型频率的个数。 [0062] 与现有技术相比,本发明实施例具有以下有益效果: [0063] (1)本发明提出了一种电磁频谱信号能量异常检测方法,实现了对电磁频谱信号数据的客观分析,以及信号能量异常的自动检测,提高了导致信号变化的重要事件发现效率; [0064] (2)为目标异常行为分析从能域维度提供了新的佐证手段,增加了目标意图分析的维度; [0066] 为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。 [0067] 图1是本发明实施例公开的一种电磁频谱信号能量异常检测方法的流程示意图; [0068] 图2是本发明实施例公开的另一种电磁频谱信号能量异常检测方法的流程示意图。 具体实施方式[0069] 为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。 [0070] 本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别不同对象,而不是用于描述特定顺序。此外,术语“包括”和“具有”以及它们任何变形,意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、装置、产品或设备没有限定于已列出的步骤或单元,而是可选地还包括没有列出的步骤或单元,或可选地还包括对于这些过程、方法、产品或设备固有的其他步骤或单元。 [0071] 在本文中提及“实施例”意味着,结合实施例描述的特定特征、结构或特性可以包含在本发明的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例,也不是与其它实施例互斥的独立的或备选的实施例。本领域技术人员显式地和隐式地理解的是,本文所描述的实施例可以与其它实施例相结合。 [0072] 本发明公开了一种电磁频谱信号能量异常检测方法,在执行重要监测任务前,获取电磁频谱信号;对电磁频谱信号进行处理,得到k个典型频率的频谱,对k个典型频率的频谱进行处理,得到信号能量波动门限;在执行重要监测任务时,实时计算电磁频谱信号的能量波动方差;将能量波动方差与信号能量波动门限进行频谱能量波动变化对比检测,得到信号能量异常检测结果;在信号能量异常时,对能量波动方差进行处理,得到能量异常程度信息;根据能量异常程度信息进行告警处理。本发明实现了对电磁频谱信号数据的客观分析,以及信号能量异常的自动检测,提高了导致信号变化的重要事件发现效率,增加了目标意图分析的维度。以下分别进行详细说明。 [0073] 实施例一 [0074] 请参阅图1,图1是本发明实施例公开的一种电磁频谱信号能量异常检测方法的流程示意图。其中,图1所描述的电磁频谱信号能量异常检测方法应用于频谱态势系统,完成实时电磁频谱信号能量异常自动检测,为电磁频谱信号事件关联分析、目标用频行为和意图分析提供有力支持,本发明实施例不做限定。如图1所示,该出电磁频谱信号能量异常检测方法可以包括以下操作: [0075] S1,在执行重要监测任务前,获取电磁频谱信号; [0076] S2,对所述电磁频谱信号进行处理,得到k个典型频率的频谱,k为整数; [0077] S3,对所述k个典型频率的频谱进行处理,得到信号能量波动门限; [0078] S4,在执行重要监测任务时,实时计算电磁频谱信号的能量波动方差; [0079] S5,将所述能量波动方差与所述信号能量波动门限进行频谱能量波动变化对比检测,得到信号能量异常检测结果; [0080] S6,在信号能量异常时,对所述能量波动方差进行处理,得到能量异常程度信息; [0081] S7,根据所述能量异常程度信息进行告警处理。 [0082] 可选的,所述对所述k个典型频率的频谱进行处理,得到信号能量波动门限,包括: [0083] S31,对所述k个典型频率的频谱进行处理,得到频谱能量方差矩阵; [0084] S32,对所述频谱能量方差矩阵进行映射,得到频谱能量方差图像信息; [0085] S33,对所述频谱能量方差图像信息进行处理,得到图像总的模糊熵; [0086] S34,对所述图像总的模糊熵进行处理,得到信号能量波动门限。 [0087] 可选的,所述对所述k个典型频率的频谱进行处理,得到频谱能量方差矩阵,包括: [0088] 利用方差计算模型,对所述k个典型频率的频谱进行处理,得到频谱能量方差矩阵; [0089] 所述方差计算模型为: [0090] [0091] 其中,i=1,2,…,k,M为采集的电磁频谱信号的帧数, 为第mi帧频谱能量方差矩阵,Smi(E)为第i个典型频率的第mi帧频谱, 为k个典型频率的频谱均值。 [0092] 可选的,所述对所述频谱能量方差矩阵进行映射,得到频谱能量方差图像信息,包括: [0093] 对所述频谱能量方差矩阵进行处理,得到图像灰度值; [0094] 所述图像灰度值为: [0095] [0096] 其中,mi=1,2,…,M,M为采集的电磁频谱信号的帧数, 为第mi帧频谱能量方差矩阵,L=256,*表示相乘,ximi为第i个典型频率第mi帧的图像灰度值; [0097] k个典型频率图像灰度值构成频谱能量方差图像信息; [0098] 所述频谱能量方差图像信息为: [0099] {xi1,xi2,…,xiM} [0100] 其中i=1,2,…,k。 [0101] 可选的,所述对所述频谱能量方差图像信息进行处理,得到图像总的模糊熵,包括: [0102] S331,对所述频谱能量方差图像信息进行处理,得到模糊集矩阵; [0103] S332,对所述频谱能量方差图像信息进行处理,得到图像直方图; [0104] S333,预设阈值,根据所述阈值对图像直方图进行划分,得到目标图像和背景图像; [0105] 可以,以直方图双峰之间的谷底处灰度值作为阈值,对图像直方图进行划分,得到目标图像和背景图像; [0106] S334,根据所述模糊集矩阵,对所述目标图像和背景图像进行处理,得到像素属于目标图像的条件概率和像素属于背景图像的条件概率; [0107] S335,对所述像素属于目标图像的条件概率和像素属于背景图像的条件概率进行处理,得到图像总的模糊熵。 [0109] 步骤一:飞蛾扑火算法初始化。 [0110] 设定飞蛾种群搜索规模,初始火焰数量,最大迭代次数,当前迭代次数,对数螺旋常数,待求变量个数的初始值。 [0111] 步骤二:生成初始种群。 [0112] 脉冲耦合神经网络模型图像分割需要确定的三个参数为连接系数、幅度系数、衰减系数,连接系数的取值范围为[0.0001,1],衰减系数的取值范围为[0.1,1],幅度系数的取值范围设定为被分割图像灰度值的上下限。设定种群规模为20。 [0113] 随机生成初始化种群 其中M为初始时飞蛾数目,设M=20。 [0114] 每个个体的空间位置用向量为: [0115] [0116] 其中,i=1,2,…,M,t为迭代的代数,d为待求变量的个数,即待求的问题的解的维度,d=3,三个变量分别代表连接系数、幅度系数、衰减系数。 [0117] 步骤三:初始化火焰空间位置。 [0118] 在搜索空间中随机生成飞蛾的位置,将得到的每只飞蛾对应的参数代入脉冲耦合神经网络模型对图像进行分割。对飞蛾所处的空间位置对应的适应度值排序,作为第1代火焰的空间位置。 [0119] 步骤四:更新飞蛾位置。 [0120] 飞蛾空间位置的更新依赖于其所环绕的火焰,更新当前代飞蛾的空间位置: [0121] Mi=S(Mi,Fj) [0122] 其中,S代表螺旋函数,Mi代表第i只飞蛾,Fj代表第j个火焰。 [0123] 步骤五:分别计算每个个体的适应度值。 [0124] 将得到的每只飞蛾对应的参数带入PCNN模型对图像进行分割,求出种群所有飞蛾相应的交叉熵,把它记为适应度值。交叉熵函数定义如下: [0125] [0126] [0127] [0128] 其中,D是交叉熵值,z表示的是图像像素的灰度值,Z是灰度值的最大值,S表示分割图的区域,S1表示目标区域,S2表示背景区域,μ1和μ2分别表示目标区域和背景区域的像素灰度值的平均值。th表示 取得最小值时的像素的最小灰度,h(z)表示图像被分割后像素灰度值为z的个数。 [0129] 步骤六:更新火焰位置。 [0130] 将更新后的飞蛾位置对应的代入参数脉冲耦合神经网络模型对图像进行分割,得到适应度值后,按照适应度值从大到小排序排列,然后与上一代火焰的适应度值进行比较,选取更优的适应度值更新火焰的位置。 [0131] 步骤七:自适应减少火焰的数量。 [0132] 把当前火焰的数量依据种群最大个体数数和最大迭代次数按照自适应减少规则逐渐减少,火焰数量的减少既可以加快飞蛾对最优解的搜索,也可以让因为火焰数目减少未与火焰匹配的飞蛾参与全局范围的最优解的搜索。火焰数量自适应减少的公式如下: [0133] [0134] 其中,N表示种群中飞蛾的数量,T表示最大迭代次数,w代表当前的迭代次数,round表示四舍五入。 [0135] 步骤八:判断算法是否达到初始时设定的最大迭代次数,若达到最大迭代次数,则结束算法,输出分割图像;若未达到最大迭代次数,则返回步骤四,进入下一次迭代。 [0136] 可选的,所述对所述图像总的模糊熵进行处理,得到信号能量波动门限,包括: [0137] S341,对所述图像总的模糊熵进行处理,得到最大模糊熵和平均模糊熵; [0138] S342,对所述最大模糊熵和平均模糊熵进行处理,得到启发函数; [0139] S343,对所述启发函数进行处理,得到最优模糊熵; [0140] S344,对所述最优模糊熵进行处理,得到图像分割阈值; [0141] S345,对所述图像分割阈值进行处理,得到信号能量波动门限。 [0142] 可选的,所述实时计算电磁频谱信号的能量波动方差,包括: [0143] 利用实时滑窗模型,计算电磁频谱信号的能量波动方差; [0144] 所述实时滑窗模型为: [0145] [0146] 其中,ni=t0,t0+1,t0+2,…,N+t0‑1,窗的宽度为N,步长表示为1, 表示第i个典型频率t0窗内电磁频谱信号的能量波动方差,Sni(E)为第i个典型频率t0窗内的频谱,为k个典型频率的频谱均值。 [0147] 可选的,所述将所述能量波动方差与所述信号能量波动门限进行频谱能量波动变化对比检测,得到信号能量异常检测信息,包括: [0148] S51,记录当前窗口内的能量波动方差和当前窗口1个步长前窗口内的能量波动方差,得到第一能量波动方差和第二能量波动方差; [0149] S52,将所述第一能量波动方差和第二能量波动方差与所述信号能量波动门限进行频谱能量波动变化对比检测,得到信号能量异常检测信息。 [0150] 可选的,所述将所述第一能量波动方差和第二能量波动方差与所述信号能量波动门限进行频谱能量波动变化对比检测,得到信号能量异常检测结果,包括: [0151] S521,如果所述信号能量波动门限大于等于第一能量波动方差,且小于第二能量波动方差,则当前时间窗口信号能量大幅度增加,信号能量异常检测结果为信号能量异常; [0152] S522,如果所述信号能量波动门限大于等于第二能量波动方差,且小于第一能量波动方差,则当前时间窗口信号能量大幅度下降,信号能量异常检测结果为信号能量异常。 [0153] 可选的,所述在信号能量异常时,对所述能量波动方差进行处理,得到能量异常程度信息,包括: [0154] 对所述能量波动方差进行处理,得到第一能量波动方差和第二能量波动方差; [0155] 利用能量异常程度计算模型,对所述第一能量波动方差和第二能量波动方差进行处理,得到能量异常程度信息; [0156] 所述能量异常程度计算模型为: [0157] [0158] 其中,E1为能量异常程度信息, 为第二能量波动方差, 为第一能量波动方差,i=1,2,…,k,ni=t0,t0+1,t0+2,…,N+t0‑1,窗的宽度为N,t0为窗的位置,k为典型频率的个数。 [0159] 实施例二 [0160] 请参阅图2,图2是本发明实施例公开的另一种电磁频谱信号能量异常检测方法的流程示意图。其中,图2所描述的电磁频谱信号能量异常检测方法应用于频谱态势系统,完成实时电磁频谱信号能量异常自动检测,为电磁频谱信号事件关联分析、目标用频行为和意图分析提供有力支持,本发明实施例不做限定。如图2所示,该出电磁频谱信号能量异常检测方法可以包括以下操作(以f1为例): [0161] 步骤1:计算各频率信号能量波动门限 [0162] 根据监测设备周围复杂的电磁环境,针对不同频率选取合理的能量波动门限是信号能量异常自动检测的先决条件。具体方法是在一定时间内持续监测电磁频谱信号并记录各典型频率(f1,f2,…,fk)的频谱,k为典型频率个数,并利用蚁群算法从中推断计算出的方差 作为信号能量波动门限。 [0163] 步骤1.1:计算频谱能量方差矩阵并映射成图像 [0164] 以f1为例,先采集M帧的频谱信号Sm1(E),计算出每帧频谱信号方差 得到方差矩阵: [0165] [0166] E为能量, 为信号强度均值。 [0167] 接着将方差矩阵映射到L为256级的图像灰度空间: [0168] {x11,x12,…,x1m} (2)式中, [0169] 步骤1.2:构建模糊集及隶属度函数 [0170] 根据形成的图像,构建出模糊集矩阵,公式如下所示: [0171] X=[μ1m1(x1m1)]1×M,m1=1,2,3,…,M (3)[0172] 其中,μ1m1(x1m1)表示矩阵中第(1,m1)个像素的灰度模糊集的隶属度为μ1m1,或图像第(1,m1)个像素x1m1具有某种特征的程度为μ1m1,且μ1m1∈[0,1]。 [0173] 模糊隶属度函数采用S函数和Z函数(Z=1‑S),如下所示: [0174] [0175] 式中,μa(k)和μb(k)分别表示灰度值为k的像素属于目标和背景的模糊隶属度函数值,a,b,c是隶属度函数的可变参量,μ0(k)为隶属度函数。 [0176] 设H={p0,p1,…,pL‑1}为图像的直方图,其中 nk为图像X中灰度值为k的像素个数,M×N表示灰度图像矩阵的大小。设阈值T将图像分为目标Ea和背景Eb两类,则其概率分布分别为: [0177] pa=p(Ea) (4) [0178] pb=p(Eb) (5) [0179] 用pa/k和pb/k分别表示灰度值为k的像素属于目标和背景的条件概率,则目标和背景的概率分布为: [0180] pka=pk×pa/k (6) [0181] pkb=pk×pb/k (7) [0182] 将灰度值为k的pa/k和pb/k与其属于目标和背景的隶属函数值相关联,则有: [0183] pa/k=μa(k) (8) [0184] pb/7=μb(k) (9) [0185] [0186] [0187] 因此,各类的类内模糊熵为: [0188] [0189] [0190] 图像总的模糊熵为: [0191] H(a,b,c)=Ha+Hb [0192] 可见,图像总的模糊熵是参数a,b,c的函数,根据最大熵准则确定参数a,b,c的值,在最佳阀值T处,目标和背景的模糊隶属度函数值相等,即应满足如下条件: [0193] μa(T)=μb(T)=0.5 (15) [0194] 由此可见,寻找最佳分割阀值相当于寻找a,b,c的最优组合。 [0195] 步骤1.3:蚁群算法求解方差门限 [0196] 设置蚁群目标函数L(H)如下所示: [0197] L(H)=H(a,b,c) (16) [0198] 令Lmax(H)和Lavg(H)分别表示最大模糊熵和平均模糊熵,则启发函数η可表示为: [0199] [0200] 经过迭代求解获取(a,b,c)的最优组合值,并利用该结果推导出图像分割阀值x10,即T,接着根据式(18)求得频率为f1时能量谱的方差阀值 [0201] [0202] 步骤2:频谱能量波动门限异常检测 [0203] 在对重点频率信号的频谱能量进行实时监测的过程中,通过滑动窗口的方式计算积累形成的频谱能量波动方差,并与对应信号能量波动门限比较,实现即时电磁频谱信号能量异常检测,其原理框图如图2所示。 [0204] 步骤2.1:实时滑窗计算信号能量波动方差 [0205] 以典型频率f1为例,采用滑动窗口的形式对收集的频谱信号能量数据Sn1(E)进行分块,窗宽表示为N,步长表示为t,并实时计算t0窗内数据的方差: [0206] [0207] 步骤2.2:频谱能量波动变化对比检测 [0208] 记录当前窗口内数据的方差为 1个步长前窗口内数据的方差为 将两个方差与方差门限进行大小对比,如满足以下任意一个公式,并保持到下一个窗口步长,则判断当前时间窗口里发生了能量异常变化: [0209] σn12’≤σ012<σn12 (20) [0210] σn12<σ012≤σn12’ (21) [0211] 其中 的计算公式如下所示: [0212] [0213] 满足式(20)说明在当前时间窗口出现信号能量大幅度增加,而满足式(21)说明在当前时间窗口出现信号能量大幅度下降。 [0214] 步骤3:能量异常严重程度分析与告警 [0215] 当发现信号能量异常后,通过计算偏离正常门限的程度E1划分异常严重程度,并向频谱监测人员提示告警,为辐射源信号的用频行为和意图分析提供支持。 [0216] 为了保证告警等级划分方法的规范性,避免重点频率间特有属性差异对等级划分的影响,衡量偏离正常门限程度E1的计算方法采用与当前频率对应波动门限的比值形式,如下所示: [0217] [0218] 其中,E1的幅值范围为[0,∞),而能量异常变化的严重程度和E1的值具有正相关性,具体异常严重程度等级划分方式见表1。 [0219] 表1异常严重程度划分表 [0220] [0221] 以上所描述的装置实施例仅是示意性的,其中作为分离部件说明的模块可以是或者也可以不是物理上分开的,作为模块显示的部件可以是或者也可以不是物理模块,即可以位于一个地方,或者也可以分布到多个网络模块上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。 [0222] 通过以上的实施例的具体描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,存储介质包括只读存储器(Read‑Only Memory,ROM)、随机存储器(Random Access Memory,RAM)、可编程只读存储器(Programmable Read‑only Memory,PROM)、可擦除可编程只读存储器(Erasable Programmable Read Only Memory,EPROM)、一次可编程只读存储器(One‑time Programmable Read‑Only Memory,OTPROM)、电子抹除式可复写只读存储器(Electrically‑Erasable Programmable Read‑Only Memory,EEPROM)、只读光盘(Compact Disc Read‑Only Memory,CD‑ROM)或其他光盘存储器、磁盘存储器、磁带存储器、或者能够用于携带或存储数据的计算机可读的任何其他介质。 [0223] 最后应说明的是:本发明实施例公开的一种电磁频谱信号能量异常检测方法所揭露的仅为本发明较佳实施例而已,仅用于说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解;其依然可以对前述各项实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或替换,并不使相应的技术方案的本质脱离本发明各项实施例技术方案的精神和范围。 |