检查系统、检查方法以及程序

技术领域

本发明涉及检查系统、检查方法以及程序，能够特别良好地用于检查铁道车辆的轨道。本申请基于2018年12月10日向日本提交的特愿2018-230834号并主张优先权，将其内容全部援引于此。

背景技术

当铁道车辆在轨道上行驶时，轨道的位置由于来自铁道车辆的载荷而变化。当产生这种轨道变化时，铁道车辆有可能显现异常的举动。因此，在专利文献1中公开有如下技术：向描述轮轴的偏转的运动方程式中代入轮轴的偏转方向上的角位移、通过进行数据同化的滤波器求出的状态变量、以及在用于支承轴箱的部件上产生的前后方向的力即前后方向力的测定值，而导出铁道车辆的轨道不规则(通端不规则量等)。

现有技术文献

专利文献

专利文献1：国际公开第2017/164133号

发明内容

发明要解决的课题

但是，在专利文献1所记载的技术中，在进行数据同化时，使用前后方向力的测定值、轮轴以及转向架(根据需要还包括车身)各自的左右方向的加速度的测定值。这些测定值虽然不使用特殊的传感器也能够得到，但配置于铁道车辆的传感器的数量优选较少。

本发明是鉴于以上那样的问题点而完成的，其目的在于减少为了检测铁道车辆的轨道不规则而使用的传感器的数量。

用于解决课题的手段

本发明的检查系统的特征在于，具有：数据取得单元，作为通过使具有车身、转向架以及轮轴的铁道车辆在轨道上行驶而测定出的测定值的数据，取得前后方向力的测定值的数据；状态变量导出单元，使用上述前后方向力的测定值，导出在使用描述了上述铁道车辆的运动的运动方程式而构成的状态方程式中应当决定的变量即状态变量；以及轨道状态导出单元，导出反映上述轨道的状态的信息，上述前后方向力是在配置于上述轮轴与设置有该轮轴的上述转向架之间的部件中产生的前后方向的力，且是根据上述轮轴的偏转方向的角位移与设置有该轮轴的上述转向架的偏转方向的角位移之差来确定的力，上述部件是用于支承轴箱的部件，上述前后方向是沿着上述铁道车辆的行驶方向的方向，上述偏转方向是将与上述轨道垂直的方向即上下方向作为转动轴的转动方向，上述状态方程式是使用上述状态变量、上述前后方向力以及转换变量来描述的方程式，上述状态变量包含上述转向架的左右方向的位移及速度、上述转向架的偏转方向的角位移及角速度、上述转向架的横摇方向的角位移及角速度、上述轮轴的左右方向的位移及速度、以及安装于上述铁道车辆的空气弹簧的横摇方向的角位移，不包含上述轮轴的偏转方向的角位移及角速度，上述横摇方向是将上述前后方向作为转动轴的转动方向，上述转换变量是将上述轮轴的偏转方向的角位移与上述转向架的偏转方向的角位移相互进行转换的变量，上述轨道状态导出单元使用作为由上述状态变量导出单元导出的上述状态变量之一的上述转向架的偏转方向的角位移、以及上述转换变量的实际值，导出上述轮轴的偏转方向的角位移的推定值，并使用所导出的上述轮轴的偏转方向的角位移的推定值来导出反映上述轨道的状态的信息，上述转换变量的实际值使用上述前后方向力的测定值来导出，上述状态变量导出单元不使用得到上述前后方向力的测定值的期间内的、上述转向架、上述轮轴以及上述车身的左右方向的加速度的测定值，而导出上述状态变量。

本发明的检查方法的特征在于，具有：数据取得工序，作为通过使具有车身、转向架以及轮轴的铁道车辆在轨道上行驶而测定出的测定值的数据，取得前后方向力的测定值的数据；状态变量导出工序，使用上述前后方向力的测定值，导出在使用描述了上述铁道车辆的运动的运动方程式而构成的状态方程式中应当决定的变量即状态变量；以及轨道状态导出工序，导出反映上述轨道的状态的信息，上述前后方向力是在配置于上述轮轴与设置有该轮轴的上述转向架之间的部件中产生的前后方向的力，且是根据上述轮轴的偏转方向的角位移与设置有该轮轴的上述转向架的偏转方向的角位移之差来确定的力，上述部件是用于支承轴箱的部件，上述前后方向是沿着上述铁道车辆的行驶方向的方向，上述偏转方向是将与上述轨道垂直的方向即上下方向作为转动轴的转动方向，上述状态方程式是使用上述状态变量、上述前后方向力以及转换变量来描述的方程式，上述状态变量包含上述转向架的左右方向的位移及速度、上述转向架的偏转方向的角位移及角速度、上述转向架的横摇方向的角位移及角速度、上述轮轴的左右方向的位移及速度、以及安装于上述铁道车辆的空气弹簧的横摇方向的角位移，不包含上述轮轴的偏转方向的角位移及角速度，上述横摇方向是将上述前后方向作为转动轴的转动方向，上述转换变量是将上述轮轴的偏转方向的角位移与上述转向架的偏转方向的角位移相互进行转换的变量，上述轨道状态导出工序使用作为由上述状态变量导出工序导出的上述状态变量之一的上述转向架的偏转方向的角位移、以及上述转换变量的实际值，导出上述轮轴的偏转方向的角位移的推定值，并使用所导出的上述轮轴的偏转方向的角位移的推定值来导出反映上述轨道的状态的信息，上述转换变量的实际值使用上述前后方向力的测定值来导出，上述状态变量导出工序不使用得到上述前后方向力的测定值的期间内的、上述转向架、上述轮轴以及上述车身的左右方向的加速度的测定值，而导出上述状态变量。

本发明的程序的特征在于，使计算机执行：数据取得工序，作为通过使具有车身、转向架以及轮轴的铁道车辆在轨道上行驶而测定出的测定值的数据，取得前后方向力的测定值的数据；状态变量导出工序，使用上述前后方向力的测定值，导出在使用描述了上述铁道车辆的运动的运动方程式而构成的状态方程式中应当决定的变量即状态变量；以及轨道状态导出工序，导出反映上述轨道的状态的信息，上述前后方向力是在配置于上述轮轴与设置有该轮轴的上述转向架之间的部件中产生的前后方向的力，且是根据上述轮轴的偏转方向的角位移与设置有该轮轴的上述转向架的偏转方向的角位移之差来确定的力，上述部件是用于支承轴箱的部件，上述前后方向是沿着上述铁道车辆的行驶方向的方向，上述偏转方向是将与上述轨道垂直的方向即上下方向作为转动轴的转动方向，上述状态方程式是使用上述状态变量、上述前后方向力以及转换变量来描述的方程式，上述状态变量包含上述转向架的左右方向的位移及速度、上述转向架的偏转方向的角位移及角速度、上述转向架的横摇方向的角位移及角速度、上述轮轴的左右方向的位移及速度、以及安装于上述铁道车辆的空气弹簧的横摇方向的角位移，不包含上述轮轴的偏转方向的角位移及角速度，上述横摇方向是将上述前后方向作为转动轴的转动方向，上述转换变量是将上述轮轴的偏转方向的角位移与上述转向架的偏转方向的角位移相互进行转换的变量，上述轨道状态导出工序使用作为由上述状态变量导出工序导出的上述状态变量之一的上述转向架的偏转方向的角位移、以及上述转换变量的实际值，导出上述轮轴的偏转方向的角位移的推定值，使用导出的上述轮轴的偏转方向的角位移的推定值导出反映上述轨道的状态的信息，上述转换变量的实际值使用上述前后方向力的测定值来导出，上述状态变量导出单元不使用得到上述前后方向力的测定值的期间内的、上述转向架、上述轮轴以及上述车身的左右方向的加速度的测定值，而导出上述状态变量。

附图说明

图1是表示铁道车辆的概要的一例的图。

图2是概念性地表示铁道车辆的构成要素的主要运动的方向的图。

图3是表示转向架的左右方向上的加速度、轮轴的左右方向上的加速度各自的测定值以及计算值的图。

图4A是表示直线轨道的通端不规则量的一例的图。

图4B是表示曲线轨道的通端不规则量的一例的图。

图5是表示检查装置的功能性构成的一例的图。

图6是表示检查装置的硬件构成的一例的图。

图7是检查装置的处理的一例的流程图。

图8是表示自相关矩阵的特征值的分布的一例的图。

图9是表示前后方向力的测定值的时间序列数据(测定值)、以及前后方向力的预测值的时间序列数据(计算值)的一例的图。

图10是表示前后方向力的高频分量的时间序列数据的一例的图。

图11是表示检查系统的构成的一例的图。

图12是表示计算例的图，且是表示通端不规则量的导出对象的轨道的曲率1/R、以及铁道车辆的行驶速度的图。

图13A是表示计算例的图，且是表示自相关矩阵R的特征值的分布的第1例的图。

图13B是表示计算例的图，且是表示自相关矩阵R的特征值的分布的第2例的图。

图14是表示计算例的图，且是表示前后方向力的测定值的时间序列数据、以及前后方向力的预测值的时间序列数据的图。

图15是表示计算例的图，且是表示前后方向力的高频分量的时间序列数据的图。

图16A是表示基于第1实施方式的方法的计算例的图，且是表示通端不规则量y

图16B是表示基于第1实施方式的方法的计算例的图，且是表示通端不规则量y

图17A是表示基于第2实施方式的方法的计算例的图，且是表示通端不规则量y

图17B是表示基于第2实施方式的方法的计算例的图，且是表示通端不规则量y

具体实施方式

以下，参照附图对本发明的实施方式进行说明。

(构思)

首先，对在实现本发明的实施方式时本发明人得到的构思进行说明。

图1是表示铁道车辆的概要的一例的图。另外，在图1中，铁道车辆向x轴的正方向前进(x轴是沿着铁道车辆的行驶方向的轴)。此外，z轴是与轨道16(地面)垂直的方向(铁道车辆的高度方向)。y轴是与铁道车辆的行驶方向垂直的水平方向(与铁道车辆的行驶方向和高度方向的双方垂直的方向)。此外，铁道车辆是营业车辆。另外，在各图中，在○中附加●的情况表示从纸面的里侧朝向跟前侧的方向，在○中附加×的情况表示从纸面的跟前侧朝向里侧的方向。

如图1所示，在本实施方式中，铁道车辆具有车身11、转向架12a、12b以及轮轴13a～13d。如此，在本实施方式中，以一个车身11具备两个转向架12a、12b以及4组轮轴13a～13d的铁道车辆为例进行说明。轮轴13a～13d具有车轴15a～15d以及设置于其两端的车轮14a～14d。在本实施方式中，以转向架12a、12b为无轴梁转向架的情况为例进行说明。另外，在图1中，为了便于表述，仅示出轮轴13a～13d的一方的车轮14a～14d，但在轮轴13a～13d的另一方也设置有车轮(在图1所示的例子中，车轮合计为8个)。此外，铁道车辆具有图1所示的构成要素以外的构成要素(在后述的运动方程式中说明的构成要素等)，但为了便于表述，在图1中省略该构成要素的图示。例如，转向架12a、12b具有转向构架以及枕簧等。此外，在各轮轴13a～13d的沿着y轴的方向的两侧配置有轴箱。此外，转向构架与轴箱通过轴箱支承装置相互结合。轴箱支承装置是配置在轴箱与转向构架之间的装置(悬架)。轴箱支承装置吸收从轨道16向铁道车辆传递的振动。此外，轴箱支承装置在限制了轴箱相对于转向构架的位置的状态下支承轴箱，以便抑制轴箱相对于转向构架向沿着x轴的方向以及沿着y轴的方向移动(优选不产生该移动)。轴箱支承装置配置在各轮轴13a～13d的沿着y轴的方向的两侧。另外，铁道车辆本身能够通过公知的技术实现，因此在此处省略其详细说明。

当铁道车辆在轨道16上行驶时，车轮14a～14d与轨道16之间的作用力(蠕变力)成为振动源，振动依次传递到轮轴13a～13d、转向架12a、12b、车身11。图2是概念性地表示铁道车辆的构成要素(轮轴13a～13d、转向架12a、12b、车身11)的主要运动的方向的图。图2所示的x轴、y轴、z轴分别与图1所示的x轴、y轴、z轴对应。

如图2所示，在本实施方式中，以轮轴13a～13d、转向架12a、12b以及车身11进行以x轴为转动轴转动的运动、以z轴为转动轴转动的运动以及沿着y轴的方向的运动的情况为例进行说明。在以下的说明中，根据需要将以x轴为转动轴转动的运动称作横摇，根据需要将以x轴为转动轴的转动方向称作横摇方向，根据需要将沿着x轴的方向称作前后方向。另外，前后方向是铁道车辆的行驶方向。在本实施方式中，沿着x轴的方向是铁道车辆的行驶方向。此外，根据需要将以z轴作为转动轴转动的运动称作偏转，根据需要将以z轴为转动轴的转动方向称作偏转方向，根据需要将沿着z轴的方向称作上下方向。另外，上下方向是与轨道16垂直的方向。此外，根据需要将沿着y轴的方向的运动称作横向振动，根据需要将沿着y轴的方向称作左右方向。另外，左右方向是与前后方向(铁道车辆的行驶方向)和上下方向(与轨道16垂直的方向)的双方垂直的方向。此外，铁道车辆也进行其他运动，但在各实施方式中为了简化说明而不考虑这些运动。但是，也可以考虑这些运动。

[第1构思]

在专利文献1所记载的技术中，将轮轴13a、13b、13c、13d的左右方向上的加速度y

图3表示转向架12a的左右方向上的加速度y

如图3所示，在轮轴13a、13b的左右方向上的加速度y

[第2构思]

如专利文献1所记载的那样，本发明人想到了如下方法：使用在配置于轮轴13a～13b(13c～13d)与设置有该轮轴13a～13b(13c～13d)的转向架12a(12b)之间的部件中产生的前后方向的力的测定值，来计算通端不规则量。在以下的说明中，根据需要将在该部件中产生的前后方向的力称作前后方向力。

使用基于对铁道车辆在直线轨道上行驶时的运动进行描述的运动方程式的公式、且是表示通端不规则量与前后方向力之间的关系的公式，来计算通端不规则量。轨道16包括直线部以及曲线部。在以下的说明中，根据需要将轨道16的直线部称作直线轨道，根据需要将轨道16的曲线部称作曲线轨道。

在通过进行数据同化的滤波器(卡尔曼滤波器)进行滤波的情况下，当使用对在曲线轨道上行驶的铁道车辆的运动进行描述的运动方程式来构成状态方程式时，状态变量有可能发散。因此，使用对在直线轨道上行驶的铁道车辆的运动进行描述的运动方程式，来构成通过进行数据同化的滤波器(卡尔曼滤波器)进行滤波的情况下的状态方程式。

在对在曲线轨道上行驶的铁道车辆的运动进行描述的运动方程式中，需要考虑在行驶时铁道车辆所受到的离心力等。因而，在对在曲线轨道上行驶的铁道车辆的运动进行描述的运动方程式中包括含有轨条(钢轨)的曲率半径的项。因此，当铁道车辆在曲线轨道上行驶时，当使用通过利用对在直线轨道上行驶的铁道车辆的运动进行描述的运动方程式而构成的进行数据同化的滤波器(卡尔曼滤波器)来导出状态变量时，有可能无法高精度地导出状态变量。

本发明人着眼于如下情况：在铁道车辆在曲线轨道上行驶的情况下，相对于在直线轨道上行驶时，前后方向力的测定值具有一定的偏置。由于通端不规则而引起的前后方向力的分量本身，不论是在曲线轨道中还是在直线轨道中都同样地产生。因此，本发明人认为通端不规则量本身与上述偏置的量无关，通过从前后方向力的测定值的时间序列数据中降低低频分量(上述偏置的举动)，由此即便使用基于对铁道车辆在直线轨道上行驶时的运动进行描述的运动方程式的公式来构成进行数据同化的滤波器(卡尔曼滤波器)，也能够从状态变量的推定值中降低由于铁道车辆在曲线轨道上行驶而引起的低频分量。根据该情况，本发明人想到了：使用降低了低频分量的前后方向力的值的时间序列数据来计算通端不规则量。通过如此地计算通端不规则量，即便使用基于对铁道车辆在直线轨道上行驶时的运动进行描述的运动方程式的公式，也能够计算出曲线轨道的通端不规则量。此外，不论是在曲线轨道还是在直线轨道，通端不规则量的计算式都成为相同的计算式。另外，即便在设计上是直线轨道，在实际上有时也具有对通端不规则量的推定精度造成影响的程度的曲率。因而，不仅在曲线轨道中，在直线轨道中，从前后方向力的测定值的时间序列数据中降低低频分量(上述偏置的举动)，也有助于通端不规则量的推定精度的提高。以下，对于即便在设计上是直线轨道、但在实际上具有对通端不规则量的推定精度造成影响的程度的曲率的轨道，也作为曲线轨道来进行说明。

(运动方程式)

接着，说明对铁道车辆的运动进行描述的运动方程式的一例。在本实施方式中，以专利文献1所记载的运动方程式为例，且以铁道车辆具有21个自由度的情况为例进行说明。即，轮轴13a～13d进行左右方向上的运动(横向振动)以及偏转方向上的运动(偏转)(2×4组＝8个自由度)。此外，转向架12a、12b进行左右方向上的运动(横向振动)、偏转方向上的运动(偏转)以及横摇方向上的运动(横摇)(3×2组＝6个自由度)。此外，车身11进行左右方向上的运动(横向振动)、偏转方向上的运动(偏转)以及横摇方向上的运动(横摇)(3×1组＝3个自由度)。此外，对于转向架12a、12b分别设置的空气弹簧(枕簧)进行横摇方向上的运动(横摇)(1×2组＝2个自由度)。此外，对于转向架12a、12b分别设置的偏转阻尼器进行偏转方向上的运动(偏转)(1×2组＝2个自由度)。

另外，自由度并不限定于21个自由度。如果增大自由度，则虽然计算精度提高，但计算负荷变高。此外，后述的卡尔曼滤波器的动作有可能变得不稳定。能够考虑这些方面而适当地决定自由度。此外，例如，通过基于非专利文献1的记载来表示各个构成要素(车身11、转向架12a、12b、轮轴13a～13d)的各个方向(左右方向、偏转方向、横摇方向)的动作，由此能够实现以下的运动方程式。因而，在此处，说明各个运动方程式的概要，省略详细的说明。另外，在以下的各式中不存在含有轨道16(轨条)的曲率半径(曲率)的项。即，以下的各式是表现铁道车辆在直线轨道上行驶的公式。在表现铁道车辆在曲线轨道上行驶的公式中，通过将轨道16(轨条)的曲率半径设为无限大(曲率为0(零))，能够得到表现铁道车辆在直线轨道上行驶的公式。

在以下的各式中，下标w表示轮轴13a～13d。(仅)附加有下标w的变量表示在轮轴13a～13d中是共同的。下标w1、w2、w3、w4分别表示轮轴13a、13b、13c、13d。

下标t、T表示转向架12a、12b。(仅)附加有下标t、T的变量表示在转向架12a、12b中是共同的。下标t1、t2分别表示转向架12a、12b。

下标b、B表示车身11。

下标x表示前后方向或者横摇方向，下标y表示左右方向，下标z表示上下方向或者偏转方向。

此外，在变量上附加的“··”、“·”分别表示2阶时间微分、1阶时间微分。

另外，在说明以下的运动方程式时，根据需要，省略出现过的变量的说明。此外，运动方程式本身与专利文献1所记载的运动方程式相同。

[轮轴的横向振动]

描述轮轴13a～13d的横向振动(左右方向上的运动)的运动方程式由以下的(1)式～(4)式表示。

[数式1]

m

[轮轴的偏转]

描述轮轴13a～13d的偏转的运动方程式由以下的(5)式～(8)式表示。

[数式2]

I

此处，所谓通端不规则，如日本工业标准(JIS E 1001:2001)所记载的那样，是指轨条的长度方向的左右的位移。通端不规则量是该位移的量。图4A以及图4B表示轮轴13a的位置处的通端不规则量y

[转向架的横向振动]

描述转向架12a、12b的横向振动(左右方向上的运动)的运动方程式由以下的(9)式、(10)式表示。

[数式3]

m

[转向架的偏转]

描述转向架12a、12b的偏转的运动方程式由以下的(11)式、(12)式表示。

[数式4]

I

[转向架的横摇]

描述转向架12a、12b的横摇的运动方程式由以下的(13)式、(14)式表示。

[数式5]

I

[车身的横向振动]

描述车身11的横向振动(左右方向上的运动)的运动方程式由以下的(15)式表示。

[数式6]

m

[车身的偏转]

描述车身11的偏转的运动方程式由以下的(16)式表示。

[数式7]

I

[车身的横摇]

描述车身11的横摇的运动方程式由以下的(17)式表示。

[数式8]

I

[偏转阻尼器的偏转]

描述配置于转向架12a的偏转阻尼器、配置于转向架12b的偏转阻尼器的偏转的运动方程式分别由以下的(18)式、(19)式表示。

[数式9]

ψ

[空气弹簧(枕簧)的横摇]

描述配置于转向架12a的空气弹簧(枕簧)、配置于转向架12b的空气弹簧(枕簧)的横摇的运动方程式分别由以下的(20)式、(21)式表示。

[数式10]

(前后方向力)

接着，对前后方向力进行说明。另外，前后方向力本身与专利文献1所记载的前后方向力相同。

一个轮轴上的左右车轮中的一方车轮的纵向蠕变力与另一方车轮的纵向蠕变力的同相的分量，是与制动力、驱动力对应的分量。因而，优选以与纵向蠕变力的反相分量对应的方式确定前后方向力。所谓纵向蠕变力的反相分量，是指一个轮轴上的左右车轮中的一方车轮的纵向蠕变力与另一方车轮的纵向蠕变力的相互成为相反相位的分量。即，所谓纵向蠕变力的反相分量，是指纵向蠕变力的使车轴扭转的方向的分量。在该情况下，前后方向力是在安装于一个轮轴的左右方向两侧的两个上述部件中产生的力的前后方向的分量中、相互成为相反相位的分量。

以下，说明以与纵向蠕变力的反相分量对应的方式确定前后方向力的情况下的前后方向力的具体例。

在轴箱支承装置是单连杆式的轴箱支承装置的情况下，轴箱支承装置具备连杆，轴箱与转向构架通过连杆连结。在该连杆的两端安装有橡胶衬套。在该情况下，前后方向力是在一个轮轴的左右方向端部分别各安装有一个的两个连杆各自承受的载荷的前后方向的分量中、相互成为相反相位的分量。此外，根据连杆的配置以及构成，连杆在前后方向、左右方向、上下方向的载荷中主要承受前后方向的载荷。因而，例如，只要对各连杆安装一个应变仪即可。通过使用该应变仪的测定值来导出该连杆所承受的载荷的前后方向的分量，由此得到前后方向力的测定值。此外，也可以代替这种方式，而通过位移仪来测定安装于连杆的橡胶衬套的前后方向的位移。在该情况下，将测定出的位移与该橡胶衬套的弹簧常数之积设为前后方向力的测定值。在轴箱支承装置是单连杆式的轴箱支承装置的情况下，上述的用于支承轴箱的部件为连杆或者橡胶衬套。

另外，在由安装于连杆的应变仪测定出的载荷中，有时不仅包含前后方向的分量，而且包含左右方向的分量和上下方向的分量中的至少任一方的分量。但是，即使在这种情况下，在轴箱支承装置的构造上，连杆所承受的左右方向的分量的载荷以及上下方向的分量的载荷与前后方向的分量的载荷相比足够小。因而，仅通过对各连杆安装一个应变仪，就能够得到具有实用上所要求的精度的前后方向力的测定值。如此，在前后方向力的测定值中有时会包含前后方向的分量以外的分量。因而，也可以对各连杆安装3个以上的应变仪，以便消除上下方向以及左右方向的应变。如此，能够提高前后方向力的测定值的精度。

在轴箱支承装置是轴梁式的轴箱支承装置的情况下，轴箱支承装置具备轴梁，轴箱与转向构架通过轴梁连结。轴梁也可以与轴箱一体地构成。在该轴梁的转向构架侧的端部安装有橡胶衬套。在该情况下，前后方向力成为在一个轮轴的左右方向的端部分别各安装有一个的两个轴梁分别承受的载荷的前后方向的分量中、相互成为相反相位的分量。此外，根据轴梁的配置构成，轴梁除了承受前后方向、左右方向、上下方向的载荷中的前后方向的载荷之外，还容易承受左右方向的载荷。因而，例如，对各轴梁安装两个以上的应变仪，以便消除左右方向的应变。通过使用这些应变仪的测定值来导出轴梁所承受的载荷的前后方向的分量，由此得到前后方向力的测定值。此外，也可以代替这种方式，而通过位移仪来测定安装于轴梁的橡胶衬套的前后方向的位移。在该情况下，将测定出的位移与该橡胶衬套的弹簧常数之积设为前后方向力的测定值。在轴箱支承装置是轴梁式的轴箱支承装置的情况下，上述的用于支承轴箱的部件为轴梁或者橡胶衬套。

另外，在由安装于轴梁的应变仪测定的载荷中，有时不仅包含前后方向以及左右方向的分量，而且包含上下方向的分量。但是，即使在这种情况下，在轴箱支承装置的构造上，轴梁所承受的上下方向的分量的载荷与前后方向的分量的载荷以及左右方向的分量的载荷相比足够小。因而，即使不以消除轴梁所承受的上下方向的分量的载荷的方式安装应变仪，也能够得到具有实用上所要求的精度的前后方向力的测定值。如此，在计测出的前后方向力中有时会包含前后方向的分量以外的分量，也可以对各轴梁安装3个以上的应变仪，以便除了消除左右方向的应变以外还消除上下方向的应变。如此，能够提高前后方向力的测定值的精度。

在轴箱支承装置是板簧式的轴箱支承装置的情况下，轴箱支承装置具备板簧，轴箱与转向构架通过板簧连结。在该板簧的端部安装有橡胶衬套。在该情况下，前后方向力成为在一个轮轴的左右方向的端部分别各安装有一个的两个板簧分别承受的载荷的前后方向的分量中、相互成为相反相位的分量。此外，根据板簧的配置构成，板簧除了承受前后方向、左右方向、上下方向的载荷中的前后方向的载荷之外，还容易受到左右方向的载荷以及上下方向的载荷。因而，例如，对各板簧安装3个以上的应变仪，以便消除左右方向以及上下方向的应变。通过使用这些应变仪的测定值来导出板簧所承受的载荷的前后方向的分量，由此得到前后方向力的测定值。此外，也可以代替这种方式，而通过位移仪来测定安装于板簧的橡胶衬套的前后方向的位移。在该情况下，将测定出的位移与该橡胶衬套的弹簧常数之积设为前后方向力的测定值。在轴箱支承装置是板簧式的轴箱支承装置的情况下，上述的用于支承轴箱的部件为板簧或者橡胶衬套。

另外，作为上述位移仪，能够使用公知的激光位移仪、涡流式的位移仪。

此外，在此处，以轴箱支承装置的方式为单连杆式、轴梁式以及板簧式的情况为例，对前后方向力进行了说明。但是，轴箱支承装置的方式并不限定于单连杆式、轴梁式以及板簧式。根据轴箱支承装置的方式，能够与单连杆式、轴梁式以及板簧式同样地确定前后方向力。

此外，在以下，为了简化说明，以对于一个轮轴得到一个前后方向力的测定值的情况为例进行说明。即，图1所示的铁道车辆具有4个轮轴13a～13d。因而，得到4个前后方向力T

(第1实施方式)

接着，对本发明的第1实施方式进行说明。

<检查装置500>

图5是表示检查装置500的功能性构成的一例的图。图6是表示检查装置500的硬件构成的一例的图。图7是表示检查装置500的处理的一例的流程图。在本实施方式中，如图1所示，以检查装置500搭载于铁道车辆的情况为例进行表示。

在图5中，检查装置500作为其功能而具有存储部501、数据取得部502、第1频率调整部503、状态变量导出部504、第2频率调整部505、轨道状态导出部506以及输出部507。

在图6中，检查装置500具有CPU601、主存储装置602、辅助存储装置603、通信电路604、信号处理电路605、图像处理电路606、I/F电路607、用户接口608、显示器609以及总线610。

CPU601统括控制检查装置500的整体。CPU601将主存储装置602用作为工作区域，执行存储于辅助存储装置603的程序。主存储装置602暂时保存数据。辅助存储装置603除了存储由CPU601执行的程序之外，还存储各种数据。辅助存储装置603存储后述状态方程式以及观测方程式。存储部501例如通过使用CPU601以及辅助存储装置603来实现。

通信电路604是用于与检查装置500的外部进行通信的电路。通信电路604例如接收前后方向力的测定值的信息。通信电路604与检查装置500的外部可以进行无线通信、也可以进行有线通信。通信电路604在进行无线通信的情况下，与设置于铁道车辆的天线连接。

信号处理电路605对于由通信电路604接收到的信号、根据CPU601进行的控制而输入的信号进行各种信号处理。数据取得部502例如通过使用CPU601、通信电路604以及信号处理电路605来实现。此外，第1频率调整部503、状态变量导出部504、第2频率调整部505以及轨道状态导出部506例如通过使用CPU601以及信号处理电路605来实现。

图像处理电路606对于根据CPU601进行的控制而输入的信号进行各种图像处理。进行了该图像处理的信号被输出至显示器609。

用户接口608是操作人员对检查装置500进行指示的部分。用户接口608例如具有按钮、开关以及转盘等。此外，用户接口608也可以具有使用了显示器609的图形用户接口。

显示器609显示基于从图像处理电路606输出的信号的图像。I/F电路607在与I/F电路607连接的装置之间进行数据的交换。在图6中，作为与I/F电路607连接的装置，示出用户接口608以及显示器609。但是，与I/F电路607连接的装置并不限定于这些。例如，便携式的存储介质也可以与I/F电路607连接。此外，用户接口608的至少一部分以及显示器609也可以处于检查装置500的外部。

输出部507例如通过使用通信电路604以及信号处理电路605、图像处理电路606、I/F电路607以及显示器609中的至少任一方来实现。

另外，CPU601、主存储装置602、辅助存储装置603、信号处理电路605、图像处理电路606以及I/F电路607与总线610连接。这些构成要素之间的通信经由总线610来进行。此外，检查装置500的硬件只要能够实现后述检查装置500的功能，则不限定于图6所示的硬件。

[存储部501]

存储部501存储后述状态变量导出部504在导出状态变量时使用的方程式。

在本实施方式中，存储部501存储状态方程式以及观测方程式。

在本实施方式中，以使用专利文献1所记载的状态方程式以及观测方程式的情况为例进行说明。

首先，对状态方程式进行说明。

在本实施方式中，在状态方程式中不包含(5)式～(8)式(描述轮轴13a～13d的偏转的运动方程式)，而如以下那样构成状态方程式。

首先，对于(9)式、(10)式(描述转向架12a、12b的横向振动(左右方向上的运动)的运动方程式)、(13)式、(14)式(描述转向架12a、12b的横摇的运动方程式)、(15)式(描述车身11的横向振动(左右方向上的运动)的运动方程式)、(16)式(描述车身11的偏转的运动方程式)、(17)式(描述车身11的横摇的运动方程式)、(18)式、(19)式(描述配置于转向架12a的偏转阻尼器、配置于转向架12b的偏转阻尼器的偏转的运动方程式)、(20)式、(21)式(描述配置于转向架12a的空气弹簧(枕簧)、配置于转向架12b的空气弹簧(枕簧)的横摇的运动方程式)，直接使用它们来构成状态方程式。

另一方面，在(1)式～(4)式(描述轮轴13a～13d的横向振动(左右方向上的运动)的运动方程式)、(11)式、(12)式(描述转向架12a、12b的偏转的运动方程式)中，包含轮轴13a～13d的偏转方向上的转动量(角位移)ψ

首先，轮轴13a～13d的前后方向力T

[数式11]

如以下的(26)式～(29)式那样定义转换变量e

[数式12]

e

e

e

e

当对(26)式～(29)式进行式变形时，得到以下的(30)式～(33)式。

[数式13]

ψ

ψ

ψ

ψ

当将(30)式～(33)式代入(1)式～(4)式的描述轮轴13a～13d的横向振动(左右方向上的运动)的运动方程式时，得到以下的(34)式～(37)式。

[数式14]

如此，通过使用转换变量e

当将(22)式～(25)式代入(11)式、(12)式(描述转向架12a、12b的偏转的运动方程式)时，得到以下的(38)式、(39)式。

[数式15]

如此，通过使用前后方向力T

此外，当将(26)式～(29)式代入(22)式～(25)式时，得到以下的(40)式～(43)式。

[数式16]

如以上那样，在本实施方式中，如(34)式～(37)式那样表示描述轮轴13a～13d的横向振动(左右方向上的运动)的运动方程式，并且如(38)式、(39)式那样表示描述转向架12a、12b的偏转的运动方程式。使用(34)式～(39)式来构成状态方程式。此外，(40)式～(43)式是常微分方程式。作为该常微分方程式的解的转换变量e

将如此求出的转换变量e

在本实施方式中，将以下的(44)式所示的变量设为状态变量，使用(9)式、(10)式、(13)式～(21)式、(34)式～(39)式的运动方程式来构成状态方程式。

[数式17]

存储部501例如基于操作人员对用户接口608的操作来输入并存储如以上那样构成的状态方程式。

接着，对观测方程式进行说明。

在本实施方式中，将车身11的左右方向上的加速度、转向架12a、12b的左右方向上的加速度以及轮轴13a～13d的左右方向上的加速度设为观测变量。该观测变量是由后述卡尔曼滤波器进行滤波的观测变量。在本实施方式中，使用(34)式～(37)式、(9)式、(10)式以及(15)式(描述横向振动的运动方程式)来构成观测方程式。

存储部501例如基于操作人员对用户接口608的操作来输入并存储如此构成的观测方程式。

在如以上那样将状态方程式以及观测方程式存储于检查装置500之后，数据取得部502、第1频率调整部503、状态变量导出部504、第2频率调整部505、轨道状态导出部506以及输出部507启动。即，在状态方程式以及观测方程式被存储于检查装置500之后开始基于图7的流程图的处理。

[数据取得部502、S701]

数据取得部502取得前后方向力的测定值的时间序列数据。前后方向力的测定方法如上所述。数据取得部502例如通过与使用用于测定前后方向力的应变仪的测定值来运算前后方向力的运算装置进行通信，由此能够取得前后方向力的测定值的时间序列数据。另外，数据取得部502不取得车身11的左右方向上的加速度的测定值的时间序列数据、转向架12a、12b的左右方向上的加速度的测定值的时间序列数据以及轮轴13a～13d的左右方向上的加速度的测定值的时间序列数据。

[第1频率调整部503、S702]

第1频率调整部503降低(优选除去)由数据取得部502取得的前后方向力的测定值的时间序列数据中所包含的低频分量的信号强度。该低频分量的信号是铁道车辆在曲率为0(零)的直线轨道上行驶的情况下计测不到、但铁道车辆在曲线轨道上行驶的情况下计测到的信号。即，铁道车辆在曲线轨道上行驶的情况下计测到的信号，能够视为对铁道车辆在曲率为0(零)的直线轨道上行驶的情况下计测到的信号叠加了该低频分量的信号而得到的信号。

本发明人研究出对自回归模型(AR(Auto-regressive)模型)进行了修正的模型。而且，本发明人想到了使用该模型来降低前后方向力的测定值的时间序列数据中所包含的低频分量的信号强度。在以下的说明中，将本发明人研究出的模型称作修正自回归模型。与此相对，将公知的自回归模型简称为自回归模型。以下，对修正自回归模型的一例进行说明。

将时刻k(1≤k≤M)的物理量y的时间序列数据的值设为y

[数式18]

(45)式中的α是自回归模型的系数。m是在自回归模型中用于对时刻k的数据y的值y

接着，使用最小二乘法，求出用于使基于自回归模型的时刻k的物理量的预测值y^

[数式19]

根据(46)式，以下的(47)式的关系成立。

[数式20]

此外，通过对(47)式进行变形(矩阵表述)，得到以下的(48)式。

[数式21]

(48)式中的R

[数式22]

基于(48)式，考虑到以下的(50)式。(50)式是根据使基于自回归模型的时刻k的物理量的预测值y^

[数式23]

此外，将(50)式中右边的自相关矩阵(由R

[数式24]

一般情况下，在求出自回归模型的系数时，使用针对系数α求解(50)式这样的方法。在(50)式中，以使通过自回归模型导出的时刻k的物理量的预测值y^

因此，本发明人着眼于与自回归模型的系数α相乘的自相关矩阵R，并进行了锐意研究。其结果，本发明人发现，使用自相关矩阵R的特征值的一部分，能够降低数据y所包含的高频分量的影响。即，本发明人发现，能够以使低频分量被强调的方式改写自相关矩阵R。

以下，对该情况的具体例子进行说明。

对自相关矩阵R进行奇异值分解。自相关矩阵R的要素是对称的。因而，当对自相关矩阵R进行奇异值分解时，如以下的(52)式那样，成为正交矩阵U、对角矩阵Σ、以及正交矩阵U的转置矩阵之积。

[数式25]

R＝U∑U

如以下的(53)式所示，(52)式的对角矩阵Σ是对角分量为自相关矩阵R的特征值的矩阵。将对角矩阵Σ的对角分量设为σ

[数式26]

根据自相关矩阵R的奇异值分解的结果得到的对角矩阵Σ的对角分量即σ

[数式27]

(54)式中的矩阵U

如果用矩阵要素来表现矩阵Σ

[数式28]

通过代替自相关矩阵R而使用矩阵R’，由此将(50)式的关系式改写为以下的(56)式那样。

[数式29]

通过对(56)式进行变形，由此作为求出系数α的公式，得到以下的(57)式。使用通过(57)式求出的系数α，根据(45)式来计算时刻k的物理量的预测值y^

[数式30]

此处，以将对角矩阵Σ的对角分量即σ

(57)式是在修正自回归模型的系数决定中利用的方程式。(57)式的矩阵U

通过对(57)式的右边进行计算，由此求出修正自回归模型的系数α。以上，对修正自回归模型的系数α的导出方法的一例进行了说明。此处，为了便于直观地理解，将成为修正自回归模型的基础的自回归模型的系数的导出方法，设为对于时刻k的物理量的预测值y^

例如，能够根据自相关矩阵R的特征值的分布，来决定从(53)式所示的自相关矩阵R中提取的特征值的数量s。

此处，上述修正自回归模型的说明中的物理量为前后方向力。前后方向力的值根据铁道车辆的状态等而变动。

因此，首先，使铁道车辆在轨道16上行驶，而得到与前后方向力的测定值相关的数据y。针对所得到的每个数据y，使用(49)式和(51)式来求出自相关矩阵R。通过对该自相关矩阵R进行由(52)式表示的奇异值分解，由此求出自相关矩阵R的特征值。图8是表示自相关矩阵R的特征值的分布的一例的图。在图8中，将对与轮轴13a的前后方向力T

在图8所示的例子中，具有比其他显著高的值的特征值存在一个。此外，虽然达不到上述具有显著高的值的特征值的程度，但与其他相比具有比较大的值且不被视为0(零)的特征值存在两个。根据该情况，作为从(53)式所示的自相关矩阵R提取的特征值的数量s，例如能够采用2或者3。无论采用哪个，结果都不会产生显著的差异。另外，根据铁道车辆的构成、轨道的构成等，具有比其他显著高的值的特征值的数量可以改变。因而，从自相关矩阵R提取的特征值的数量s只要为1个以上，则并不限定于这些值。

每当由数据取得部502以规定的取样周期取得前后方向力的测定值y的时间序列数据的时刻k的值y

首先，第1频率调整部503基于前后方向力的测定值y的时间序列数据、以及预先设定的数M、m，使用(49)式和(51)式来生成自相关矩阵R。

接着，第1频率调整部503对自相关矩阵R进行奇异值分解，由此导出(52)式的正交矩阵U以及对角矩阵Σ，并根据对角矩阵Σ导出自相关矩阵R的特征值σ

接着，第1频率调整部503将自相关矩阵R的多个特征值σ

接着，第1频率调整部503基于前后方向力的测定值y的时间序列数据、特征值σ

然后，第1频率调整部503基于修正自回归模型的系数α以及前后方向力的测定值y的时间序列数据，通过(45)式来导出前后方向力的测定值y的时间序列数据的时刻k的预测值y^

图9是表示前后方向力的测定值的时间序列数据(测定值)以及前后方向力的预测值的时间序列数据(计算值)的一例的图。另外，在本实施方式中，得到4个前后方向力T

在图9中，轮轴13a的前后方向力T

因而，在图9中，如果从轮轴13a～13d的前后方向力T

因此，第1频率调整部503从前后方向力的测定值y

图10是表示前后方向力的高频分量的时间序列数据的一例的图。图10的纵轴表示前后方向力T

第1频率调整部503如以上那样导出前后方向力T

[状态变量导出部504，S703]

状态变量导出部504将观测方程式设为由存储部501存储的观测方程式，将状态方程式设为由存储部501存储的状态方程式，通过卡尔曼滤波器来决定(44)式所示的状态变量的推定值。此时，状态变量导出部504使用由第1频率调整部503生成的前后方向力T

卡尔曼滤波器是进行数据同化的方法之一。即，卡尔曼滤波器是以使能够观测的变量(观测变量)的测定值与推测值之间的差异变小(最小)的方式、决定未观测的变量(状态变量)的推测值的方法的一例。状态变量导出部504求出观测变量的测定值与推测值之间的差异变小(最小)的卡尔曼增益，并求出此时的未观测的变量(状态变量)的推测值。在卡尔曼滤波器中，使用以下的(58)式的观测方程式、以及以下的(59)式的状态方程式。

Y＝HX+V …(58)

X·＝ΦX+W …(59)

在(58)式中，Y是保存观测变量的测定值的向量。H是观测模型。X是保存状态变量的向量。V是观测噪声。在(59)式中，X·表示X的时间微分。Φ是线性模型。W是系统噪声。另外，卡尔曼滤波器本身能够通过公知的技术来实现，因此省略其详细说明。

在专利文献1所记载的技术中，将测定值(车身11的左右方向上的加速度的测定值、转向架12a、12b的左右方向上的加速度的测定值以及轮轴13a～13d的左右方向上的加速度的测定值)，直接用作为被作为观测变量的测定值而赋予的值。与此相对，在本实施方式中，如在[第1构思]的部分中说明的那样，在数据同化时，作为本来被作为观测变量的测定值而赋予的值，不赋予测定值，而赋予预先设定的恒定值。在本实施方式中，设为加速度的时间序列数据的平均值为0(零)，而将被作为观测变量赋予的恒定值全部设为0(零)。因而，在本实施方式中，状态变量导出部504在进行数据同化时，以观测变量的推定值相对于恒定值(此处为0(零))的误差成为最小或者该误差的期待值成为最小的方式导出状态变量的推定值。

状态变量导出部504通过以规定的取样周期决定(44)式所示的状态变量的推定值，由此生成(44)式所示的状态变量的推定值的时间序列数据。

[第2频率调整部505、S704]

当通过第1频率调整部503未充分除去前后方向力的测定值的时间序列数据中所包含的低频分量的信号强度时，在由状态变量导出部504生成的状态变量的推定值的时间序列数据中有可能残留由于铁道车辆在曲线轨道上行驶而引起的低频分量的信号。因此，第2频率调整部505降低(优选除去)由状态变量导出部504生成的状态变量的推定值的时间序列数据中所包含的低频分量的信号强度。另外，在能够以通过第1频率调整部503充分除去前后方向力的测定值的时间序列数据中所包含的低频分量的信号强度的方式、确定从(53)式所示的自相关矩阵R中提取的特征值的数量s的情况下，无需第2频率调整部505的处理。

在本实施方式中，第2频率调整部505与第1频率调整部503相同，使用修正自回归模型来降低状态变量的推定值的时间序列数据中所包含的低频分量的信号强度。

第2频率调整部505以规定的取样周期对每个状态变量进行以下的处理。

此处，上述修正自回归模型的说明中的物理量成为状态变量。即，状态变量的数据y成为由状态变量导出部504生成的状态变量的推定值的时间序列数据。状态变量的推定值均根据铁道车辆的状态而变动。

首先，第2频率调整部505基于状态变量的推定值的数据y、预先设定的数M、m，使用(49)式和(51)式来生成自相关矩阵R。

接着，第2频率调整部505通过对自相关矩阵R进行奇异值分解，由此导出(52)式的正交矩阵U以及对角矩阵Σ，并根据对角矩阵Σ导出自相关矩阵R的特征值σ

接着，第2频率调整部505将自相关矩阵R的多个特征值σ

接着，第2频率调整部505基于状态变量的推测值的数据y、特征值σ

然后，第2频率调整部505基于修正自回归模型的系数α以及状态变量的推测值的数据y，通过(45)式来导出状态变量的推测值的数据y的时刻k的预测值y^

然后，第2频率调整部505从状态变量的推测值的数据y中减去状态变量的预测值y^

[轨道状态导出部506、S705]

当向(5)式～(8)式的描述轮轴13a～13d的偏转的运动方程式中代入(22)式～(25)式时，得到以下的(60)式～(63)式。

[数式31]

在本实施方式中，如(60)式～(63)式所示那样确定表示前后方向力T

轨道状态导出部506根据(30)式～(33)式来计算轮轴13a～13d的偏转方向上的转动量(角位移)ψ

然后，轨道状态导出部506根据通端不规则量y

轨道状态导出部506计算出使相位一致后的通端不规则量y

另外，轨道状态导出部506也可以对于使相位一致后的通端不规则量y

此外，轨道状态导出部506也可以计算使相位一致后的通端不规则量y

检查装置500使用铁道车辆在通端不规则量的导出对象的行驶区间行驶的期间由数据取得部502取得的各取样时刻的前后方向力的测定值的时间序列数据，执行第1频率调整部503、状态变量导出部504、第2频率调整部505以及轨道状态导出部506的处理。

如此，轨道状态导出部506能够得到铁道车辆在通端不规则量的导出对象的行驶区间行驶的期间的各取样时刻的通端不规则量y

另外，轨道状态导出部506不一定需要如上述那样计算出各取样时刻的铁道车辆的行驶位置。例如，轨道状态导出部506也可以使用GPS(Global Positioning System：全球定位系统)来求出各取样时刻的铁道车辆的行驶位置。

[输出部507、S706]

输出部507输出由轨道状态导出部506计算出的最终的通端不规则量y

[总结]

如以上那样，在本实施方式中，检查装置500将前后方向力T

此外，在本实施方式中，检查装置500根据前后方向力的测定值y的时间序列数据生成自相关矩阵R，并使用对自相关矩阵R进行奇异值分解而得到的特征值中的、从最大的特征值起的s个特征值，来决定对前后方向力的测定值y的时间序列数据进行近似的修正自回归模型的系数α。因而，能够以前后方向力的测定值y的时间序列数据中所包含的低频分量的信号残留、高频分量不残留的方式决定系数α。检查装置500通过将比时刻k靠前的时刻k-l(1≤l≤m)的前后方向力的测定值y的时间序列数据赋予到如此确定了系数α的修正自回归模型中，由此计算出该时刻k的前后方向力的预测值y^

[变形例]

在本实施方式中，在数据同化时，作为原本被作为观测变量的测定值而赋予的值，赋予了预先设定的恒定值。该恒定值并不限定于0(零)。例如，也可以得到搭载有检查装置500的铁道车辆或者与该铁道车辆等同的铁道车辆(构造与该铁道车辆相同的铁道车辆)在通端不规则量y

在本实施方式中，以使用修正自回归模型的情况为例进行了说明。但是，并不一定需要使用修正自回归模型来从前后方向力的测定值y的时间序列数据中降低由于铁道车辆在曲线轨道上行驶而产生的低频分量的信号。例如，在能够确定由于铁道车辆在曲线轨道上行驶而产生的频带的情况下，也可以使用高通滤波器来从前后方向力的测定值y的时间序列数据中降低由于铁道车辆在曲线轨道上行驶而产生的低频分量的信号。此外，在铁道车辆所行驶的轨道是(曲率为0(零)的理想的)直线轨道或者虽然在设计上是直线轨道但具有不对通端不规则量的推定精度造成影响的程度的曲率的轨道的情况下，无需第1频率调整部503以及第2频率调整部505的处理。

此外，在本实施方式中，以成为使相位一致时的基准的轮轴是轮轴13a的情况为例进行了说明。但是，成为基准的轮轴也可以是轮轴13a以外的轮轴13b、13c或者13d。

此外，在本实施方式中，以使用卡尔曼滤波器的情况为例进行了说明。但是，只要使用以观测变量的推定值相对于恒定值的误差成为最小或者该误差的期待值成为最小的方式导出状态变量的推定值的滤波器(即，进行数据同化的滤波器)即可，不一定需要使用卡尔曼滤波器。例如，也可以使用粒子滤波器。另外，作为观测变量的推定值相对于恒定值的误差，例如可列举观测变量的推定值与恒定值之间的平方误差。

此外，在本实施方式中，以导出通端不规则量的情况为例进行了说明。但是，作为反映轨道16的状态的信息，只要导出反映轨道不规则(轨道16外观上的不良)的信息即可，不一定需要导出通端不规则量。例如，也可以通过在通端不规则量的基础上或者代替通端不规则量而进行以下的(64)式～(67)式的计算，由此导出铁道车辆在直线轨道上行驶时产生的横向压力(车轮与导轨之间的左右方向的应力)。其中，Q

[数式32]

此外，在本实施方式中，以包含表示车身11的状态的状态变量的情况为例进行了说明。但是，车身11是由于车轮14a～14d与轨道16之间的作用力(蠕变力)而产生的振动最后传播到的部分。因而，例如，在判断为在车身11中由于该传播而产生的影响较小的情况下，也可以不包含表示车身11的状态的状态变量。在这种情况下，不需要(1)式～(21)式的运动方程式中的(15)式～(17)式(描述车身11的横向振动、偏转、横摇的运动方程式)、以及(18)式、(19)式(描述配置于转向架12a的偏转阻尼器、配置于转向架12b的偏转阻尼器的偏转的运动方程式)。此外，在(1)式～(21)式的运动方程式中，将与车身相关的状态量(包含下标b的状态量)、以及包含与车身相关的状态量(包含下标b的状态量)的{}内的值(例如(21)式的左边第3项的

此外，在本实施方式中，以转向架12a、12b为无轴梁转向架的情况为例进行了说明。但是，转向架12a、12b并不限定于无轴梁转向架。此外，根据铁道车辆的构成要素、铁道车辆所承受的力、以及铁道车辆的运动方向等，能够适当改写运动方程式。即，运动方程式并不限定于本实施方式中例示的运动方程式。在运动方程式中表示铁道车辆承受不依赖于状态变量的外力的情况下，在状态方程式中包含表示该外力的项。

(第2实施方式)

接着，对第2实施方式进行说明。在第1实施方式中，以如下情况为例进行了说明：在数据同化时将原本被作为观测变量的测定值而赋予的值(车身11的左右方向上的加速度、转向架12a、12b的左右方向上的加速度以及轮轴13a～13d的左右方向上的加速度)设为恒定值(0(零))，并使用进行数据同化的滤波器(卡尔曼滤波器)来导出状态变量。与此相对，在本实施方式中，对不进行数据同化地导出状态变量的情况进行说明。如此，本实施方式与第1实施方式主要的不同点在于导出状态变量的方法(状态变量导出部504所具有的功能)。因而，在本实施方式的说明中，对于与第1实施方式相同的部分赋予与在图1～图10中赋予的符号相同的符号等而省略详细的说明。

在本实施方式中，存储部501不存储状态方程式((58)式)以及观测方程式((59)式)，而存储以下的(68)式的运动方程式。

X·＝cΦX …(68)

(68)式是将使用(44)式所示的状态变量表现了(9)式、(10)式、(13)式～(21)式、(34)式～(39)式的运动方程式而得到的公式(在(68)式中将c设为1的公式)，变更为与该公式相比状态变量的时间变化变小的公式的一例。具体而言，(68)式是在使用(44)式所示的状态变量表现了(9)式、(10)式、(13)式～(21)式、(34)式～(39)式的运动方程式而得到的公式中，对与状态变量的一阶时间微分(X·)的项用等号连结的项乘以遗忘系数(forgetting factor)c而得到的公式。即，(68)式是通过在(59)式的状态方程式中导入遗忘系数c而使系统噪声W为0(零)的公式。

遗忘系数c是预先确定的值，(理论上)是超过0且为1以下的值(0

但是，遗忘系数c必须被选择为，使求解(68)式而得到的状态变量的推定值(解)不发散。只要求解(68)式而得到的状态变量的推定值(解)不发散，则遗忘系数c的值为1.0时的该状态变量的推定值(解)就成为精度最高的解。但是，当遗忘系数c的值为1.0时，求解(68)式而得到的状态变量的推定值(解)发散(求不出解)的可能性较高。

从这样的观点出发，也可以使遗忘系数c的上限值小于1.0而进行遗忘系数c的选择。即，遗忘系数c的值例如能够从超过0.0且小于1.0(0.0

另外，在遗忘系数c的值为1.0的情况下，(68)式成为(9)式、(10)式、(13)式～(21)式、(34)式～(39)式的运动方程式本身(仅是使用状态变量来表现了该运动方程式)。

状态变量导出部504使用由第1频率调整部503生成的前后方向力T

如以上那样，在本实施方式中，检查装置500向在将系统噪声W设为0(零)的状态方程式中对状态变量的时间微分项X·以外的项乘以遗忘系数c而得到的方程式，赋予前后方向力T

在本实施方式中也能够采用在第1实施方式中说明过的各种变形例。当在运动方程式中含有不依赖于状态方程式X的外力等的情况下，(68)式如以下的(69)式那样表示。

X·＝c(ΦX+Gf) …(69)

G是在运动方程式中保存不依赖于状态方程式的项的向量。F是与向量G对应的矩阵。

(第3实施方式)

接着，对第3实施方式进行说明。

在第1、第2实施方式中，以搭载于铁道车辆的检查装置500对最终的通端不规则量y

图11是表示检查系统的构成的一例的图。在图11中，检查系统具有数据收集装置1110a、1110b以及数据处理装置1120。图11还表示数据收集装置1110a、1110b以及数据处理装置1120的功能性构成的一例。另外，数据收集装置1110a、1110b以及数据处理装置1120的硬件例如能够通过图6所示的硬件来实现。因而，省略数据收集装置1110a、1110b以及数据处理装置1120的硬件构成的详细的说明。

在铁道车辆的各自中各搭载有1个数据收集装置1110a、1110b。数据处理装置1120配置在指挥中心。指挥中心例如对多个铁道车辆的运行进行集中管理。

[数据收集装置1110a、1110b]

数据收集装置1110a、1110b能够通过相同的装置来实现。数据收集装置1110a、1110b具有数据取得部1111a、1111b以及数据发送部1112a、1112b。

<数据取得部1111a、1111b>

数据取得部1111a、1111b具有与数据取得部502相同的功能。即，数据取得部1111a、1111b与数据取得部502同样地取得前后方向力的测定值的时间序列数据。用于得到前后方向力的测定值的构成与在第1实施方式中说明过的构成相同。

[数据发送部1112a、1112b]

数据发送部1112a、1112b将由数据取得部1111a、1111b取得的前后方向力的测定值的时间序列数据发送至数据处理装置1120。在本实施方式中，数据发送部1112a、1112b将由数据取得部1111a、1111b取得的前后方向力的测定值的时间序列数据通过无线通信发送至数据处理装置1120。此时，数据发送部1112a、1112b将搭载有数据收集装置1110a、1110b的铁道车辆的识别编号附加到由数据取得部1111a、1111b取得的前后方向力的测定值的时间序列数据中。如此，数据发送部1112a、1112b发送附加了铁道车辆的识别编号的前后方向力的测定值的时间序列数据。

<数据处理装置1120>

[数据接收部1121]

数据接收部1121接收由数据发送部1112a、1112b发送的前后方向力的测定值的时间序列数据。在该前后方向力的测定值的时间序列数据中附加有该前后方向力的测定值的时间序列数据的发送源即铁道车辆的识别编号。

[数据存储部1122]

数据存储部1122存储由数据接收部1121接收到的前后方向力的测定值的时间序列数据。数据存储部1122针对每个铁道车辆的识别编号来存储前后方向力的测定值的时间序列数据。数据存储部1122基于铁道车辆当前的运行状況、以及前后方向力的测定值的时间序列数据的接收时刻，确定该前后方向力的测定值的时间序列数据的接收时刻的铁道车辆的位置，并将确定出的位置的信息与该前后方向力的测定值的时间序列数据相互建立关联地存储。另外，数据收集装置1110a、1110b也可以收集铁道车辆当前位置的信息，并将所收集到的信息包含在前后方向力的测定值的时间序列数据中。

[数据读出部1123]

数据读出部1123读出由数据存储部1122存储的前后方向力的测定值的时间序列数据。数据读出部1123能够读出由数据存储部1122存储的前后方向力的测定值的时间序列数据中、由操作人员指定的数据。此外，数据读出部1123也能够在预先确定的时刻读出符合预先确定的条件的前后方向力的测定值的时间序列数据。在本实施方式中，例如基于铁道车辆的识别编号以及位置中的至少任一个来决定由数据读出部1123读出的前后方向力的测定值的时间序列数据。

存储部501、第1频率调整部503、状态变量导出部504、第2频率调整部505、轨道状态导出部506以及输出部507与在第1实施方式中说明过的相同。因而，在此处省略其详细说明。另外，第1频率调整部503代替由数据取得部502取得的前后方向力的测定值的时间序列数据而使用由数据读出部1123读出的前后方向力的测定值的时间序列数据，来决定前后方向力T

<总结>

如以上那样，在本实施方式中，搭载于铁道车辆的数据收集装置1110a、1110b收集前后方向力的测定值的时间序列数据而发送至数据处理装置1120。配置于指挥中心的数据处理装置1120存储从数据收集装置1110a、1110b接收到的前后方向力的测定值的时间序列数据，并使用所存储的前后方向力的测定值的时间序列数据来计算最终的通端不规则量y

<变形例>

在本实施方式中，以从数据收集装置1110a、1110b向数据处理装置1120直接发送计测数据的情况为例进行了说明。但是，不一定需要如此进行。例如，也可以利用云计算来构建检查系统。

此外，在本实施方式中，也能够采用在第1实施方式以及第2实施方式中说明过的各种变形例。

此外，在第1、第2实施方式中，以存储部501、数据取得部502、第1频率调整部503、状态变量导出部504、第2频率调整部505、轨道状态导出部506以及输出部507包含于一个装置的情况为例进行了说明。但是，不一定需要如此。也可以通过多个装置来实现存储部501、数据取得部502、第1频率调整部503、状态变量导出部504、第2频率调整部505、轨道状态导出部506以及输出部507的功能。在该情况下，使用该多个装置来构成检查系统。

(计算例)

接着，对计算例进行说明。在本计算例中，通过第1实施方式的方法导出最终的通端不规则量y

此外，对于第1实施方式的方法，作为被作为观测变量的测定值而赋予的值(车身11、转向架12a、12b以及轮轴13a～13d的左右方向上的加速度的测定值)，不赋予预先设定的恒定值，而直接赋予测定值(即，专利文献1所记载的方法)，由此导出最终的通端不规则量y

图12是表示本计算例，且表示通端不规则量的导出对象的轨道16的曲率1/R以及铁道车辆的行驶速度v的图。在图12中，曲线1201表示铁道车辆的行驶速度，曲线1202表示轨道16的曲率1/R。另外，图12的横轴是在将基准的时刻设为0(零)的情况下从该基准的时刻起的经过时间(秒)。

图13A以及图13B是表示本计算例，且表示自相关矩阵R的特征值的分布的图。图13A表示针对轮轴13a的前后方向力T

图14是表示本计算例，且表示前后方向力T

图15是表示本计算例，且表示前后方向力T

图16A以及图16B是表示使用图15所示的前后方向力T

图17A是表示使用图15所示的前后方向力T

当将图16A的计算值与图16B的计算值进行比较时，可知通过第1实施方式的方法导出的通端不规则量y

(其他实施方式)

以上说明的本发明的实施方式能够通过计算机执行程序来实现。此外，记录有上述程序的计算机可读取的记录介质以及上述程序等计算机程序产品也能够作为本发明的实施方式加以应用。作为记录介质，例如能够使用软盘、硬盘、光盘、光磁盘、CD-ROM、磁带、非易失性存储卡、ROM等。

此外，以上说明的本发明的实施方式都只是实施本发明时的具体化的例子，不能通过这些来限定性地解释本发明的技术范围。即，本发明在不脱离其技术思想或者主要特征的情况下，能够以各种方式实施。

另外，能够将专利文献1的说明书以及附图的内容全部援用于此。

产业上的可利用性

本发明例如能够用于检查铁道车辆。