发明名称 交通流短时预测方法
摘要 本发明公开了一种交通流短时预测方法,包括如下步骤:在选定的路段上设置检测器,按照预设的时间周期采集交通流数据;预处理获得的交通流数据,判断车流量和速度是否处于预期范围;建立交通流短时预测模型;检验上述模型是否符合平稳性要求,如果不符合,在进行差分处理,直到其符合平稳性要求;对符合平稳性要求的模型进行参数估计;采用上述模型预测交通流,并评价采用相关评价指标对其进行评价。通过线性ARIMA模型和非线性EGARCH-M模型的结合,本发明能够更好地追踪交通流的数据特征,控制异常数据带来的不利影响,具有更高的预测精度和可靠性,在各项评价指标方面优于现有方法。
申请公布号 CN103903452B 申请公布日期 2015.10.07
申请号 CN201410088566.4 申请日期 2014.03.11
申请人 东南大学 发明人 叶智锐;王超
分类号 G08G1/065(2006.01)I;G08G1/01(2006.01)I 主分类号 G08G1/065(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 柏尚春
主权项 一种交通流短时预测方法,其特征在于,包括如下步骤:S1、在选定的路段上设置检测器,按照预设的时间周期采集交通流数据;S2、预处理获得的交通流数据,判断车流量和速度是否处于预期范围;S3、建立交通流短时预测模型;S31、建立ARIMA模型:<img file="FDA0000736183740000011.GIF" wi="513" he="138" />其中,X<sub>t</sub>为时间序列;X<sub>t‑i</sub>表示相邻为i的时间序列;p和q分别为AR和MA项;<img file="FDA0000736183740000012.GIF" wi="61" he="59" />和θ<sub>j</sub>是未知系数;ε<sub>t‑j</sub>为t‑j时的随机误差;i=1,2,…,p,j=0,1,…,q;S32:建立ARIMA‑EGARCH‑M复合模型:<img file="FDA0000736183740000013.GIF" wi="848" he="143" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mo>&dtri;</mo><msub><mi>X</mi><mi>t</mi></msub><mo>=</mo><msub><mi>X</mi><mi>t</mi></msub><mo>-</mo><msub><mi>X</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo></mrow>]]></math><img file="FDA0000736183740000014.GIF" wi="350" he="76" /></maths>a<sub>t</sub>=σ<sub>t</sub>ε<sub>t</sub><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>l</mi><mi>n</mi><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mi>t</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>=</mo><msub><mi>&alpha;</mi><mn>0</mn></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>&mu;</mi></munderover><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>&lsqb;</mo><mo>|</mo><mfrac><msub><mi>a</mi><mrow><mi>t</mi><mo>-</mo><mi>i</mi></mrow></msub><msub><mi>&sigma;</mi><mrow><mi>t</mi><mo>-</mo><mi>i</mi></mrow></msub></mfrac><mo>|</mo><mo>-</mo><mi>u</mi><mo>&rsqb;</mo><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>&mu;</mi></munderover><msub><mi>&gamma;</mi><mi>i</mi></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>v</mi></munderover><msub><mi>&beta;</mi><mi>j</mi></msub><mi>l</mi><mi>n</mi><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mrow><mi>t</mi><mo>-</mo><mi>j</mi></mrow><mn>2</mn></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000736183740000015.GIF" wi="1208" he="194" /></maths>其中,<img file="FDA0000736183740000016.GIF" wi="626" he="170" />a<sub>t</sub>=σ<sub>t</sub>ε<sub>t</sub>,σ<sub>t</sub>、ε<sub>t</sub>分别为t时所对应的方差和随机误差;<img file="FDA0000736183740000017.GIF" wi="52" he="60" />为差分符号;d为差分阶数;γ<sub>i</sub>表示非对称效应;a<sub>t‑i</sub>=σ<sub>t‑i</sub>ε<sub>t‑i</sub>,σ<sub>t‑i</sub>、ε<sub>t‑i</sub>分别为t‑i所对应的方差和随机误差;<img file="FDA0000736183740000018.GIF" wi="62" he="61" />为待定参数;μ和ν分别为GARCH和ARCH项;X<sub>t‑1</sub>表示与X<sub>t</sub>相邻的时间序列;a<sub>t‑j</sub>=σ<sub>t‑j</sub>ε<sub>t‑j</sub>,σ<sub>t‑j</sub>、ε<sub>t‑j</sub>分别为t‑j所对应的方差和随机误差;α<sub>0</sub>,α<sub>i</sub>和β<sub>j</sub>均为未知系数;i=1,2,…,μ,j=1,2,…,ν;S33:令ε<sub>t</sub>服从广义误差分布,其密度函数为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>f</mi><mrow><mo>(</mo><msub><mi>&epsiv;</mi><mi>t</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>n</mi><mo>&CenterDot;</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mrow><mo>|</mo><mrow><msub><mi>&epsiv;</mi><mi>t</mi></msub><mo>/</mo><mi>&lambda;</mi></mrow><mo>|</mo></mrow><mi>n</mi></msup><mo>}</mo></mrow><mrow><mi>&lambda;</mi><mo>&CenterDot;</mo><msup><mn>2</mn><mrow><mn>1</mn><mo>+</mo><mfrac><mn>1</mn><mi>n</mi></mfrac></mrow></msup><mo>&CenterDot;</mo><mi>&Gamma;</mi><mrow><mo>(</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000736183740000021.GIF" wi="571" he="277" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>&lambda;</mi><mo>=</mo><msup><mn>2</mn><mrow><mo>-</mo><mfrac><mn>1</mn><mi>n</mi></mfrac></mrow></msup><msqrt><mfrac><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mn>1</mn><mo>/</mo><mi>n</mi><mo>)</mo></mrow></mrow><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mn>3</mn><mo>/</mo><mi>n</mi><mo>)</mo></mrow></mrow></mfrac></msqrt></mrow>]]></math><img file="FDA0000736183740000022.GIF" wi="352" he="160" /></maths>其中,Γ为伽马函数;n为描述尾厚的分布参数;λ为中间变量;S34:建立ARIMA‑EGARCH‑M‑GED模型:<img file="FDA0000736183740000023.GIF" wi="844" he="144" /><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mo>&dtri;</mo><msub><mi>X</mi><mi>t</mi></msub><mo>=</mo><msub><mi>X</mi><mi>t</mi></msub><mo>-</mo><msub><mi>X</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo></mrow>]]></math><img file="FDA0000736183740000024.GIF" wi="350" he="71" /></maths>a<sub>t</sub>=σ<sub>t</sub>ε<sub>t</sub><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>ln</mi><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mi>t</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>=</mo><msub><mi>&alpha;</mi><mn>0</mn></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>&mu;</mi></munderover><mrow><msub><mi>&alpha;</mi><mi>i</mi></msub><mrow><mo>&lsqb;</mo><mrow><mrow><mo>(</mo><mrow><mrow><mo>|</mo><mfrac><msub><mi>a</mi><mrow><mi>t</mi><mo>-</mo><mi>i</mi></mrow></msub><msub><mi>&sigma;</mi><mrow><mi>t</mi><mo>-</mo><mi>i</mi></mrow></msub></mfrac><mo>|</mo></mrow><mo>-</mo><msqrt><mfrac><mn>2</mn><mi>&pi;</mi></mfrac></msqrt></mrow><mo>)</mo></mrow><mo>+</mo><msub><mi>&gamma;</mi><mi>i</mi></msub><mfrac><msub><mi>a</mi><mrow><mi>t</mi><mo>-</mo><mi>i</mi></mrow></msub><msub><mi>&sigma;</mi><mrow><mi>t</mi><mo>-</mo><mi>i</mi></mrow></msub></mfrac></mrow><mo>&rsqb;</mo></mrow><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>v</mi></munderover><mrow><msub><mi>&beta;</mi><mi>j</mi></msub><mi>ln</mi><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mrow><mi>t</mi><mo>-</mo><mi>j</mi></mrow><mn>2</mn></msubsup><mo>)</mo></mrow></mrow></mrow></mrow>]]></math><img file="FDA0000736183740000025.GIF" wi="1370" he="195" /></maths>其中,<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>E</mi><mrow><mo>|</mo><msub><mi>&epsiv;</mi><mi>t</mi></msub><mo>|</mo></mrow><mo>=</mo><msqrt><mfrac><mn>2</mn><mi>&pi;</mi></mfrac></msqrt><mo>;</mo></mrow>]]></math><img file="FDA0000736183740000026.GIF" wi="270" he="150" /></maths>S4、检验上述模型是否符合平稳性要求,如果不符合,在进行差分处理,直到其符合平稳性要求;S5、对符合平稳性要求的模型进行参数估计;S6、采用上述模型预测交通流,并评价采用相关评价指标对其进行评价。
地址 210096 江苏省南京市四牌楼2号