交叉口信号配时

⑥各个相位的有效绿灯时间和显示绿灯时间,计算式所示：gej=Ge*max[yi,yi……]/YGe—总有效绿灯时间,就是C0减去L.计算各个相位的显示绿灯时间,公式如下所示：gj=gei-aj+lj3.2程序计算结果表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输入的数据和输出结果4延误分析与服务水平测定4.1延误估算方法延误是一个影响因素十分复杂的指标.设计交叉口时各车道延误用下式估算：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—分析时段的持续时长〔h〕,取0.25h；e—定时信号取0.5；各进口道的平均信控延误按该进口道中的各车道延误的加权平均数估算：d1=Σdiqi/Σqi式中：dA—进口道A的平均信控延误〔s/pcu〕；di—进口道A中第i车道的平均信控延误〔s/pcu〕；qi—进口道A中第i车道的小时交通量换算为其中高峰15min的交通流量〔辆/15min〕;整个交叉口的平均信控延误按交叉口各进口道延误的加权平均数估计：dA=ΣdAqA/ΣqA式中：d1—交叉口每车的平均信控延误〔s/pcu〕；qa—进口道A的高峰15min交通流率〔辆/15min〕；4.2服务水平表4延误服务水平表服务水平每车信控延误〔s〕服务水平每车信控延误〔s〕A<=10B11~20C21~35D36~55E56~80F>805结果分析5.1结果对比WEBSTER法的配时仿真实验结果和实测值的比较如下表：表相位配时对比分析相位灯色显示配时值实测值第一相位红灯13797绿灯6360黄灯33第二相位红灯129127绿灯7130黄灯33第三相位红灯155119绿灯4538黄灯33第四相位红灯188149绿灯129黄灯32表5.1.2延误对比分析现有延误设计延误延误值15．79.5服务等级BA由上表可以看出,实际采用的配时方案相关参数与我们通过调查得到的优化配时方案得到的配时方案相关参数都大致接近,说明实际情况与理论分析基本符合,说明配时方案可行.参考文献[1]吴兵《交通管理与控制》第四版人民交通2009.1[2]徐吉谦陈学武《交通工程总论》第三版人民交通2008.6[3]全永燊《城市交通控制》人民交通1989附录本次课程设计是我们小组几名同学共同完成的,起初是调查交叉口的交通流量,早高峰同学们要起的很早的,之后是两个人一小组测量车速和红绿灯时间并且在此基础上算出绿信比和信号周期.之后以班级为单位分上下午在机房进行信号配时,用C语言或其它形式将自己的数据编程序,算出所需要的数据.附录一：交通流量调查表交叉口:西南路/五一路方向:南进口时间:调查者:张博宇,姚博时间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交叉口:西南路/五一路方向:南进口时间:调查者:郭佟,左晓晨14,15直行15右转小中大和小中大和274543596374253539.5465563.5194637512258243742.5583674.5236542573367.5273439.5644478262437575268.5287854.5583368.5254337592470293441.5573469.5313341.5593369.5283440.5562365294647563468.5323444.5562365312340513463.5293339.5553569.5282439573469.5273337.5613371.5交叉口:西南路/五一路方向:南进口时间:调查者:张博宇17专左小中大和2103.520144028300345011.55118.590111420771110.55106.5111012.5110215101011.51120146119.5801108201181111.5交叉口:西南路/五一路方向:北进口时间:调查者:杨宁,柳继波时间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交叉口:西南路/五一路方向:北进口时间:调查者:马心怡16专左小中大和70211819.52217421971110.56211191112.550571212.56118718.5617.571110.551753213.571314.510314.591112.5交叉口:西南路/五一路方向:东进口时间:调查者:马心怡,柳继波时间4,5直行加左转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交叉口:西南路/五一路方向:东进口时间:调查者:杨宁,周贺7直行8直行小中大和小中大和18180101.521122.52014124201114.5201225.52115.5191122.55118.517174028202050291221680087720141111100114115.54117.5172211013151118.51001181810131313301517118.500121414200212121013交叉口:西南路/五一路方向:西进口时间:调查者:张博宇,姚博时间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交叉口:西南路/五一路方向:西进口时间:调查者:左晓晨3直行小中大和183430.5212432242435253233.5194229213331.5222331234437253437.5261333.5242231231330.5242333232434193329.5172326222229232332其他数据调查表初始积余车辆调查者:马亮,杨丽红,白俏,金曦车道编号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<"请输入停止线到冲突点的距离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]为各车道对应的饱和流量*/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]为各车道对应的饱和流量*/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]为各车道对应的饱和流量*/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]为各车道对应的饱和流量*/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中积余车辆的持续时间*/{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<>/*绿灯期车流到达率校正系数*/{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<"服务水平为A级\n">;elseif<dl>10&&dl<=20>printf<"服务水平为B级\n">;elseif<dl>20&&dl<=35>printf<"服务水平为C级\n">;elseif<dl>35&&dl<=55>printf<"服务水平为D级\n">;elseif<dl>55&&dl<=80>printf<"服务水平为E级\n">;elseif<dl>80>printf<"服务水平为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<"请输入一个周期内的相位数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<无实测数据时可取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<"请输入停止线到冲突点的距离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]为各车道对应的饱和流量*/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]为各车道对应的饱和流量*/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]为各车道对应的饱和流量*/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]为各车道对应的饱和流量*/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中积余车辆的持续时间*/{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<>/*绿灯期车流到达率校正系数*/{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>;