[发明专利]基于傅里叶变换的数据驱动交叉口信号控制实际执行周期估计方法

摘要

一种基于傅里叶变换的数据驱动交叉口信号控制实际执行周期估计方法，包括以下步骤：1)从数据库中获取某个路口某个车道某一整天的过车数据2)对该车道当天一整天的过车数据进行快速傅里叶变换；3)找到频域离散幅频序列Y(f)的波峰所对应的频率；4)找出信号机实际执行周期所对应频率，转化成周期。与现有技术相比，本发明提出了基于快速傅里叶变换的交叉口信号灯周期确定方法，将时域的问题转换到频域以达到准确计算信号灯周期的目的。

主权项

1.一种基于傅里叶变换的数据驱动交叉口信号控制实际执行周期估计方法，其特征在于，所述方法包括以下步骤：1)从数据库中获取某个路口某个车道某一整天的过车数据{pn}，其中，pn为第n辆车通过停车线的时刻，单位为一天中的整数秒，n＝1，2，......N；N为该交叉口该车道该天通过停车线的总车辆数；2)对该车道当天一整天的过车数据进行快速傅里叶变换，过程如下：2.1)将过车数据{pn}转换成以时间为标号的0/1序列y(t)，时间标号t的单位为一天中的整数秒，计算采用如下公式：2.2)利用快速傅里叶变换近似计算时域离散序列y(t)的频域离散幅频序列Y(f)，在进行快速傅里叶变换时，离散序列的长度选择为217＝131072，并且时域离散序列y(t)的采样频率设置为1Hz，因此频域离散幅频序列Y(f)的频域标号f的采样间隔可以精确到0.00001Hz；3)分析频域离散幅频序列Y(f)的波峰，估计一整天中所有周期值其中M为该天的分段配时中使用的不同周期的个数，过程如下：3.1)找到频域离散幅频序列Y(f)的波峰所对应的频率{f′i}，i＝1，2，…，n，n为实际找到的波峰个数，因为在频域离散幅频序列Y(f)中的绝大多数频率f所对应的幅值Y(f)不为0且数值较小，若将这部分频率中的波峰作为之后分析的样本会严重影响实际估计周期的准确性，因此首先要设定一个阈值K，筛选出频域离散幅频序列Y(f)中大于K的那部分离散散幅频序列Y′(f)；经过多次试验，阈值K选取为频域离散幅频序列Y(f)中从大到小排序第300个点所对应的值；筛选后的离散散幅频序列Y′(f)为锯齿波，往往在一个波峰临近位置还有数个波峰，因此寻找波峰时依据该点峰值大于前后相邻5个频率所对应的Y(f)来代表这个区域的波峰，即{f′i}中的频率f满足以下公式：Y(f‑i)＜Y(f)，i＝1，2，3，4，5Y(f)＞Y(f+i)，i＝1，2，3，4，53.2)分析{f′i}，找出信号灯实际执行周期所对应的频率{f″i}，i＝1，2，…，n，n为实际执行的周期个数，因为过车数据序列具有一定的周期特征，在经过傅里叶变换以后，原周期对应频率以及其2倍频率位置都会在频域上出现波峰，利用这个特征对{f′i}中的频率进行筛选；若{f′i}中存在一个频率fa和另一个频率fb，其中fa近似等于fb的2倍，误差±0.5％，则认为fa是实际执行的周期所对应的频率，即fa和fb之间满足以下关系：1.995＜fa/fb＜2.005找到所有满足条件的频率f_a得到实际执行周期所对应的频率{f″_i}，i＝1，2，…，n，将其中的频率取倒数即是估计得到的实际执行的周期