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;
交通平均流量模拟(matlab)
最新推荐文章于 2023-01-07 21:24:33 发布