发明名称 一种基于客流量的快速公交到站时刻预测方法
摘要 本发明公开了一种基于客流量的快速公交到站时刻预测方法。实现本发明所述方法的系统由站点RFID采集单元、客流信息视频采集单元、交通信息存储单元、信息处理单元和信息发布单元组成。利用RFID采集单元采集车辆到站时刻信息,以及视频采集系统采集站点客流信息,由信息处理单元计算快速公交车辆停站时间,并结合车辆历次到站信息预测快速公交车辆到达交叉口时刻,根据交叉路口的红绿灯信息判断快速公交车辆是否能够一次通过交叉口,从而计算到达下一个站点时刻,并通过信息发布系统发布。本发明考虑了客流量对快速公交到站时刻的影响,剔除了拥堵等显著影响常规公交而不存在于快速公交的影响因素,提高了预测的精度。
申请公布号 CN103606272B 申请公布日期 2015.12.09
申请号 CN201310618153.8 申请日期 2013.11.27
申请人 北京工业大学 发明人 高超;刘小明;杨孝宽;曹静;张统洋
分类号 G08G1/00(2006.01)I;G08G1/123(2006.01)I 主分类号 G08G1/00(2006.01)I
代理机构 北京思海天达知识产权代理有限公司 11203 代理人 张慧
主权项 一种基于客流量的快速公交到站时刻预测方法,实现该方法的系统包括:站点RFID采集单元,客流信息视频采集单元,交通信息存储单元,信息处理单元和信息发布单元;其特征在于,所述方法包括以下步骤:步骤一,由站点RFID采集单元获得快速公交车辆到站时刻信息,并将数据送到信息处理单元;到站信息用符号T<sub>k</sub>(i)表示,表示编号为k的快速公交车辆到达第i个站点的时刻,i=1,2,……,I,I为站点总数;k=1,2,……,K,K为快速公交车辆当天发车总数,当天发出的第1辆车编号为1,第2辆车编号为2,第k辆车编号为k;步骤二,由客流信息视频采集单元获得快速公交车辆客流信息,并将数据送到信息处理单元;N<sub>k</sub>(i)表示编号为k的快速公交车辆到达第i个站点时该站点的上车乘客数量;K<sub>k</sub>(i)表示编号为k的快速公交车辆到达第i个站点时车内乘客密度,等于车内乘客数量N<sub>k</sub>(i)与车内面积的比值;步骤三,由交通信息存储单元获得绿灯每个周期的起、止时刻信息和快速公交车辆历次到站时刻信息;交通信息存储单元提供交叉口一天里所有的绿灯起始时刻和绿灯结束时刻;快速公交车辆历次到站时刻信息用于计算站点间行程时间;客流信息包括历次车辆在每个停靠站的上下车乘客数;步骤四,预测快速公交车辆到站时刻;假设编号为k的快速公交车辆刚刚通过编号为i的站点,编号为k‑1的快速公交车辆已通过编号为i+1的站点,编号为k和k‑1的快速公交车辆之间的站点数为N<sub>k,k‑1</sub>,N<sub>k,k‑1</sub>等于RFID采集单元获得的编号为k和k‑1的快速公交车刚刚经过的站点数的差;假设编号为i和i+1的站点之间交叉口的编号为依次为j,j+1,……,J,J为最后一个交叉口的编号;(1)计算快速公交车辆进出编号为i的站点所用时间;车辆进出站点所用的时间主要由上车乘客数、下车乘客数、车内乘客数决定,可以近似看作是上车乘客数、车内乘客密度、下车乘客率的一个非线性函数,经非线性拟合得到编号为k的快速公交车辆进出编号为i的站点用时ΔT<sub>k</sub>(i)的表达式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&Delta;T</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&rho;</mi><mi>b</mi></msub><mo>&CenterDot;</mo><msub><mi>N</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&rho;</mi><mi>l</mi></msub><mo>&CenterDot;</mo><msub><mi>K</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>N</mi><mo>&CenterDot;</mo><mover><mrow><mi>P</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow><mo>&OverBar;</mo></mover><mo>+</mo><msub><mi>&rho;</mi><mi>k</mi></msub><mo>&CenterDot;</mo><msub><mi>K</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>+</mo><mi>L</mi><mo>+</mo><mover><msub><mi>&Delta;T</mi><mi>e</mi></msub><mo>&OverBar;</mo></mover><mo>+</mo><mover><msub><mi>&Delta;T</mi><mi>o</mi></msub><mo>&OverBar;</mo></mover></mrow>]]></math><img file="FDA0000732964040000011.GIF" wi="1389" he="91" /></maths>式中,<img file="FDA0000732964040000012.GIF" wi="104" he="81" />为乘客在第i站点下车乘客率;N为车内核定乘客数;<img file="FDA0000732964040000013.GIF" wi="98" he="86" />为从减速到停稳车辆平均用时,<img file="FDA0000732964040000021.GIF" wi="100" he="79" />为从加速到正常行驶车辆平均用时;ρ<sub>b</sub>为上车用时系数,ρ<sub>l</sub>为下车用时系数,ρ<sub>k</sub>为车内乘客密度对下车时间延迟的影响系数;L为由于少数乘客上下车时长较长导致的损失时间;ρ<sub>b</sub>、ρ<sub>l</sub>、ρ<sub>k</sub>、L利用交通信息存储单元中停靠时间与上车乘客数、下车乘客数、车内乘客数的关系进行数据非线性拟合得到;<img file="FDA0000732964040000022.GIF" wi="100" he="84" />通过下述公式进行计算:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mrow><mi>P</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>1</mn></msub></munderover><msub><mi>P</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow><msub><mi>M</mi><mn>1</mn></msub></mfrac></mrow>]]></math><img file="FDA0000732964040000023.GIF" wi="317" he="218" /></maths>式中,P<sub>k</sub>(i)为编号为k的快速公交车辆在i站点的下车乘客数与车内核定乘客数N的比值,下车乘客数由交通信息存储单元提供;M<sub>1</sub>为存储单元中数据样本量;<img file="FDA0000732964040000024.GIF" wi="96" he="83" />与<img file="FDA0000732964040000025.GIF" wi="95" he="83" />通过下述公式进行计算:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mover><msub><mi>&Delta;T</mi><mi>o</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mover><mi>V</mi><mo>&OverBar;</mo></mover><mover><msub><mi>a</mi><mi>o</mi></msub><mo>&OverBar;</mo></mover></mfrac></mrow>]]></math><img file="FDA0000732964040000026.GIF" wi="204" he="157" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mover><msub><mi>&Delta;T</mi><mi>e</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mover><mi>V</mi><mo>&OverBar;</mo></mover><mover><msub><mi>a</mi><mi>e</mi></msub><mo>&OverBar;</mo></mover></mfrac></mrow>]]></math><img file="FDA0000732964040000027.GIF" wi="198" he="155" /></maths>式中,<img file="FDA0000732964040000028.GIF" wi="51" he="72" />为快速公交车辆匀速行驶速度,<img file="FDA0000732964040000029.GIF" wi="66" he="83" />为快速公交车辆平均加速度;<img file="FDA00007329640400000210.GIF" wi="64" he="83" />为快速公交车辆平均减速度;<img file="FDA00007329640400000211.GIF" wi="165" he="83" />和<img file="FDA00007329640400000212.GIF" wi="64" he="83" />3个参数通过雷达测速调查得到;(2)计算快速公交车辆到达编号为J的交叉口的时刻T<sub>ck</sub>(J);编号为k的快速公交车辆进入编号为j的交叉口的时刻T<sub>ck</sub>(j)为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>T</mi><mi>ck</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>T</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&Delta;T</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><mover><msub><mi>&Delta;T</mi><mi>e</mi></msub><mo>&OverBar;</mo></mover><mo>+</mo><mover><msub><mi>&Delta;T</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover></mrow>]]></math><img file="FDA00007329640400000213.GIF" wi="721" he="91" /></maths>式中,<img file="FDA00007329640400000214.GIF" wi="124" he="90" />为从i站点到j交叉口的平均用时,计算公式为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mover><msub><mi>&Delta;T</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>2</mn></msub></munderover><msub><mi>&Delta;T</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><msub><mi>M</mi><mn>2</mn></msub></mfrac></mrow>]]></math><img file="FDA00007329640400000215.GIF" wi="404" he="216" /></maths>式中,M<sub>2</sub>为样本数量,ΔT<sub>k</sub>(i,j)为样本数据中编号为k的快速公交车辆从i站点到j交叉口的用时,由交通信息存储单元提供数据;计算到达编号为j+1的交叉口的时刻T<sub>ck</sub>(j+1):<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>T</mi><mi>ck</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>T</mi><mi>ck</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><mover><msub><mi>&Delta;T</mi><mrow><mi>j</mi><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&OverBar;</mo></mover></mtd><mtd><msub><mi>T</mi><mi>ck</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>&Element;</mo><mo>[</mo><msub><mi>T</mi><mi>bj</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>T</mi><mi>ej</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo></mtd></mtr><mtr><mtd><msub><mi>T</mi><mi>bj</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mover><msub><mi>&Delta;T</mi><mrow><mi>j</mi><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&OverBar;</mo></mover><mo>+</mo><mover><msub><mi>&Delta;T</mi><mi>o</mi></msub><mo>&OverBar;</mo></mover><mo>/</mo><mn>2</mn></mtd><mtd><msub><mi>T</mi><mi>ck</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>&NotElement;</mo><mo>[</mo><msub><mi>T</mi><mi>bj</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>T</mi><mi>ej</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000732964040000031.GIF" wi="1301" he="177" /></maths>式中,T<sub>bj</sub>(t)和T<sub>ej</sub>(t)分别为当时的绿灯起始和终止时刻,由交通信息存储单元提供数据;<img file="FDA0000732964040000032.GIF" wi="147" he="91" />为绿灯信号情况下快速公交车辆到达编号为j的交叉口进口时刻与进入第j+1交叉口时刻的时间差的算术平均值;<img file="FDA0000732964040000033.GIF" wi="152" he="92" />由下述公式进行计算:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mover><msub><mi>&Delta;T</mi><mrow><mi>j</mi><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>3</mn></msub></munderover><msub><mi>&Delta;T</mi><mrow><mi>k</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow><msub><mi>M</mi><mn>3</mn></msub></mfrac></mrow>]]></math><img file="FDA0000732964040000034.GIF" wi="425" he="218" /></maths>式中,M<sub>3</sub>为样本数量,ΔT<sub>k,j,j+1</sub>为绿灯信号情况下快速公交车辆到达编号为j的交叉口进口时刻与进入编号为j+1的交叉口时刻的时间差,由交通信息存储单元提供数据;按照上面的求解方法,依次求出T<sub>ck</sub>(j+1),T<sub>ck</sub>(j+2),……,T<sub>ck</sub>(J);(3)计算快速公交车辆到达编号为i+1的站点的时刻;编号为k的快速公交车辆进入第i+1站点的时刻T<sub>k</sub>(i+1)为:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>T</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>T</mi><mi>ck</mi></msub><mrow><mo>(</mo><mi>J</mi><mo>)</mo></mrow><mo>+</mo><mover><msub><mi>&Delta;T</mi><mrow><mi>J</mi><mo>,</mo><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&OverBar;</mo></mover><mo>+</mo><mover><msub><mi>&Delta;T</mi><mi>e</mi></msub><mo>&OverBar;</mo></mover></mtd><mtd><msub><mi>T</mi><mi>ck</mi></msub><mrow><mo>(</mo><mi>J</mi><mo>)</mo></mrow><mo>&Element;</mo><mo>[</mo><msub><mi>T</mi><mi>bJ</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>T</mi><mi>eJ</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo></mtd></mtr><mtr><mtd><msub><mi>T</mi><mi>bJ</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mover><msub><mi>&Delta;T</mi><mrow><mi>J</mi><mo>,</mo><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&OverBar;</mo></mover><mo>+</mo><mover><msub><mi>&Delta;T</mi><mi>o</mi></msub><mo>&OverBar;</mo></mover><mo>/</mo><mn>2</mn><mo>+</mo><mover><msub><mi>&Delta;T</mi><mi>e</mi></msub><mo>&OverBar;</mo></mover></mtd><mtd><msub><mi>T</mi><mi>ck</mi></msub><mrow><mo>(</mo><mi>J</mi><mo>)</mo></mrow><mo>&NotElement;</mo><mo>[</mo><msub><mi>T</mi><mi>bJ</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>T</mi><mi>eJ</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000732964040000035.GIF" wi="1365" he="181" /></maths>式中,T<sub>bJ</sub>(t)和T<sub>eJ</sub>(t)分别为当时的绿灯起始和终止时刻,由交通信息存储单元提供数据;<img file="FDA0000732964040000036.GIF" wi="140" he="89" />为绿灯信号情况下快速公交车辆到达编号为J的交叉口进口时刻与即将开始减速进入第i+1站点时刻的时间差的算术平均值;<img file="FDA0000732964040000037.GIF" wi="146" he="96" />由下述公式进行计算:<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mover><msub><mi>&Delta;T</mi><mrow><mi>J</mi><mo>,</mo><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>4</mn></msub></munderover><msub><mi>&Delta;T</mi><mrow><mi>k</mi><mo>,</mo><mi>J</mi><mo>,</mo><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow><msub><mi>M</mi><mn>4</mn></msub></mfrac></mrow>]]></math><img file="FDA0000732964040000038.GIF" wi="417" he="215" /></maths>式中,M<sub>4</sub>为样本数量,ΔT<sub>k,J,i+1</sub>为绿灯信号情况下快速公交车辆到达编号为J的交叉口进口时刻与即将开始减速进入第i+1站点时刻的时间差,由交通信息存储单元提供数据;(4)重复(1)~(3),计算编号为k的快速公交车辆到达第i+2,i+3,……,i+N<sub>k,k‑1</sub>站点的时刻;(5)计算快速公交车辆到达第i+n站点预计需要的时间和所需站点数;编号为k的公交车辆到达第i+n站点预计需要的时间为:ΔT=T<sub>k</sub>(i+n)‑t式中,ΔT为编号为k的公交车辆到达第i+n站点预计需要的时间,n=1,2,……,N<sub>k,k‑1</sub>,t为现在时刻;编号为k的快速公交车辆到达i+n站点的所需站数为n‑1;步骤五,将车辆到站信息发布在可变信息板上。
地址 100124 北京市朝阳区平乐园100号
您可能感兴趣的专利