技術(shù)領(lǐng)域：
：本發(fā)明涉及一種尾號(hào)限行條件下的交通分配修正方法，屬于城市道路網(wǎng)絡(luò)交通分配
技術(shù)領(lǐng)域：
。
背景技術(shù)：
：：機(jī)動(dòng)車尾號(hào)限行措施是我國(guó)特有的一種交通制度，起初作為一種交通管理控制策略誕生時(shí)的目的是緩解城市交通擁堵，最早為2007年北京奧運(yùn)會(huì)測(cè)試賽期間北京市實(shí)施的機(jī)動(dòng)車尾號(hào)單雙號(hào)限行。近年來(lái)隨著我國(guó)各大城市機(jī)動(dòng)車保有量的迅猛增長(zhǎng)，交通需求逐漸接近甚至大于交通供給，最直觀的表現(xiàn)就是城市道路的交通擁堵，很多城市除了出臺(tái)私家車限購(gòu)政策外，也希望通過(guò)尾號(hào)限行措施更直接地緩解交通擁堵，而尾號(hào)限行措施的有效性估計(jì)是必不可少的環(huán)節(jié)。目前國(guó)際通用的城市道路網(wǎng)絡(luò)交通分配方法已經(jīng)十分成熟，然而針對(duì)我國(guó)特有的機(jī)動(dòng)車尾號(hào)限行條件下的交通分配修正方法，還是研究空白。由于尾號(hào)限行條件下駕駛員的路徑選擇結(jié)果與未實(shí)施尾號(hào)限行措施時(shí)是截然不同的，同時(shí)尾號(hào)限行之后分方式od分布矩陣會(huì)發(fā)生相應(yīng)變化，依舊采用未被修正的交通分配方法進(jìn)行尾號(hào)限行措施有效性測(cè)試是不合適的，因此目前已經(jīng)實(shí)施過(guò)尾號(hào)限行措施的十余個(gè)城市，僅能通過(guò)實(shí)際實(shí)施尾號(hào)限行方案后的路段交通量情況進(jìn)行方案評(píng)估，進(jìn)而對(duì)尾號(hào)限行方案進(jìn)行相應(yīng)調(diào)整，而不能提前對(duì)所要實(shí)施的尾號(hào)限行方案進(jìn)行有效性測(cè)試。在智慧城市、智慧交通建設(shè)背景下，許多城市都將建立城市虛擬交通系統(tǒng)測(cè)試平臺(tái)，而尾號(hào)限行方案測(cè)試將會(huì)成為城市交通管理控制方案量化測(cè)試平臺(tái)的功能之一。根據(jù)最新的城市基礎(chǔ)路網(wǎng)信息、交通小區(qū)信息、分方式od分布矩陣數(shù)據(jù)和用于測(cè)試的尾號(hào)限行方案的實(shí)施對(duì)象、作用時(shí)間和作用范圍，對(duì)od分布矩陣和交通分配的具體流程進(jìn)行修正，得到尾號(hào)限行條件下的交通分配結(jié)果，即可對(duì)所測(cè)試的尾號(hào)限行方案實(shí)施效果進(jìn)行量化評(píng)估，避免直接實(shí)施不恰當(dāng)?shù)奈蔡?hào)限行方案而給城市居民的交通出行造成不便。為此，本發(fā)明提出了一種尾號(hào)限行條件下的交通分配修正方法。技術(shù)實(shí)現(xiàn)要素：本發(fā)明的目的是針對(duì)現(xiàn)有尾號(hào)限行條件下駕駛員的路徑選擇結(jié)果與未實(shí)施尾號(hào)限行措施時(shí)是截然不同的，然而目前還沒(méi)有針對(duì)尾號(hào)限行條件下的交通分配修正方法的問(wèn)題，提出一種尾號(hào)限行條件下的交通分配修正方法，充分考慮了我國(guó)各大城市現(xiàn)有的機(jī)動(dòng)車尾號(hào)限行方案，對(duì)交通分配方法進(jìn)行修正，彌補(bǔ)了尾號(hào)限行條件下交通分配修正方法的空白，對(duì)城市虛擬交通測(cè)試平臺(tái)中準(zhǔn)確、快速實(shí)現(xiàn)機(jī)動(dòng)車尾號(hào)限行策略實(shí)施效果的量化測(cè)試具有較大意義。上述的目的通過(guò)以下技術(shù)方案實(shí)現(xiàn)：一種尾號(hào)限行條件下的交通分配修正方法，該方法包括如下步驟：一種尾號(hào)限行條件下的交通分配修正方法，該方法包括如下步驟：(h)獲取最新的城市基礎(chǔ)路網(wǎng)信息，包括交叉口編號(hào)node_id(無(wú)量綱)、路段編號(hào)link_id(無(wú)量綱)、路段長(zhǎng)度ll(單位為m)和路段實(shí)用通行能力cl(單位為pcu/h)；(i)獲取用于交通分配的初始路權(quán)表w，其方法為先獲取步驟(a)中各路段編號(hào)link_id所對(duì)應(yīng)的設(shè)計(jì)速度sd(單位為km/h)和最高限速sl(單位為km/h)，然后采用公式(1)計(jì)算得到各路段編號(hào)link_id所對(duì)應(yīng)的自由流行程時(shí)間tl0(單位為s)，這些路段編號(hào)link_id所對(duì)應(yīng)的行程時(shí)間即為初始路權(quán)表w：tl0＝3600×{0.001×ll/min{sd,sl}}(1)；(j)獲取用于城市交通分配的交通小區(qū)信息，包含每一個(gè)交通小區(qū)的編號(hào)cell_id(無(wú)量綱)、每一個(gè)交通小區(qū)所包含的節(jié)點(diǎn)編號(hào)node_id和路段編號(hào)link_id；(k)獲取最新的城市各出行方式od分布矩陣；(l)獲取用于測(cè)試的區(qū)域機(jī)動(dòng)車尾號(hào)限行方案，包含工作日中每日限行時(shí)間內(nèi)受限制的機(jī)動(dòng)車尾號(hào)個(gè)數(shù)nl(單位為個(gè)，取值為1、2、……、9)、尾號(hào)限行方案所對(duì)應(yīng)的步驟(d)中的出行方式(受尾號(hào)限行的出行方式)、尾號(hào)限行區(qū)域所完全覆蓋的步驟(c)中交通小區(qū)(受尾號(hào)限行小區(qū))的編號(hào)cell_id；或者：獲取用于測(cè)試的城市部分主干道機(jī)動(dòng)車尾號(hào)限行方案，包含工作日中每日限行時(shí)間內(nèi)受限制的機(jī)動(dòng)車尾號(hào)個(gè)數(shù)nl(單位為個(gè)，取值為1、2、……、9)、尾號(hào)限行方案所對(duì)應(yīng)的步驟(d)中的受尾號(hào)限行的出行方式、尾號(hào)限行方案所覆蓋的受尾號(hào)限行路段的編號(hào)link_id”；(m)根據(jù)步驟(e)中確定的受限制機(jī)動(dòng)車尾號(hào)個(gè)數(shù)nl確定第一次和第二次交通分配分別所用od分布矩陣的比例p1和p2，取值規(guī)則如公式(2)：(n)采用尾號(hào)限行條件下的多路徑-容量限制交通分配方法進(jìn)行迭代交通分配；(i)輸出各路段所分配的機(jī)動(dòng)車交通量vl。所述的一種尾號(hào)限行條件下的交通分配修正方法，步驟(d)中所述的城市各出行方式od分布矩陣，包括：私人小客車出行方式od分布矩陣carod(單位為pcu/h)、常規(guī)公交出行方式od分布矩陣busod(單位為pcu/h)、軌道交通出行方式od分布矩陣railod(單位為人/h)、摩托車出行方式od分布矩陣motorod(單位為pcu/h)、出租車出行方式od分布矩陣taxiod(單位為pcu/h)、非機(jī)動(dòng)車出行方式od分布矩陣bikeod(單位為人/h)和步行出行方式od分布矩陣walkod(單位為人/h)中的幾種。所述的一種尾號(hào)限行條件下的交通分配修正方法，所述步驟(f)中將受尾號(hào)限行的出行方式od分布矩陣拆分為不受限尾號(hào)部分和受限尾號(hào)部分兩份，并進(jìn)行兩次交通分配，第一次分配不受限尾號(hào)部分od分布矩陣，第二次分配受限尾號(hào)od分布矩陣。所述的一種尾號(hào)限行條件下的交通分配修正方法，所述步驟(g)中采用的尾號(hào)限行條件下的多路徑-容量限制交通分配方法，其迭代交通分配過(guò)程包含以下步驟：(g12)交通分配迭代次數(shù)記為k，第k(k＝1,2,……k)次迭代加載的od量比例為qk，且令k＝1；(g13)將步驟(d)中各出行方式的od分布矩陣分別乘以qk，得到用于第k次迭代交通分配的od分布矩陣，開始第k次迭代；(g14)對(duì)于步驟(e)中獲取的受尾號(hào)限行的出行方式，將其對(duì)應(yīng)的步驟(g2)中得到的od分布矩陣乘以p1得到用于第一次交通分配的od分布矩陣(不受限尾號(hào)部分)，將其對(duì)應(yīng)的步驟(g2)中得到的od分布矩陣乘以p2之后再進(jìn)行od分布矩陣轉(zhuǎn)移修正，即得到用于第二次交通分配的od分布矩陣(受限尾號(hào)部分)；或者：對(duì)于步驟(e)中獲取的受尾號(hào)限行的出行方式，將其對(duì)應(yīng)的步驟(g2)中得到的od分布矩陣乘以p1得到用于第一次交通分配的od分布矩陣(不受限尾號(hào)部分)，將其對(duì)應(yīng)的步驟(g2)中得到的od分布矩陣乘以p2得到用于第二次交通分配的od分布矩陣(受限尾號(hào)部分)；(g15)對(duì)于步驟(d)中不受尾號(hào)限行的出行方式，利用步驟(g3)中od分布矩陣轉(zhuǎn)移修正時(shí)產(chǎn)生的轉(zhuǎn)移交通量更新不受尾號(hào)限行的出行方式od分布矩陣，結(jié)合路權(quán)表w和更新后不受尾號(hào)限行的出行方式od分布矩陣進(jìn)行第k次迭代交通分配，并將分配后的交通量累計(jì)加載到對(duì)應(yīng)路段機(jī)動(dòng)車交通量vl(單位為pcu/h)上；(g16)根據(jù)步驟(g4)中得到的各路段機(jī)動(dòng)車交通量vl，采用美國(guó)聯(lián)邦公路局路阻函數(shù)模型更新路權(quán)表w中各路段的行程時(shí)間tl(單位為s)，即公式(3)：tl＝tl0[1+α(vl/cl)β](3)其中α＝0.15，β＝4；(g17)采用步驟(g3)中得到的不受限尾號(hào)部分od分布矩陣和步驟(g5)更新過(guò)后的路權(quán)表w進(jìn)行第一次交通分配，并將分配后的交通量累計(jì)加載到對(duì)應(yīng)路段機(jī)動(dòng)車交通量vl上；(g18)根據(jù)步驟(g6)中得到的各路段機(jī)動(dòng)車交通量vl，采用公式(3)更新路權(quán)表w中各路段的行程時(shí)間tl；(g19)將受尾號(hào)限行小區(qū)的編號(hào)cell_id與步驟(c)中得到的交通小區(qū)信息匹配，得到這些交通小區(qū)所包含的所有路段的路段編號(hào)link_id，并將步驟(g7)中得到的路權(quán)表w中這些路段編號(hào)對(duì)應(yīng)的行程時(shí)間tl的值更改為不小于99999的任意數(shù)值，得到尾號(hào)限行條件下的路權(quán)表wl；或者：將步驟(g7)中得到的路權(quán)表w中受尾號(hào)限行路段編號(hào)link_id對(duì)應(yīng)的行程時(shí)間tl的值更改為不小于99999的任意數(shù)值，得到尾號(hào)限行條件下的路權(quán)表wl；(g20)采用步驟(g3)中得到的受限尾號(hào)部分od分布矩陣和步驟(g8)中得到的尾號(hào)限行條件下的路權(quán)表wl進(jìn)行第二次交通分配，并將分配后的交通量累計(jì)加載到對(duì)應(yīng)路段機(jī)動(dòng)車交通量vl上；(g21)根據(jù)步驟(g9)中得到的各路段機(jī)動(dòng)車交通量vl，采用公式(3)更新路權(quán)表w中各路段的行程時(shí)間tl；(g22)令k＝k+1，如果k≤k，則返回步驟(g2)，否則迭代交通分配結(jié)束。所述的一種尾號(hào)限行條件下的交通分配修正方法，所述步驟(g3)中進(jìn)行od分布矩陣轉(zhuǎn)移修正的過(guò)程包含以下步驟：(g3d)記交通發(fā)生小區(qū)i與交通吸引小區(qū)j之間的分布交通量為qij(i∈{cell_id}，j∈{cell_id})，記qij的轉(zhuǎn)移量為qtij，對(duì)于受尾號(hào)限行的小區(qū)編號(hào)cell_id所對(duì)應(yīng)的所有交通發(fā)生量和吸引量，先令qtij＝qij，再令qij＝0；(g3e)對(duì)于od分布矩陣中的其余qij，計(jì)算尾號(hào)限行條件下的繞行率γij；(g3f)結(jié)合od量轉(zhuǎn)移函數(shù)qtij＝f(γij,qij)，計(jì)算其余qij的轉(zhuǎn)移量qtij，再令qij＝qij-qtij，即完成od分布矩陣轉(zhuǎn)移修正。所述的一種尾號(hào)限行條件下的交通分配修正方法，所述步驟(g3b)中計(jì)算尾號(hào)限行條件下的繞行率γij的步驟如下：(g3b1)將受尾號(hào)限行小區(qū)的編號(hào)cell_id與權(quán)利要求1所述步驟(c)中得到的交通小區(qū)信息匹配，得到這些交通小區(qū)所包含的所有路段的路段編號(hào)link_id，并將權(quán)利要求4中當(dāng)前步驟下的最新路權(quán)表w中這些路段編號(hào)對(duì)應(yīng)的行程時(shí)間tl的值更改為不小于99999的任意數(shù)值，得到尾號(hào)限行條件下的路權(quán)表wl；(g3b2)分別根據(jù)路權(quán)表w和尾號(hào)限行條件下的路權(quán)表wl計(jì)算od對(duì)ij之間的最短路權(quán)dij和dlij；(g3b3)采用公式(4)計(jì)算繞行率γij：γij＝dlij/dij(4)。所述的一種尾號(hào)限行條件下的交通分配修正方法，所述步驟(g3c)中所述od量轉(zhuǎn)移函數(shù)qtij＝f(γij,qij)的表達(dá)式為公式(5)：其中β1、β0為待標(biāo)定系數(shù)，γ0為指定置信水平下所有交通量均不更改出行方式的最大轉(zhuǎn)移率，γ1為指定置信水平下所有交通量均更改出行方式的最小轉(zhuǎn)移率，γ0和γ1可通過(guò)問(wèn)卷調(diào)查直接獲取，采用公式(6)即可標(biāo)定β1和β0：所述的一種尾號(hào)限行條件下的交通分配修正方法，所述步驟(g4)中更新不受尾號(hào)限行的出行方式od分布矩陣的方法是：將得到的所有轉(zhuǎn)移交通量qtij按照不受尾號(hào)限行的出行方式od分布矩陣中同一od對(duì)ij的人次交通量比例進(jìn)行拆分并各自疊加進(jìn)去，再轉(zhuǎn)換為各自的原單位交通量即得到更新后不受尾號(hào)限行的出行方式od分布矩陣。所述的一種尾號(hào)限行條件下的交通分配修正方法，私人小客車、常規(guī)公交、摩托車和出租車這四種出行方式的人次交通量與原單位交通量相互轉(zhuǎn)換的轉(zhuǎn)換系數(shù)分別為θcar＝1.3、θbus＝40.0、θmotor＝1.0、θtaxi＝1.3，原單位交通量乘以轉(zhuǎn)換系數(shù)即得到人次交通量，人次交通量除以轉(zhuǎn)換系數(shù)即變?yōu)樵瓎挝唤煌?；軌道交通、非機(jī)動(dòng)車和步行這三種出行方式的人次交通量即為原單位交通量，無(wú)需轉(zhuǎn)換。有益效果：本發(fā)明提出的一種尾號(hào)限行條件下的交通分配修正方法，通過(guò)將受尾號(hào)限制的出行方式od分布矩陣拆分為不受限尾號(hào)部分和受限尾號(hào)部分，并分別進(jìn)行交通分配，將區(qū)域或者部分城市主干道的尾號(hào)限行方案從od分布矩陣、路權(quán)表和分配過(guò)程三方面很好地整合到交通分配的相關(guān)流程中，既保證了交通分配方法的易于計(jì)算性，又切實(shí)、準(zhǔn)確地考慮了尾號(hào)限行方案的作用對(duì)象、作用時(shí)間和作用空間。本發(fā)明充分考慮了我國(guó)各大城市現(xiàn)有的機(jī)動(dòng)車尾號(hào)限行方案，對(duì)交通分配方法進(jìn)行修正，彌補(bǔ)了尾號(hào)限行條件下交通分配修正方法的空白，對(duì)城市虛擬交通測(cè)試平臺(tái)中準(zhǔn)確、快速實(shí)現(xiàn)機(jī)動(dòng)車尾號(hào)限行策略實(shí)施效果的量化測(cè)試具有較大意義。附圖說(shuō)明圖1為本發(fā)明方法的總體流程圖。具體實(shí)施方式下面結(jié)合具體實(shí)施方式，進(jìn)一步闡明本發(fā)明，應(yīng)理解下述具體實(shí)施方式僅用于說(shuō)明本發(fā)明而不用于限制本發(fā)明的范圍。實(shí)施例1：本實(shí)施例中所述的方法可以實(shí)現(xiàn)區(qū)域機(jī)動(dòng)車尾號(hào)限行方案實(shí)施，具體步驟是：如圖1所示為本發(fā)明的總體流程圖，下面結(jié)合實(shí)施例進(jìn)一步說(shuō)明。(a)獲取最新的城市基礎(chǔ)路網(wǎng)信息，包括交叉口編號(hào)node_id(無(wú)量綱)、路段編號(hào)link_id(無(wú)量綱)、路段長(zhǎng)度ll(單位為m)和路段實(shí)用通行能力cl(單位為pcu/h)；實(shí)施例選取一個(gè)簡(jiǎn)化的城市，其基礎(chǔ)路網(wǎng)信息如下：交叉口信息:node_id＝1,2，……，24；路段信息：見表1。表1實(shí)施例中所獲取的路段信息(b)獲取用于交通分配的初始路權(quán)表w，其方法為先獲取步驟(a)中各路段編號(hào)link_id所對(duì)應(yīng)的設(shè)計(jì)速度sd(單位為km/h)和最高限速sl(單位為km/h)，然后采用公式(1)計(jì)算得到各路段編號(hào)link_id所對(duì)應(yīng)的自由流行程時(shí)間tl0(單位為s)，這些路段編號(hào)link_id所對(duì)應(yīng)的行程時(shí)間即為初始路權(quán)表w。本實(shí)施例中本步驟所需的數(shù)據(jù)及采用公式(1)計(jì)算所得的自由流行程時(shí)間tl0見表1，則實(shí)施例的初始路權(quán)表w見表2。表2本實(shí)施例的初始路權(quán)表w路段編號(hào)link_id路權(quán)(s)128218314424…………833…………1146…………7634(c)獲取用于城市交通分配的交通小區(qū)信息，包含每一個(gè)交通小區(qū)的編號(hào)cell_id(無(wú)量綱)、每一個(gè)交通小區(qū)所包含的節(jié)點(diǎn)編號(hào)node_id和路段編號(hào)link_id；本實(shí)施例中的交通小區(qū)信息見表3。表3本實(shí)施例的交通小區(qū)信息(d)獲取最新的城市各出行方式od分布矩陣；本實(shí)施例中出行方式選擇私人小客車、常規(guī)公交和非機(jī)動(dòng)車三種，最新的城市各出行方式od分布矩陣見表4-6。表4本實(shí)施例中的carod(pcu/h)cell_id12345交通發(fā)生量合計(jì)1300010005000100010001100021000300050001000100011000350005000100050005000210004100010005000300010001100051000100050001000300011000交通吸引量合計(jì)110001100021000110001050065000表5本實(shí)施例中的busod(pcu/h)cell_id12345交通發(fā)生量合計(jì)17030703020220230707020302203707010070703804302070703022052030703070220交通吸引量合計(jì)2202203802202201260表6本實(shí)施例中的bikeod(人/h)(e)獲取用于測(cè)試的區(qū)域機(jī)動(dòng)車尾號(hào)限行方案，包含工作日中每日限行時(shí)間內(nèi)受限制的機(jī)動(dòng)車尾號(hào)個(gè)數(shù)nl(單位為個(gè)，取值為1、2、……、9)、尾號(hào)限行方案所對(duì)應(yīng)的步驟(d)中的出行方式(受尾號(hào)限行的出行方式)、尾號(hào)限行區(qū)域所完全覆蓋的步驟(c)中交通小區(qū)(受尾號(hào)限行小區(qū))的編號(hào)cell_id；本實(shí)施例中的區(qū)域機(jī)動(dòng)車尾號(hào)限行方案為交通小區(qū)3在工作日07:00～09:00、17:00～19:00針對(duì)私人小客車限行2個(gè)尾號(hào)，即nl＝2、受尾號(hào)限行的出行方式為私人小客車、受尾號(hào)限行小區(qū)編號(hào)cell_id為3。(f)根據(jù)步驟(e)中確定的受限制機(jī)動(dòng)車尾號(hào)個(gè)數(shù)nl確定第一次和第二次交通分配分別所用od分布矩陣的比例p1和p2，取值規(guī)則如公式(2)；本實(shí)施例中p1＝80％、p2＝20％。(g)采用多路徑-容量限制交通分配方法進(jìn)行交通分配，具體步驟如下：(g1)交通分配迭代次數(shù)記為k，第k(k＝1,2,……k)次迭代加載的od量比例為qk，且令k＝1；本實(shí)施例中k＝5，q1＝30％、q2＝25％、q3＝20％、q4＝15％、q5＝10％。(g2)將步驟(d)中各出行方式的od分布矩陣分別乘以qk，得到用于第k次迭代交通分配的od分布矩陣，開始第k次迭代；本實(shí)施例用于第1次迭代的od分布矩陣見表7-9。表7第1次迭代的carod(pcu/h)cell_id12345交通發(fā)生量合計(jì)190030015003003003300230090015003003003300315001500300150015006300430030015009003003300530030015003009003300交通吸引量合計(jì)3300330063003300330019500表8第1次迭代的busod(pcu/h)cell_id12345交通發(fā)生量合計(jì)121921966629212169663212130212111449621219665692192166交通吸引量合計(jì)66661146666378表9第1次迭代的bikeod(人/h)cell_id12345交通發(fā)生量合計(jì)16001501501501810682150600150181501068315015015015015075041501815060015010685181501501506001068交通吸引量合計(jì)10681068750106810685022(g3)對(duì)于步驟(e)中獲取的受尾號(hào)限行的出行方式，將其對(duì)應(yīng)的步驟(g2)中得到的od分布矩陣乘以p1得到用于第一次交通分配的od分布矩陣(不受限尾號(hào)部分)，將其對(duì)應(yīng)的步驟(g2)中得到的od分布矩陣乘以p2之后再進(jìn)行od分布矩陣轉(zhuǎn)移修正，即得到用于第二次交通分配的od分布矩陣(受限尾號(hào)部分)；本實(shí)施例中對(duì)于carod,不受限尾號(hào)部分od分布矩陣見表10,乘以p2之后得到的用于轉(zhuǎn)移修正的od分布矩陣分別見表11；表10不受限尾號(hào)部分carod(pcu/h)cell_id12345交通發(fā)生量合計(jì)172024012002402402640224072012002402402640312001200240120012005040424024012007202402640524024012002407202640交通吸引量合計(jì)2640264050402640264015600表11用于轉(zhuǎn)移修正的carod(pcu/h)cell_id12345交通發(fā)生量合計(jì)11806030060606602601803006060660330030060300300126046060300180606605606030060180660交通吸引量合計(jì)66066012606606603900(g3a)記交通發(fā)生小區(qū)i與交通吸引小區(qū)j之間的分布交通量為qij(i∈{cell_id}，j∈{cell_id})，記qij的轉(zhuǎn)移量為qtij，對(duì)于受尾號(hào)限行的小區(qū)編號(hào)cell_id所對(duì)應(yīng)的所有交通發(fā)生量和吸引量，先令qtij＝qij，再令qij＝0；本實(shí)施例中，對(duì)于表11中的交通小區(qū)3，先令qt1*＝q1*、qt2*＝q2*、qt*1＝q*1、qt*2＝q*2，再令q1*＝0、q2*＝0、q*1＝0、q*1＝0、，此處下標(biāo)中“*”依次取1、2、……、5；(g3b)對(duì)于od分布矩陣中的其余qij，計(jì)算尾號(hào)限行條件下的繞行率γij；(g3b1)將受尾號(hào)限行小區(qū)的編號(hào)cell_id與權(quán)利要求1所述步驟(c)中得到的交通小區(qū)信息匹配，得到這些交通小區(qū)所包含的所有路段的路段編號(hào)link_id，并將權(quán)利要求3中當(dāng)前步驟下的最新路權(quán)表w中這些路段編號(hào)對(duì)應(yīng)的行程時(shí)間tl的值更改為不小于99999的任意數(shù)值，得到尾號(hào)限行條件下的路權(quán)表wl；本實(shí)施例第1次迭代過(guò)程中，本步驟的當(dāng)前最新路權(quán)表w即為表2，將表2中小區(qū)3所包含的link_id的路權(quán)改為99999，得到wl。(g3b2)分別根據(jù)路權(quán)表w和尾號(hào)限行條件下的路權(quán)表wl計(jì)算od對(duì)ij之間的最短路權(quán)dij和dlij；本實(shí)施例中當(dāng)前步驟下，根據(jù)表2計(jì)算od對(duì)之間的最短路權(quán)見表12，根據(jù)路權(quán)表wl計(jì)算od對(duì)之間的最短路權(quán)見表13：表12最短路權(quán)dij(s)cell_id123451236418——4185682418236——5684183——————————4418568——2364185568418——418236表13最短路權(quán)dlij(s)cell_id123451236418——41814122418236——14124183——————————44181412——23641851412418——418236(g3b3)采用公式(4)計(jì)算繞行率γij；本實(shí)施例中當(dāng)前步驟下得到的繞行率見表14：表14繞行率γijcell_id1234511.01.0——1.02.521.01.0——2.51.03——————————41.02.5——1.01.052.51.0——1.01.0(g3c)od量轉(zhuǎn)移函數(shù)qtij＝f(γij,qij)，計(jì)算其余qij的轉(zhuǎn)移量qtij，再令qij＝qij-qtij，即完成od分布矩陣轉(zhuǎn)移修正；步驟(g3c)中所述od量轉(zhuǎn)移函數(shù)qtij＝f(γij,qij)的表達(dá)式為公式(5)：其中β1、β0為待標(biāo)定系數(shù)，γ0為指定置信水平下所有交通量均不更改出行方式的最大轉(zhuǎn)移率，γ1為指定置信水平下所有交通量均更改出行方式的最小轉(zhuǎn)移率，γ0和γ1可通過(guò)問(wèn)卷調(diào)查直接獲取，采用公式(6)即可標(biāo)定β1和β0：本實(shí)施例本步驟中，在0.95的置信水平下獲取到γ0＝1.2、γ1＝3.0，根據(jù)公式(6)和(5)進(jìn)行計(jì)算，結(jié)果見表15、16：表15受限尾號(hào)部分carod(pcu/h)cell_id1234511806006017260180017603000004601701806051760060180表16carod的轉(zhuǎn)移交通量(pcu/h)cell_id12345100600432003043033003006030030040431800054306000(g4)對(duì)于步驟(d)中不受尾號(hào)限行的出行方式，利用步驟(g3)中od分布矩陣轉(zhuǎn)移修正時(shí)產(chǎn)生的轉(zhuǎn)移交通量更新不受尾號(hào)限行的出行方式od分布矩陣，結(jié)合路權(quán)表w和更新后不受尾號(hào)限行的出行方式od分布矩陣進(jìn)行第k次迭代交通分配，并將分配后的交通量累計(jì)加載到對(duì)應(yīng)路段機(jī)動(dòng)車交通量vl(單位為pcu/h)上；本實(shí)施例中將表16中各od對(duì)之間的轉(zhuǎn)移交通量，按照表8和表9同一od對(duì)之間的人次交通量比例進(jìn)行拆分并各自疊加進(jìn)去，再轉(zhuǎn)換為各自的原單位交通量，例如對(duì)于qt15＝43pcu/h，表8和表9中同一od對(duì)之間的人次交通量比例為(6×40＝240)：18＝40：3，則將40人次的轉(zhuǎn)移交通量轉(zhuǎn)換為1pcu/h的常規(guī)公交出行交通量疊加至表8的q15中，將3人次的轉(zhuǎn)移交通量即3人次的非機(jī)動(dòng)車出行交通量疊加至表9的q15中。更新后的busod和bikeod分別見表17和表18。表17更新后的busodcell_id12345交通發(fā)生量合計(jì)121921976729212179673212130212111449721219675792192167交通吸引量合計(jì)67671146767382表18更新后的bikeodcell_id12345交通發(fā)生量合計(jì)16001501501502110712150600150211501071315015015015015075041502115060015010715211501501506001071交通吸引量合計(jì)10711071750107110715034當(dāng)前步驟中，利用表2和表17進(jìn)行機(jī)動(dòng)車交通量的第1次迭代分配，并將分配后的交通量累計(jì)加載到對(duì)應(yīng)路段機(jī)動(dòng)車交通量vl上。(g5)根據(jù)(g4)中得到的各路段機(jī)動(dòng)車交通量vl，采用公式(3)更新路權(quán)表w中各路段的行程時(shí)間tl；本實(shí)例本步驟得到的更新后路權(quán)表w見表19。表19路權(quán)表w(g6)采用步驟(g3)中得到的不受限尾號(hào)部分od分布矩陣和步驟(g5)更新過(guò)后的路權(quán)表w進(jìn)行第一次交通分配，并將分配后的交通量累計(jì)加載到對(duì)應(yīng)路段機(jī)動(dòng)車交通量vl上；本實(shí)施例本步驟中，采用表10和表19進(jìn)行第一次交通分配，并將分配后的交通量累計(jì)加載到對(duì)應(yīng)路段機(jī)動(dòng)車交通量vl上。(g7)根據(jù)步驟(g6)中得到的各路段機(jī)動(dòng)車交通量vl，采用公式(3)更新路權(quán)表w中各路段的行程時(shí)間tl；(g8)將受尾號(hào)限行小區(qū)的編號(hào)cell_id與步驟(c)中得到的交通小區(qū)信息匹配，得到這些交通小區(qū)所包含的所有路段的路段編號(hào)link_id，并將步驟(g7)中得到的路權(quán)表w中這些路段編號(hào)對(duì)應(yīng)的行程時(shí)間tl的值更改為不小于99999的任意數(shù)值，得到尾號(hào)限行條件下的路權(quán)表wl；當(dāng)前步驟得到的尾號(hào)限行條件下的路權(quán)表wl見表20。表20路權(quán)表wl路段編號(hào)link_id路權(quán)(s)130219315426…………899999…………1199999…………7637(g9)采用步驟(g3)中得到的受限尾號(hào)部分od分布矩陣和步驟(g8)中得到的尾號(hào)限行條件下的路權(quán)表wl進(jìn)行第二次交通分配，并將分配后的交通量累計(jì)加載到對(duì)應(yīng)路段機(jī)動(dòng)車交通量vl上；本實(shí)施例當(dāng)前步驟中采用表15和表20進(jìn)行交通分配，并將分配后的交通量累計(jì)加載到對(duì)應(yīng)路段機(jī)動(dòng)車交通量vl上。(g10)根據(jù)步驟(g9)中得到的各路段機(jī)動(dòng)車交通量vl，采用公式(3)更新路權(quán)表w中各路段的行程時(shí)間tl；本實(shí)施例當(dāng)前步驟中得到的更新后路權(quán)表w見表21。表21路權(quán)表w路段編號(hào)link_id路權(quán)(s)135222317430…………842…………1158…………7641(g11)令k＝k+1，如果k≤k，則返回步驟(g2)，否則迭代交通分配結(jié)束；令k依次取2、3、4、5，返回步驟(g2)繼續(xù)第2、3、4、5次迭代交通分配，完成第5次迭代交通分配后，迭代結(jié)束。(h)輸出各路段所分配的機(jī)動(dòng)車交通量vl。本實(shí)施例中通過(guò)5次迭代交通分配，得到最終各路段所分配的機(jī)動(dòng)車交通量vl。本實(shí)例中，可以看出，通過(guò)本發(fā)明方法的八個(gè)步驟，可以簡(jiǎn)單快速地完成區(qū)域機(jī)動(dòng)車尾號(hào)限行條件下的分方式od分布矩陣修正及交通分配過(guò)程，該修正方法的邏輯性強(qiáng)、易于程序化實(shí)現(xiàn)。實(shí)施例2：本實(shí)施例與實(shí)施例1的不同之處在于，本實(shí)施例中所述的方法，可以實(shí)現(xiàn)僅在城市部分主干道實(shí)施的機(jī)動(dòng)車尾號(hào)限行方案，具體是：所述步驟(e)為獲取用于測(cè)試的城市部分主干道機(jī)動(dòng)車尾號(hào)限行方案，包含工作日中每日限行時(shí)間內(nèi)受限制的機(jī)動(dòng)車尾號(hào)個(gè)數(shù)nl(單位為個(gè)，取值為1、2、……、9)、尾號(hào)限行方案所對(duì)應(yīng)的步驟(d)中的出行方式(受尾號(hào)限行的出行方式)、尾號(hào)限行方案所覆蓋的路段(受尾號(hào)限行路段)的編號(hào)link_id；本實(shí)施例中的城市部分主干道機(jī)動(dòng)車尾號(hào)限行方案為路段4、8、42、46在工作日07:00～09:00、17:00～19:00針對(duì)私人小客車限行2個(gè)尾號(hào)，即nl＝2，受尾號(hào)限行的出行方式為私人小客車，受尾號(hào)限行路段的編號(hào)link_id取值為4、8、42、46。所述步驟(g3)為“對(duì)于步驟(e)中獲取的受尾號(hào)限行的出行方式，將其對(duì)應(yīng)的步驟(g2)中得到的od分布矩陣乘以p1得到用于第一次交通分配的od分布矩陣(不受限尾號(hào)部分)，將其對(duì)應(yīng)的步驟(g2)中得到的od分布矩陣乘以p2得到用于第二次交通分配的od分布矩陣(受限尾號(hào)部分)”；本實(shí)施例中對(duì)于carod,不受限尾號(hào)部分od分布矩陣同表10,受限尾號(hào)部分od分布矩陣同表11。所述步驟(g4)為“對(duì)于步驟(d)中不受尾號(hào)限行的出行方式，結(jié)合路權(quán)表w和步驟(g2)中獲得的不受尾號(hào)限行的出行方式od分布矩陣進(jìn)行第k次迭代交通分配，并將分配后的交通量累計(jì)加載到對(duì)應(yīng)路段機(jī)動(dòng)車交通量vl上”；本實(shí)施例當(dāng)前步驟中，利用表2和表8進(jìn)行機(jī)動(dòng)車交通量的第1次迭代分配，并將分配后的交通量累計(jì)加載到對(duì)應(yīng)路段機(jī)動(dòng)車交通量vl上。所述步驟(g8)為“將步驟(g7)中得到的路權(quán)表w中受尾號(hào)限行路段編號(hào)link_id對(duì)應(yīng)的行程時(shí)間tl的值更改為不小于99999的任意數(shù)值，得到尾號(hào)限行條件下的路權(quán)表wl”。當(dāng)前步驟得到的尾號(hào)限行條件下的路權(quán)表wl見表21。表21路權(quán)表wl路段編號(hào)link_id路權(quán)(s)131220315499999…………899999…………1148…………7636其余步驟均與實(shí)施例1相同。本實(shí)例中，可以看出，通過(guò)本發(fā)明方法同樣可以簡(jiǎn)單快速地完成僅在城市部分主干道實(shí)施機(jī)動(dòng)車尾號(hào)限行方案的條件下的交通分配過(guò)程，此時(shí)無(wú)需對(duì)分方式od分布矩陣進(jìn)行修正，該修正方法的邏輯性強(qiáng)、易于程序化實(shí)現(xiàn)。應(yīng)當(dāng)指出，上述實(shí)施實(shí)例僅僅是為清楚地說(shuō)明所作的舉例，而并非對(duì)實(shí)施方式的限定，這里無(wú)需也無(wú)法對(duì)所有的實(shí)施方式予以窮舉。本實(shí)施例中未明確的各組成部分均可用現(xiàn)有技術(shù)加以實(shí)現(xiàn)。對(duì)于本
技術(shù)領(lǐng)域：
的普通技術(shù)人員來(lái)說(shuō)，在不脫離本發(fā)明原理的前提下，還可以做出若干改進(jìn)和潤(rùn)飾，這些改進(jìn)和潤(rùn)飾也應(yīng)視為本發(fā)明的保護(hù)范圍。當(dāng)前第1頁(yè)12