一种区域电磁环境的动态特性表征方法
【技术领域】
[0001] 本发明涉及一种动态特性表征方法，特别是涉及一种区域电磁环境的动态特性表 征方法。
【背景技术】
[0002] 区域电磁环境中各种电磁干扰辐射信号种类繁多、调制体制各异，且在传播空间 的叠加，信号参数的动态变化，传播路径及其相应影响的不确定性，表明电磁信号在时间、 能量上相互干扰发生异常变化造成电磁环境动态性异常激烈。对区域中的接收设备造成严 重危害。
[0003] 文献"AnalysisofComplexityinBattlefieldElectromagnetic Environment4thIEEE.ConferenceoIndustrialElectronicsandApplications Xi'an,China:IEEE, 2009"公开了一种对区域战场电磁环境存在的异动信号特征分析方法。 该方法利用基于信号密度分布的方法对区域电磁环境复杂度进行估计与表征，为战场电磁 环境分析奠定基础。但是，该方法只是估计电磁环境复杂度特征，而电磁环境具有很多特 征，比如动态性、涌现性、非线性等。因此，难以表征真实区域电磁环境特性。
【发明内容】
[0004] 为了克服现有区域电磁环境的动态特性表征方法实用性差的不足，本发明提供一 种区域电磁环境的动态特性表征方法。该方法首先基于信号基础理论对区域接收点处电磁 信号集合的进行建模；其次，利用数学形态学滤波原理对接收的电磁信号集合进行交替滤 波，滤除电磁信号中的脉冲噪声和白噪声；最后，运用分形理论求解电磁环境接收点集合的 关联维数。对区域电磁环境动态性量化表征。本发明根据区域电磁信号表示和数学形态学 理论，不仅对区域L个接收点N个辐射源产生的电磁信号进行了建模，并且滤除了环境中突 发噪声，实现了对接收信号的动态性表征。
[0005] 本发明解决其技术问题所采用的技术方案是：一种区域电磁环境的动态特性表征 方法，其特点是采用以下步骤：
[0006] 步骤一、将N个辐射源在某空间区域形成的电磁信号环境表示为L个点处电磁信 号的集合：
[0007] 取⑴上⑴，···土⑴} (1)其中，ΕΛ接收点i处合成电磁信号，
[00081
[0009] 式中，Ajt)为电场的幅值，匕⑴为频率，$·(?)为相位，&为场强方向或极化方向， Fi为福射天线方向图，H为福射源i的j径信号的传递函数，Σ表不传输信号多径合成，Λ仁(t)为辐射源i相对接收点运动速度为Vl (t)时产生的频移。
[0010] 步骤二、基于数学形态学的电磁环境信号处理。
[0011] 通过形态学混合开、闭运算交替滤波，滤除电磁信号中的脉冲噪声和白噪声。具体 公式如下：
[0012] (3):
[0013] <4) \ f- ./. f-
[0014] 式中，Ejt)是接收点一维离散电磁信号，取值范围为Ejt) = {0，1，2,… ，N-1} ;g(m)是结构元素，取值范围为G= {0, 1，2，"'M-l} ; ?表示膨胀算 子，表示为
?表示腐蚀算子，表示为 mζζ\j
[0015] 步骤三、基于分形理论求解电磁环境关联维数D(m)。
[0016]
[0017]
[0018]
[0019] 式中，Η是Heaviside函数，当r>0，H(r)为 1，r< 0 是H(r)为 0。 {Ei，E2，…，Ε1; ···}为区域形成的电磁合成信号的集合，重构延时后得到新集合记为：
τ为延迟时间，τ=kΛt，Λt是时 间序列的采样时间间隔，N^N-On-l)τ为的总点数。
[0020] 本发明的有益效果是：该方法首先基于信号基础理论对区域接收点处电磁信号集 合的进行建模；其次，利用数学形态学滤波原理对接收的电磁信号集合进行交替滤波，滤除 电磁信号中的脉冲噪声和白噪声；最后，运用分形理论求解电磁环境接收点集合的关联维 数。对区域电磁环境动态性量化表征。本发明根据区域电磁信号表示和数学形态学理论， 不仅对区域L个接收点Ν个辐射源产生的电磁信号进行了建模，并且滤除了环境中突发噪 声，实现了对接收信号的动态性表征。
[0021] 下面结合附图和【具体实施方式】对本发明作详细说明。
【附图说明】
[0022] 图1是本发明方法中QPSK信号关联维数D= 1. 7的曲线。
[0023] 图2是本发明方法中QPSK和16QAM合成关联维数D= 1. 8的曲线。
[0024] 图3是本发明方法中DPSK和ASK合成关联维数D= 1. 6的曲线。
[0025] 图4是本发明方法中DPSKASK16QAMQPSK合成关联维数D= 2. 2的曲线。
【具体实施方式】
[0026] 参照图1-4。本发明区域电磁环境的动态特性表征方法具体步骤如下：
[0027] 步骤一、区域电磁环境信号表示。
[0028] N个福射源在某空间区域形成的电磁信号环境表不为L个点处电磁信号的集合：
[0029] {EjthEjt)，…上⑴} (1)
[0030] 其中，EiS接收点i处合成电磁信号，
[0031]
[0032]Ai⑴为电场的幅值，匕⑴为频率，％·(〇为相位，叾.为场强方向或极化方向， 福射天线方向图，为福射源i的j径信号的传递函数，Σ表不传输信号多径合成，Λf; (t) 为辐射源i相对接收点运动速度为Vi(t)时产生的频移。
[0033] 步骤二、基于数学形态学的电磁环境信号处理。
[0034] 通过形态学混合开、闭运算交替滤波，滤除电磁信号中的脉冲噪声和白噪声。具体
公式如下，
[0035] (3)
[0036] (4)
[0037] 式中屯⑴是接收点一维离散电磁信号，取值范围为Ejt)= {0，1，2，···，Ν-1} ;g(m)是结构元素，取值范围为G= {0，1，2，···，Μ-1};十表示膨 胀算子，表示戈
；?表示腐蚀算子，表示为
[0038] 步骤三、基于分形理论求解电磁环境关联维数D(m)。
[0039]
[0040]
[0041]
[0042]式中Η是Heaviside函数，当r>0，H(r)为 1，r< 0 是H(r)为 0。 {Ei，E2，…，Ε1; ···}为区域形成的电磁合成信号的集合，重构延时后得到新集合记为：
:，τ为延迟时间，τ=kΛt，Λt是时 间序列的采样时间间隔，N^N-On-l)τ为的总点数。
[0043] 参照表1首先运用Matlab模拟产生了四种不同调制方式的通信信号作为构成复 杂区域电磁环境的主要要素，其次对于其余的信号和杂波一并作为噪声信号模拟产生。运 用基于数学形态学和分型理论对该区域四中不同调制的通信信号及其相互作用后信号关 联维数的计算，解算其分行维数对区域电磁环境动态性进行特征。
[0044]表1仿真信号参数设置
[0045]
[0046]图1所示曲线表示区域中是1. 4GHz的QPSK信号时接收点集合的关联维数；直线 表示最小二乘拟合后，该区域的关联维数。
[0047] 图2所示曲线表示区域中是1. 4GHz的QPSK信号和1. 2GHz的16QAM合成时接收 点集合的关联维数；直线表示最小二乘拟合后，该区域的关联维数。
[0048] 图3所示曲线表示区域中是1. 1GHz的DPSK信号和1. 2GHz的ASK合成时接收点 集合的关联维数；直线表示最小二乘拟合后，该区域的关联维数。
[0049] 图4所示曲线表示区域中是1. 1GHz的DPSK信号、1. 2GHz的ASK、1. 4GHz的QPSK 信号和1. 2GHz的16QAM合成时接收点集合的关联维数；直线表示最小二乘拟合后，该区域 的关联维数。
【主权项】
1. 一种区域电磁环境的动态特性表征方法，其特征在于包括以下步骤： 步骤一、将N个辐射源在某空间区域形成的电磁信号环境表示为L个点处电磁信号的 集合： -,EL(t)} (1) 其中，EiS接收点i处合成电磁信号，式中，Ajt)为电场的幅值，&(〇为频率，&(?)为相位，g为场强方向或极化方向，FiS 福射天线方向图，为福射源i的j径信号的传递函数，Σ表不传输信号多径合成，Λf; (t) 为辐射源i相对接收点运动速度为Vi(t)时产生的频移； 步骤二、基于数学形态学的电磁环境信号处理； 通过形态学混合开、闭运算交替滤波，滤除电磁信号中的脉冲噪声和白噪声；具体公式 如下： 开运算：(Ei〇g) (n) =EiΘg?g (3) 闭运算：(Ei·g) (n) =Ei?gΘg (4) 式中，Ejt)是接收点一维离散电磁信号，取值范围为Ejt) = {0, 1，2,… ，N-1} ;g(m)是结构元素，取值范围为G= {0, 1，2, "·，Μ-1} ; ?表示膨胀算子，?Θ表示腐蚀算子，表示为步骤三、基于分形理论求解电磁环境关联维数D(m);式中，Η是Heaviside函数，当r>0,H(r)为 1，r< 0 是H(r)为 0 ; {Ε1;E2,…，…} 为区域形成的电磁合成信号的集合，重构延时后得到新集合记为：Yn(m，τ) = (En，En+1，… En+〇n1)τ)，η= 1，…，Nm，τ为延迟时间，τ=kΛt，Λt是时间序列的采样时间间隔，Nm = N-(m-l)τ为的总点数。
【专利摘要】本发明公开了一种区域电磁环境的动态特性表征方法，用于解决现有区域电磁环境的动态特性表征方法实用性差的技术问题。技术方案是：首先，基于信号基础理论对区域接收点处电磁信号集合的进行建模；其次，利用数学形态学滤波原理对接收的电磁信号集合进行交替滤波，滤除电磁信号中的脉冲噪声和白噪声；最后，运用分形理论求解电磁环境接收点集合的关联维数。对区域电磁环境动态性量化表征。本发明根据区域电磁信号表示和数学形态学理论，不仅对区域L个接收点N个辐射源产生的电磁信号进行了建模，并且滤除了环境中突发噪声，实现了对接收信号的动态性表征。
【IPC分类】G06F17/50
【公开号】CN105389431
【申请号】CN201510750250
【发明人】郭淑霞, 张磊, 高颖, 王亚锋, 陈旭, 于学伟, 董文华, 刘瑞兵, 李瑛
【申请人】西北工业大学
【公开日】2016年3月9日
【申请日】2015年11月6日