交通平均流量模拟（matlab）

function gaps = gaplength(x, L, ncar);
gaps = zeros(1, ncar);
gaps = x([2:end 1]) - x;
gaps(gaps <= 0) = gaps(gaps <= 0) + L;



function flux = ns(rho, p, L, tmax)
flux = 0;
ncar = round(rho * L);
x = sort(randperm(L, ncar));
vmax = 5;
v = vmax * ones(1, ncar);
for i = 1 : tmax
    v = min(v+1, vmax);
    gaps = gaplength(x, L, ncar);
    v = min(gaps - 1, v);
    v = max(v - rand(1, ncar) < p, 0);
    x = x + v;
    x(x > L) = x(x > L) - L;
    flux = flux + sum(v)/L;
end
avgflux = flux / tmax;