交通分布預(yù)測(cè)模型

1、復(fù)旦大學(xué)第五章 交通分布預(yù)測(cè)交通分布預(yù)測(cè)模型模型 5.1 定義與符號(hào)定義與符號(hào)5.2 增長(zhǎng)系數(shù)法增長(zhǎng)系數(shù)法5.3 最大熵模型最大熵模型5.4 重力模型重力模型教材第五章教材第五章 175205頁(yè)頁(yè)本章目錄本章目錄5.5 重力模型的標(biāo)定重力模型的標(biāo)定5.6 三段比例法三段比例法5.7 其他數(shù)學(xué)模型其他數(shù)學(xué)模型5.8 實(shí)際考慮實(shí)際考慮1 研究?jī)?nèi)容研究?jī)?nèi)容2 增長(zhǎng)系數(shù)法增長(zhǎng)系數(shù)法及評(píng)價(jià)及評(píng)價(jià)3 重力重力模型模型及評(píng)價(jià)及評(píng)價(jià)4 最大熵模型最大熵模型簡(jiǎn)介簡(jiǎn)介課堂講解思路課堂講解思路交通發(fā)生與吸引交通分布方式劃分交通流分配如何得到OD矩陣交通分布預(yù)測(cè)交通分布預(yù)測(cè) DO123合計(jì) DO123合計(jì)117.0

2、7.0 4.0 28.0 1 38.6 27.0 38.0 6.0 51.0 2 91.9 34.0 5.0 17.0 26.0 3 36.0 合計(jì)28.0 50.0 27.0 105.0 合計(jì)39.3 90.3 36.9 166.5 表表 現(xiàn)狀現(xiàn)狀OD表（單位：萬(wàn)次）表（單位：萬(wàn)次） 表表 將來(lái)的發(fā)生與吸引交通量將來(lái)的發(fā)生與吸引交通量如何通過(guò)現(xiàn)狀OD表、交通需求預(yù)測(cè)結(jié)果求得每個(gè)OD小區(qū)的交通量？2. 增長(zhǎng)系數(shù)法在現(xiàn)狀交通分布量給定的情況下，預(yù)測(cè)將來(lái)的交通分布量。2.1 單增長(zhǎng)系數(shù)法每格同步擴(kuò)大若只有整個(gè)規(guī)劃區(qū)域的平均增長(zhǎng)系數(shù)，則可以假設(shè)它適用于矩陣的每一單元。表 基年出行矩陣12341550

3、100200355250510030045535010051002554100200250205702053554556201635i表 未來(lái)年出行矩陣（ =1.2）12341660120240426260612036054636012061203064120240300246842464265467441962i2. 增長(zhǎng)系數(shù)法2.2 單約束增長(zhǎng)系數(shù)法知行求列、知列求行已知每一小區(qū)出行產(chǎn)生（或吸引的值）的期望增長(zhǎng)值信息，可利用由發(fā)點(diǎn)決定的出行產(chǎn)生率來(lái)估計(jì)矩陣中相應(yīng)行（列）的元素值。表 未來(lái)年出行矩陣DO123合計(jì)合計(jì)預(yù)測(cè)值預(yù)測(cè)值117.07.04.028.038.627.038.06.051

4、.091.934.05.017.026.036.0合計(jì)合計(jì)28.050.027.0105.0166.5表 基年出行矩陣DO123合計(jì)合計(jì)目標(biāo)值目標(biāo)值123.436 9.650 5.51438.638.6212.614 68.475 10.81291.991.935.538 6.923 23.53836.036.0合計(jì)合計(jì)41.588 85.048 39.865166.5166.52. 增長(zhǎng)系數(shù)法2.3雙約束增長(zhǎng)率法知行知列求平均每個(gè)小區(qū)的出行產(chǎn)生量與吸引量均已知，即每個(gè)小區(qū)的發(fā)生量與吸引量的增長(zhǎng)率將不同，可取兩者平均迭代。表表 基年出行矩陣基年出行矩陣 123總計(jì)預(yù)測(cè)系數(shù)117.0 7.0 4

5、.0 28.0 38.6 1.378627.0 38.0 6.0 51.0 91.9 1.802034.0 5.0 17.0 26.0 36.0 1.3846總計(jì)28.0 50.0 27.0 105預(yù)測(cè)39.3 90.3 36.9 90.3 系數(shù)1.40361.8060 1.366717*1/2(1.3786+1.4036)2. 增長(zhǎng)系數(shù)法表表 未來(lái)未來(lái)年出年出行矩陣行矩陣 第二次迭代結(jié)果第二次迭代結(jié)果 DO123合計(jì)122.819 11.080 5.270 39.169 211.226 70.585 9.462 91.273 35.427 7.995 22.637 36.058 合計(jì)39.4

6、71 89.660 37.369 166.500 收斂判定：行、列所有系數(shù)誤差小于3%，即所得系數(shù)在（0.97，1.03）詳見(jiàn)例詳見(jiàn)例12.4 增長(zhǎng)系數(shù)法的優(yōu)點(diǎn)增長(zhǎng)系數(shù)法的優(yōu)點(diǎn)（1）結(jié)構(gòu)簡(jiǎn)單、實(shí)用的比較多，不需要交通小區(qū)之間的距離和時(shí)間。（2）可以適用于小時(shí)交通量或日交通量等的預(yù)測(cè)，也可以獲得各種交通目的的OD交通量。（3）對(duì)于變化較小的OD表預(yù)測(cè)非常有效。（鐵路車站間）2.4 增長(zhǎng)系數(shù)法的缺點(diǎn)增長(zhǎng)系數(shù)法的缺點(diǎn)（1）必須有所有小區(qū)的OD交通量。（2）對(duì)象地區(qū)發(fā)生如下大規(guī)模變化時(shí)，該方法不適用： 將來(lái)的交通小區(qū)分區(qū)發(fā)生變化（有新開(kāi)發(fā)區(qū)時(shí)）； 交通小區(qū)之間的行駛時(shí)間發(fā)生變化時(shí)； 土地利用發(fā)生較大變

7、化時(shí)。（3）交通小區(qū)之間的交通量值較小時(shí)，存在如下問(wèn)題： 若現(xiàn)狀交通量為零，那么將來(lái)預(yù)測(cè)值也為零； 對(duì)于可靠性較低的OD交通量，將來(lái)的預(yù)測(cè)誤差將被擴(kuò)大。（4）將來(lái)交通量?jī)H用一個(gè)增長(zhǎng)系數(shù)表示缺乏合理性。 1234目標(biāo)15501002003554002050005046035010051002554004100200250205707021554003553201230 目標(biāo)年 2604005008021962 1234目標(biāo)13.40.761.0355.342040020388.200388460365.52.85.9345.74204004191.28.3433.1101.07347022604

8、005008021962 目標(biāo)260400500802 1962增長(zhǎng)系數(shù)法增長(zhǎng)系數(shù)法的的困境困境小區(qū)交通量為小區(qū)交通量為0迭代后的矩陣迭代后的矩陣1、目標(biāo)年與用模型算出的總量之間有很大的差別。2、這一行中的和不可能是460，因?yàn)槲ㄒ坏姆橇阍氐淖畲笾祪H為400。3.重力模型法（Gravity Model）3.1模型的提出 考慮到交通阻抗，1955年Casey提出重力模型2ijjidPPijq 重力模型重力模型 iPjP分別表示i小區(qū)和j小區(qū)的人口；ijd表示i，j小區(qū)之間的距離；改進(jìn)后的改進(jìn)后的重力模型重力模型 i、j分區(qū)之間的出行量（i為產(chǎn)生區(qū)，j為吸引區(qū)）預(yù)測(cè)值； 兩分區(qū)之間的交通阻抗；

9、分別為分區(qū)i的出行產(chǎn)生量、分區(qū)j的吸引量； 待定系數(shù)。ijijijP AqKRijqijRijP A、K 、 、阻抗 可以是出行時(shí)間、距離、油耗等因素的綜合。城市交通：時(shí)間地區(qū)交通：距離ijR3.重力模型法（Gravity Model）采用線性回歸方法標(biāo)定模型中的參數(shù)，兩邊取自然對(duì)數(shù)：ijijijP AqKR參數(shù)的標(biāo)定lnlnlnlnlnijijijqKPAR此處 可從現(xiàn)狀調(diào)查數(shù)據(jù)中取若干個(gè)分區(qū)作為樣本，待標(biāo)定的參數(shù)有：,ijijijP ARq故要用多元線性回歸方法，設(shè)：123ln ,(1,) (1,ln ,ln ,ln)ijijijYq XX X XPARln, , ,K 3.重力模型法（G

10、ravity Model）根據(jù)經(jīng)驗(yàn)， 取值范圍0.51.0，多數(shù)情況下，可取 ，則問(wèn)題簡(jiǎn)化成了：0123ln,bK bbb01123() ()TTbbBX XX Ybb 、1ijijijPAqKR這只需要用一元線性回歸方法就能標(biāo)定參數(shù)K和 了。 例題13.重力模型法（Gravity Model）問(wèn)題：誤差太大原因（1）模型系數(shù)無(wú)法保證： 和 ，既對(duì)系數(shù)K沒(méi)有約束范圍； （2）內(nèi)內(nèi)出行的阻抗值問(wèn)題，既一個(gè)分區(qū)內(nèi)的出行 值不好確定。解決方法（1）增加約束條件 （2）修改“阻抗函數(shù)” ijjiqAijRijijqP3.重力模型法（Gravity Model）阻抗函數(shù)ijijijP AqKRijR 就

11、是阻抗的函數(shù)，并且可以有多種形式。一般地，用來(lái)表示阻抗函數(shù)，常用的有：冪型：指數(shù)型：（冪與指數(shù)）復(fù)合型：半鐘型：離散型：()ijf R()ijijf RR()exp()ijijf Rb R ()ijbRijijf ReR1()ijijf RabR()ijijmmmf Rr在公式 中，3.重力模型法（Gravity Model）3.2 單約束重力模型單約束重力模型約束條件：ijijqPijiiqA或ijijijPAqKR把重力模型公式 帶入上面第一個(gè)等式，得：ijjijiijjjijijPAAqKKPPRR從而得：1/jijjKAR對(duì)于一般的阻抗函數(shù) ，就可以寫成： ()ijf R1()jijj

12、KA f R3.重力模型法（Gravity Model）K叫做“行約束系數(shù)”，如果從列的角度進(jìn)行約束，同理可以得出“列約束系數(shù)”：用行約束系數(shù)后，重力模型就變成：1()jijjKA f R1()iijiKP f R()()ijijijjijP A fRqA fR此模型的參數(shù)比原來(lái)少一個(gè)，無(wú)需單獨(dú)標(biāo)定K，只要標(biāo)定 中的參數(shù)，就可以算出K。()ijf R3.重力模型法（Gravity Model）參數(shù)的標(biāo)定：以 為例，采用“試算法”來(lái)標(biāo)定b。()ijbRijf Re算法：步1：給b一個(gè)初值，如b=1.步2：從模型 算得現(xiàn)狀的出行量“理論值” ，得現(xiàn)狀理論分布表。步3：計(jì)算現(xiàn)狀實(shí)際PA分布表的平均交

13、通阻抗 ；再計(jì)算理論分布表的平均阻抗 ，求兩者之間的相對(duì)誤差：ijijbRijijbRjjP A eqA eijq1ijijijRqRQ 1ijijijRq RQ3.重力模型法（Gravity Model）當(dāng) 接受關(guān)于b的假設(shè)，否則執(zhí)行下一步。步4：當(dāng) ，既 ，應(yīng)減少b值，令b=b/2；反之增加b值，令b=2b。返回第2步。算法結(jié)束。 100%RRR3%0RR3.3 雙約束重力模型：同時(shí)引進(jìn)行約束系數(shù)和列約束系數(shù)的模型叫雙約束重力模型，起形式是：()ijijijijqK K PA f R1()ijjijjKK A f R1()jiiijjKK P f R 分別為行約束系數(shù)、列約束系數(shù)ijKK、

14、3.重力模型法（Gravity Model）雙約束重力模型雙約束重力模型參數(shù)標(biāo)定：參數(shù)標(biāo)定：以以 為例標(biāo)定參數(shù)，仍然采用試算法。為例標(biāo)定參數(shù)，仍然采用試算法。()ijijf RR步步1：給參數(shù)：給參數(shù) 取初值，可取取初值，可取 。步步2：用迭代法求約束系數(shù)：用迭代法求約束系數(shù) ：2-1.首先令各個(gè)列約束系數(shù)首先令各個(gè)列約束系數(shù)2-2.將各列約束系數(shù)代入式將各列約束系數(shù)代入式求各個(gè)行約束系數(shù)求各個(gè)行約束系數(shù)2-3.再將求得的各個(gè)行約束系數(shù)代入再將求得的各個(gè)行約束系數(shù)代入求各個(gè)列約束系數(shù)求各個(gè)列約束系數(shù)2-4.比較前后兩批列約束系數(shù)，考察它們的相對(duì)誤差是比較前后兩批列約束系數(shù)，考察它們的相對(duì)誤差是

15、否小于否小于3%，若是，轉(zhuǎn)至第，若是，轉(zhuǎn)至第3步，否則返回步，否則返回2-2步。步。1ijKK、1(1,., )jKjn1()ijjijjKK A f RiK1()jiiijjKK P f RjK3.重力模型法（Gravity Model）3：將求得的約束系數(shù)代入 中，用現(xiàn)狀值求現(xiàn)狀理論分布表步4：計(jì)算現(xiàn)狀實(shí)際PA分布表的平均交通阻抗 ；再計(jì)算理論分布表的平均阻抗 ，求兩者之間的相對(duì)誤差 。當(dāng) 接受關(guān)于 的假設(shè)，否則執(zhí)行下一步。步5：當(dāng) ，既 ，應(yīng)減少 值，令 = /2；反之增加 值，令 =2 。返回第2步。算法結(jié)束。例題2 ()ijijijijqK K PA f RijPA、i jq1ijijijRq RQ1ijijijRqRQ 3%0RR 3.4 重力模型的優(yōu)點(diǎn)的優(yōu)點(diǎn)（1）直觀上容易理解。（2）能考慮路網(wǎng)的變化和土地利用對(duì)人們的出行產(chǎn)生的影響。（3）能比較敏感地反映交通小區(qū)之間行駛時(shí)間變化的情況。3.5 重力模型的的缺點(diǎn)缺點(diǎn)（1）缺乏對(duì)人出行行