1.一种基于扩展卡尔曼滤波测速方法，其特征在于，包括以下步骤：S1建立机车运行模型方程 函数f是关于机车轮对转速ω与机车车速x的状态方程，函数h是关于x的函数，wk、νk+1为独立同分布的零均值白噪声；

S2求解函数f、h关于xk的偏导Ak、H(k)；

S3基于扩展卡尔曼滤波方程进行机车车速预测：前一状态预测方程： 前一状态的协方差预测方程： 滤波器增益计算方程：其中Qk，R分别是高斯白噪声的协方差矩阵；

S4判断轨面粘着系数是否发生跳变：

当 时，轨面未发生跳变，进入S5；

当 时，轨面发生跳变，进入S6；

S5维持原扩展卡尔曼滤波状态矩阵完成测速：状态校正方程： 协方差误差校正方程： 最终输出预测速度

S6计算扩展卡尔曼滤波估计偏差，并更新卡尔曼滤波状态矩阵完成测速：经过扩展卡尔曼滤波器得到速度估计值 xk为机车运行理想速度，则卡尔曼滤波估计偏差Δxk为 那么状态更新方程： Δx偏差协方差为： 协方差误差校正方程： 最终输出预测速度

2.根据权利要求1所述的一种基于扩展卡尔曼滤波测速方法，其特征在于，包括以下步骤：将机车机车轮对转速ω与机车车速x的状态方程离散化得到函数f具体为：r为车轮半径，W为等效轴重，g为重力加速度，μ为轮轨间粘着系数，B为转动摩擦力矩系数，Rg为齿轮箱传动比，J为机车车轮转动惯量，uk为机车控制力矩，M为机车质量，a0、a1、a2为阻力系数；

h具体为[μ(k+1) x(k+1)]T，T为采样周期，μ为c·e-a(ω·r-x)-d·e-b(ω·r-x)，a、b、c、d为轨面参数。

3.根据权利要求1所述的一种基于扩展卡尔曼滤波测速方法，其特征在于，Ak具体为：H(k)具体为：

4.根据权利要求1所述的一种基于扩展卡尔曼滤波测速方法，其特征在于，步骤S4中的阈值具体为 xs为蠕滑速度具体为xs＝ω·r-x。

5.一种基于权利要求1-4任一一项所述的扩展卡尔曼滤波测速方法的测速系统，其特征在于，包括机车运行模块、扩展卡尔曼滤波器、轨面检测模块、滤波偏差估计模块、状态更新模块；

机车运行模块接收机车控制力矩uk，用于根据机车动力学模型与轮轨粘着模型建立状态方程，得到关于机车速度x和轮对速度ω的非线性状态方程。

扩展卡尔曼滤波器与机车运行模块连接，用于对机车车速进行实时预测；

轨面检测模块分别与扩展卡尔曼滤波器、滤波偏差估计模块连接，用于根据机车行驶过程中轮轨接触面粘着系数的变化，判断轨面是否发生跳变；

滤波偏差估计模块，用于当轨面发生跳变后，根据机车运行理想速度与经扩展卡尔曼滤波器估计的速度比较获取扩展卡尔曼滤波器的估计偏差值；

状态更新模块分别与滤波偏差估计模块、扩展卡尔曼滤波器连接，根据扩展卡尔曼滤波器的估计偏差更新状态矩阵并反馈至扩展卡尔曼滤波器。