交叉口信號配時

⑥各個相位的有效綠燈時間和顯示綠燈時間,計算式所示：gej=Ge*max[yi,yi……]/YGe—總有效綠燈時間,就是C0減去L.計算各個相位的顯示綠燈時間,公式如下所示：gj=gei-aj+lj3.2程序計算結(jié)果表4各車道交通流量道路出口車道1車道2車道3車道4車道5車道6車道7車道8北進口2139013190791037033南進口300117126109881104631西進口17988817644555119東進口16089842791564517道路出口車道10車道11車道12車道13車道14車道15車道16車道17北進口7179120144693333433南進口6567112153682373430周期為203第一相位綠燈時間63s第二相位綠燈時間71s第三相位綠燈時間45s第四相位綠燈時間12s第一相位綠信比0.375第二相位綠信比0.1875第三相位綠信比0.74375第四相位綠信比0.05625圖4輸入的數(shù)據(jù)和輸出結(jié)果4延誤分析與服務(wù)水平測定4.1延誤估算方法延誤是一個影響因素十分復雜的指標.設(shè)計交叉口時各車道延誤用下式估算：d=d1+d2d1=0.5C<1-λ>2/<1-min[1,x]λ>d2=900T[<x-1>+√[<x-1>2+8ex/<CAP*T>]式中：C—周期時長〔s〕λ—所計算車道的綠信比；x—所計算車道的飽和度；CAP—所計算車道的通行能力；T—分析時段的持續(xù)時長〔h〕,取0.25h；e—定時信號取0.5；各進口道的平均信控延誤按該進口道中的各車道延誤的加權(quán)平均數(shù)估算：d1=Σdiqi/Σqi式中：dA—進口道A的平均信控延誤〔s/pcu〕；di—進口道A中第i車道的平均信控延誤〔s/pcu〕；qi—進口道A中第i車道的小時交通量換算為其中高峰15min的交通流量〔輛/15min〕;整個交叉口的平均信控延誤按交叉口各進口道延誤的加權(quán)平均數(shù)估計：dA=ΣdAqA/ΣqA式中：d1—交叉口每車的平均信控延誤〔s/pcu〕；qa—進口道A的高峰15min交通流率〔輛/15min〕；4.2服務(wù)水平表4延誤服務(wù)水平表服務(wù)水平每車信控延誤〔s〕服務(wù)水平每車信控延誤〔s〕A<=10B11~20C21~35D36~55E56~80F>805結(jié)果分析5.1結(jié)果對比WEBSTER法的配時仿真實驗結(jié)果和實測值的比較如下表：表相位配時對比分析相位燈色顯示配時值實測值第一相位紅燈13797綠燈6360黃燈33第二相位紅燈129127綠燈7130黃燈33第三相位紅燈155119綠燈4538黃燈33第四相位紅燈188149綠燈129黃燈32表5.1.2延誤對比分析現(xiàn)有延誤設(shè)計延誤延誤值15．79.5服務(wù)等級BA由上表可以看出,實際采用的配時方案相關(guān)參數(shù)與我們通過調(diào)查得到的優(yōu)化配時方案得到的配時方案相關(guān)參數(shù)都大致接近,說明實際情況與理論分析基本符合,說明配時方案可行.參考文獻[1]吳兵《交通管理與控制》第四版人民交通2009.1[2]徐吉謙陳學武《交通工程總論》第三版人民交通2008.6[3]全永燊《城市交通控制》人民交通1989附錄本次課程設(shè)計是我們小組幾名同學共同完成的,起初是調(diào)查交叉口的交通流量,早高峰同學們要起的很早的,之后是兩個人一小組測量車速和紅綠燈時間并且在此基礎(chǔ)上算出綠信比和信號周期.之后以班級為單位分上下午在機房進行信號配時,用C語言或其它形式將自己的數(shù)據(jù)編程序,算出所需要的數(shù)據(jù).附錄一：交通流量調(diào)查表交叉口:西南路/五一路方向:南進口時間:調(diào)查者:張博宇,姚博時間12直行13直行小中大和小中大和6:30-6:35191020.5273337.56:35-6:40210021293441.56:40-6:452141292826436:45-6:50183532.5343446.56:50-6:55263234.5313545.56:55-7:002740332745437:00-7:05223026.52944437:05-7:102740333245487:10-7:151940253446527:15-7:20171018.53025437:20-7:25101113.52743397:25-7:30202023243436.57:30-7:352020232645427:35-7:40233027.52824397:40-7:45221125.5293441.57:45-7:50202125303340.57:50-7:551900193223417:55-8:00201021.5283440.5交叉口:西南路/五一路方向:南進口時間:調(diào)查者:郭佟,左曉晨14,15直行15右轉(zhuǎn)小中大和小中大和274543596374253539.5465563.5194637512258243742.5583674.5236542573367.5273439.5644478262437575268.5287854.5583368.5254337592470293441.5573469.5313341.5593369.5283440.5562365294647563468.5323444.5562365312340513463.5293339.5553569.5282439573469.5273337.5613371.5交叉口:西南路/五一路方向:南進口時間:調(diào)查者:張博宇17專左小中大和2103.520144028300345011.55118.590111420771110.55106.5111012.5110215101011.51120146119.5801108201181111.5交叉口:西南路/五一路方向:北進口時間:調(diào)查者:楊寧,柳繼波時間9直行10,11直行小中大和小中大和6:30-6:3512113.517176:35-6:4015117133196:40-6:452020255356:45-6:5018119.5208366:50-6:551212145246:55-7:0016163227497:00-7:05182212626417:05-7:1024125.52928487:10-7:1511112.5222267:15-7:2026127.5323387:20-7:2519192126367:25-7:302121126247:30-7:351818123187:35-7:4015117185287:40-7:452020227367:45-7:501414166287:50-7:551414231432.57:55-8:0018119.57313交叉口:西南路/五一路方向:北進口時間:調(diào)查者:馬心怡16專左小中大和70211819.52217421971110.56211191112.550571212.56118718.5617.571110.551753213.571314.510314.591112.5交叉口:西南路/五一路方向:東進口時間:調(diào)查者:馬心怡,柳繼波時間4,5直行加左轉(zhuǎn)6直行小中大和小中大和6:30-6:351724281421196:35-6:4018222163020.56:40-6:4514422191224.56:45-6:50291336.51721226:50-6:5525251702216:55-7:00472502142317:00-7:055444681922267:05-7:1033331722247:10-7:15451148.51701197:15-7:20301439.51602207:20-7:25421655.5211022.57:25-7:30252291900197:30-7:35141219.5171120.57:35-7:40301337.5181121.57:40-7:452824391603227:45-7:50211226.51721227:50-7:5532133.5151118.57:55-8:0028232141117.5交叉口:西南路/五一路方向:東進口時間:調(diào)查者:楊寧,周賀7直行8直行小中大和小中大和18180101.521122.52014124201114.5201225.52115.5191122.55118.517174028202050291221680087720141111100114115.54117.5172211013151118.51001181810131313301517118.500121414200212121013交叉口:西南路/五一路方向:西進口時間:調(diào)查者:張博宇,姚博時間1專左2直行小中大和小中大和6:00-6:05455460.5223434.56:05-6:10487670.5233333.56:10-6:15535672.52424356:15-6:20544264233435.56:20-6:25549067.52423336:25-6:30598173213331.56:30-6:35549475.52024316:35-6:40459364.5213433.56:40-6:45235030.52623356:45-6:50494055253335.56:50-6:555121562124326:55-7:00537165.5203330.57:00-7:05576372213433.57:05-7:10593267.52442347:10-7:15567168.5203330.57:15-7:205341612125347:20-7:25504260223434.57:25-7:30485259.5203330.5交叉口:西南路/五一路方向:西進口時間:調(diào)查者:左曉晨3直行小中大和183430.5212432242435253233.5194229213331.5222331234437253437.5261333.5242231231330.5242333232434193329.5172326222229232332其他數(shù)據(jù)調(diào)查表初始積余車輛調(diào)查者:馬亮,楊麗紅,白俏,金曦車道編號1245678暫停車輛20132312222115離開車輛251829172726219101112131416179131020152532311420202521313837信號周期信號周期綠燈時間紅燈時間黃燈時間信號周期北進口631373203南進口711293203東進口451553203西進口121883203附錄二：#include"stdio.h"#include"stdlib.h"#include"math.h"floatS[18],l[5],Ge[5],maxA,maxB,maxC,maxD,D[5],g[5],max,A;floatdi,dk=0,dj=0,lkk;floatda[5],cap[18];floatP[18],Ql[18],Fa[18];floatTu[18],Qb[18],Ku[18],Ju[17],min;floatdA[5],dab[5],da1[17],T=0.25,x[18],Ds[5],Tu[18],Dd[5],Hu[18],e,d[18],Co,D[5],Q[18];floatQb[18],CAPP[18],v,o,Fs[18],y[18],Du[18],j[5],Ql[18],Qb[18],d[18],d1[18],d2[18],d3[18];inti;floatfun1<>/*綠燈間隔時間計算公式*/{floatz,u,t,h,I;printf<"請輸入停止線到?jīng)_突點的距離z<m>:">;scanf<"%f",&z>;printf<"請輸入車輛在進口道上的行駛車速u<m/s>:">;scanf<"%f",&u>;printf<"請輸入車輛制動時間t<s>:">;scanf<"%f",&t>;h=z/u;/*此處需插入一個界面*/I=h+t;printf<"計算得綠燈間隔時間I為%fs",I>;if<I<3>printf<"黃燈時間A=3s">;elseprintf<"<其中3s配以黃燈,其余時間配以紅燈>\n">;return<I>;}floatfun2<>/*最大流量比計算公式*/{inti,c,e;printf<"請輸入第一相位各車道的交通流量〔1-3〕:">;for<i=1;i<4;i++>scanf<"%f",&Q[i]>;for<c=1;c<4;c++>y[c]=Q[c]/S[c];/*S[c]為各車道對應(yīng)的飽和流量*/maxA=y[1];/*第一相位最大流量比*/for<e=1;e<=3;e++>if<y[e]<=maxA>maxA=y[e];printf<"第一相位最大流量比為%f\n",maxA>;printf<"請輸入第二相位各車道的交通流量〔4-8〕:">;for<i=4;i<9;i++>scanf<"%f",&Q[i]>;for<c=4;c<9;c++>y[c]=Q[c]/S[c];/*S[c]為各車道對應(yīng)的飽和流量*/maxB=y[4];/*第二相位最大流量比*/for<e=4;e<=9;e++>if<y[e]>=maxB>maxB=y[e];printf<"第二相位最大流量比為%f\n",maxB>;printf<"請輸入第三相位各車道的交通流量〔9-15〕:">;for<i=9;i<16;i++>scanf<"%f",&Q[i]>;for<c=9;c<16;c++>y[c]=Q[c]/S[c];/*S[c]為各車道對應(yīng)的飽和流量*/maxC=y[9];/*第二相位最大流量比*/for<e=9;e<=15;e++>if<y[e]>=maxC>maxC=y[e];printf<"第三相位最大流量比為%f\n:",maxC>;printf<"請輸入第四相位各車道的交通流量〔16-17〕:">;for<i=16;i<18;i++>scanf<"%f",&Q[i]>,y[c]=Q[c]/S[c];/*S[c]為各車道對應(yīng)的飽和流量*/maxD=y[16];/*第二相位最大流量比*/for<e=16;e<=17;e++>if<y[e]>=maxD>maxD=y[e];printf<"第四相位最大流量比為%f\n:",maxD>;max=maxA+maxB+maxC+maxD;printf<"此交叉口的流量比為%f:",max>;return<max>;}fun31<>{inti;for<i=1;i<5;i++>g[i]=Ge[i]-A+l[i],printf<"第%d相位的顯示綠燈時間為%fs\n",i,g[i]>;}fun32<>{inti;for<i=1;i<3;i++> Ge[i]=g[i]+A-l[i], printf<"第%d相位的有效綠燈時間為%fs\n",i,Ge[i]>;}intfun15<>{for<i=1;i<18;i++> {printf<"請輸入第%d車道的分析期初初始積余車輛<輛>:",i>;scanf<"%f",&Qb[i]>;printf<"請輸入第%d車道的分析期離開車輛<輛>:",i>;scanf<"%f",&Ql[i]>;}}fun3<>{for<i=1;i<18;i++>{if<i<4>CAPP[1]+=S[i]*D[1];elseif<i>=4&&i<9>CAPP[2]+=S[i]*D[2];elseif<i>=9&&i<16>CAPP[3]+=S[i]*D[3];elseCAPP[4]+=S[i]*D[4];}}fun4<>{for<i=1;i<5;i++>Ds[i]=<float>1.5*Co*<1-D[i]>;}fun5<>/*不飽和延誤的計算公式*/{floatj[5],min,Hu[18];for<i=1;i<5;i++>j[i]=0.5*Co*<1-D[i]>*<1-D[i]>;for<i=1;i<18;i++>{if<i<4>cap[i]=S[i]*D[1];elseif<i>=4&&i<9>cap[i]=S[i]*D[2];elseif<i>=9&&i<16>cap[i]=S[i]*D[3];elsecap[i]=S[i]*D[4];}for<i=1;i<18;i++>x[i]=<float>Q[i]/cap[i];for<i=1;i<18;i++>{min=x[i];if<min>1>min=1;lkk=<float>1-min;if<i<4>Hu[i]=<float>1-min*D[1];elseif<i>=4&&i<9>Hu[i]=<float>1-min*D[2];elseif<i>=9&&i<16>Hu[i]=<float>1-min*D[3];elseHu[i]=<float>1-min*D[4];}for<i=1;i<18;i++>{if<i<4>Du[i]=j[1]/Hu[i];elseif<i>=4&&i<9>Du[i]=j[2]/Hu[i];elseif<i>=9&&i<16>Du[i]=j[3]/Hu[i];elseDu[i]=j[4]/Hu[i];}}fun6<>/*在T中積余車輛的持續(xù)時間*/{min=T;for<i=1;i<18;i++>{{if<i<4> Ku[i]=cap[i]*lkk,Ju[i]=Qb[i]/Ku[i];elseif<i>=4&&i<9>Ku[i]=cap[i]*lkk,Ju[i]=Qb[i]/Ku[i];elseif<i>=9&&i<16>Ku[i]=cap[i]*lkk,Ju[i]=Qb[i]/Ku[i];elseKu[i]=cap[i]*lkk,Ju[i]=Qb[i]/Ku[i];}{if<min<Ju[i]>Tu[i]=min;elseTu[i]=Ju[i];}}}fun7<>/*綠燈期車流到達率校正系數(shù)*/{for<i=1;i<18;i++>P[i]=<Ql[i]-Qb[i]>/Ql[i];for<i=1;i<18;i++>{if<i<4>Fa[i]=<1-P[i]>/<<float>1-D[1]>;elseif<i>=4&&i<9>Fa[i]=<1-P[i]>/<1-D[2]>;elseif<i>=9&&i<16>Fa[i]=<1-P[i]>/<1-D[1]>;elseFa[i]=<1-P[i]>/<1-D[2]>;}}fun8<>/*均勻延誤計算公式*/{floatar[18],an[18],am[18],al[18],az[18];for<i=1;i<18;i++>if<i<4>ar[i]=Ds[1]*Tu[i]/T,an[i]=Fs[i]*Du[i]*<T-Tu[i]>,d1[i]=ar[i]+an[i]/T;elseif<"i>=4&&i<9">d1[i]=<Ds[2]*Tu[i]/T>+<Fs[i]*Du[i]*<T-Tu[i]>/T>;elseif<"i>=9&&i<16">d1[i]=<Ds[3]*Tu[i]/T>+<Fs[i]*Du[i]*<T-Tu[i]>/T>;elsed1[i]=<Ds[4]*Tu[i]/T>+<Fs[i]*Du[i]*<T-Tu[i]>/T>;}floatfun13<floatx>{if<x>=0.000000&&x<=0.500000>e=0.13;elseif<x>0.500000&&x<=0.600000>e=0.2;elseif<x>0.600000&&x<=0.700000>e=0.28;elseif<x>0.700000&&x<=0.800000>e=0.35;elseif<x>0.800000&&x<=1.000000>e=0.43;elseif<x>1>e=0.5;elseprintf<"輸入錯誤\n">;return<e>;}fun9<>/*隨機附加延誤計算公式*/{floatc[18],ad,ac,ae,ar;for<i=1;i<18;i++>{e=fun13<x[i]>;ac=pow<<x[i]-1>,2>;ae=<float>8*e*x[i];if<i<4>ar=cap[i]*T;elseif<i>=4&&i<9>ar=cap[i]*T;elseif<i>=9&&i<16>ar=cap[i]*T;elsear=cap[i]*T;c[i]=ac+ae/ar;ad=sqrt<c[i]>;d2[i]=<float>900*T*<<x[i]-1>+ad>;}}fun10<>/*初始排隊延誤計算公式*/{floatmin,ji[18],jk[18],jh[18];for<i=1;i<18;i++>{min=x[i];{if<min>1>min=1;}if<Tu[i]==T>{ji[i]=<float>3600*Qb[i];jk[i]=<float>1800*T*<1-min>;if<i<4>d3[i]=ji[i]/cap[i]-jk[i];elseif<i>=4&&i<9>d3[i]=ji[i]/cap[i]-jk[i];elseif<i>=9&&i<16>d3[i]=ji[i]/cap[i]-jk[i];elsed3[i]=ji[i]/cap[i]-jk[i];}else{if<i<4>d3[i]=<float>1800*Qb[i]*Tu[i]/<T*cap[i]>;elseif<i>=4&&i<9>d3[i]=<float>1800*Qb[i]*Tu[i]/<T*cap[i]>;elseif<i>=9&&i<16>d3[i]=<float>1800*Qb[i]*Tu[i]/<T*cap[i]>;elsed3[i]=<float>1800*Qb[i]*Tu[i]/<T*cap[i]>;}}}fun16<>{for<i=1;i<18;i++>d[i]=d1[i]+d2[i]+d3[i];}fun17<>{for<i=1;i<18;i++>d[i]=Du[i]+d2[i];}fun11<>/*各進口道平均信控延誤計算公式*/{intk;for<i=1;i<4;i++>dab[1]+=<d[i]*<Q[i]/4>>,da1[1]+=<Q[i]/4>;dA[1]=<float>dab[1]/da1[1];for<i=4;i<9;i++>dab[2]+=<d[i]*<Q[i]/4>>,da1[2]+=<Q[i]/4>;dA[2]=dab[2]/da1[2];for<i=9;i<12;i++>dab[3]+=<d[i]*<Q[i]/4>>,da1[3]+=<Q[i]/4>;dA[3]=dab[3]/da1[3];for<i=12;i<16;i++>dab[4]+=<d[i]*<Q[i]/4>>,da1[4]+=<Q[i]/4>;dA[4]=dab[4]/da1[4];i=16;dab[3]+=<d[i]*<Q[i]/4>>,da1[3]+=<Q[i]/4>;dA[3]=dab[3]/da1[3];i=17;dab[4]+=<d[i]*<Q[i]/4>>,da1[4]+=<Q[i]/4>;dA[4]=dab[4]/da1[4];}fun14<floatdl>{if<dl>=0&&dl<=10>printf<"服務(wù)水平為A級\n">;elseif<dl>10&&dl<=20>printf<"服務(wù)水平為B級\n">;elseif<dl>20&&dl<=35>printf<"服務(wù)水平為C級\n">;elseif<dl>35&&dl<=55>printf<"服務(wù)水平為D級\n">;elseif<dl>55&&dl<=80>printf<"服務(wù)水平為E級\n">;elseif<dl>80>printf<"服務(wù)水平為F級\n">;elseprintf<"error">;}fun12<>/*整個交叉口*/{for<i=1;i<5;i++>{dk+=da1[i];dj+=<float><dA[i]*dab[i]>;}di=<float>dj/dk;printf<"此交叉口的信控延誤為%f\n",di>;fun14<di>;}voidmain<>{intL,i,I,k,n;floata,b,Y,ge,h[17];printf<"請輸入黃燈時間A<一般取值3s>:">;scanf<"%f",&A>;/*此處需插入一個界面*/printf<"請輸入一個周期內(nèi)的相位數(shù)k:">;scanf<"%d",&k>;for<i=1;i<4;i++>printf<"請輸入第一相位第%d車道的車頭時距:",i>,scanf<"%f",h+i>,S[i]=3600/h[i];for<i=4;i<9;i++>printf<"請輸入第二相位第%d車道的車頭時距:",i>,scanf<"%f",h+i>,S[i]=3600/h[i];for<i=9;i<16;i++>printf<"請輸入第三相位第%d車道的車頭時距:",i>,scanf<"%f",h+i>,S[i]=3600/h[i];for<i=16;i<18;i++>printf<"請輸入第四相位第%d車道的車頭時距:",i>,scanf<"%f",h+i>,S[i]=3600/h[i];I=fun1<>;Y=fun2<>;for<i=1;i<5;i++> printf<"\n請輸入第%d相位的啟動損失時間l<無實測數(shù)據(jù)時可取3s>:",i>,scanf<"%f",&l[i]>;for<i=1,L=0;i<=k;i++>L+=<l[i]>+I-A;/*I為綠燈間隔時間*/a=1.5*L+5;b=1-Y;Co=a/b;printf<"計算得周期時長Co為%ds\n",Co>;ge=Co-L;printf<"總有效綠燈時間為%d\n",ge>;Y=maxA+maxB+maxC+maxD;a=ge*<maxA>;b=ge*maxB;v=ge*maxC;o=ge*<maxD>;Ge[1]=a/Y;Ge[2]=b/Y;Ge[3]=v/Y;Ge[4]=o/Y;for<n=1;n<5;n++>D[n]=Ge[n]/Co,printf<"第%d相位的綠信比為%f\n",n,D[n]>;for<i=1;i<5;i++>printf<"第%d相位有效綠燈時間為%f\n",i,Ge[i]>;fun31<>;fun15<>;fun3<>;fun4<>;fun5<>;fun6<>;fun7<>;fun8<>;fun9<>;fun10<>;fun17<>;fun11<>;fun12<>; }附錄二：#include"stdio.h"#include"stdlib.h"#include"math.h"floatS[18],l[5],Ge[5],maxA,maxB,maxC,maxD,D[5],g[5],max,A;floatdi,dk,dj,lkk;floatda[5],cap[18];floatP[18],Ql[18],Fa[18];floatTu[18],Qb[18],Ku[18],Ju[17],min;floatdA[5],dab[5],da1[17],T=0.25,x[18],Ds[5],Tu[18],Dd[5],Hu[18],e,d[18],Co,D[5],Q[18];floatQb[18],CAPP[18],v,o,Fs[18],y[18],Du[18],j[5],Ql[18],Qb[18],d[18],d1[18],d2[18],d3[18];inti;floatfun1<>/*綠燈間隔時間計算公式*/{floatz,u,t,h,I;printf<"請輸入停止線到?jīng)_突點的距離z<m>:">;scanf<"%f",&z>;printf<"請輸入車輛在進口道上的行駛車速u<m/s>:">;scanf<"%f",&u>;printf<"請輸入車輛制動時間t<s>:">;scanf<"%f",&t>;h=z/u;/*此處需插入一個界面*/I=h+t;printf<"計算得綠燈間隔時間I為%fs",I>;if<I<3>printf<"黃燈時間A=3s">;elseprintf<"<其中3s配以黃燈,其余時間配以紅燈>\n">;return<I>;}floatfun2<>/*最大流量比計算公式*/{inti,c,e;printf<"請輸入第一相位各車道的交通流量〔1-3〕:">;for<i=1;i<4;i++>scanf<"%f",&Q[i]>;for<c=1;c<4;c++>y[c]=Q[c]/S[c];/*S[c]為各車道對應(yīng)的飽和流量*/maxA=y[1];/*第一相位最大流量比*/for<e=1;e<=3;e++>if<y[e]<=maxA>maxA=y[e];printf<"第一相位最大流量比為%f\n",maxA>;printf<"請輸入第二相位各車道的交通流量〔4-8〕:">;for<i=4;i<9;i++>scanf<"%f",&Q[i]>;for<c=4;c<9;c++>y[c]=Q[c]/S[c];/*S[c]為各車道對應(yīng)的飽和流量*/maxB=y[4];/*第二相位最大流量比*/for<e=4;e<=9;e++>if<y[e]>=maxB>maxB=y[e];printf<"第二相位最大流量比為%f\n",maxB>;printf<"請輸入第三相位各車道的交通流量〔9-15〕:">;for<i=9;i<16;i++>scanf<"%f",&Q[i]>;for<c=9;c<16;c++>y[c]=Q[c]/S[c];/*S[c]為各車道對應(yīng)的飽和流量*/maxC=y[9];/*第二相位最大流量比*/for<e=9;e<=15;e++>if<y[e]>=maxC>maxC=y[e];printf<"第三相位最大流量比為%f\n:",maxC>;printf<"請輸入第四相位各車道的交通流量〔16-17〕:">;for<i=16;i<18;i++>scanf<"%f",&Q[i]>,y[c]=Q[c]/S[c];/*S[c]為各車道對應(yīng)的飽和流量*/maxD=y[16];/*第二相位最大流量比*/for<e=16;e<=17;e++>if<y[e]>=maxD>maxD=y[e];printf<"第四相位最大流量比為%f\n:",maxD>;max=maxA+maxB+maxC+maxD;printf<"此交叉口的流量比為%f:",max>;return<max>;}fun31<>{inti;for<i=1;i<5;i++>g[i]=Ge[i]-A+l[i],printf<"第%d相位的顯示綠燈時間為%fs\n",i,g[i]>;}fun32<>{inti;for<i=1;i<5;i++> Ge[i]=g[i]+A-l[i], printf<"第%d相位的有效綠燈時間為%fs\n",i,Ge[i]>;}intfun15<>{for<i=1;i<18;i++> {printf<"請輸入第%d車道的分析期初初始積余車輛<輛>:",i>;scanf<"%f",&Qb[i]>;printf<"請輸入第%d車道的分析期離開車輛<輛>:",i>;scanf<"%f",&Ql[i]>;}}fun3<>{for<i=1;i<18;i++>{if<i<4>CAPP[1]+=S[i]*D[1];elseif<i>=4&&i<9>CAPP[2]+=S[i]*D[2];elseif<i>=9&&i<16>CAPP[3]+=S[i]*D[3];elseCAPP[4]+=S[i]*D[4];}}fun4<>{for<i=1;i<5;i++>Ds[i]=<float>1.5*Co*<1-D[i]>;}fun5<>/*不飽和延誤的計算公式*/{floatj[5],min,Hu[18];for<i=1;i<5;i++>j[i]=0.5*Co*<1-D[i]>*<1-D[i]>;for<i=1;i<18;i++>{if<i<4>cap[i]=S[i]*D[1];elseif<i>=4&&i<9>cap[i]=S[i]*D[2];elseif<i>=9&&i<16>cap[i]=S[i]*D[3];elsecap[i]=S[i]*D[4];}for<i=1;i<18;i++>x[i]=<float>Q[i]/cap[i];for<i=1;i<18;i++>{min=x[i];if<min>1>min=1;lkk=<float>1-min;if<i<4>Hu[i]=<float>1-min*D[1];elseif<i>=4&&i<9>Hu[i]=<float>1-min*D[2];elseif<i>=9&&i<16>Hu[i]=<float>1-min*D[3];elseHu[i]=<float>1-min*D[4];}for<i=1;i<18;i++>{if<i<4>Du[i]=j[1]/Hu[i];elseif<i>=4&&i<9>Du[i]=j[2]/Hu[i];elseif<i>=9&&i<16>Du[i]=j[3]/Hu[i];elseDu[i]=j[4]/Hu[i];}}fun6<>/*在T中積余車輛的持續(xù)時間*/{min=T;for<i=1;i<18;i++>{{if<i<4> Ku[i]=cap[i]*lkk,Ju[i]=Qb[i]/Ku[i];elseif<i>=4&&i<9>Ku[i]=cap[i]*lkk,Ju[i]=Qb[i]/Ku[i];elseif<i>=9&&i<16>Ku[i]=cap[i]*lkk,Ju[i]=Qb[i]/Ku[i];elseKu[i]=cap[i]*lkk,Ju[i]=Qb[i]/Ku[i];}{if<min<Ju[i]>Tu[i]=min;elseTu[i]=Ju[i];}}}fun7<>/*綠燈期車流到達率校正系數(shù)*/{for<i=1;i<18;i++>P[i]=<Ql[i]-Qb[i]>/Ql[i];for<i=1;i<18;i++>{if<i<4>Fa[i]=<1-P[i]>/<<float>1-D[1]>;elseif<i>=4&&i<9>Fa[i]=<1-P[i]>/<1-D[2]>;elseif<i>=9&&i<16>Fa[i]=<1-P[i]>/<1-D[1]>;elseFa[i]=<1-P[i]>/<1-D[2]>;}}fun8<>/*均勻延誤計算公式*/{floatar[18],an[18],am[18],al[18],az[18];for<i=1;i<18;i++>if<i<4>ar[i]=Ds[1]*Tu[i]/T,an[i]=Fs[i]*Du[i]*<T-Tu[i]>,d1[i]=ar[i]+an[i]/T;elseif<"i>=4&&i<9">d1[i]=<Ds[2]*Tu[i]/T>+<Fs[i]*Du[i]*<T-Tu[i]>/T>;elseif<"i>=9&&i<16">d1[i]=<Ds[3]*Tu[i]/T>+<Fs[i]*Du[i]*<T-Tu[i]>/T>;elsed1[i]=<Ds[4]*Tu[i]/T>+<Fs[i]*Du[i]*<T-Tu[i]>/T>;}floatfun13<floatx>{if<x>=0.000000&&x<=0.500000>e=0.13;elseif<x>0.500000&&x<=0.600000>e=0.2;elseif<x>0.600000&&x<=0.700000>e=0.28;elseif<x>0.700000&&x<=0.800000>e=0.35;elseif<x>0.800000&&x<=1.000000>e=0.43;elseif<x>1>e=0.5;elseprintf<"輸入錯誤\n">;return<e>;}fun9<>/*隨機附加延誤計算公式*/{floatc[18],ad,ac,ae,ar;for<i=1;i<18;i++>{e=fun13<x[i]>;ac=pow<<x[i]-1>,2>;ae=<float>8*e*x[i];if<i<4>ar=cap[i]*T;elseif<i>=4&&i<9>ar=cap[i]*T;elseif<i>=9&&i<16>ar=cap[i]*T;elsear=cap[i]*T;c[i]=ac+ae/ar;ad=sqrt<c[i]>;d2[i]=<float>900*T*<<x[i]-1>+ad>;}}fun10<>/*初始排隊延誤計算公式*/{floatmin,ji[18],jk[18],jh[18];for<i=1;i<18;i++>{min=x[i];{if<min>1>min=1;}if<Tu[i]==T>{ji[i]=<float>3600*Qb[i];jk[i]=<float>1800*T*<1-min>;if<i<4>d3[i]=ji[i]/cap[i]-jk[i];elseif<i>=4&&i<9>d3[i]=ji[i]/cap[i]-jk[i];elseif<i>=9&&i<16>d3[i]=ji[i]/cap[i]-jk[i];elsed3[i]=ji[i]/cap[i]-jk[i];}else{if<i<4>d3[i]=<float>1800*Qb[i]*Tu[i]/<T*cap[i]>;elseif<i>=4&&i<9>d3[i]=<float>1800*Qb[i]*Tu[i]/<T*cap[i]>;elseif<i>=9&&i<16>d3[i]=<float>1800*Qb[i]*Tu[i]/<T*cap[i]>;elsed3[i]=<float>1800*Qb[i]*Tu[i]/<T*cap[i]>;}}}fun16<>{for<i=1;i<18;i++>d[i]=d1[i]+d2[i]+d3[i];}fun17<>{for<i=1;i<18;i++>d[i]=Du[i]+d2[i];}fun11<>/*各進口道平均信控延誤計算公式*/{intk;for<i=1;i<4;i++>dab[1]+=<d[i]*<Q[i]/4>>,da1[1]+=<Q[i]/4>;