交叉口信号配时

摘要道路交叉口是指两条或两条以上道路的相交处.车辆、行人汇集、转向和疏散的必经之地,为交通的咽喉.因此,正确设计道路交叉口,合理组织、管理交叉口交通,是提高道路通行能力和保障交通安全的重要方面.此次交叉口信号灯控制配时的调查地点是西南路和五一路交叉口.该交叉口地处市区西南部,属于平面十字型交叉口.西南路方向路段为双向五车道；五一路方向由东向西黄线以北是五车道,黄线以南是五车道,五一路由西向东黄线以北是两车道,黄线以南是三车道.周围分布饭店、居民住宅区、净水厂等,是一个非常重要的交叉口,并且西南路是主干道.本组通过实际观测的方法测得了道路交叉口的交通流量等信息.西南路车流量比五一路车流量大很多,在五一路方向均有左转车流,西南路只在南进口存在左转车流,另外在五一路西路口和西南路南路口均有直行加右转相位.且西南路南进口的左转仅限公交且车流量极少.到目前为止,定时信号的配时方法在国际上主要有英国的WEBSTER法,澳大利亚ARRB法与美国HCM法等.我们在《交通管理与控制》课本中已经学会了webster法和HCM法,我国有停车线法和冲突点法等方法.随着研究不断深入,定时信号的配时方法也在进一步的改进.本设计采用的方法以英国的WEBSTER法为主.本次设计本小组分工合作,共采集了车道宽、交通流量、车头时距、信号灯信号显示与周期等数据.并且对数据作出了运算整理.摒弃了有问题的数据,保证使用严谨的数据进行运算.关键字道路交叉口,信号配时,WEBSTER法,相位,课程设计.目录第一章现状交通调查1.1西南路与五一路交叉口现状概况 (1)1.2交通流量调整 (2)1.3交叉口几何尺寸调查 (2)第二章信号相位分析2.1实地观测 (3)2.2理论依据 (5)2.3具体算法步骤 (5)2.4必要性分析结果 (6)第三章制定配时方案3.1信号配时方案原理 (7)3.2程序计算结果 (8)第四章延误分析与服务水平测定4.1延误估算方法 (10)4.2服务水平 (10)第五章结果分析5.1结果对比 (12)第一章现状交通调查1.1 西南路/五一路交叉口现状概况道路交叉口是指两条或两条以上道路的相交处.车辆、行人汇集、转向和疏散的必经之地,为交通的咽喉.因此,正确设计道路交叉口,合理组织、管理交叉口交通,是提高道路通行能力和保障交通安全的重要方面.此次交叉口信号灯控制配时的调查地点是西南路和五一路交叉口.该交叉口地处市区西南部,属于平面十字型交叉口.西南路方向路段为双向五车道；五一路由东向西方向为双向五车道,包括自行车道,五一路由西向东黄线以北是两辆车道、黄线以南是三车道,周围分布饭店、居民住宅区、净水厂和大连沙河口区的部分政府部门等,是一个非常重要的交叉口,并且西南路是主干道.根据实地观察测量和分析讨论,本组对整个交叉口形状、车道划分与交通流运行轨迹进行了绘制,如下1.2交通流量调查本组通过实际观测的方法测得了道路交叉口的交通流量等信息.具体观测时间为.具体测量方法为,对各车道分配人员进行定时计数以获得交通流量,对直行和左转的车道在红灯结束时按车计时获得车头时距,对各相位红绿灯分别计时等.交通流量通过在交叉口的高峰时间观测2小时获得,整理后的高峰小时流量见下表.1.3 交叉口几何尺寸调查交叉口几何尺寸通过实际测量,整理后的信息见下表2信号相位分析2.1实地观测由本次设计实地观测的西南路/五一路交叉口情况可知：西南路车流量比五一路车流量大很多,且只有西南路南进口和五一西进口存在左转车流.且西南路南进口的左转仅限公交且车流量极少.而我们想知道是否有必要设置专有的左转相位,就必须求出各相位的左转通行能力与实测值进行比较.2.2理论依据可插车间隙理论虽然是用于计算无信号机的平面交叉口通行能力,而无信号交叉口遵循的是优先规则控制,即其通行顺序依次为：主要道路直行流向、主要道路左转流向、次要道路直行流向、次要道路左转流向.两相位交叉口车流运行特征跟交叉口优先规则控制有一些相似之处,因此可以根据可穿插间隙理论来分析其左转车道的通行能力.2.3具体算法步骤根据可插入间隙理论,假设主要道路车流按泊松分布到达,可建立计算左转通行能力Q ’max 的公式如下：)1/(maxqh q e Qe Q ---='τ式中：Q ’max ---次要道路可通过的最大交通量Q---主要道路交通量q--- Q/3600t c ---次要道路可以穿过主要道路车流的临界间隙〔s 〕 h---次要道路车辆连续通过时的车头时距〔s 〕2.4必要性分析结果由计算可以得出：可见每个相位都能满足其冲突左转交通量.故没必要增设一个专用左转相位.所以我们决定对西南路/五一路交叉口采用四相位配时方案是合理的.3制定配时方案到目前为止,定时信号的配时方法在国际上主要有英国的WEBSTER法,澳大利亚ARRB法与美国HCM法等.我国有停车线法和冲突点法等方法.随着研究不断深入,定时信号的配时方法也在进一步的改进.本设计采用的方法以英国的WEBSTER法为主.其信号配时设计流程图和信号相位基本方案如下：3.1信号配时方案原理①首先计算每车道的饱和流量S,使用下式进行计算：S＝3600/hh——饱和车头时距②然后计算流量比.公式如下.y i=q i/s i③计算流量比的总和,公式如下式：y=Σmax[y j,y j……]= Σmax[〔q d/s d〕j, 〔q d/s d〕j……]④启动损失时间L=Σ<l+I-A>⑤信号周期时长的计算,公式如下所示：C=<1.5l+5>/<1-y>C0—周期时长,Y—流量比总和,L—信号总损失时间,⑥各个相位的有效绿灯时间和显示绿灯时间,计算式所示：ge j=Ge*max[y i,y i……]/YGe—总有效绿灯时间,就是C0减去L.计算各个相位的显示绿灯时间,公式如下所示：g j= ge i - a j + l j3.2程序计算结果表4 各车道交通流量周期为203第一相位绿灯时间63s第二相位绿灯时间71s第三相位绿灯时间45s第四相位绿灯时间12s第一相位绿信比0.375第二相位绿信比0.1875第三相位绿信比0.74375第四相位绿信比0.05625图4 输入的数据和输出结果4延误分析与服务水平测定4.1延误估算方法延误是一个影响因素十分复杂的指标.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；2.各进口道的平均信控延误按该进口道中的各车道延误的加权平均数估算：d1=Σd i q i/Σq i式中：dA—进口道A的平均信控延误〔s/pcu〕；di—进口道A中第i车道的平均信控延误〔s/pcu〕；qi—进口道A中第i车道的小时交通量换算为其中高峰15min 的交通流量〔辆/15min〕;3.整个交叉口的平均信控延误按交叉口各进口道延误的加权平均数估计：d A=Σd A q A/Σq A式中：d1—交叉口每车的平均信控延误〔s/pcu〕；qa—进口道A的高峰15min交通流率〔辆/15min〕；4.2服务水平表4 延误服务水平表5结果分析5.1 结果对比WEBSTER法的配时仿真实验结果和实测值的比较如下表：由上表可以看出,实际采用的配时方案相关参数与我们通过调查得到的优化配时方案得到的配时方案相关参数都大致接近,说明实际情况与理论分析基本符合,说明配时方案可行.参考文献[1]吴兵《交通管理与控制》第四版人民交通 2009.1[2]徐吉谦陈学武《交通工程总论》第三版人民交通2008.6[3]全永燊《城市交通控制》人民交通 1989附录本次课程设计是我们小组几名同学共同完成的,起初是调查交叉口的交通流量,早高峰同学们要起的很早的,之后是两个人一小组测量车速和红绿灯时间并且在此基础上算出绿信比和信号周期.之后以班级为单位分上下午在机房进行信号配时,用C语言或其它形式将自己的数据编程序,算出所需要的数据.附录一：交通流量调查表其他数据调查表附录二：#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;float di,dk=0,dj=0,lkk;float da[5],cap[18];float P[18],Ql[18],Fa[18];float Tu[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];int i;float fun1 <> /*绿灯间隔时间计算公式*/{ float z,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为%f s",I>;if<I<3>printf<"黄灯时间A=3s">;elseprintf<"<其中3s配以黄灯,其余时间配以红灯>\n">;return <I>;}float fun2 <> /*最大流量比计算公式*/{ int i,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<>{ int i;for<i=1;i<5;i++>g[i]=Ge[i]-A+l[i],printf<"第%d相位的显示绿灯时间为%f s\n",i,g[i]>;}fun32<>{int i;for<i=1;i<3;i++>Ge[i]=g[i]+A-l[i],printf<"第%d相位的有效绿灯时间为%f s\n",i,Ge[i]>;}int fun15<>{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];else if<i>=4&&i<9>CAPP[2]+=S[i]*D[2];else if<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<> /*不饱和延误的计算公式*/{ float j[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];else if<i>=4&&i<9>cap[i]=S[i]*D[2];else if<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];else if<i>=4&&i<9>Hu[i]=<float>1-min*D[2];else if<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];else if<i>=4&&i<9>Du[i]=j[2]/Hu[i];else if<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];else if<i>=4&&i<9>Ku[i]=cap[i]*lkk,Ju[i]=Qb[i]/Ku[i];else if<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]>;else if<i>=4&&i<9>Fa[i]=<1-P[i]>/<1-D[2]>;else if<i>=9&&i<16>Fa[i]=<1-P[i]>/<1-D[1]>;elseFa[i]=<1-P[i]>/<1-D[2]>;}}fun8<> /*均匀延误计算公式*/ { float ar[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;else if<"i>=4&&i<9">d1[i]=<Ds[2]*Tu[i]/T>+<Fs[i]*Du[i]*<T-Tu[i]>/T> ;else if<"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> ;}float fun13< float x>{if <x>=0.000000&&x<=0.500000>e=0.13;else if<x>0.500000&&x<=0.600000>e=0.2;else if<x>0.600000&&x<=0.700000>e=0.28;else if<x>0.700000&&x<=0.800000>e=0.35;else if<x>0.800000&&x<=1.000000>e=0.43;else if<x>1>e=0.5;elseprintf<"输入错误\n">;return <e>;}fun9<> /*随机附加延误计算公式*/{ float c[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;else if<i>=4&&i<9>ar=cap[i]*T;else if<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<> /*初始排队延误计算公式*/{ float min,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];else if<i>=4&&i<9>d3[i]=ji[i]/cap[i]-jk[i];else if<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]>;else if<i>=4&&i<9>d3[i]=<float>1800*Qb[i]*Tu[i]/<T*cap[i]>;else if<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<> /*各进口道平均信控延误计算公式*/{int k;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<float dl>{if <dl>=0&&dl<=10>printf<"服务水平为A级\n">; else if<dl>10&&dl<=20>printf<"服务水平为B级\n">; else if<dl>20&&dl<=35>printf<"服务水平为C级\n">;else if<dl>35&&dl<=55>printf<"服务水平为D级\n">;else if<dl>55&&dl<=80>printf<"服务水平为E级\n">;else if<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>;}void main<>{int L,i,I,k,n;float a,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为%d s\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;float di,dk,dj,lkk;float da[5],cap[18];float P[18],Ql[18],Fa[18];float Tu[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];int i;float fun1 <> /*绿灯间隔时间计算公式*/{ float z,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为%f s",I>;if<I<3>printf<"黄灯时间A=3s">;elseprintf<"<其中3s配以黄灯,其余时间配以红灯>\n">;return <I>;}float fun2 <> /*最大流量比计算公式*/{ int i,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<>{ int i;for<i=1;i<5;i++>g[i]=Ge[i]-A+l[i],printf<"第%d相位的显示绿灯时间为%f s\n",i,g[i]>;}fun32<>{int i;for<i=1;i<5;i++>Ge[i]=g[i]+A-l[i],printf<"第%d相位的有效绿灯时间为%f s\n",i,Ge[i]>;}int fun15<>{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];else if<i>=4&&i<9>CAPP[2]+=S[i]*D[2];else if<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<> /*不饱和延误的计算公式*/{ float j[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];else if<i>=4&&i<9>cap[i]=S[i]*D[2];else if<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];else if<i>=4&&i<9>Hu[i]=<float>1-min*D[2];else if<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];else if<i>=4&&i<9>Du[i]=j[2]/Hu[i];else if<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];else if<i>=4&&i<9>Ku[i]=cap[i]*lkk,Ju[i]=Qb[i]/Ku[i];else if<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]>;else if<i>=4&&i<9>Fa[i]=<1-P[i]>/<1-D[2]>;else if<i>=9&&i<16>Fa[i]=<1-P[i]>/<1-D[1]>;elseFa[i]=<1-P[i]>/<1-D[2]>;}}fun8<> /*均匀延误计算公式*/ { float ar[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;else if<"i>=4&&i<9">d1[i]=<Ds[2]*Tu[i]/T>+<Fs[i]*Du[i]*<T-Tu[i]>/T> ;else if<"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> ;}float fun13< float x>{if <x>=0.000000&&x<=0.500000>e=0.13;else if<x>0.500000&&x<=0.600000>e=0.2;else if<x>0.600000&&x<=0.700000>e=0.28;else if<x>0.700000&&x<=0.800000>e=0.35;else if<x>0.800000&&x<=1.000000>e=0.43;else if<x>1>e=0.5;elseprintf<"输入错误\n">;return <e>;}fun9<> /*随机附加延误计算公式*/{ float c[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;else if<i>=4&&i<9>ar=cap[i]*T;else if<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<> /*初始排队延误计算公式*/{ float min,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];else if<i>=4&&i<9>d3[i]=ji[i]/cap[i]-jk[i];else if<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]>;else if<i>=4&&i<9>d3[i]=<float>1800*Qb[i]*Tu[i]/<T*cap[i]>;else if<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<> /*各进口道平均信控延误计算公式*/{int k;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[3]=dab[3]/da1[3];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[3]/da1[3];}fun14<float dl>{if <dl>=0&&dl<=10>printf<"服务水平为A级\n">; else if<dl>10&&dl<=20>printf<"服务水平为B级\n">;else if<dl>20&&dl<=35>printf<"服务水平为C级\n">;else if<dl>35&&dl<=55>printf<"服务水平为D级\n">;else if<dl>55&&dl<=80>printf<"服务水平为E级\n">;else if<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>;}void main<>{int L,i,I,k,ge,n;float a,b,Y,h[18];printf<"请输入黄灯时间A<一般取值3s>:">;scanf<"%f",&A>; /*此处需插入一个界面*/printf<"请输入一个周期内的绿灯间隔数k:">;scanf<"%d",&k>;for<i=1;i<9;i++>{if<i<=3>{printf<"请输入第一相位第%d车道的车头时距:",<0+i>>;scanf<"%f",h+i>;S[i]=3600/h[i];}else{printf<"请输入第二相位第%d车道的车头时距:",<i>>;scanf<"%f",&h[i]>;S[i]=3600/h[i];}}for<i=1;i<10;i++>{if<i<=7>printf<"请输入第三相位第%d车道的车头时距:",8+i>,scanf<"%f",&h[8+i]>,S[8+i]=3600/h[8+i];elseprintf<"请输入第四相位第%d车道的车头时距:",8+i>,scanf<"%f",&h[8+i]>,S[8+i]=3600/h[8+i];}printf<"请输入绿灯间隔时间I:">;scanf<"%f",&I>;Y=fun2<>;printf<"请输入现有交叉口的信号周期Co:">;scanf<"%f",&Co>;for<i=1;i<5;i++>printf<"请输入第%d相位的启动损失时间l<无实测数据时可取3s>:",i>,scanf<"%f",&l[i]>;for<i=1;i<=4;i++>L+=<l[i]>+I-A,printf<"请输入第%d相位的实际显示绿灯时间:",i>,scanf<"%f",&g[i]>;fun32<>;for<n=1;n<5;n++>D[n]=Ge[n]/Co,printf<"第%d相位的绿信比为%.2f\n",n,D[n]>;for<i=1;i<3;i++>printf<"第%d相位有效绿灯时间为%.1f\n",i,Ge[i]>;fun15<>;fun3<>;fun4<>;fun5<>;fun6<>;fun7<>;fun8<>;fun9<>;fun10<>;fun16<>;fun11<>;fun12<>;}。