车辆感测方法与流程

本发明涉及车辆感测方法。

背景技术：

以往，例如如专利文献1记载，公开有一种系统，当车辆通过振动传感器附近时，感测来自前轮的影响带来的振动，接着感测来自后轮的振动，并且与来自离开一定距离设置的其它振动传感器的信息合并，从而计算行驶速度，对前轮和后轮通过时的振动进行计数从而估计车辆的通过数量。

专利文献1：日本特开2004-252520号公报

但是，根据专利文献1记载的处理，由于噪声或谐振等扰动影响，有时难以正确地感测车轴，所以存在出现车辆的感测遗漏或误检测的问题。

技术实现要素：

本车辆感测方法的特征在于，包括：获取检测结构物的振动的第一传感器的第一输出和检测所述结构物中与所述第一传感器不同的部位的振动的第二传感器的第二输出的工序；检测所述第一输出的多个峰值和所述第一输出的多个峰值的检测时刻的工序；检测所述第二输出的多个峰值和所述第二输出的多个峰值的检测时刻的工序；从所述第一输出的多个峰值中提取在所述结构物上移动的一台车辆的量的第一多个峰值的工序；从所述第二输出的多个峰值中提取在所述结构物上移动的所述一台车辆的量的第二多个峰值的工序；在将所述第一多个峰值的数量设为n，将所述第二多个峰值的数量设为m时，在所述n大于所述m的情况下，根据所述第一多个峰值生成种的峰值的第一组合，对于所述第一组合中的各组合，计算在所述检测时刻顺序彼此之间所述第一多个峰值的检测时刻与所述第二多个峰值的检测时刻的第一差量，对于所述第一组合中的各组合，计算所述第一差量的第一波动指标，从所述第一组合中选择所述第一波动指标最小的峰值的组合的工序；以及在所述n比所述m小的情况下，根据所述第二多个峰值生成种的峰值的第二组合，对于所述第二组合中的各组合，计算在所述检测时刻顺序彼此之间所述第一多个峰值的检测时刻与所述第二多个峰值的检测时刻的第二差量，对于所述第二组合中的各组合，计算所述第二差量的第二波动指标，从所述第二组合中选择所述第二波动指标最小的峰值的组合的工序。

优选地，在上述的车辆感测方法中，包括：计算所述第一多个峰值以及所述第二多个峰值的能量大小的工序；在所述n比所述m大的情况下，对于所述第一组合中的各组合，计算在所述检测时刻顺序彼此之间所述第一多个峰值的能量大小与所述第二多个峰值的能量大小的第三差量的工序；对于所述第一组合中的各组合，计算所述第三差量的第三波动指标的工序；在所述n比所述m小的情况下，对于所述第二组合中的各组合，计算在所述检测时刻顺序彼此之间所述第一多个峰值的能量大小与所述第二多个峰值的能量大小的第四差量的工序；以及对于所述第二组合中的各组合，计算所述第四差量的第四波动指标的工序。

优选地，在上述的车辆感测方法中，所述车辆具有三个以上的车轴，提取所述第一多个峰值的工序以及提取所述第二多个峰值的工序是提取所述车轴中的多个车轴通过所述第一传感器或者所述第二传感器所需要的时间中在最短时间以上且最长时间以下的时间内检测到的所述第一多个峰值以及所述第二多个峰值的工序。

优选地，在上述的车辆感测方法中，所述第一传感器以及所述第二传感器是加速度传感器。

附图说明

图1是说明车辆感测系统的图。

图2是实施方式一涉及的从车轴检测到车辆感测为止的流程图。

图3是说明车辆通过时的传感器输出例子的图。

图4是根据传感器输出说明车轴通过时刻的图。

图5是说明车辆通过时的车轴通过时刻的图。

图6是说明多个峰值的组合的生成的图。

图7是说明时刻差值的导出及其方差值的图。

图8是实施方式二涉及的从车轴检测到车辆感测为止的流程图。

图9是说明车辆通过时的能量的计算的图。

图10是根据传感器输出说明能量的图。

图11是说明多个峰值的组合的图。

图12是说明能量差值的导出和方差值的图。

附图标记说明

100…车辆感测系统；1…测量装置；2a…作为第一传感器的第一加速度传感器；2b…作为第二传感器的第二加速度传感器；3…结构物；4…车辆；5、6、7、8…车轴；l1…区间；p1…最大峰值；p5a、p6a、p7a、p8a、p5b、p6b、p7b、p8b…峰值；t5a、t6a、t7a、t8a、t5b、t6b、t7b、t8b…检测时刻；w1a…作为第一输出的第一检测数据；w1b…作为第二输出的第二检测数据。

具体实施方式

下面，参照附图说明本发明的实施方式。需要说明的是，为了将各层和各部件示出为可识别的程度的大小，下面的附图中的各层和各部件的尺寸与实际不同。

实施方式一

车辆感测系统

首先，参照图1说明实施方式一涉及的车辆感测系统100的简要构成。需要说明的是，在本实施方式中，例举加速度传感器作为第一传感器以及第二传感器的理想的一例进行说明，但是除了加速度传感器之外，还可以应用电磁式振动传感器、光学振动传感器等。

图1是说明车辆感测系统的简要图。

如图1所示，车辆感测系统100具有测量装置1、作为第一传感器的第一加速度传感器2a以及作为第二传感器的第二加速度传感器2b。测量装置1与第一加速度传感器2a以及测量装置1与第二加速度传感器2b通过有线或者无线连接，能够进行高速的数据通信。在测量装置1中，接收作为来自第一加速度传感器2a的第一输出的第一检测数据以及作为来自第二加速度传感器2b的第二输出的第二检测数据，按照后述的感测方法，处理各输出数据，从而能够感测在桥梁等的结构物3上移动或者通过的车辆4。

第一加速度传感器2a和第二加速度传感器2b隔开一定的距离安装在结构物3。换言之，第二加速度传感器2b设置在与设置有第一加速度传感器2a的部位不同的部位。在行驶的车辆4进入结构物3上的一侧配置有第一加速度传感器2a，在车辆4退出结构物3上的一侧配置有第二加速度传感器2b。在结构物3上行驶的车辆4具有四个车轴5、6、7、8，每当车辆4的车轴5、6、7、8通过第一加速度传感器2a以及第二加速度传感器2b上时发生振动。第一加速度传感器2a以及第二加速度传感器2b感测各车轴5、6、7、8的振动，比较第一加速度传感器2a的第一检测数据和第二加速度传感器2b的第二检测数据，从而能够感测在结构物3上行驶的车辆4。需要说明的是，在本实施方式中，例举车辆4的车轴为四个的情况进行说明，但是，并不限定于此，还可以是车轴为三个以上。需要说明的是，在本实施方式中，最短时间表示最小车轴通过时间，最长时间表示最大车轴通过时间。

车辆感测方法

其次，对于实施方式一涉及的车辆感测方法，参照图2～图7进行说明。

图2是实施方式一涉及的到车辆感测为止的流程图。图3是说明车辆通过时的传感器输出例子的图。图4是根据传感器输出说明车轴通过时刻的图。图5是说明车辆通过例子的车轴通过时刻的图。图6是说明多个峰值的组合的生成的图。图7是说明时刻差值的导出及其方差值的图。

如图2所示，首先，在步骤s1中，实施方式一涉及的车辆感测方法通过测量装置1接收在结构物3上行驶的车辆4通过第一加速度传感器2a上时检测到的作为第一输出的第一检测数据和相同地通过第二加速度传感器2b上时检测到的作为第二输出的第二检测数据，从而获取第一检测数据和第二检测数据。

需要说明的是，如图3所示，第一检测数据w1a以及第二检测数据w1b包括车辆4通过第一加速度传感器2a以及第二加速度传感器2b上时的与四个车轴5、6、7、8对应的峰值和检测到该峰值时的检测时刻。

此外，为了去除噪声或多余的振动，对于获取到的第一检测数据w1a和第二检测数据w1b应用滤波处理。作为例子，可以例举40hz～100hz的带通滤波处理。

其次，在步骤s2中，根据第一加速度传感器2a检测到的第一检测数据w1a，检测通过第一加速度传感器2a上的车轴5、6、7、8的振动峰值和那时的检测时刻。在图3示出的例子中，检测超过检测数据w1的阈值th1和阈值th2的大小的对应于车轴5、6、7、8的多个检测峰值p5a、p6a、p7a、p8a以及那时的检测时刻t5a、t6a、t7a、t8a。此外，作为随着车轴通过带来的振动的峰值的检测方法，如图4所示，将在超过阈值th1和阈值th2的大小的区间l1中传感器输出的大小最大峰值p1作为检测峰值，并将检测到该峰值p1的时刻t1作为检测时刻。

其次，在步骤s3中，根据第二加速度传感器2b检测到的第二检测数据w1b，检测通过第二加速度传感器2b上的车轴5、6、7、8的振动的峰值和那时的检测时刻。在图3示出的例子中，检测除了未超过阈值th1和阈值th2的大小的车轴7之外的超过阈值th1和阈值th2的大小的对应于车轴5、6、8的多个检测峰值p5b、p6b、p8b和那时的检测时刻t5b、t6b、t8b。作为随着车轴通过带来的振动的峰值的检测方法，与根据第一加速度传感器2a的检测数据检测峰值的情况相同。

其次，在步骤s4中，从在步骤s2中检测到的第一检测数据w1a的多个峰值及其检测时刻中，提取一台车辆的量的第一检测数据w1a的多个峰值及其检测时刻。在这里，将在步骤s4中提取的峰值的数量设为n。一台车辆的量的提取方法为：规定车辆通过中的最低速度vmin和最高速度vmax、车辆中的最短车轴长度lmin和最长车轴长度lmax，计算“最小车轴通过时间tmin＝lmax/vmin”和“最大车轴通过时间tmax＝lmin/vmax”。

将x设为正整数时，从某第x个的峰值的检测时刻起到下一个第x+1个峰值的检测时刻为止的时间差δt在最小车轴通过时间≤δt≤最大车轴通过时间的情况下，将第x个的峰值添加至一台量的多个峰值。在时间差δt＜最小车轴通过时间的情况下，将两个车轴视为同一轴，不添加至一台量的多个峰值。在时间差δt＞最大车轴通过时间的情况下，将一台量的多个峰值确定为一台量，并进行提取。在图5示出的例子中，例举了提取第一检测数据w1a中的一台车辆的量的振动的峰值p5a、p6a、p7a、p8a及其检测时刻t5a、t6a、t7a、t8a，提取到的峰值的数量为4。

其次，在步骤s5中，与步骤s4相同地，从在步骤s3中检测到的第二检测数据w1b的多个峰值及其检测时刻中提取一台车辆的量的第二检测数据的多个峰值及其检测时刻。在这里，将在步骤s5中提取的峰值的数量设为m。在图5示出的例子中，例举了提取了第二检测数据w1b中的一台车辆的量的振动的峰值p5b、p6b、p8b及其检测时刻t5b、t6b、t8b，提取到的峰值的数量为3。

其次，在步骤s6中，比较在步骤s4中提取到的峰值的数量n和在步骤s5中提取到的峰值的数量m，在n与m相等的情况下，进入步骤s17的处理。此外，在n＞m的情况下，进入步骤s7的处理。此外，在n＜m的情况下，进入步骤s12的处理。

在n＝m的情况下，在步骤s17中，将第一检测数据w1a的多个峰值和第二检测数据w1b的多个峰值视为同一车辆的峰值，作为一台车辆进行提取。

通过以上处理，结束车辆感测方法的流程。

在n＞m的情况下，在步骤s7中，根据第一检测数据w1a的多个峰值，生成种的第一检测数据w1a的多个峰值的组合、即第一组合，以使峰值数量成为m。在这里，是表示组合(combination)的数量的符号，从不同的a个中提取r个的组合可以以数学式(1)表示。

【数学式1】

需要说明的是，在图5示出的例子中，n＝4、m＝3的情况，所以组合数量是4种。

在图6示出的例子中，是4种组合的多个峰值。生成出4种组合，组合a是p5a、p6a、p7a，组合b是p5a、p6a、p8a，组合c是p5a、p7a、p8a，组合d是p6a、p7a、p8a。

其次，在步骤s8中，从在步骤s7中生成的种的第一检测数据w1a的多个峰值的组合a～d中各选择一组，并且反复执行到步骤s10为止的步骤，直到选择种的组合全部为止。

其次，在步骤s9中，在步骤s8中取出的一组中，根据第一检测数据w1a的多个峰值的检测时刻，对于第二检测数据w1b的多个峰值的检测时刻计算检测时刻顺序彼此之间的第一差量。在图7示出的例子中，对于图6示出的第一检测数据w1a的多个峰值的组合a～d每一个的检测时刻，计算图5示出的一台车辆的量的第二检测数据w1b的多个峰值的检测时刻的差量。以图6示出的组合a～d为例，则在组合a中，针对第一轴计算t5a与t5b的差值，针对第二轴计算t6a与t6b的差值，针对第三轴计算t7a与t8b的差值。在组合b中，针对第一轴计算t5a与t5b的差值，针对第二轴计算t6a与t6b的差值，针对第三轴计算t8a与t8b的差值。在组合c中，针对第一轴计算t5a与t5b的差值，针对第二轴计算t7a与t6b的差值，针对第三轴计算t8a与t8b的差值。在组合d中，针对第一轴计算t6a与t5b的差值，针对第二轴计算t7a与t6b的差值，针对第三轴计算t8a与t8b的差值。

其次，在步骤s10中，计算方差来作为在步骤s9中计算的多个检测时刻的差量的第一波动指标。需要说明的是，各组合的检测时刻的差量的方差可以按照数学式(2)计算。

【数学式2】

在这里，σ²是方差，i是正整数，xi是各差量值，μ是各差量值的平均。

需要说明的是，在本实施方式中作为波动指标的一例计算了方差，但是，只要能求出波动情形则可以是除了方差之外的指标。例如，还可以计算标准偏差。

在图7示出的例子中，计算了第一检测数据w1a的多个峰值的组合a～d的各轴的差量值的方差值。组合a的方差值是0.026247，组合b的方差值是0.016997，组合c的方差值是0.222161，组合d的方差值是0.201163。

其次，在步骤s11中，从在从步骤s7到步骤s10中求出的多个方差值选择方差值最小的组合。在图7示出的例子中，选择方差值为0.016997的最小的组合b作为一台车辆的量的第一检测数据w1a的多个峰值。

其次，在步骤s17中，将第一检测数据w1a的多个峰值和第二检测数据w1b的多个峰值视为同一车辆的峰值，作为1台车辆进行提取。

通过以上处理结束车辆感测方法的流程。

在n＜m的情况下，在步骤s12中，根据第二检测数据w1b的多个峰值，生成种的第二检测数据w1b的多个峰值的组合、即第二组合，以使峰值数量成为n。

其次，在步骤s13中，从在步骤s12中生成的种的第二检测数据w1b的多个峰值的组合各选择一组，并且反复进行到步骤s15为止的步骤，直到选择种的组合全部为止。

其次，在步骤s14中，在步骤s13中取出的一组中，根据第二检测数据w1b的多个峰值的检测时刻，对于第一检测数据w1a的多个峰值的检测时刻计算检测时刻顺序彼此之间的第二差量。

其次，在步骤s15中，作为在步骤s14中计算的多个检测时刻的差量的第二波动指标计算方差。

其次，在步骤s16中，从在从步骤s13到步骤s15中求出的多个方差值选择方差值最小的组合。

其次，在步骤s17中，将第一检测数据w1a的多个峰值和第二检测数据w1b的多个峰值视为同一车辆的峰值，作为一台车辆进行提取。

通过以上处理结束车辆感测方法的流程。

根据如上所述的车辆感测方法，在由于噪声或谐振的扰动影响而使第一多个峰值的数量n与第二多个峰值的数量m不同的情况下，展开为或者种组合，针对峰值的每个组合，计算在检测时刻顺序彼此之间第一多个峰值的检测时刻与第二多个峰值的检测时刻的差量，并计算差量的方差值，从而可以选择认为最适合的组合，能够检测出通过第一加速度传感器2a检测的车辆4和通过第二加速度传感器2b检测的车辆4相同。因此，能够降低车辆4的误检测以及检测遗漏。从而，能够得到高精度的车辆检测效果。

此外，在提取第一多个峰值的工序以及提取第二多个峰值的工序中，根据最小车轴通过时间和最大车轴通过时间的限制，可以提取一台车辆分的第一多个峰值以及第二多个峰值。从而，缩小组合处理的峰值数量的误差，能够得到更高精度的车辆检测效果。

此外，通过利用加速度传感器2a、2b，灵敏度比其它的传感器出色，可以以高速的采样频率获取数据，所以即使在车辆4的行驶速度较快的情况下，也可以感测车轴，能够得到更高精度的车辆检测效果。

实施方式二

车辆感测方法

其次，对于实施方式二涉及的车辆感测方法，参照图8～图12进行说明。图8是根据实施方式二的从车轴检测到车辆感测为止的流程图。图9是说明车辆通过例子的能量的计算的图。图10是根据传感器输出说明能量的图。图11是说明多个峰值的组合的图。图12是说明能量差值的导出和方差值的图。

根据实施方式二的车辆感测方法，如图8所示，首先，在步骤s21中，与实施方式一的步骤s1相同地，在测量装置1接收在结构物3上行驶的车辆4通过第一加速度传感器2a上时检测到的第一检测数据w1a和相同地通过第二加速度传感器2b上时检测到的第二检测数据w1b，从而获取第一检测数据w1a和第二检测数据w1b。

此外，为了去除噪声和多余的振动，对于获取的第一检测数据w1a和第二检测数据w1b应用滤波处理。作为例子，可以例举40hz～100hz的带通滤波处理。

其次，在步骤s22中，根据第一加速度传感器2a检测到的第一检测数据w1a，检测通过第一加速度传感器2a上的车轴5、6、7、8的振动的峰值和那时的检测时刻和那时的能量值。需要说明的是，与实施方式一的步骤s2相同地，车轴通过带来的振动的峰值、检测时刻的检测方法是将超过阈值th1和阈值th2的大小的区间l1中传感器输出的大小最大峰值p1作为检测峰值，将检测到该峰值p1的时刻t1作为检测时刻。此外，如图9所示，车轴通过带来的振动的峰值的能量值的检测方法是将检测数据w1的区间l1的波形数据的绝对值的总和作为能量e1。

其次，在步骤s23中，根据第二加速度传感器2b检测到的第二检测数据w1b，检测通过第二加速度传感器2b上的车轴5、6、7、8的振动的峰值和那时的检测时刻和那时的能量值。与实施方式一的步骤s3相同地，检测除了未超过阈值th1和阈值th2的大小的车轴7的检测峰值p7b和检测时刻t7b之外的超过阈值th1和阈值th2的大小的对应车轴5、6、8的多个检测峰值p5b、p6b、p8b和那时的检测时刻t5b、t6b、t8b和能量e5b、e6b、e8b。随着车轴通过带来的振动的峰值的检测方法与根据第一加速度传感器2a的检测数据检测峰值的情况相同。

其次，在步骤s24中，在步骤s22中检测到的第一检测数据w1a的多个峰值及其检测时刻中，提取一台车辆的量的第一检测数据w1a的多个峰值及其检测时刻和那时的能量值。在这里，将在步骤s24提取的峰值的数量设为n。一台车辆的量的提取方法与实施方式一的步骤s4相同。在图10示出的例子中，例举了提取第一检测数据w1a中的一台车辆的量的振动的峰值p5a、p6a、p7a、p8a和那时的检测时刻t5a、t6a、t7a、t8a和能量e5a、e6a、e7a、e8a，提取的峰值的数量为4。

其次，在步骤s25中，在步骤s23中检测到的第二检测数据w1b的多个峰值及其检测时刻中，与步骤s24相同地，提取一台车辆的量的第二检测数据w1b的多个峰值及其检测时刻和那时的能量值。在这里，将在步骤s25中提取的峰值的数量设为m。在图10示出的例子中，例举了提取第二检测数据w1b中的一台车辆的量的振动的峰值p5b、p6b、p8b和那时的检测时刻t5b、t6b、t8b和能量e5b、e6b、e8b，提取到的峰值的数量为3。

其次，在步骤s26中，比较在步骤s24中提取到的峰值的数量n和在步骤s25中提取到的峰值的数量m，在n与m相等的情况下，进入步骤s37的处理。在n＞m的情况下，进入步骤s27的处理。在n＜m的情况下，进入步骤s32的处理。

在n与m相等的情况下，在步骤s37中，将第一检测数据w1a的多个峰值和第二检测数据w1b的多个峰值视为同一车辆的峰值，作为一台车辆进行提取，结束车辆感测的流程。

在n＞m的情况下，在步骤s27中，根据第一检测数据w1a的多个峰值，生成种的第一检测数据w1a的多个峰值的组合、即第一组合，以使峰值数量变为m。

需要说明的是，在图10示出的例子是n＝4、m＝3的情况，所以组合数量是4种。

在图11示出的例子中表示出4种组合的多个峰值。生成了4种组合，组合a是p5a、p6a、p7a，组合b是p5a、p6a、p8a，组合c是p5a、p7a、p8a，组合d是p6a、p7a、p8a。

其次，在步骤s28中，从在步骤s27中生成的种的第一检测数据w1a的多个峰值的组合各选择一组，并且反复到步骤s30为止的步骤，直到选择种的组合全部。

其次，在步骤s29中，根据作为在步骤s28中取出的组合中的一组的第一检测数据w1a的多个峰值的检测时刻和能量，对于第二检测数据w1b的多个峰值的检测时刻和能量，计算在检测时刻顺序彼此之间各检测时刻的差量和作为各能量的差量的第三差量。各检测时刻的差量的计算方法在实施方式一中已经说明，所以省略说明。在图12示出的例子中，对于图11示出的第一检测数据w1a的多个峰值的组合a～d每一个的能量，计算图10示出的一台车辆的量的第二检测数据w1b的多个峰值的能量的差量。以图12示出的组合a～d为例，则在组合a中，针对第一轴计算e5a与e5b的差值，针对第二轴计算e6a与e6b的差值，针对第三轴计算e7a与e8b的差值。在组合b中，针对第一轴计算e5a与e5b的差值，针对第二轴计算e6a与e6b的差值，针对第三轴计算e8a与e8b的差值。在组合c中，针对第一轴计算e5a与e5b的差值，针对第二轴计算e7a与e6b的差值，针对第三轴计算e8a与e8b的差值。在组合d中，针对第一轴计算e6a与e5b的差值，针对第二轴计算e7a与e6b的差值，针对第三轴计算e8a与e8b的差值。

其次，在步骤s30中，计算在步骤s29计算的多个检测时刻的差量的方差和作为多个能量的差量的第三波动指标的方差。多个检测时刻的差量的方差在实施方式一中已经说明，所以省略说明。在图12示出的例子中，计算第一检测数据w1a的多个峰值的组合a～d的各轴的能量差量值的方差值。组合a的方差值是0.370075，组合b的方差值是0.012472，组合c的方差值是0.111455，组合d的方差值是0.108012。

需要说明的是，在本实施方式中作为波动指标的一例计算了方差，但是，只要能够求出波动情形则可以是除了方差之外的指标。例如，还可以是标准偏差。

其次，在步骤s31中，从在从步骤s27起到步骤s30中求出的检测时刻的多个方差值和能量的多个方差值中选择方差值最小的组合。如图7所示，与实施方式一的步骤s11相同地，检测时刻的方差值的最小的是组合b，在图12示出的例子中，能量的方差值的最小的是方差值为0.012472的组合b，作为一台车辆的量的第一检测数据的多个峰值选择。

在n＜m的情况下，在步骤s32中，根据第二检测数据w1b的多个峰值，生成种的第二检测数据w1b的多个峰值的组合、即第二组合，以使峰值数量变为n。

其次，在步骤s33中，从在步骤s32中生成的种的第二检测数据w1b的多个峰值的组合中各选择一个，并且反复直到步骤s35为止的步骤，直到选择种的组合全部为止。

其次，在步骤s34中，根据作为在步骤s33中取出的组合中的一个的第二检测数据的多个峰值的检测时刻和能量，对于第一检测数据w1a的多个峰值的检测时刻和能量，计算在检测时刻顺序彼此之间各检测时刻的差量和各能量的第四差量。

其次，在步骤s35中，计算在步骤s34计算的多个检测时刻的差量的方差和作为多个能量的差量的第四波动指标的方差。

其次，在步骤s36中，从步骤s33到步骤s35中求出的检测时刻的多个方差值和能量的多个方差值中选择方差值最小的组合。

最后，在步骤s37中，将第一检测数据w1a的多个峰值和第二检测数据w1b的多个峰值视为同一车辆的峰值，作为一台车辆进行提取。

通过以上处理，结束车辆感测方法的流程。

根据如上所述的车辆感测方法，除了实施方式一的检测时刻的组合的选择处理之外，还针对峰值的每个组合，计算在检测时刻顺序彼此之间第一多个峰值的能量与第二多个峰值的能量的差量，通过计算差量的方差值，从而能够选择基于能量的最合适的组合。由此，即使在车辆4的通过速度在第一加速度传感器2a和第二加速度传感器2b发生变化的情况下，也能够得到更高精度的车辆检测效果。

下面，记载从实施方式导出的内容。

车辆感测方法的特征在于，包括：获取检测结构物的振动的第一传感器的第一输出和检测所述结构物中与所述第一传感器不同的部位的振动的第二传感器的第二输出的工序；检测所述第一输出的多个峰值和所述第一输出的多个峰值的检测时刻的工序；检测所述第二输出的多个峰值和所述第二输出的多个峰值的检测时刻的工序；从所述第一输出的多个峰值中提取在所述结构物上移动的一台车辆的量的第一多个峰值的工序；从所述第二输出的多个峰值中提取在所述结构物上移动的所述一台车辆的量的第二多个峰值的工序；在将所述第一多个峰值的数量设为n，将所述第二多个峰值的数量设为m时，在所述n大于所述m的情况下，根据所述第一多个峰值生成种的峰值的第一组合，对于所述第一组合中的各组合，计算在所述检测时刻顺序彼此之间所述第一多个峰值的检测时刻与所述第二多个峰值的检测时刻的第一差量，对于所述第一组合中的各组合，计算所述第一差量的第一波动指标，从所述第一组合中选择所述第一波动指标最小的峰值的组合的工序；以及，在所述n比所述m小的情况下，根据所述第二多个峰值生成种的峰值的第二组合，对于所述第二组合中的各组合，计算在所述检测时刻顺序彼此之间所述第一多个峰值的检测时刻与所述第二多个峰值的检测时刻的第二差量，对于所述第二组合中的各组合，计算所述第二差量的第二波动指标，从所述第二组合中选择所述第二波动指标最小的峰值的组合的工序。

根据该方法，即使在由于噪声和谐振的扰动影响而使第一多个峰值的数量n与第二多个峰值的数量m不同的情况下，展开成或者种的组合，并针对峰值的每个组合，计算在检测时刻顺序彼此之间第一多个峰值的检测时刻与第二多个峰值的检测时刻的差量，并计算差量的方差值，从而可以选择最合适的组合，能够检测出在第一传感器检测的车辆与在第二传感器检测的车辆为相同。因此，能够降低车辆的误检测以及检测遗漏。从而，能够得到高精度的车辆检测效果。

优选地，在上述的车辆感测方法中，包括：计算所述第一多个峰值以及所述第二多个峰值的能量大小的工序；在所述n比所述m小的情况下，对于所述第一组合中的各组合，计算在所述检测时刻顺序彼此之间所述第一多个峰值的能量大小与所述第二多个峰值的能量大小的第三差量的工序；对于所述第一组合中的各组合，计算所述第三差量的第三波动指标的工序；在所述n比所述m小的情况下，对于所述第二组合中的各组合，计算在所述检测时刻顺序彼此之间所述第一多个峰值的能量大小与所述第二多个峰值的能量大小的第四差量的工序；以及，对于所述第二组合中的各组合，计算所述第四差量的第四波动指标的工序。

根据该方法，除了上述的检测时刻的组合的选择处理之外，还针对峰值的每个组合，计算在检测时刻顺序彼此之间第一多个峰值的能量与第二多个峰值的能量的差量，通过计算差量的方差值，从而能够选择基于能量的最合适的组合。由此，即使在移动体的通过速度在第一加速度传感器和第二加速度传感器发生变化的情况下，也能够得到更高精度的车辆检测效果。

优选地，在上述的车辆感测方法中，所述车辆具有三个以上的车轴，提取所述第一多个峰值的工序以及提取所述第二多个峰值的工序是提取所述车轴中的多个车轴通过所述第一传感器或者所述第二传感器所需要的时间中在最短时间以上且最长时间以下的时间内检测到的所述第一多个峰值以及所述第二多个峰值的工序。

根据该方法，在提取第一多个峰值的工序以及提取第二多个峰值的工序中，根据最短时间和最长时间的限制，可以提取一台车辆分的第一多个峰值以及第二多个峰值。从而，缩小组合处理的峰值数量的误差，能够得到更高精度的车辆检测效果。

优选地，在上述的车辆感测方法中，所述第一传感器以及所述第二传感器是加速度传感器。

根据该方法，通过利用加速度传感器，灵敏度比其它的传感器出色，可以以高速的采样频率来获取数据，所以即使在车辆的行驶速度较快的情况下，也可以感测车轴，从而能够得到更高精度的车辆检测效果。