1.一种基于证据推理规则的轨道高低不平顺幅值估计方法，其特征在于该方法包括以下步骤：(1)设定GJ-4型轨检车安装在车轴与车体上的加速度计获得车轴和车厢位置的时域振动加速度信号为r1(t)和r2(t)，其单位为G(重力加速度，9.8m/s2)，轨检车以100-120km/h的时速，每隔0.15-0.3m采样一次加速度振动信号，共采集TS次，一般TS>5000，则采样时刻t＝

1,2,…,TS；设轨检车在采样时刻t利用惯性测量方法测得垂直位移为d(t)，其单位为mm；

(2)将步骤(1)中每个采样时刻的时域振动加速度信号r1(t)和r2(t)分别以5.25m的窗口长度进行短时傅里叶变换，然后求得各频域幅值绝对值的平均值作为车轴和车体的频域特征信号f1(t)和f2(t)， 其中 分别为输入特征信号f1(t)和f2(t)的最小和最大值；将垂直位移d(t)取绝对值记为Ir(t)，Ir(t)∈[l1，l2]，其中l1和l2分别是Ir(t)的最小和最大值；将f1(t)、f2(t)和Ir(t)表示成样本集合S＝{[f1(t),f2(t),Ir(t)]|t＝1,2,…,TS}，其中[f1(t),f2(t),Ir(t)]为一个样本向量；

(3)设定轨道高低不平顺幅值Ir的结果参考值集合D＝{Dn|n＝1,…,N}，频域振动信号fi的输入参考值集合 N为轨道高低不平顺幅值的结果参考值Ir的个数，Ji为频域振动信号fi的参考值个数；

(4)将T个样本向量[f1(t),f2(t),Ir(t)]中的样本对(fi(t),Ir(t))分别用定性信息转换方法变化为关于参考值相似度的形式，具体步骤如下：(4-1)样本对(fi(t),Ir(t))的输入值fi(t)匹配参考值 的相似度分布为其中

αi,j'＝0 j'＝1,...,Ji,j'≠j,j+1    (1c)αi,j表示输入值fi(t)匹配参考值 的相似度；

(4-2)样本对(fi(t),Ir(t))的结果值Ir(t)匹配参考值Dn的相似度分布为SO(Ir(t))＝{(Dn,γn)|n＝1,...,N}       (2a)其中

γn'＝0 n'＝1,...,N,n'≠n,n+1  (2c)γn表示结果值Ir(t)匹配参考值Dn的相似度；

(4-3)根据步骤(4-1)和步骤(4-2)，样本对(fi(t),Ir(t))可以被转化地表示为相似度分布的形式(αi,jγn,αi,j+1γn,αi,jγn+1,αi,j+1γn+1)，其中，αi,jγn表示样本对(fi(t),Ir(t))中输入值匹配参考值 同时结果值Ir(t)匹配参考值Dn的综合相似度；

(5)根据步骤(4)，将样本集S中的所有样本对转化为综合相似度的形式，用它们可构造结果参考值和输入参考值之间的投点统计表，如下表1所示，其中an,j表示所有输入值fi(t)匹配参考值 并且结果值Ir(t)匹配参考值Dn的样本对(fi(t),Ir(t))综合相似度的和，表示所有结果值Ir(t)匹配参考值Dn的样本对综合相似度的和， 表示所有输入值fi(t)匹配参考值 的样本对综合相似度的和，并有表1 样本对(fi(t),Ir(t))的投点统计表

(6)根据步骤(5)中的投点统计表，可获得当输入值fi(t)取参考值 时，结果值Ir(t)为参考值Dn的信度为并有 则可定义对应于参考值 的证据为

因此，可构造出如表2所示的证据矩阵表来描述输入fi和结果Ir之间的关系；

表2 输入fi的证据矩阵表

(7)定义证据的可靠性ri描述输入信息源fi评估轨道高低不平顺幅值Ir的能力，具体获取步骤如下：(7-1)定义输入值fi(t)与结果值Ir(t)的相对变化值为(7-2)根据(7-1)定义的相对变化值，可获取反映输入特征信号与不平顺幅值变化趋势的评价因子为(7-3)根据(7-2)获得的评价因子，可由下式计算输入信息源fi的可靠性(8)给定样本集中的任意一组输入样本向量F(t)＝(f1(t),f2(t))，根据步骤(6)获得的输入证据矩阵表和步骤(7)获得的输入信息源可靠性，可利用证据推理规则推理出初始的轨道高低不平顺估计值 具体步骤如下：(8-1)对于输入值fi(t)，其必然落入某两个参考值构成的区间 此时这两个参考值对应的证据 和 被激活，则输入值fi(t)的证据可由参考值证据 和 以加权和的形式获得ei＝{(Dn,pn,i),n＝1,...,N}        (9a)(8-2)利用式(9a)和式(9b)获得f1(t)和f2(t)的证据e1和e2，并设定初始证据权重wi＝ri，利用证据推理规则对它们进行融合，得到融合结果为O(F(t))＝{(Dn,pn,e(2)),n＝1,...,N}    (10a)(8-3)根据步骤(8-2)得到融合结果O(F(t))，估计的高低不平顺幅值 可由下式推理得到

(9)基于均方误差构建参数优化模型，具体步骤如下：

(9-1)确定优化参数集合 wi表示证据的权

重，其他参数分别设定为D1＝l1,DN＝l2,

(9-2)将最小化均方误差作为优化目标函数

s.t. 0≤wi≤1,i＝1,2       (12b)D2

式(12b)-(12d)表示优化参数需满足的约束条件；

(9-3)利用基于梯度的方法或者非线性优化软件包，获得最优的参数集合P，从现行列车车轴和车体上安装的加速度计采集输入特征信号，将其利用步骤(2)处理，再一次重复步骤(4)～步骤(8)即可得到更为精确的轨道高低不平顺估计值