发明名称 基于时间序列的高速公路交通流量预测方法
摘要 一种基于时间序列的高速公路交通流量预测方法,包括选择一个时间标度,统计建立交通流量时间序列Q=(x<sub>i</sub>)的步骤;根据所选择的时间标度,自回归项个数p和移动平均项个数q的取值范围的步骤;求取自回归项个数p和移动平均项个数q的步骤;将最优自回归项个数p和移动平均项个数q,结合交通流量时间序列Q,采用极大似然估计进行拟合得到最优ARMA模型,得到历史测量值的权值参数和历史误差值的权值参数的步骤;以及求取不同时间标度下的交通流量的预测序列<img file="DDA0000576400470000011.GIF" wi="168" he="77" />的步骤。本发明所的时间序列模型能更好的满足高速公路各类型流量的预测需求,提高了预测的普适性。运算简单,提高了预测的效率和速度,满足了高速公路交通预测的工程需求。
申请公布号 CN104269055A 申请公布日期 2015.01.07
申请号 CN201410495385.3 申请日期 2014.09.24
申请人 四川省交通科学研究所 发明人 戴元;盛鹏;代雨婷;柏吉琼;谢强;王俊峰
分类号 G08G1/01(2006.01)I;G06Q10/04(2012.01)I 主分类号 G08G1/01(2006.01)I
代理机构 成都信博专利代理有限责任公司 51200 代理人 卓仲阳
主权项 一种基于时间序列的高速公路交通流量预测方法,其特征在于,包括步骤1:选择一个时间标度,统计建立交通流量时间序列Q=(x<sub>i</sub>),其中i=1,2,...,n;所述时间标度为月、日或小时;步骤2:根据所选择的时间标度,设定自回归项个数p和移动平均项个数q的取值范围:当时间标度为月,p、q的取值范围为1至3;当时间标度为日,p、q的取值范围为1至5;当时间标度为小时,p、q的取值范围为1至10;步骤3:求取最优自回归项个数p和移动平均项个数q:任意选取自回归项个数p和移动平均项个数q的一种组合,结合交通流量时间序列Q,采用极大似然估计进行拟合得到ARMA模型,根据ARMA模型计算出赤池信息量准则AIC值与贝叶斯信息准则BIC值;以赤池信息量准则AIC值且贝叶斯信息准则BIC值为最小时所对应的自回归项个数p和移动平均项个数q为最优自回归项个数p和移动平均项个数q;步骤4:将最优自回归项个数p和移动平均项个数q,结合交通流量时间序列Q,采用极大似然估计进行拟合得到最优ARMA模型,其自回归项即历史测量值的权值参数为<img file="FDA0000576400450000011.GIF" wi="266" he="65" />移动平均项即历史误差值的权值参数为θ<sub>1</sub>,θ<sub>2</sub>,...,θ<sub>q</sub>;步骤5:求取不同时间标度下的交通流量的预测序列<img file="FDA0000576400450000012.GIF" wi="201" he="83" />其中i=1,2,3,...,n,n+1;包括(1)当时间标度为月:a.当i=1时,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>x</mi><mo>^</mo></mover><mi>i</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>x</mi><mi>k</mi></msub></mrow><mi>n</mi></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000576400450000013.GIF" wi="247" he="194" /></maths>预测误差<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&epsiv;</mi><mi>i</mi></msub><mo>=</mo><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mn>1</mn></msub><mo>;</mo></mrow>]]></math><img file="FDA0000576400450000014.GIF" wi="240" he="73" /></maths>b.当i≥2时,<img file="FDA0000576400450000015.GIF" wi="842" he="143" />其中,当i≤k时,令x<sub>i‑k</sub>=0;当i≤j时,令ε<sub>i‑j</sub>=0;当i>j时,预测误差<img file="FDA0000576400450000016.GIF" wi="243" he="72" />(2)当时间标度为日:a.当i=1时,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mover><mi>x</mi><mo>^</mo></mover><mi>i</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>x</mi><mi>k</mi></msub></mrow><mi>n</mi></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000576400450000017.GIF" wi="247" he="194" /></maths>预测误差<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&epsiv;</mi><mi>i</mi></msub><mo>=</mo><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mn>1</mn></msub><mo>;</mo></mrow>]]></math><img file="FDA0000576400450000018.GIF" wi="240" he="73" /></maths>b.当2≤i<14时,<img file="FDA0000576400450000019.GIF" wi="836" he="137" />其中,当i≤k时,令x<sub>i‑k</sub>=0;当i≤j时,令ε<sub>i‑j</sub>=0;当i>j时,预测误差<img file="FDA00005764004500000110.GIF" wi="239" he="72" />c.当i≥14时,<img file="FDA0000576400450000021.GIF" wi="916" he="144" />其中,预测误差<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&epsiv;</mi><mi>i</mi></msub><mo>=</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mi>i</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0000576400450000022.GIF" wi="241" he="72" /></maths>调节权值系数<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>&rho;</mi><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mi>i</mi><mo>-</mo><mn>6</mn></mrow><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>x</mi><mi>k</mi></msub></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>-</mo><mn>13</mn></mrow><mrow><mi>i</mi><mo>-</mo><mn>8</mn></mrow></munderover><msub><mi>x</mi><mi>j</mi></msub></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000576400450000023.GIF" wi="270" he="282" /></maths>(3)当时间标度为小时:a.当i=1时,<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mover><mi>x</mi><mo>^</mo></mover><mi>i</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>x</mi><mi>k</mi></msub></mrow><mi>n</mi></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000576400450000024.GIF" wi="247" he="195" /></maths>预测误差<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>&epsiv;</mi><mi>i</mi></msub><mo>=</mo><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mn>1</mn></msub><mo>;</mo></mrow>]]></math><img file="FDA0000576400450000025.GIF" wi="240" he="74" /></maths>b.当2≤i≤168时,<img file="FDA0000576400450000026.GIF" wi="841" he="144" />其中,当i≤k时,令x<sub>i‑k</sub>=0;当i≤j时,令ε<sub>i‑j</sub>=0;当i>j时,预测误差<img file="FDA0000576400450000027.GIF" wi="240" he="73" />c.当i>168时,<img file="FDA0000576400450000028.GIF" wi="917" he="144" />其中,预测误差<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>&epsiv;</mi><mi>i</mi></msub><mo>=</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mi>i</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0000576400450000029.GIF" wi="246" he="72" /></maths>当<img file="FDA00005764004500000211.GIF" wi="319" he="57" />时,令调节权值系数ρ=1;当<img file="FDA00005764004500000212.GIF" wi="333" he="67" />时,令调节权值系数<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mi>&rho;</mi><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mn>7</mn></munderover><msub><mi>x</mi><mrow><mi>i</mi><mo>-</mo><mn>24</mn><mo>&times;</mo><mi>k</mi></mrow></msub></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>7</mn></munderover><msub><mi>x</mi><mrow><mi>i</mi><mo>-</mo><mn>24</mn><mo>&times;</mo><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></mfrac><mo>.</mo></mrow>]]></math><img file="FDA00005764004500000210.GIF" wi="345" he="281" /></maths>
地址 610041 四川省成都市武侯区武侯祠横街4号