車道被占用對城市道路通行能力影響

故橫斷面實際通行能力，運用Gibbs故橫斷面實際通行能力，運用Gibbs抽樣仿真方法，通過excel軟件解決數(shù)據(jù)缺失問題，并用eviews進行ARMAyt=?2.0232+?t?0.7243?t?1? 虛數(shù)i2= j j j j 自然指數(shù)e 車道iP{N(t+?t)|N t時刻狀態(tài)為N(t)時，t+?t時刻狀態(tài)為N(t+?t) 額定座位≤19 額定座位≥19 載質(zhì)量≤2 載質(zhì)量≥2 額定座位≤19 額定座位≥19 載質(zhì)量≤2 載質(zhì)量≥2時間小車大車當(dāng)量轉(zhuǎn)換時間小車大車當(dāng)量轉(zhuǎn)換 GibbsX,Y等大寫字母表[X],[Y][X|Y],[Y|X]則表示條件分布。在二元的情況下，如果[X—Y]和GibbsX,Y等大寫字母表[X],[Y][X|Y],[Y|X]則表示條件分布。在二元的情況下，如果[X—Y]和[Y—X][X,Y]產(chǎn)生服從[X,Y]的隨機點列{(x,y)m}={(x1,y1),(x2,y2),...,(xn,yn),...},數(shù)列”{(x)n}={xxx(y)}={y1y2yn}分別服從[X,Y], g(x,y)=E[g(x, n→∞Step1任意選取X的一個可能取值點x1，根據(jù)[Y|X=x1]產(chǎn)生隨機數(shù)y1對(x1,y1)成為隨機點列{(x,y)n}x2為隨機點列(x,y)n中第二個點； 紅燈10 綠燈 7,8,8.5,9,9.5,10,10.5,11對應(yīng)用數(shù)字（1,2,3,4,5,6,7,8）分別表示這8個數(shù)。運用Gibbs抽樣填 交通 X 1 2 1 1 3 3 1 1 2 1 3 2 1 1 時間車輛當(dāng)量時間車輛當(dāng)量時間車輛當(dāng)量 yt=?2.0232+?t?0.7243yt=?2.0232+?t?0.7243?t?1?①首先將兩樣本數(shù)據(jù)(X1,X2,···,Xm)和(Y1,Y2,···,Yn)②分別對樣本(X1,X2,···,Xm)和(Y1,Y2,···,Yn)Wx/m和Wy/n,其③計算樣本(X1,X2,···,Xm)每個秩先于樣本(Y1,Y2,···,Yn)每個秩的個數(shù)U1,以及樣本(Y1,Y2,···,Yn)每個秩先于樣本(X1,X2,···,Xm)每個秩的個數(shù)U2U1=Wx?1/2m(m+1)U2=Wy1/2n(n1),且U1U2=mn.⑥計算曼-惠特尼U統(tǒng)計量。即U=W k(k+U? Z=12mn(m+n+車流量車流量車流量車流量車流量WilcoxonW-車流量車流量車流量車流量車流量WilcoxonW-WilcoxonW- WilcoxonW-我們以視頻1為比較基準(zhǔn)，分別計算60S為周期，綠燈和紅燈時候視頻2 生滅過程的定義如下：設(shè)一個系統(tǒng)具有可數(shù)個狀態(tài)（車輛數(shù)）0，1，2,···,或有限個狀態(tài)0,1,2···,m,令系統(tǒng)在時刻t的狀態(tài)記為N(t),在任一時刻生滅過程的定義如下：設(shè)一個系統(tǒng)具有可數(shù)個狀態(tài)（車輛數(shù)）0，1，2,···,或有限個狀態(tài)0,1,2···,m,令系統(tǒng)在時刻t的狀態(tài)記為N(t),在任一時刻t,若系統(tǒng)處于N(t)=j,則在時間區(qū)間(t;tt)內(nèi)，只能轉(zhuǎn)移到相鄰的狀態(tài)或者保持狀態(tài)不變，而不允許有兩次或者更多P{N(t+?t)=j+1|N(t)=j}=λj?t+(P{N(t+?t)=j?1|N(t)=j}=μj?t+P{N(t+?t)=j|N(t)=j}=1?(λj+μj)?t+N(t+?t)=j,N(t)=j;N(tj1;N(tj1.Pj(t+?t)=λj?1?tPj?1(t)+(1?(λj+μj)?t)Pj(t)+μj?tPj+1(t)+Pj(t+?t)?=(t)?(λ+μ)P(t)+j?1 j+1Pj′(t)=λj?1Pj?1(t)?(λj+μj)Pj(t)+當(dāng)j0Pn′(t)=λn?1Pn?1(t)?Pj′(t)=λPj?1(t)?(λ+μ)Pj(t)+Pj′(t)=λPj?1(t)?(λ+μ)Pj(t)+(j=1,2,···,n?>=AP ?λ? ?λ?μA=我們來求解形如P=rextx為矩陣A的特征值，r為x對應(yīng)的特征向量；則特征方程|xIA|0的解如下:|xI?Ax+x+λ+x+λ+μx+λ+x+=An+1?μAn?λAn+x+x+λ+x+λ+μx+λ+x+x+x+λ+x+λ+μx+λ+x+是一個n×n所以An+1=(xλμ)An所以|xIA|=xλμ)AnλμAn?1μAnλAnλμAn?1顯然x=0是其中一個特征值，由|An|=0(x+λ+μ+√??i)n+1?(x+λ+μ?An其中√a24bc,axλμ,b?μ,c設(shè)復(fù)角θ為a+√4bc?a2i由An0得：校園咖客收集整理(cosθ+isinθ)n+1=[cos(2Π?θ)+isin(2Π?(cosθ+isinθ)n+1=(cosθ?isincos((n+1)θ)+isin((n+1)θ)=cos((n+1)θ)?isin((n+所以sin((n+1)θ)=0,則θ ,k=1,2,···,由cos(kΠx= λμ—λ?μ,k=1,2,···,n+所以A的n+1個特征值x0x1···xn和對應(yīng)的n+1特征向量r0r1···n;其中x0=x=2√λμ—λ?(k>n+λ?=μ時，λ+μ>2√λμ,所以xk<=xt 初始狀況下有k輛車已經(jīng)在排隊，則初始條件Pk(0)=1,P=xt 初始狀況下有k輛車已經(jīng)在排隊，則初始條件Pk(0)=1,Pj(0)=0(j?=k時)0=c0r(0,0)+c1r(1,0)+···+1=c0r(0,k)+c1r(1,k)+···+0=c0r(0,n)+c1r(1,n)+···+xP(t) r(k, ∑xm(t) iP(t) r(k,i=0s(t)=l0m(t)?max{Pi}=λ=λ0+a1? (0.21l0m(t)<λ0? (0.21l0m(t)≥這時初始條件變?yōu)镻k(t0)=1,Pj(t0)=0(j?=k時)路口流入量為λ′2,小區(qū)路口流入量a1,小區(qū)路口流出量a2,道路通行能力考慮到紅燈和綠燈的持續(xù)時間都是30秒，即?t30秒，隊長的變化量由于紅綠燈不s(t+?t)?s(t)=λ′1?μ+a1? (t∈(60k,60k+λ′2?μ+a1? (t∈(60k+30,60k+其中事故持續(xù)時間t的單位是秒，k12··令w1λ′1μa1a2,w2λ′2μa1s(t)=s(0)+k(w1+w2)+w1(t? (t∈(60k,60k+s(0)+k(w1+w2)+w1+w2(t?60k? (t∈(60k+30,60k+ ∑xm( ∑xm(t) iP(t) r(k,i=0s(t)=s(0)+k(w1+w2)+w1(t? (t∈(60k,60k+s(0)+k(w1+w2)+w1+w2(t?60k? (t∈(60k+30,60k+ ------------------達率分別為0.89輛/0.5min和-0.67輛/0.5min ------------------達率分別為0.89輛/0.5min和-0.67輛/0.5min 占比例最重，中間道路往往最長，所以S(T)=0.44Ml0=140,其中l(wèi)0=7（米）M=46，考慮到系統(tǒng)的最大容量Mmax=60,化，當(dāng)概率P46(t1)=maxP46(t)于是用 x r(k,是到達時間。其中：λ=25,μ=x0=x=2√λμ—λ?(k>n+=··0=c0r(0,0)+c1r(1,0)+=··0=c0r(0,0)+c1r(1,0)+···+··1=c0r(0,k)+c1r(1,k)+···+··0=c0r(0,n)+c1r(1,n)+···+ =λ=1500 λ=1500=s(t+?t)?s(t)=λ′1?μ? (t∈(60k,60k+ (t∈(60k+30,60k+s(t)=s(0)+k(w1+w2)+w1(t? (t∈(60k,60k+s(0)+k(w1+w2)+w1+w2(t?60k? (t∈(60k+30,60k+s(t+?t)?s(t)= (t∈(60k,60k+ (t∈(60k+30,60k+s(t)= (t∈(60k,60k+ (t∈(60k+30,60k+ matlab matlab程序見附錄46輛車時的最大概率在6.6566.656478.1467.97分matlab（程序見附錄由下圖可知在386.2,,王習(xí)娟，M|M|1|C排隊論系統(tǒng)非穩(wěn)態(tài)解分析，連云港師范高等?？茖W(xué)校學(xué)報,第二期：46頁functionm=f(t)y=0;forift<(30+60*k)&&t>=(60*k)elseift>=(30+60*k)&&t<(60*k+60)elsey=y;functionuu=solve1(u,w,n,t)fork=1:i-s=a(k,f)*sifi<nforp=(a(o,f)-foro=1:i-p=(a(o,f)-endrr(n,:)=-endrr(n,:)=-Q=zeros(n+1,n+1);Q(1,1)=-w;fori=1:n-x=eig(Q);B(1)=1;fori=1:nendforforforfunctionuu=solve2(u,w,n,t)functionuu=solve2(u,w,n,t)fork=1:i-s=a(k,f)*sifi<nforp=(a(o,f)-foro=1:i-p=(a(o,f)-Q=zeros(n+1,n+1);Q(1,1)=-w;Q=zeros(n+1,n+