发明名称 一种基于运行车速预测的城市干道车队离散模型的建模方法
摘要 一种基于运行车速预测的城市干道车队离散模型的建模方法,它涉及一种城市干道车队离散模型的建模方法,具体涉及一种基于运行车速预测的城市干道车队离散模型的建模方法。本发明为了解决现有车队离散模型中最低车速、最高车速、平均车速设置依靠工程师经验,无法准确预测下游端面车队到达图式的问题。本发明通过在城市干道上交叉口出口下游的每条车道上布设感应线圈检测器,利用感应线圈检测器测量每个车队中每辆机动车经过感应线圈检测器前端和后端的时间,进而建立城市干道车队离散模型。本发明用于建立城市干道车队离散模型。
申请公布号 CN102855755A 申请公布日期 2013.01.02
申请号 CN201210327184.3 申请日期 2012.09.06
申请人 哈尔滨工业大学 发明人 别一鸣;裴玉龙;王琳虹;孙烨
分类号 G08G1/00(2006.01)I;G08G1/042(2006.01)I;G06F19/00(2006.01)I 主分类号 G08G1/00(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 1.一种基于运行车速预测的城市干道车队离散模型的建模方法,其特征在于:所述一种基于运行车速预测的城市干道车队离散模型的建模方法的具体步骤如下:步骤一、在城市干道(1)上交叉口出口下游的每条车道上布设感应线圈检测器(2),且每个感应线圈检测器(2)与其最近的交通信号机(3)连接,信号机(3)采集每一辆机动车经过感应线圈检测器(2)后端边缘与前端边缘时产生的脉冲信号,然后记录每个脉冲信号产生的时刻;步骤二、识别经过感应线圈检测器(2)的机动车是否属于同一个车队;步骤三、计算同一个车队中的最高速度、最低速度、平均速度,并预测当前经过感应线圈检测器(2)的车队的速度值;具体步骤为:步骤三(一)、设当前经过感应线圈检测器(2)的为第M个车队,该车队中共有M<sub>N</sub>辆机动车,第n辆车触发感应线圈检测器(2)后端边缘的时刻为t<sub>Mn</sub>,其中1≤n≤M<sub>n</sub>,触发感应线圈检测器(2)前端边缘的时刻为t′<sub>Mn</sub>,那么第n辆车经过感应线圈检测器(2)时的速度V<sub>Mn</sub>为:<maths num="0001"><![CDATA[<math><mrow><msub><mi>V</mi><mi>Mn</mi></msub><mo>=</mo><mfrac><mn>2</mn><mrow><msubsup><mi>t</mi><mi>Mn</mi><mo>&prime;</mo></msubsup><mo>-</mo><msub><mi>t</mi><mi>Mn</mi></msub></mrow></mfrac></mrow></math>]]></maths>①;步骤三(二)、计算第M个车队中机动车的最大车速V<sub>M max</sub>:<maths num="0002"><![CDATA[<math><mrow><msub><mi>V</mi><mrow><mi>M</mi><mi>max</mi></mrow></msub><mo>=</mo><mi>max</mi><mrow><mo>(</mo><msub><mi>V</mi><mrow><mi>M</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>V</mi><mrow><mi>M</mi><mn>2</mn></mrow></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><msub><mi>V</mi><msub><mi>MM</mi><mi>N</mi></msub></msub><mo>)</mo></mrow></mrow></math>]]></maths>②,公式②中max为取最大函数;步骤三(三)、计算第M个车队中机动车的最小车速V<sub>M min</sub>:<maths num="0003"><![CDATA[<math><mrow><msub><mi>V</mi><mrow><mi>M</mi><mi>min</mi></mrow></msub><mo>=</mo><mi>min</mi><mrow><mo>(</mo><msub><mi>V</mi><mrow><mi>M</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>V</mi><mrow><mi>M</mi><mn>2</mn></mrow></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><msub><mi>V</mi><msub><mi>MM</mi><mi>N</mi></msub></msub><mo>)</mo></mrow></mrow></math>]]></maths>③,公式③中min为取最小函数;步骤三(四)、计算第M个车队的机动车平均车速<img file="FDA00002107612600014.GIF" wi="79" he="60" /><maths num="0004"><![CDATA[<math><mrow><msub><mover><mi>V</mi><mo>&OverBar;</mo></mover><mi>M</mi></msub><mo>=</mo><mi>mean</mi><mrow><mo>(</mo><msub><mi>V</mi><mrow><mi>M</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>V</mi><mrow><mi>M</mi><mn>2</mn></mrow></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><msub><mi>V</mi><msub><mi>MM</mi><mi>N</mi></msub></msub><mo>)</mo></mrow></mrow></math>]]></maths>④,公式④中mean为取均值函数;步骤三(五)、预测当前经过感应线圈检测器(2)的车队的速度值,设当前经过感应线圈检测器(2)的为第K个车队,采用三阶移动平均方法预测该车队的最高车速V′<sub>K max</sub>、最低车速V′<sub>K min</sub>、平均车速<img file="FDA00002107612600021.GIF" wi="82" he="67" /><maths num="0005"><![CDATA[<math><mrow><msubsup><mi>V</mi><mrow><mi>K</mi><mi>max</mi></mrow><mo>&prime;</mo></msubsup><mo>=</mo><mfrac><mrow><msub><mi>V</mi><mrow><mrow><mo>(</mo><mi>K</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>max</mi></mrow></msub><mo>+</mo><msub><mi>V</mi><mrow><mrow><mo>(</mo><mi>K</mi><mo>-</mo><mn>2</mn><mo>)</mo></mrow><mi>max</mi></mrow></msub><mo>+</mo><msub><mi>V</mi><mrow><mrow><mo>(</mo><mi>K</mi><mo>-</mo><mn>3</mn><mo>)</mo></mrow><mi>max</mi></mrow></msub></mrow><mn>3</mn></mfrac></mrow></math>]]></maths>⑤<maths num="0006"><![CDATA[<math><mrow><msubsup><mi>V</mi><mrow><mi>K</mi><mi>min</mi></mrow><mo>&prime;</mo></msubsup><mo>=</mo><mfrac><mrow><msub><mi>V</mi><mrow><mrow><mo>(</mo><mi>K</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>min</mi></mrow></msub><mo>+</mo><msub><mi>V</mi><mrow><mrow><mo>(</mo><mi>K</mi><mo>-</mo><mn>2</mn><mo>)</mo></mrow><mi>min</mi></mrow></msub><mo>+</mo><msub><mi>V</mi><mrow><mrow><mo>(</mo><mi>K</mi><mo>-</mo><mn>3</mn><mo>)</mo></mrow><mi>min</mi></mrow></msub></mrow><mn>3</mn></mfrac></mrow></math>]]></maths>⑥<maths num="0007"><![CDATA[<math><mrow><msubsup><mover><mi>V</mi><mo>&OverBar;</mo></mover><mi>K</mi><mo>&prime;</mo></msubsup><mo>=</mo><mfrac><mrow><msub><mover><mi>V</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>K</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msub><mo>+</mo><msub><mover><mi>V</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>K</mi><mo>-</mo><mn>2</mn><mo>)</mo></mrow></msub><mo>+</mo><msub><mover><mi>V</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>K</mi><mo>-</mo><mn>3</mn><mo>)</mo></mrow></msub></mrow><mn>3</mn></mfrac></mrow></math>]]></maths>⑦公式⑤⑥⑦中,V<sub>(K-1)max</sub>表示第K-1个车队的最高车速,V<sub>(K-2)max</sub>表示第K-2个车队的最高车速,V<sub>(K-3)max</sub>表示第K-3个车队的最高车速,V<sub>(K-1)min</sub>表示第K-1个车队的最低车速,V<sub>(K-2)min</sub>表示第K-2个车队的最低车速,V<sub>(K-3)min</sub>表示第K-3个车队的最低车速,<img file="FDA00002107612600025.GIF" wi="108" he="74" />表示第K-1个车队的平均车速,<img file="FDA00002107612600026.GIF" wi="114" he="74" />表示第K-2个车队的平均车速,<img file="FDA00002107612600027.GIF" wi="112" he="74" />表示第K-3个车队的平均车速;步骤四、计算行程时间系数β:<maths num="0008"><![CDATA[<math><mrow><mi>&beta;</mi><mo>=</mo><mfrac><msubsup><mover><mi>V</mi><mo>&OverBar;</mo></mover><mi>K</mi><mo>&prime;</mo></msubsup><msubsup><mi>V</mi><mrow><mi>k</mi><mi>max</mi></mrow><mo>&prime;</mo></msubsup></mfrac></mrow></math>]]></maths>⑧;步骤五、计算平均行程时间T<sub>α</sub>,用时段数表示<maths num="0009"><![CDATA[<math><mrow><msub><mi>T</mi><mi>&alpha;</mi></msub><mo>=</mo><mi>int</mi><mrow><mo>(</mo><mfrac><msup><mi>L</mi><mo>&prime;</mo></msup><mrow><msubsup><mover><mi>V</mi><mo>&OverBar;</mo></mover><mi>K</mi><mo>&prime;</mo></msubsup><mo>&times;</mo><mi>&Delta;t</mi></mrow></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>⑨,公式⑨中,L′表示感应线圈检测器(2)所在断面与下游观测断面之间的距离,Δt表示单位时段长度,int表示取整函数;步骤六、计算车队离散系数α,车队离散系数α受城市干道(1)上车流运行方向的车道数的影响,设路段上车道数为X,则车队离散系数α的计算公式为:<maths num="0010"><![CDATA[<math><mrow><mi>&alpha;</mi><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1.126</mn><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>-</mo><msup><mrow><mo>(</mo><mfrac><mrow><msup><mi>q</mi><mo>&prime;</mo></msup><mo>/</mo><mi>s</mi><mo>-</mo><mn>0.786</mn></mrow><mn>0.107</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow></msup><mo>+</mo><mn>0.793</mn><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>-</mo><msup><mrow><mo>(</mo><mfrac><mrow><msup><mi>q</mi><mo>&prime;</mo></msup><mo>/</mo><mi>s</mi><mo>-</mo><mn>0.809</mn></mrow><mn>1.312</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow></msup><mo>,</mo><mi>x</mi><mo>=</mo><mn>2</mn></mtd></mtr><mtr><mtd><mn>0.16</mn><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>-</mo><msup><mrow><mo>(</mo><mfrac><mrow><msup><mi>q</mi><mo>&prime;</mo></msup><mo>/</mo><mi>s</mi><mo>-</mo><mn>0.741</mn></mrow><mn>0.146</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow></msup><mo>+</mo><mn>0.771</mn><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>-</mo><msup><mrow><mo>(</mo><mfrac><mrow><msup><mi>q</mi><mo>&prime;</mo></msup><mo>/</mo><mi>s</mi><mo>-</mo><mn>0.706</mn></mrow><mn>0.778</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow></msup><mo>,</mo><mi>x</mi><mo>=</mo><mn>3</mn></mtd></mtr><mtr><mtd><mn>0.141</mn><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>-</mo><msup><mrow><mo>(</mo><mfrac><mrow><msup><mi>q</mi><mo>&prime;</mo></msup><mo>/</mo><mi>s</mi><mo>-</mo><mn>0.766</mn></mrow><mn>0.062</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow></msup><mo>+</mo><mn>0.771</mn><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>-</mo><msup><mrow><mo>(</mo><mfrac><mrow><msup><mi>q</mi><mo>&prime;</mo></msup><mo>/</mo><mi>s</mi><mo>-</mo><mn>0.701</mn></mrow><mn>0.553</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow></msup><mo>,</mo><mi>x</mi><mo>=</mo><mn>4</mn></mtd></mtr><mtr><mtd><mn>0.084</mn><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>-</mo><msup><mrow><mo>(</mo><mfrac><mrow><msup><mi>q</mi><mo>&prime;</mo></msup><mo>/</mo><mi>s</mi><mo>-</mo><mn>0.744</mn></mrow><mn>0.06</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow></msup><mo>+</mo><mn>0.815</mn><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>-</mo><msup><mrow><mo>(</mo><mfrac><mrow><msup><mi>q</mi><mo>&prime;</mo></msup><mo>/</mo><mi>s</mi><mo>-</mo><mn>0.742</mn></mrow><mn>0.416</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow></msup><mo>,</mo><mi>x</mi><mo>=</mo><mn>5</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>⑩,公式⑩中,S表示感应线圈检测器(2)所在断面下游交叉口一条进口车道的饱和流率,q′表示预测第K个车队经过感应线圈检测器(2)时的流率,q′的计算公式为:<maths num="0011"><![CDATA[<math><mrow><msup><mi>q</mi><mo>&prime;</mo></msup><mo>=</mo><mfrac><mrow><mo>[</mo><msub><mrow><mo>(</mo><mi>K</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>N</mi></msub><mo>+</mo><msub><mrow><mo>(</mo><mi>K</mi><mo>-</mo><mn>2</mn><mo>)</mo></mrow><mi>N</mi></msub><mo>+</mo><msub><mrow><mo>(</mo><mi>K</mi><mo>-</mo><mn>3</mn><mo>)</mo></mrow><mi>N</mi></msub><mo>]</mo><mo>/</mo><mn>3</mn></mrow><mrow><mrow><mo>(</mo><msub><mi>T</mi><mrow><mi>K</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>T</mi><mrow><mi>K</mi><mo>-</mo><mn>2</mn></mrow></msub><mo>+</mo><msub><mi>T</mi><mrow><mi>K</mi><mo>-</mo><mn>3</mn></mrow></msub><mo>)</mo></mrow><mo>/</mo><mn>3</mn></mrow></mfrac></mrow></math>]]></maths>其中(K-1)<sub>N</sub>表示第K-1个车队所包含的车辆数,T<sub>K-1</sub>表示第K-1个车队头车触发感应线圈检测器(2)与车队尾车触发感应线圈检测器(2)之间的时间差,(K-2)<sub>N</sub>表示第K-2个车队所包含的车辆数,T<sub>K-2</sub>表示第K-2个车队头车触发感应线圈检测器(2)与车队尾车触发感应线圈检测器(2)之间的时间差,(K-3)<sub>N</sub>表示第K-3个车队所包含的车辆数,T<sub>K-3</sub>表示第K-3个车队头车触发感应线圈检测器(2)与车队尾车触发感应线圈检测器(2)之间的时间差;步骤七、计算平滑系数F:<maths num="0012"><![CDATA[<math><mrow><mi>F</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><mi>&alpha;</mi><mo>&times;</mo><mi>&beta;</mi><mo>&times;</mo><msub><mi>T</mi><mi>&alpha;</mi></msub></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>步骤八、将经过感应线圈检测器(2)的第K个车队的机动车划分至不同时段,具体步骤为:步骤八(一)、以第K个车队头车触发感应线圈检测器(2)后端边缘的时刻记为0时刻,车队尾车触发感应线圈检测器(2)后端边缘的时刻记为T<sub>K</sub>;步骤八(二)、以Δt为单位时段长度,将[0,T<sub>K</sub>]时段内经过感应线圈检测器(2)的机动车划分至不同时段;步骤九、建立车队离散模型,预测下游断面的车流到达图式:Q′(m)=F·Q(m-T<sub>α</sub>)+(1-F)·Q(m-1),m>T<sub>α</sub>其中Q′(m)表示第K个车队在第m个时段到达下游断面的车辆数,Q(m-T<sub>α</sub>)表示第K个车队在m-T<sub>α</sub>个时段通过感应线圈检测器(2)所在断面的车辆数。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号