发明名称 基于EEMD特征提取和PNN的非平稳电磁干扰信号模式识别方法
摘要 本发明公开了基于EEMD特征提取和PNN的非平稳电磁干扰信号模式识别方法,属于电磁信号识别与分析研究领域;具体步骤为:首先,利用频谱仪分别测量不同类电磁设备的电磁干扰信号,并划分为训练样本和测试样本;然后,对所有训练样本进行层次EEMD分解,得到每个训练样本每层的固有模态函数IMF;接着,分别计算每个训练样本所有层IMF的能量分布和峭度,作为该训练样本的特征向量;并将特征向量和电磁设备的类别标签,进行训练得到PNN分类器;最后,利用PNN分类器,识别不同测试样本的频谱数据所属的类型;优点在于:利用EEMD方法对信号提取特征,具有更好的自适应性;将峭度和IMF能量分布作为特征向量,有效区分不同类型的电磁干扰信号。
申请公布号 CN106485209A 申请公布日期 2017.03.08
申请号 CN201610846687.X 申请日期 2016.09.23
申请人 北京航空航天大学 发明人 李红裔;赵迪;苏东林;宁博明;赵连坤;黄子晏
分类号 G06K9/00(2006.01)I;G06K9/62(2006.01)I 主分类号 G06K9/00(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 赵文利
主权项 一种基于EEMD特征提取和PNN的非平稳电磁干扰信号模式识别方法,其特征在于,具体步骤如下:步骤一、利用频谱仪分别测量实际电磁设备,获取不同类设备的电磁干扰信号,并将频谱数据划分为训练样本和测试样本;步骤二、针对所有训练样本,分别进行层次EEMD分解,得到每个训练样本每层的固有模态函数IMF;具体步骤为:步骤201、针对频谱数据中的某个训练样本X,初始化EEMD的数目,构造不同幅值的高斯白噪声;训练样本X长度为n,X=[x<sub>1</sub>,x<sub>2</sub>,...,x<sub>j</sub>,...x<sub>n</sub>];EEMD的数目为M;不同幅值的高斯白噪声为M组,第m组幅值的高斯白噪声N<sub>m</sub>=[n<sub>1</sub><sup>(m)</sup>,n<sub>2</sub><sup>(m)</sup>,...n<sub>j</sub><sup>(m)</sup>,...n<sub>n</sub><sup>(m)</sup>];m=1,2,...M;步骤202、在训练样本X上依次添加给定幅值的高斯白噪声,得到M个干扰信号;第m个被高斯白噪声干扰的信号为:X<sub>m</sub>=X‑N<sub>m</sub>;M个干扰信号集合为:[X<sub>1</sub>,X<sub>2</sub>,...X<sub>m</sub>,...X<sub>M</sub>];步骤203、将每个干扰信号利用EMD方法,分解为I层固有模态函数IMF,共得到M×I个固有模态函数;M×I个固有模态函数集合表示如下:<maths num="0001"><math><![CDATA[<mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>c</mi><mrow><mn>1</mn><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>c</mi><mrow><mn>1</mn><mo>,</mo><mn>2</mn></mrow></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>c</mi><mrow><mn>1</mn><mo>,</mo><mi>m</mi></mrow></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>c</mi><mrow><mn>1</mn><mo>,</mo><mi>M</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>c</mi><mrow><mn>2</mn><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>c</mi><mrow><mn>2</mn><mo>,</mo><mn>2</mn></mrow></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>c</mi><mrow><mn>2</mn><mo>,</mo><mi>m</mi></mrow></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>c</mi><mrow><mn>2</mn><mo>,</mo><mi>M</mi></mrow></msub></mtd></mtr><mtr><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd><mtd><mo>...</mo></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><msub><mi>c</mi><mrow><mi>i</mi><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>c</mi><mrow><mi>i</mi><mo>,</mo><mn>2</mn></mrow></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>c</mi><mrow><mi>i</mi><mo>,</mo><mi>m</mi></mrow></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>c</mi><mrow><mi>i</mi><mo>,</mo><mi>M</mi></mrow></msub></mtd></mtr><mtr><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd><mtd><mo>...</mo></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><msub><mi>c</mi><mrow><mi>I</mi><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>c</mi><mrow><mi>I</mi><mo>,</mo><mn>2</mn></mrow></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>c</mi><mrow><mi>I</mi><mo>,</mo><mi>m</mi></mrow></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>c</mi><mrow><mi>I</mi><mo>,</mo><mi>M</mi></mrow></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001119294040000011.GIF" wi="622" he="451" /></maths>c<sub>i,m</sub>表示第m个干扰信号X<sub>m</sub>分解的第i层固有模态函数;i=1,2,3,...I;m=1,2,3,...M;步骤204、针对每层固有模态函数IMF,分别计算该层M个固有模态函数的平均值作为该层最终的IMF值,得到训练样本X的每层固有模态函数IMF;训练样本X的第i层固有模态函数IMF<sub>i</sub>值为该层所有IMF的平均值,公式如下:<maths num="0002"><math><![CDATA[<mrow><msub><mi>IMF</mi><mi>i</mi></msub><mo>=</mo><mover><msub><mi>c</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>M</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>c</mi><mrow><mi>i</mi><mo>,</mo><mi>m</mi></mrow></msub></mrow>]]></math><img file="FDA0001119294040000012.GIF" wi="429" he="131" /></maths>步骤205、对所有训练样本重复上述步骤,得到每个训练样本每层的固有模态函数IMF。步骤三、针对每个训练样本,分别计算所有层IMF的能量分布和峭度,作为该训练样本的特征向量;步骤四、将所有训练样本的特征向量和电磁设备的类别标签,输入未训练的PNN进行训练得到PNN分类器;步骤五、对于每个测试样本,输入PNN分类器,识别频谱数据所属的类型。
地址 100191 北京市海淀区学院路37号