主权项 |
一种基于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>‾</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>M</mi></mfrac><munderover><mo>Σ</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分类器,识别频谱数据所属的类型。 |