1.一种陵区入侵振动信号特征提取与分类识别方法，其特征在于：该方法包括以下步骤：

步骤1：采集不高于200Hz频率的入侵低频振动信号：探测系统将通过分布式设置的多光栅传感器在同一时段采集到的入侵低频振动信号转化为电信号x(t)；

步骤2：EEMD分解：将转化之后的电信号x(t)加入一定幅值的高斯分布白噪声na(t)，合成后的信号表示为xa(t)＝na(t)+x(t)；再将已加入白噪声的合成信号xa(t)通过经验模态分解成m个本征模态函数(IMF)分量cj(t)和一个余项ra(t) ，表示为式中m表示IMF分量个数，再次针对各个本征模态信号分别加入

均方根相等的不同高斯白噪声，经由经验模态再次分解为若干个IMF分量与一个余项，至此提取得到N组不同的IMF分量，N满足 εn为待分析信号与分解得到IMF分量之和的标准偏差，ε为添加的白噪声的幅值；N为迭代次数，最后将N组同阶IMF进行总体平均计算得到：其中的cb(t)是通过EEMD分解得出的第b个IMF分量；综上，EEMD分解的最终结果为：步骤3：计算EEMD能量熵：信号经EEMD分解得到m个IMF分量c1(t)，c2(t)，c3(t)…cm(t)和余项r(t)，m个IMF分量的能量分别为E1,E2,E3...Em，其中， 余项r(t)的影响忽略不计，振动信号的能量分布为：E＝{E1,E2,E3...Em}，定义EEMD的能量熵值为式子中pk＝Ek/Q表示第k个IMF的能量比， H表征信号能量分布的不确定度；

步骤4：入侵信号的特征提取：振动信号经过EEMD算法分解之后得到不同尺度下的平稳信号，将敏感性强的峭度作为不同尺度平稳信号的特征向量进行提取，具体过程如下：首先计算经过EEMD方法分解得到的m个IMF分量的各个峭度：Tb表示第b个本征模态函数分量的峭度；p表示第b个分量的样本点数；再将峭度进行归一化处理，得到特征向量[T1',T2',T3',T4'......T’m]；

步骤5：采用粒子群算法优化后的SVM对入侵振动信号进行分类识别：将多个传感器中测得的若干振波中振幅最高的振波作为主振频率，其他作为分量，对应于粒子群算法中的个体最优解Pbest–主振频率分析和群体最优解Gbest–平均频率分析，以此得出的最优解作为SVM的最优决策函数，若SVM输出为1，则有入侵行为发生，进而继续判断是爆破行为还是盗掘行为，若输出为0则为非入侵行为。

2.根据权利要求1所述的一种陵区入侵振动信号特征提取与分类识别方法，其特征在于：步骤2中所述的高斯分布白噪声na(t)的幅值取采集信号幅值标准差的0.1～0.4倍。

3.根据权利要求1所述的一种陵区入侵振动信号特征提取与分类识别方法，其特征在于：步骤2中所述的迭代次数N取100～400次。