//?2?Roads?Merging?NaSch?CA?model
//?By?David?Clarke?（?dave@maths.tcd.ie?）
//?Original?source?available?at?www.maths.tcd.ie/~dave/traffic/
//Compile?with?cc?-o?traffic?traffic.c

#include
#include
#include

#define?leadin_d?500?
#define?overlap_d?10
#define?leadout_d?500
#define?max_length?2100??//?greater?then?leadin+overlap+leadout
#define?speed_limit1_d?5
#define?speed_limit2_d?5
#define?steps_d?20000?		//Number?of?Steps?updates
#define?cars1_d?0		//?cars?to?start?with?（0?is?empty?road）
#define?cars2_d?0
#define?cars0_d?0?//?leadout?cars
#define?prob_d?0.5		//?probability?of?stopping
#define?slow_start_d?0.0
#define?relax_d?600		//?steps?ignore?before?measuring
#define?flux_time_d?2000		//steps?flux?measured?over
#define?print_d?0
#define?split?10000000
#define?once?0?			//?if?1?run?for?only?one?density?for?scatter?plot.
#define?max_density?1?	//?limiter?0->1?can?stop?program?for?a?lower?density?

typedef?struct?{
??int?car;
??int?vel;
??int?moved;
??int?pre_vel;			//previous?velocity?used?in?slow?to?start
??int?jtime;
??int?origin;
??int?num;
}?road_t;
FILE?*out1;
FILE?*outScatter1?*outScatter2;
FILE?*out_jtime1?*out_jtime2;

//Function?definitions
void?print_road?（int）;
void?print_jams?（）;
void?print_flow_density（）;
void?print_scatter?（int?c?int?tstep）;
void?setup?（int）;
void?setup_leadout?（）;
void?setspeed?（int）;
void?update?（int?c?int?tstep）;
void?measure_density（）;
void?variables?（int?argc?char?*argv[]）;

//global?variables
road_t?Road1[max_length]Road2[max_length];
int?leadinoverlapleadout?speed_limit[3]?steps?cars[3]?relax?density_width?flux_timetot_lengthprint;
double?prob?slow_start;
int?count[3]density_start[3]density_stop[3]?time_count?density[3]?jtimeMax[3]jtimeMin[3];

int?main?（int?argc?char?*argv[]）
{		/*?Start?up?*/
??int?ij?n?dpr;
??road_t?*roadP;
??variables?（argc?argv）;	//take?in?variables?from?command?line
??tot_length=leadin+overlap+leadout;
//density?measurement?parameters.
??density_start[1]=density_start[2]=leadin/2;
??density_start[0]=leadin+overlap+leadout/2;
??density_stop[1]=density_stop[2]=leadin;
??density_stop[0]=leadin+overlap+leadout;
??out1?=?fopen?（“density_flow.dat“?“w“）;
??outScatter1?=?fopen?（“trafficYScatter1.dat“?“w“）;
??outScatter2?=?fopen?（“trafficYScatter2.dat“?“w“）;
??out_jtime1?=?fopen?（“jtime1.dat“?“w“）;
??out_jtime2?=?fopen?（“jtime2.dat“?“w“）;
??srand?（time?（0））;		//?random?seed


while（cars[1]		/*?Set?Up?*/
	n?=?0;
	setup?（1）;			//?initial?positioning?of?cars
	setup?（2）;			//?initial?positioning?of?cars
	setup_leadout（）;
	count[1]=count[2]=1;
	if（print）{	print_road（1）;???????//?initial?print
		print_road（2）;???????//?initial?print
	}
		/*?Relax?*/	
	while（n?		setspeed（1）;?
		update（1-1）;?
		setspeed（2）;?
		update（2-1）;}???

			//	Give?each?car?a?number:
	for（i=0;i		if（Road1[i].c

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????14553??2008-10-16?17:05??tarc\traffic.c

?????目錄??????????0??2009-03-04?21:49??tarc

-----------?---------??----------?-----??----

????????????????14553????????????????????2