1.一种基于地磁信号的车型识别方法，其特征在于：步骤一、在地面下设置三分量磁场传感器；取m种不同车型的车辆；每个车型均对应W辆试验车辆；将1赋值给i和j；

步骤二、用第i种车型对应的第j辆试验车辆经过三分量磁场传感器的上方；三分量磁场传感器进行Nij次磁场采集，得到X方向磁场数据集S1、Y方向磁场数据集S2、Z方向磁场数据集S3；

其中，s1b为三分量磁场传感器在第b次采集中得到的磁场x轴方向分量，s2b为三分量磁场传感器在第b次采集中得到的磁场y轴方向分量，s3b为三分量磁场传感器在第b次采集中得到的磁场z轴方向分量，b＝1,2,…,Nij；

步骤三、对X方向磁场数据集S1、Y方向磁场数据集S2、Z方向磁场数据集S3内的元素分别进行归一化处理，得到X方向归一化数据集S′1、Y方向归一化数据集S′2、Z方向归一化数据集S′3；

其中， |sa，max|为

中最大值的绝对值，|sa，min|为 中最小值的绝对值，max(|sa，max|，|sa，min|)等于|sa，max|、|sa，min|中较大的那个数；

步骤四、求取X方向归一化数据集S′1的均值 最大值s′1，max、最小值s′1,min、极差R1、最大值的相对位置E1,max、最小值的相对位置E1,min、正负能量比PG1、极大值个数、极小值个数，Y方向归一化数据集S′2的均值 最大值s′2，max、最小值s′2,min、极差R2、最大值的相对位置E2,max、最小值的相对位置E2,min、正负能量比PG2、极大值个数、极小值个数，Z方向归一化数据集S′3的均值 最大值s′3,max、最小值s′3,min、极差R3、最大值的相对位置E3,max、最小值的相对位置E3,min、正负能量比PG3、极大值个数、极小值个数；

s′a,max为 中的最大值；s′a，min为 中的最小值；Ra＝sa，max-sa，min；Ea,max＝σa,max/Nij，σa,max为s′a，max在归一化数据集S′a中的序数；Ea,min＝σa,min/Nij，σa,min为s′a,min在归一化数据集S′a中的序数；Pga＝Aa/Ba；Aa为归一化数据集S′a中所有大于0的元素的平方和；Ba为归一化数据集S′a中所有小于0的元素的平方和，a＝1,2,

3；

计算X方向归一化数据集S′1、Y方向归一化数据集S′2、Z方向归一化数据集S′3的平均能量步骤五、a＝1,2,3，依次执行步骤六；

步骤六、对归一化数据集S′a进行二尺度小波包变换，得到四个二尺度分解系数集；对四个二尺度分解系数集分别取平均值、最大值、近似熵；并求取所得的第一个二尺度分解系数集的平均能量 其中，da，l为第一个二尺度分解系数集的第l个元素；Ma为第一个二尺度分解系数集的元素个数；

步骤七、将步骤四、五、六所得的X方向归一化数据集S′1的均值、最大值、最小值、极差、最大值的相对位置、最小值的相对位置、正负能量比、极大值个数、极小值个数，Y方向归一化数据集S′2的均值、最大值、最小值、极差、最大值的相对位置、最小值的相对位置、正负能量比、极大值个数、极小值个数，Z方向归一化数据集S′3的均值、最大值、最小值、极差、最大值的相对位置、最小值的相对位置、正负能量比、极大值个数、极小值个数，X方向归一化数据集S′1对应的第一个二尺度分解系数集的平均能量、四个二尺度分解系数集的平均值、最大值、近似熵，Y方向归一化数据集S′2对应的第一个二尺度分解系数集的平均能量、四个二尺度分解系数集的平均值、最大值、近似熵，Z方向归一化数据集S′3对应的第一个二尺度分解系数集的平均能量、四个二尺度分解系数集的平均值、最大值、近似熵，依次排序作为对应第i个车型的第j辆试验车辆的64个特征值；进入步骤八；

步骤八、若j＜W，则将j增大1，并执行步骤二至七；若i＜m，且j＝W，则将i增大1，将1赋值给j增大1，并执行步骤二至七；若i＝m，且j＝W，则进入步骤九；

步骤九、建立初始特征数据集F如下：

其中，fk,i,j为第i个车型的第j辆试验车辆对应的第k个特征值，i＝1,2,…,m，j＝1,

2,…,W，k＝1,2,…,64；

步骤十、采用混合式的特征选择方法进行特征筛选，并得到分类模型；

10.1、采用ReliefF算法计算初始特征数据集F内64个二维数组的ReliefF权重；将初始特征数据集F中ReliefF权重低于权重阈值Wt的二维数组删除，得到一级优化特征集合Fs；

其中，0.05≤Wt≤0.15；将1赋值给p，将2赋值给q；

10.2、计算一级优化特征集合Fs内第p个二维数组与第q个二维数组的相关系数其中， 若COR(p,

q)＞Wc，则将配对数组(p,q)放入初始为空集的Pos数组；其中，0.7≤Wc≤0.8；进入步骤

10.3；

10.3、若q＜r，则将q增大1，并执行步骤10.2；若q＝r，且p＜r-1，则将p增大1后，将p+1赋值给q，并执行步骤10.2；若q＝r，且p＝r-1，则进入步骤10.4；r为一级优化特征集合Fs内二维数组的个数；

10.4、若Pos数组不是空数组，则执行步骤10.5及10.6；否则，进入步骤10.7；

10.5、将Pos数组中重复次数最多的元素作为初始目标序数；若仅存在一个初始目标序数，则将该初始目标序数作为最终目标序数ξ；若同时存在多个初始目标序数，则将各初始目标序数对应的二维数组作为目标二维数组，取各目标二维数组中ReliefF权重最小的那个目标二维数组对应的初始目标序数作为最终目标序数ξ；

10.6、将最终目标序数ξ对应的那个二维数组标记为无效数组，并将Pos数组中含有最终目标序数ξ的配对数组均删除；进入步骤10.4；

10.7、将一级优化特征集合Fs中所有被标记为无效数组的二维数组均删除，得到二级优化特征集合Fb；

10.8、用包装算法将对二级优化特征集合Fb进行筛选，得到三级优化特征集合F'；三级优化特征集合F'内二维数组的个数为v；三级优化特征集合F'内v个二维数组对应的特征作为目标特征；

10.9、根据三级优化特征集合F'内的各个元素的大小及其对应的车型，通过C4.5算法获取三级优化特征集合F'内v个二维数组对应的特征的目标阈值；

步骤十一、当被测车辆经过三分量磁场传感器的上方时，三分量磁场传感器根据自身采样频率进行N'次磁场采集，得到X方向被测数据集C1、Y方向被测数据集C2、Z方向被测数据集C3；

步骤十二、根据X方向被测数据集C1、Y方向被测数据集C2、Z方向被测数据集C3求取被测车辆的v个目标特征值；将被测车辆的v个目标特征值与v个目标阈值进行大小比较，进而确定被测车辆的车型；

进一步地，步骤四中，磁场数据集的极大值、极小值的判定方法如下：若s′ab满足以下两个条件之一，则s′ab为极大值：(1)s′ab＞s′a,b-1，s′ab＞s′a,b+1，并且满足s′ab-s′a,b-1＞0.2或s′ab-s′a,b+1＞0.2；

(2)s′ab-s′a,b-1＞0.1，s′a,b-1-s′a,b-2＞0.1，s′ab-s′a,b+1＞0.1，s′a,b+1-s′a,b+2＞0.1；

若s′ab满足以下两个条件之一，则s′ab为极小值：(1)s′ab＜s′a,b-1，s′ab＜s′a,b+1，并且满足s′a,b-1-s′ab＞0.2或s′a,b+1-s′ab＞0.2；

(2)s′a,b-1-s′ab＞0.1，s′a,b-2-s′a,b-1＞0.1，s′a,b+1-s′ab＞0.1，s′a,b+2-s′a,b+1＞0.1。

2.根据权利要求1所述的一种基于地磁信号的车型识别方法，其特征在于：所述的二尺度小波包变换采用的小波函数为db4小波。

3.根据权利要求1所述的一种基于地磁信号的车型识别方法，其特征在于：所述的包装算法采用论文《Wrappers for Feature Subset Selection》中的描述的wrapper算法框架。