双圆道交通流CA模型Matlab程序.docVIP

%function [rho, flux, vmean] = ns2(rho, p, L, tmax) %注意：主函数一定不能像子函数一样以function开头。此程序为含控制按钮NS模型改编双车道，无换道。 % NS: This script implements the Nagel Schreckenberg cellular automata based % traffic model. Car move forward governed by NS algorithm: % % 1. Acceleration. If the vehicle can speed up without hitting the speed % limit vmax it will add one to its velocity, vn - vn + 1. Otherwise, % the vehicle has constant speed, vn 鈫? vn. % % 2. Collision prevention. If the distance between the vehicle and the car % ahead of it, dn, is less than or equal to vn , i.e. the nth vehicle % will collide if it doesnt slow down, then vn 鈫? dn 鈭? 1. % % 3. Random slowing. Vehicles often slow for non-traffic reasons (cell % phones, coffee mugs, even laptops) and drivers occasionally make % irrational choices. With some probability pbrake, vn =? vn -1, % presuming vn 0. % % 4. Vehicle movement. The vehicles are deterministically moved by their % velocities, xn - xn + vn . % % USAGE: flux = ns(rho, p, L, tmax, isdraw) % rho = density of the traffic % p = probability of random braking % L = length of the load % tmax = number of the iterations % animation = if show the animation （动画）of the traffic % spacetime = if plot the space-time after the simuation ended. % flux = flux of the traffic % % zhou lvwen: zhou.lv.wen@ rho = 0.25; p = 0.25; L = 100; tmax = 100; i=0; %============================================= %build the GUI %define the plot button plotbutton=uicontrol(style,pushbutton,... string,Run, ... fontsize,12, ... position,[60,400,50,20], ... callback, run=1;); %define the stop button erasebutton=uicontrol(style,pushbutton,... string,Stop, ... fontsize,12, ... position,[120,400,50,20], ... callback,freeze=1;); %============================================= %define the Quit button quitbutton=uicontrol(style,pushbutton,...