[发明专利]一种应用于轨检仪基本弦轨向数据的滤波方法

摘要

本发明公开了一种应用于轨检仪基本弦轨向数据的滤波方法，包括以下步骤：1)首先对基本弦轨向原始数据的一维数组X求一阶前向差分得到X1；然后对X1求均值得到数Y1；再分别求取X1中每个数据点的偏差，得到偏差数组Y2；最后对Y2求方差，得到数Y3；2)对X求均值得到数Z1；然后分别求取X中每个数据点的偏差，得到偏差数组Z2；再对Z2求方差，得到数Z3；3)计算Z2中的每一个数据点的权值；4)再根据归一化处理后的权值对X进行加权计算，得到归一化值C；5)根据递推公式求出X经过滤波后的一维数组。本发明计算速度快，适合轨检仪处理数据速度的要求，且该方法能有效滤除信号中的尖噪声干扰。

主权项

1.一种应用于轨检仪基本弦轨向数据的滤波方法，其特征在于，包括以下步骤：步骤1：输入一组基本弦轨向原始数据，由此生成基本弦轨向原始数据的一维数组X，一维数组X中数据点的个数记为n+1；步骤2：首先对一维数组X求一阶前向差分，得到一阶前向差分数组X1；然后对一阶前向差分数组X1求均值，得到数Y1；再分别求取一阶前向差分数组X1中每个数据点的偏差，得到偏差数组Y2；最后对该偏差数组Y2求方差，得到数Y3；步骤3：对一维数组X求均值，得到数Z1；然后分别求取一维数组X中每个数据点的偏差，得到偏差数组Z2；再对该偏差数组Z2求方差，得到数Z3；步骤4：由步骤2和步骤3中得出的数据，对于偏差数组Z2中的每一个数据点Z2(m)，如果|Z2(m)|>2|Z3|或m＝n+1，则令该数据点的权值G(m)＝0,否则由此求得每个数据点的权值；步骤5：对步骤4中的权值G(m)(m＝1…,n+1)首先求和，得到权值和g，然后进行如下处理：再根据G1(m)对一维数组X进行加权计算，得到一个最优值，即归一化值C，步骤6：根据步骤5中的归一化值C和递推公式求出基本弦轨向原始数据一维数组X经过滤波后的一维数组X′，递推公式为：X′(i+2)＝X′(i+1)+C·[X(i+2)-X′(i+2)]，i＝0,1,…,n-2，其初始条件为X'(1)＝X(1)。