元胞自动机模型程序

%function NaSch02()
N=2000;
vmax=5;
timestep=500;
ss=1;
averagev=zeros();
flow=zeros();
avcardis=zeros();%000平均车间距
avrv=zeros();%000平均相对速度
avFF=zeros();%平均车辆作用力
for rho=0.05:0.05:1%车辆密度
cells=rand(1,N)<rho;
carposition=find(cells);%车辆位置初始化
carn=length(carposition);
carv=round(vmax*rand(1,carn));%速度初始化
cardis=zeros(1,carn);
rv=zeros(1,carn);%记录车辆先对速度
DMIN=4;D=40;FMAX=100;%
F=zeros(1,carn);M=2;%..........................................
%车间距、相对速度、作用力初始化
for mm=1:carn-1
cardis(mm)=carposition(mm+1)-carposition(mm)-1;
rv(mm)=carv(mm+1)-carv(mm);%rrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrr if DMIN<cardis(mm)<D
F(mm)=M*rv(mm)^2/cardis(mm);%..........................
else
if cardis(mm)>D
rv(mm)=0;
F(mm)=0;
else
rv(mm)=0;
F(mm)=FMAX;
end
end%..........................
end
mm=mm+1;
cardis(mm)=N+carposition(1)-carposition(mm)-1;
rv(mm)=carv(1)-carv(mm);%rrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrr rrrrrrr
if DMIN<cardis(mm)<D
F(mm)=M*rv(mm)^2/cardis(mm);%..........................
else
if cardis(mm)>D
rv(mm)=0;
F(mm)=0;
else
rv(mm)=0;
F(mm)=FMAX;
end
end%................................
%车间距、相对速度、作用力初始化完成
memor_cells=zeros(timestep,N);
memor_cellv=zeros(timestep,N);
memor_F=zeros(timestep,N);%............................
memor_cardis=zeros(timestep,N);%000
memor_rv=zeros(timestep,N);%000
sumv=0;pmin=0.1;pmax=0.3; flag=0;%
sumcardis=0;%00000000
sumrv=0;%000
for jj=1:timestep%开始更新
for ii=1:carn
if rv(ii)>=0
flag=1;
else
flag=0;
end
carv(ii)=min(carv(ii)+1,vmax);%加速更新
carv(ii)=min(carv(ii),cardis(ii));%减速
if flag==1;
pslow=pmin;
else
pslow=pmax;
end
if rand(1)<pslow
carv(ii)=max(carv(ii)-1,0);%随机慢化
end
if (jj>100)
sumv=sumv+carv(ii);%vvvvvvvvvvvvvvvvvvvvvvvvvvvv速度求和end
%carpositon变化
newposition=carposition(ii)+carv(ii);
if newposition<=N
carposition(ii)=newposition;
else
carposition(ii)=newposition-N;
end
memor_cells(jj,carposition(ii))=1;
end
%速度位置更新完，进行车间距、相对速度的更新
for mm=1:carn-1
if carposition(mm)<carposition(mm+1)
cardis(mm)=carposition(mm+1)-carposition(mm)-1;
if cardis(mm)<vmax%
rv(mm)=carv(mm+1)-carv(mm);%rrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrr rrrrrrrrrrrrr
else%
rv(mm)=0;
end%
if cardis(mm)>D||cardis(mm)<DMIN
rv(mm)=0;
end
%F(mm)=M*rv(mm)^2/cardis(mm);%。

else
cardis(mm)=N+carposition(mm+1)-carposition(mm)-1;
rv(mm)=carv(mm+1)-carv(mm);%rrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrr rrrrrrrrrrrrrrrrrrrrr
if cardis(mm)>D||cardis(mm)<DMIN
rv(mm)=0;
end
%
F(mm)=M*rv(mm)^2/cardis(mm);%。

end
end
mm=mm+1;
if carposition(mm)<carposition(1)
cardis(mm)=carposition(1)-carposition(mm)-1;
rv(mm)=carv(1)-carv(mm);%rrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrr rrrrrrrrrrrrrrrrrrrr
%F(mm)=M*rv(mm)^2/cardis(mm);%。

else
cardis(mm)=N+carposition(1)-carposition(mm)-1;
rv(mm)=carv(mm)-carv(1);%rrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrr
rrrrrrrrrrrrrrrrrrrrrr
if cardis(mm)>D||cardis(mm)<DMIN
rv(mm)=0;
end
% F(mm)=M*rv(mm)^2/cardis(mm);%。

end
%速度位置更新完，进行车间距、相对速度的更新完毕
for mm=1:carn-1
if (jj>100)
sumcardis=sumcardis+cardis(mm);%00000000车间距求和
sumrv=sumrv+rv(mm);%000
end
end
end%timestep时间步jj结束
cardistwo=sumcardis/(carn*(timestep-100));%000求平均车间距
rvtwo=sumrv/(carn*(timestep-100));%000求平均相对车速
FF=M*rvtwo^2/cardistwo;%000
vv=sumv/(carn*(timestep-100));%求平均车速
ww=rho*vv;
averagev(ss)=vv;
flow(ss)=ww;
avcardis(ss)=cardistwo;%0000
avrv(ss)= rvtwo;%
avFF(ss)= FF;%
ss=ss+1;
end%rho循环结束
rho=0.05:0.05:1;
plot(avrv,flow);
xlabel('avrv');
ylabel('flow');。