Matlab实证ns模型,matlab ns模型 流量与密度关系有问题

是单车道NS模型，最后结果应该是凸形的图，但是我的结果是随着密度的增大而增大

不知道是不是求流量的公式有问题

clf

clear all

%CA setup

n=1000; %数据初始化

z=zeros(2,n); %元胞个数

dens=0.01;%密度初值

num_car=fix(dens*n);

z=roadstart(z,num_car); %道路状态初始化，路段上车辆随机分布

vmax=1; %最大速度

%mindu=zeros(1,20);%用于画基本图时的密度数组

%van=zeros(1,20);%平均流量数组

v=speedstart(z,vmax); %速度初始化

q=1;

while dens<=1 %用于画基本图

m=1; %开始求每个密度下的时空状态

v_ave=0; %用于求平均流量

%边界条件处理，搜素首末车，控制进出，使用开口条件

a=searchleadcar(z);

b=searchlastcar(z);

[z,v]=border_control(z,a,b,v,vmax);  %？？

i=searchleadcar(z); %搜索首车位置

for j=1:i

if i-j+1==n

[z,v]=leadcarupdate(z,v);

continue;

else

%=&#