, ,

(1.五邑大學(xué) 信息工程學(xué)院, 廣東 江門(mén) 529020;2.電子科技大學(xué)中山學(xué)院 機(jī)電工程學(xué)院, 廣東 中山 528400)

0 引言

社會(huì)經(jīng)濟(jì)的快速增長(zhǎng)導(dǎo)致汽車(chē)數(shù)量劇增，高速公路擁堵時(shí)有發(fā)生，擁堵降低了高速公路主線(xiàn)的通行能力，增加了車(chē)輛延誤時(shí)間，增大了燃油消耗，加大了環(huán)境污染，造成了巨大的經(jīng)濟(jì)損失。為了解決高速公路擁堵，一方面必須投入更多的建設(shè)資金、修建更多的道路；另一方面必須研究新技術(shù)，用先進(jìn)的方法對(duì)高速公路進(jìn)行控制和管理，從而充分利用現(xiàn)有的道路資源、減少交通擁堵。高速公路入口匝道控制是解決主線(xiàn)交通擁堵的有效辦法，其主要目標(biāo)是調(diào)節(jié)進(jìn)入高速公路的車(chē)輛數(shù)目使得主線(xiàn)的交通需求不超過(guò)交通容量[1]。匝道控制方法根據(jù)控制對(duì)象不同可分為單匝道控制[2-3]和多匝道協(xié)調(diào)控制[4-5]，其中多匝道協(xié)調(diào)控制綜合考慮了高速公路不同路段之間的交互影響，成為匝道控制方法的主流發(fā)展方向。

目前已有多種多匝道協(xié)調(diào)控制方法，相關(guān)研究也取得了較好的結(jié)果，但現(xiàn)場(chǎng)應(yīng)用還需加強(qiáng)[6]。文獻(xiàn)[7]采用定量分層模型研究快速路多匝道協(xié)調(diào)控制，并根據(jù)交通狀態(tài)和臨界密度的差值比確定匝道控制次序，該方法具有一定的優(yōu)越性，但可能出現(xiàn)某個(gè)入口匝道車(chē)輛排隊(duì)長(zhǎng)度較長(zhǎng)。文獻(xiàn)[8]利用動(dòng)態(tài)模糊神經(jīng)網(wǎng)絡(luò)研究快速路匝道控制，該方法能使主線(xiàn)密度不超過(guò)臨界密度，但文中只對(duì)單匝道控制進(jìn)行了仿真，未對(duì)多匝道控制進(jìn)行相應(yīng)的分析。文獻(xiàn)[9]采用模糊控制研究快速路匝道控制，該方法對(duì)主線(xiàn)流量控制具有較好的效果，然而模糊邏輯的隸屬度函數(shù)和規(guī)則庫(kù)的選取依賴(lài)于操作者的經(jīng)驗(yàn)知識(shí)。文獻(xiàn)[10]采用最優(yōu)控制方法來(lái)減少高速公路擁擠，最優(yōu)控制在實(shí)際應(yīng)用中效果較好。文獻(xiàn)[11]研究了基于迭代學(xué)習(xí)的高速公路多匝道調(diào)節(jié)，但是迭代學(xué)習(xí)依賴(lài)于重復(fù)性交通條件，當(dāng)某天的交通異常不滿(mǎn)足重復(fù)性時(shí)，控制效果會(huì)變差。近年來(lái)隨著智能優(yōu)化算法的廣泛應(yīng)用，匝道控制算法已經(jīng)向智能優(yōu)化算法發(fā)展，智能優(yōu)化算法用于交通控制的參數(shù)優(yōu)化或優(yōu)化控制取得了較好的效果。文獻(xiàn)[12]利用遺傳算法對(duì)匝道控制多個(gè)性能指標(biāo)進(jìn)行綜合優(yōu)化，效果良好。除了多匝道協(xié)調(diào)控制外，還有多交叉口協(xié)調(diào)控制[13]。

從以上分析可以看出，目前多匝道協(xié)調(diào)控制的研究大多局限于理論層面，實(shí)用的技術(shù)較少。針對(duì)上述問(wèn)題，本文采用廣泛應(yīng)用的比例積分微分(Proportional-Integral-Derivative, PID)控制器，并將智能優(yōu)化算法與大系統(tǒng)分層控制策略相結(jié)合。大系統(tǒng)分層控制將整個(gè)控制任務(wù)分解成多層來(lái)實(shí)現(xiàn)，具有意義明確、易于實(shí)現(xiàn)的特點(diǎn)。多匝道控制系統(tǒng)分為協(xié)調(diào)控制層和直接控制層，直接控制層采用PID控制，粒子群算法用來(lái)對(duì)直接控制層的PID參數(shù)進(jìn)行優(yōu)化，最后通過(guò)仿真實(shí)例進(jìn)行了控制效果和控制性能的驗(yàn)證。

1 高速公路多匝道系統(tǒng)的原理模型

1.1 多匝道系統(tǒng)的結(jié)構(gòu)模型以及入口匝道流量控制的難點(diǎn)

圖1是高速公路多匝道系統(tǒng)的結(jié)構(gòu)模型，該高速公路有M個(gè)路段，每個(gè)路段至多含有一個(gè)入口匝道和一個(gè)出口匝道(某個(gè)路段允許不含有入口匝道或出口匝道)，rj表示從路段j入口匝道進(jìn)入到主線(xiàn)的交通流量，εj表示從路段j出口匝道離開(kāi)的交通流量，j-1/2表示路段j-1和路段j的分界處，j+1/2表示路段j和路段j+1的分界處。當(dāng)改變路段j入口匝道流量rj時(shí)，就會(huì)改變路段j的主線(xiàn)流量。由于主線(xiàn)交通流會(huì)向上游或下游傳播，即：當(dāng)交通順暢時(shí)，交通流向下游傳播，下游交通會(huì)受到上游交通的影響；反之，當(dāng)交通擁擠時(shí)，交通流向上游傳播，上游交通會(huì)受到下游交通的影響；因此，改變路段j的主線(xiàn)流量會(huì)影響其它路段的主線(xiàn)流量，從而影響其它路段的入口匝道流量。從以上分析可以看出，rj的改變會(huì)影響其它路段的入口匝道調(diào)節(jié)率，也就是說(shuō)，各個(gè)入口匝道的調(diào)節(jié)存在耦合關(guān)系，是相互關(guān)聯(lián)和相互協(xié)調(diào)的，而不是彼此孤立的。多匝道協(xié)調(diào)控制具有強(qiáng)的耦合性、非線(xiàn)性和時(shí)變性，這是入口匝道流量控制的難點(diǎn)。解決該難點(diǎn)問(wèn)題的一種有效途徑是建立宏觀交通模型，然后采用系統(tǒng)分層和粒子群優(yōu)化的PID控制器對(duì)該交通模型進(jìn)行控制。

圖1 高速公路多匝道系統(tǒng)的結(jié)構(gòu)模型

1.2 宏觀交通模型

交通控制系統(tǒng)中一般采用宏觀交通模型[6]。設(shè)交通流滿(mǎn)足一階LWR(Lighthill-Whitham-Richards)宏觀交通模型：

?q/?x+?p/?t=0

(1)

q=f(p)

(2)

式中,p和q分別代表交通流密度和流量，x是空間長(zhǎng)度，t是時(shí)間。

式(1)是守恒方程，表示道路上的車(chē)輛數(shù)守恒，式(2)稱(chēng)為交通流基本圖，反映現(xiàn)場(chǎng)的交通變量之間的函數(shù)關(guān)系。基本圖滿(mǎn)足：(a) 邊界條件f(0)=f(pjam)=0，其中pjam為阻塞密度；(b) 存在一個(gè)臨界密度pc，對(duì)應(yīng)流量最大值f(pc)。當(dāng)p∈(0,pc)時(shí)，f'(p)>0；當(dāng)p∈(pc,pjam)時(shí)，f'(p)<0，即隨著p增大，q減少，最終出現(xiàn)交通堵塞。

采用式(1)的偏微分方程描述的系統(tǒng)稱(chēng)為分布參數(shù)系統(tǒng)，由有限差分法可得到式(1)的數(shù)值式[1]：

pj(n+1)=pj(n)+

(3)

式中,pj(n),j=1,2,...,M是在時(shí)間步數(shù)n時(shí)第j個(gè)路段的交通密度，rj(n)和εj(n)分別表示在時(shí)間步數(shù)n時(shí)路段j的入口匝道和出口匝道流量。Δt是時(shí)間步長(zhǎng)，lj是第j個(gè)路段的長(zhǎng)度。fj+1/2表示數(shù)值流量，j+1/2表示路段j和路段j+1的分界處。fj+1/2計(jì)算方法如下：

當(dāng)pj+1(n)≠pj(n)時(shí)，

-|μ|j+1/2(n)[pj+1(n)-pj(n)]}

(4)

當(dāng)pj+1(n)=pj(n)時(shí)，

fj+1/2(n)=f[pj+1(n)]=f[pj(n)]

(5)

式(3)中第j個(gè)路段的出口匝道流量εj(n)為該路段流量f[pj(n)]的一部分，即：

εj(n)=sj(n)f[pj(n)],0≤sj(n)≤1

(6)

在有限差分法中，CFL(Courant-Friedichs-Lewy)條件是數(shù)值穩(wěn)定性和收斂性判別條件，當(dāng)式(7)的CFL條件成立時(shí)，式(3)的有限差分模型滿(mǎn)足數(shù)值穩(wěn)定性和收斂性[1, 3]。

(7)

因此，Δt和lj的參數(shù)選取需滿(mǎn)足CFL條件。例如：選取f′[pj(n)]=νf=97.3 km/h，sj(n)=0.25，lj=1 km時(shí)，對(duì)應(yīng)最大的時(shí)間步長(zhǎng)為Δt=29.6 s。CFL條件的含義可理解為：以自由流速度vf行駛的車(chē)輛在一個(gè)時(shí)間步長(zhǎng)內(nèi)不能通過(guò)路段j，即：Δt·vf≤Δt[1+sj(n)]·vf≤lj。

2 多匝道協(xié)調(diào)控制系統(tǒng)的實(shí)現(xiàn)

2.1 協(xié)調(diào)控制系統(tǒng)結(jié)構(gòu)

多匝道協(xié)調(diào)控制系統(tǒng)采用分層控制結(jié)構(gòu)，其系統(tǒng)框圖如圖2所示，系統(tǒng)由協(xié)調(diào)控制層和直接控制層組成，其中協(xié)調(diào)控制層負(fù)責(zé)選取合適的交通模型，并根據(jù)各路段的交通條件確定期望密度軌跡；直接控制層根據(jù)密度偏差改變各入口匝道的調(diào)節(jié)率，使得各路段的實(shí)際密度達(dá)到期望密度值。

協(xié)調(diào)控制層考慮各路段的交通信息以及被控路段與其上游、下游的相互關(guān)聯(lián)與影響，確定被控路段的期望密度。(1) 當(dāng)路段沒(méi)有出現(xiàn)交通擁擠時(shí)，在交通繁忙期間，期望密度設(shè)為臨界密度的負(fù)鄰域，這樣既不會(huì)出現(xiàn)交通擁擠，又可維持較大的道路交通流量。(2) 當(dāng)某個(gè)路段出現(xiàn)交通擁擠時(shí)，該擁擠路段的下游路段，期望密度仍設(shè)為臨界密度的負(fù)鄰域，這樣可以為擁擠路段的車(chē)輛提供較大的出口流量容量，但對(duì)于擁擠路段的上游路段，期望密度應(yīng)該減小，這可理解為通過(guò)降低該路段的上游交通流量來(lái)消除該路段的擁擠。

圖2 多匝道協(xié)調(diào)控制系統(tǒng)框圖

目前，PID控制器在工業(yè)控制領(lǐng)域被廣泛應(yīng)用，其中比例環(huán)節(jié)用于跟蹤控制系統(tǒng)的偏差信號(hào)，從而減少偏差；積分環(huán)節(jié)可以消除靜差；微分環(huán)節(jié)的作用是反映偏差信號(hào)的變化趨勢(shì)，從而提高系統(tǒng)的響應(yīng)速度。由于PID控制有諸多優(yōu)點(diǎn)，因此本文直接控制層采用PID控制。

使用PID控制器關(guān)鍵在于參數(shù)的整定，一般的經(jīng)驗(yàn)試湊法難以獲得滿(mǎn)意的結(jié)果，本文引入粒子群算法來(lái)實(shí)現(xiàn)多匝道PID參數(shù)的優(yōu)化。

2.2 直接控制層的具體實(shí)現(xiàn)

路段j直接控制層的系統(tǒng)框圖如圖3所示。圖3中的交通模型由協(xié)調(diào)控制層負(fù)責(zé)選取，協(xié)調(diào)控制層同時(shí)也決定模型參數(shù)，并確定路段j的期望密度軌跡pdj。圖3的控制器是一種非線(xiàn)性的反饋型控制器，期望密度pdj(n)與實(shí)際密度pj(n)相比較，產(chǎn)生誤差信號(hào)ej(n)，誤差ej(n)、ej(n-1)與ej(n-2)作用于PID控制器，得到uj(n)，即：

Δuj(n)=kPj[ej(n)-ej(n-1)]+kIjej(n)

+kDj[ej(n)-2ej(n-1)+ej(n-2)]

(8)

uj(n)=Δuj(n)+uj(n-1)

(9)

式中,kPj、kIj和kDj分別是第j個(gè)路段PID控制器的比例系數(shù)、積分系數(shù)和微分系數(shù)，其數(shù)值由粒子群算法尋優(yōu)來(lái)得到。PID控制器的輸出uj(n)與路段j入口匝道調(diào)節(jié)率rj(n)相等，即：uj(n)=rj(n)，rj(n)饋入到交通模型中產(chǎn)生實(shí)際的交通密度pj(n+1)，pj(n+1)經(jīng)一步延時(shí)后反饋到輸入端。以上過(guò)程是第j個(gè)路段在時(shí)間步數(shù)n時(shí)的控制過(guò)程。

本文協(xié)調(diào)控制層選擇1.2節(jié)的宏觀交通模型，第j個(gè)路段的PID控制器產(chǎn)生下一步的控制信號(hào)rj(n+1)的具體控制過(guò)程如下：

對(duì)當(dāng)前時(shí)間步數(shù)n，各路段的密度pj(n),j=1,2,...,M已知，由式(2)可求出各路段流量qj(n),j=1,2,...,M，通過(guò)式(4)和式(5)計(jì)算，可得到數(shù)值流量fj+1/2(n),j=1,2,...,M-1，再根據(jù)第n步的控制信號(hào)rj(n)由式(3)和式(6)求出pj(n+1),j=1,2,...,M，因此可求出密度誤差ej(n+1)，即：ej(n+1)=pdj(n+1)-pj(n+1)，由ej(n+1)和以前時(shí)間步數(shù)的誤差ej(n)與ej(n-1)再根據(jù)式(8)和式(9)就可通過(guò)PID控制器得到uj(n+1)，也即得到rj(n+1)。

由以上過(guò)程可知，只要已知各路段的初始密度和初始調(diào)節(jié)率，則各路段的密度曲線(xiàn)和調(diào)節(jié)率曲線(xiàn)便可由上述步驟計(jì)算得到。

圖3 路段j直接控制層的系統(tǒng)框圖

2.3 粒子群算法參數(shù)優(yōu)化

為了找到各個(gè)匝道最佳的PID參數(shù)kPj、kIj和kDj，本文采用改進(jìn)的粒子群優(yōu)化算法。兩個(gè)改進(jìn)措施為：1) 隨著迭代次數(shù)的變化，動(dòng)態(tài)改變慣性權(quán)值，從而提高粒子群優(yōu)化算法的收斂性能；2) 采用粒子群優(yōu)化的模擬退火算法，從而提高粒子群優(yōu)化算法的全局優(yōu)化能力。PID參數(shù)的優(yōu)化過(guò)程如下：

1) 設(shè)置初始參數(shù)，例如：初始慣性權(quán)值、初始退火溫度和初始退火速度。

2) 設(shè)置種群規(guī)模和初始種群范圍，隨機(jī)產(chǎn)生粒子群的初始種群。

3) 定義目標(biāo)函數(shù)和適應(yīng)度函數(shù)。目標(biāo)函數(shù)E由式(10)定義，它是密度誤差的平方和，而適應(yīng)度函數(shù)是目標(biāo)函數(shù)的倒數(shù)。

(10)

4) 計(jì)算每個(gè)粒子的適應(yīng)度函數(shù)，并選擇全局最優(yōu)粒子。

5) 更新每個(gè)粒子的速度和位置。

6) 計(jì)算更新后的粒子群的適應(yīng)度函數(shù)，選擇個(gè)體的最優(yōu)狀態(tài)，并選擇群體的全局最優(yōu)位置。

7) 采用模擬退火算法更新個(gè)體的最佳位置和更新群體的全局最佳位置。

8) 重復(fù)以上步驟直到滿(mǎn)足收斂準(zhǔn)則或者達(dá)到最大的迭代次數(shù)。

9) 輸出各個(gè)匝道最佳的比例系數(shù)、積分系數(shù)和微分系數(shù)。

3 仿真案例

對(duì)圖2所示的多匝道協(xié)調(diào)控制系統(tǒng)進(jìn)行仿真實(shí)驗(yàn)，以一段長(zhǎng)9 km、雙向4車(chē)道的高速公路為仿真對(duì)象，按每段1 km分成9段，假設(shè)第3段、第5段和第7段各含一個(gè)入口匝道，第2段和第5段各含一個(gè)出口匝道。設(shè)主線(xiàn)的上游進(jìn)口流量為0～1 500輛/小時(shí)/車(chē)道，各入口匝道交通需求范圍為0～1 000輛/小時(shí)。協(xié)調(diào)控制層選擇1.2節(jié)的宏觀交通模型，并選擇式(2)的流量-密度關(guān)系方程為Green-Shields提出的關(guān)系式[6]，即：

q=f(p)=vf(p-p2/pjam)

(11)

式中，取vf=97.3千米/小時(shí)，pjam=74.0輛/千米/車(chē)道，對(duì)應(yīng)最大交通容量qm=1 800.0輛/小時(shí)/車(chē)道，臨界密度pc=37.0輛/千米/車(chē)道。

圖4 各路段實(shí)際交通密度變化三維圖

選擇期望密度pd為34.0輛/千米/車(chē)道，9個(gè)路段的初始密度分別設(shè)為16.0，54.0，27.5，21.0，28.0，46.0，25.0，24.0和41.0輛/千米/車(chē)道，其中第2路段、第6路段和第9路段存在交通擁堵。由于第3路段、第5路段和第7路段有入口匝道控制器，因此協(xié)調(diào)控制層必須確定這3個(gè)路段的期望密度軌跡。第3路段的上游路段有初始擁擠，但第3路段的下游路段無(wú)初始擁擠，因此第3路段的期望密度從初始值上升到期望值34.0輛/千米/車(chē)道。同理，第7路段的期望密度從初始值上升到期望值34.0輛/千米/車(chē)道。但是，由于第5路段的

圖5 密度跟蹤曲線(xiàn)

圖6 入口匝道調(diào)節(jié)率變化曲線(xiàn)

下游路段有初始擁擠，因此第5路段的期望密度開(kāi)始時(shí)維持較低的水平，然后再上升到期望值34.0輛/千米/車(chē)道，這樣有利于路段6消除擁擠。

仿真時(shí)間步長(zhǎng)選為Δt=20s，式(6)的出口流量系數(shù)為sj(n)=0～0.25。通過(guò)協(xié)調(diào)控制層得到各路段的期望密度后，由直接控制層實(shí)施PID控制，PID控制器的比例系數(shù)、積分系數(shù)和微分系數(shù)由粒子群算法優(yōu)化得到。

用Matlab軟件編寫(xiě)仿真程序，粒子群算法得到最優(yōu)值為kP3=130.093 6，kI3=275.700 7，kD3=41.210 5，kP5=330.305 5，kI5=355.890 6，kD5=85.436 7，kP7=235.293 1，kI7=316.275 9，kD7=62.184 2。所有路段的實(shí)際交通密度變化三維效果圖如圖4所示，可見(jiàn)，開(kāi)始時(shí)在路段2、路段6和路段9都存在嚴(yán)重?fù)矶?，?jīng)協(xié)調(diào)控制后，路段2和路段9大約過(guò)了8分鐘，路段6大約過(guò)了13分鐘，交通狀態(tài)進(jìn)入正常流量的穩(wěn)定區(qū)域，最后達(dá)到期望密度，消除了擁堵。有入口匝道的路段密度跟蹤曲線(xiàn)如圖5所示，可見(jiàn)，實(shí)際的交通密度可以跟蹤期望的交通密度，說(shuō)明協(xié)調(diào)控制系統(tǒng)具有良好的動(dòng)態(tài)和穩(wěn)態(tài)性能。各入口匝道車(chē)輛調(diào)節(jié)率如圖6所示。

4 結(jié)論

本文采用大系統(tǒng)分層控制策略和粒子群算法對(duì)高速公路多匝道進(jìn)行了協(xié)調(diào)控制，將多匝道控制問(wèn)題按不同的任務(wù)分解成協(xié)調(diào)控制層和直接控制層，前者負(fù)責(zé)模型選擇、參數(shù)調(diào)整和確定期望密度，后者采用PID控制器實(shí)施控制，并引入粒子群算法對(duì)PID控制器的比例系數(shù)、積分系數(shù)和微分系數(shù)進(jìn)行優(yōu)化，以改善控制系統(tǒng)的動(dòng)態(tài)性能?；贛atlab的9個(gè)路段仿真結(jié)果表明，當(dāng)其中3個(gè)路段存在嚴(yán)重?fù)矶聲r(shí)，系統(tǒng)達(dá)到期望密度的平均時(shí)間約為13分鐘，反映了系統(tǒng)具有良好的動(dòng)態(tài)性能，能有效消除交通擁堵和維持主線(xiàn)交通流穩(wěn)定。

[1]李 健. 匝道交通控制理論與方法 [M]. 北京: 北京交通大學(xué)出版社, 2013.

[2]Ali Sercan Kesten, Murat Ergün, Tetsuo Yai. An analysis on efficiency and equity of fixed-time ramp metering [J]. Journal of Transportation Technologies, 2013, 3(2A): 48-56.

[3]Liang X R, Lu Q, Lin P Q, et al. Freeway ramp metering using fuzzy logic and genetic proportional control [J]. Traffic Engineering & Control, 2015, 56(2): 219-225.

[4]Ioannis Papamichail, Apostolos Kotsialos, Ioannis Margonis, et al. Coordinated ramp metering for freeway networks-A model-predictive hierarchical control approach [J]. Transportation Research Part C, 2010, 18(3): 311-331.

[5]Meshkat A, Zhi M, Vrancken J L M, et al. Coordinated ramp metering with priorities [J]. IET Intelligent Transport Systems, 2015, 9(6): 639-645.

[6]梁新榮, 劉艷艷, 滿(mǎn)國(guó)永, 等. 基于小腦模型關(guān)節(jié)控制器與PID復(fù)合的高速公路交通流密度控制 [J]. 控制理論與應(yīng)用, 2013, 30(10): 1281-1286.

[7]涂輝招, 王 穎, 謝欣睿. 基于定量分層模型的多匝道協(xié)調(diào)控制次序 [J]. 同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2017, 45(1): 39-45.

[8]張 偉, 肖日東, 鄧 晶. 基于遺傳算法的動(dòng)態(tài)模糊神經(jīng)網(wǎng)絡(luò)城市快速路入口匝道控制 [J]. 公路交通科技, 2017, 34(2): 129-134.

[9]王興宇, 吳文祥. 基于交通流預(yù)測(cè)的入口匝道模糊控制研究 [J]. 工業(yè)控制計(jì)算機(jī), 2017, 30(1): 1-2.

[10]Pasquale C, Papamichail I, Roncoli C, et al. Two-class freeway traffic regulation to reduce congestion and emissions via nonlinear optimal control [J]. Transportation Research Part C, 2015, 55(1): 85-99.

[11]Hou Z S, Xu J X, Yan J W. An iterative learning approach for density control of freeway traffic flow via ramp metering [J]. Transportation Research Part C, 2008, 16(1): 71-97.

[12]柴 干, 謝永利. 高速公路匝道多性能指標(biāo)優(yōu)化控制方法 [J]. 中國(guó)公路學(xué)報(bào), 2012, 25(6): 115-122.

[13]王 浩, 時(shí)柏營(yíng), 楊曉光. 自組織式交通信號(hào)協(xié)調(diào)控制模型 [J]. 計(jì)算機(jī)測(cè)量與控制, 2012, 20(11): 2979-2982.