NaSch模型编程及数值模拟

模型参数取值：Lroad = 1000，p = 0.3，Vmax = 5。边界条件：周期性边界。输出结果：流量密度图，时空图。

function [new_matrix_cellsnew_v]=border_control（matrix_cellsabvvmax） %边界条件，开口边界，控制车辆出入
%出口边界，若头车在道路边界，则以一定该路0.9离去
n = length（matrix_cells）;
for j=1:n
    if road（j） ==1
        if j+v（j） <= n
            index = j+v（j）;
        else
            index = j+v（j） - n; %周期性边界
        end
        %碰撞检测
        if matrix_cells（index） == 1
            disp（‘Collision detected‘）
        end
        matrix_cells（index） = 1;
        v（index） = velocities（j）;
    end
end
new_matrix_cells=matrix_cells;
new_v=v;
% % 开口性边界
% if a==n
%     rand（‘state‘sum（clock）*rand（1））;
%     p_1=rand（1）; %产生随机概率
%     if p_1<=1 %如果随机概率小于0.9，则车辆离开路段，否则不离口
%         matrix_cells（n）=0;
%         v（n）=0;
%     end
% end
% %入口边界，泊松分布到达，1s内平均到达车辆数为q，t为1s
% if b>vmax
%     t=1;
%     q=0.25;
%     x=1;
%     p=（q*t）^x*exp（-q*t）/prod（x）; %1s内有1辆车到达的概率
%     rand（‘state‘sum（clock）*rand（1））;
%     p_2=rand（1）;
%     if p_2<=p
%         m=min（b-vmaxvmax）;
%         matrix_cells（m）=1;
%         v（m）=m;
%     end
% end
%周期性边界
% s = v（a）;
% matrix_cells（a）=0;
% v（a）=0;
% if a + s > n
%     if b > a + s - n
%         matrix_cells（a + s-n） = 1;
%         v（a+s-n） = s;
%     else
%         if b-1>0
%             matrix_cells（b-1）=1;
%             v（b-1）=s;
%         else
%             matrix_cells（n）=1;
%             v（n）=s;
%         end
%     end
% else
%     matrix_cells（a+s） = 1;
%     v（a+v（a）） = s;
% end
% new_matrix_cells=matrix_cells;
% new_v=v;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1622  2018-10-05 16:27  Nasch元胞自动机模型\border_control.m
     文件         204  2018-09-29 12:38  Nasch元胞自动机模型\leadcarupdate.m
     文件        3766  2018-10-14 19:44  Nasch元胞自动机模型\NaSch.m
     文件         168  2018-09-30 12:18  Nasch元胞自动机模型\randslow.m
     文件         256  2018-09-29 10:48  Nasch元胞自动机模型\roadstart.m
     文件         347  2018-09-29 12:22  Nasch元胞自动机模型\searchfrontcar.m
     文件         309  2018-09-29 12:36  Nasch元胞自动机模型\searchlastcar.m
     文件         225  2018-09-29 12:36  Nasch元胞自动机模型\searchleadcar.m
     文件         243  2018-09-29 12:36  Nasch元胞自动机模型\speedstart.m