电脑桌面
添加小米粒文库到电脑桌面
安装后可以在桌面快捷访问

交通流元胞自动机程序

交通流元胞自动机程序_第1页
1/7
交通流元胞自动机程序_第2页
2/7
交通流元胞自动机程序_第3页
3/7
% 车流密度不变下的多车道仿真(包括单车道) % nc:车道数目(1 或2),nl:车道长度 % v:平均速度,d:换道次数(1000 次)p:车流密度 % dt:仿真步长时间,nt:仿真步长数目 % fp:车道入口处新进入车辆的概率 v = 0;p=0;d=0; nl = 100;nc = 1; dt=0.01;nt=1000; fp = 0.5; [ v d p ] = multi_driveway( nl,nc,fp,dt,nt ); function [ v d p ] = multi_driveway( nl,nc,fp,dt,nt ) % 在某一特定车流密度下的(车流密度由fp 决定)单、双车道仿真模型 % nc:车道数目(1 或2),nl:车道长度——输入参数 % v:平均速度,d:换道次数(1000 次)p:车流密度——输出参数 % dt:仿真步长时间,nt:仿真步长数目——输入参数 % fp:车道入口处新进入车辆的概率——输入参数 % test: % nl = 400;fp = 0.5; % nc = 2;dt=0.01;nt=500; %构造元胞矩阵 B=ones(2*nc+1,nl+2); %奇数行为不可行车道 B(1:2:(2*nc+1),:)=1.2; %初始化仿真元胞状态(1 为无车,0 为有车) bb=B(2:2:2*nc,:);bb(bb~=0)=1;B(2:2:2*nc,:)=bb;B(2:2:2*nc,end)=0; %显示初始交通流图 figure(1); H=imshow(B,[]); set(gcf,'position',[241 132 560 420]) ;%241 132 560 420 set(gcf,'doublebuffer','on'); %241 title('cellular-automation to traffic modeling','color','b'); %初始化化存储元胞上车辆状态的矩阵 S(1:nc,nl) = 0; Q(1:nc,1:2) = 0; Acc(1:nc,1:(nl+2))=0; %初始化换道频率、平均速度、车流密度相关变量 ad = 0; av(1:nt) = 0; ap(1:nt) = 0; c = 1; for n = 1:nt A=B(2:2:2*nc,:); %确定前n-2 个车辆的状态 S(:,:) = 0; S(A(:,1:end-2)==0&A(:,2:end-1)==1&A(:,3:end)==1)=2;%加速的车 S(A(:,1:end-2)==0&A(:,2:end-1)==0)=3;%停车的车 S(A(:,1:end-2)==0&A(:,2:end-1)==1&A(:,3:end)==0)=1;%减速行驶的车 %确定最后2 两个元胞的状态 Q(:,:) = 0; Q(A(:,end-1)==0&A(:,end)==0) = 1; Q(A(:,end-1)==0&A(:,end)==1) = 2; Q(A(:,end-1)==1&A(:,end)==0) = 2; Q(:,end) = 1; %获得所有元胞上车辆的状态 Acc = [ S Q ]; %换路规则 if(nc>1&&n>nl/2) %遍历每一个元胞 for g = 1:length(Acc(1,:)) %停车状态车辆如另一条路有2 空位则换路 if( Acc(1,g)==3&&Acc(2,g)==0&&Acc(2,g+1)==0) A(1,g)=1; A(2,g)=0; a...

1、当您付费下载文档后,您只拥有了使用权限,并不意味着购买了版权,文档只能用于自身使用,不得用于其他商业用途(如 [转卖]进行直接盈利或[编辑后售卖]进行间接盈利)。
2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。
3、如文档内容存在违规,或者侵犯商业秘密、侵犯著作权等,请点击“违规举报”。

碎片内容

交通流元胞自动机程序

确认删除?
VIP
微信客服
  • 扫码咨询
会员Q群
  • 会员专属群点击这里加入QQ群
客服邮箱
回到顶部