車道被占用對城市道路通行能力及影響(全國一等獎)-資料下載頁
2025-09-25 19:14本頁面
【正文】 [a,b]=find(randperm(L)=N)。cells_posi=b。%變化規(guī)則for i=1:T%定義車頭間距 if cells_posi(N)cells_posi(1) headways(N)=Lcells_posi(N)+cells_posi(1)1。 else headways(N)=cells_posi(1)cells_posi(N)1。 end for j=N1:1:1 if cells_posi(j+1)cells_posi(j) headways(j)=cells_posi(j+1)cells_posi(j)1。 else headways(j)=L+cells_posi(j+1)cells_posi(j)1。 end end %速度變化 cells_veloNS1=min([v_max1,cells_velo(1),max(0,headways(1)1)])。 %NS規(guī)則下第一輛車的速度估計值 cells_velo(N)=min([v_max,cells_velo(N)+1,headways(N)+cells_veloNS1])。 %the speed of the Nth vehicle under VV rules for j=N1:1:1 cells_veloNS=min([v_max1,cells_velo(j+1),max(0,headways(j+1)1)])。 %NS規(guī)則下前一輛車的速度估計值 cells_velo(j)=min([v_max,cells_velo(j)+1,headways(j)+cells_veloNS])。 %the speed under VV rules end %以概率P隨機慢化 if rand()P。 cells_velo=max(cells_velo1,0)。 end %位置更新 for j=N:1:1 cells_posi(j)=cells_posi(j)+cells_velo(j)。 if cells_posi(j)=L cells_posi(j)=cells_posi(j)L。 end end %采集數(shù)據(jù)作圖 if i10*L+10000 %采用每組的后30個變量取平均 speed(m)=sum(cells_velo)/N。%平均速度 m=m+1。 endendflow(p)=(sum(speed)/30)*dens。%不同密度下的流量數(shù)組density(p)=dens。dens=dens+。p=p+1。endmean_wait_time = waiting_time/sum(output)mean_wait_time = 507s附錄六圖8的畫法程序clear。figure(39。position39。,[50 132 600 600],39。doublebuffer39。,39。on39。) T = 1。 % hrs to simulateline_width=4。%道寬line_length=40。%道長t=10。dt=70。allow_left = 1。%allow_left 取1表示允許左拐，取0則相反crossroads = create_crossroads(line_width,line_length)。%====================================================%為了顯示出圖形來，做了一個show_matrix %show_crossroads=show_matrix(crossroads,4)。 %H=image(show_crossroads)。 %axis off %%====================================================entry_vector = create_entry(T)。waiting_time = 0。output = 0。for i = 1:T*3600 crossroads=line_move(crossroads,line_width,line_length)。 crossroads = new_cars(crossroads, entry_vector(1,i), line_length, line_width, allow_left)。 %allow new cars to enter crossroads=center_move(crossroads,line_width,line_length)。 crossroads=line_to_center(crossroads, line_length, line_width)。 [crossroads,t]=red_light(crossroads,line_width,line_length,t,dt)。 waiting_time = waiting_time + pute_wait(crossroads)。 %pute waiting time during timestep i output = output + pute_output(crossroads,line_width , line_length )。 crossroads= clear_boundary(crossroads,line_width,line_length)。 %==================================================== show_crossroads=show_matrix(crossroads,4)。 % set(H,39。CData39。,show_crossroads)。 % pause() % %====================================================end mean_wait_time = waiting_time/sum(output)29
法律信息相關推薦
鄂ICP備17016276號-1