traffic1.m

关于交通流量的仿真

function traffic1(n) 
%       n是进站前的公路的格数
%    Example:   traffic(100); 
if nargin==0; n=100; end 
B=round(5*rand(9,n+2)); B([1,3,5,7,9],:)=1.2;
bb=B([2,4,6,8],:);bb(bb~=0)=1;B([2,4,6,8],:)=bb;
Z=B;H=imshow(Z,[]); 
set(gcf,'position',[241 132 560 420]) ;set(gcf,'doublebuffer','on'); 
xlabel('Please press "space" key and stop this program!',... 
 'fontsize',12,'color','r');          k=1; 
title('cellular-automation to traffic modeling','color','b'); 
while k; 	
s=get(gcf,'currentkey'); 
  if strcmp(s,'space'); clc;k=0; end 
  A=B([2,4,6,8],:);C=A(:,2:end);
  C(A(:,2:end)==0&A(:,1:end-1)==0)=0;%规则一 
  C(A(:,2:end)==1&A(:,1:end-1)==0)=0;%规则二
  C(A(:,2:end)==0&A(:,1:end-1)==1)=1;%规则三
  C=[round(5*rand(4,1)) C];C(C~=0)=1;B([2,4,6,8],:)=C;
  set(H,'CData',B); pause(0.1); 
end 
figure(gcf);
end