求助matlab

N=100;%车辆数量
L=1000;%道路长度
vmax=120/3.6;%最大速度
rhomax=N/L;%最大密度                                                                                     
T=1.5;% 安全时间头车间
a=1.2;%加速度
b=2; % 减速度
dt =0.1;% 时间步长
tmax=100;%最大时间
t=0:dt:tmax;% 时间向量
x=zeros(N,length(t));%位置矩阵
v=zeros(N,length(t));% 速度矩阵
%初始化位置和速度
x(:,1)= L/N*(0:N-1)';
v(:,1)= vmax*ones(N,1);
% FVD模型的主循环
for i = 2:length(t)
   for j = 1:N
     if j == N
       dx = L+x(1,i-1)-x(j,i-1);
       dv = v(1,i-1)-v(j,i-1);
     else
       dx = x(j+1,i-1)-x(j,i-1);
       dv = v(j+1,i-1)-v(j,i-1);
     end
       vsafe = dx/T - dv;
       vfree = vmax*(1 - rhomax*dx/N);
       v(j,i) = v(j,i-1) + dt*(a*(vsafe - v(j,i-1)) + x(j,i)==x(j,i-1) + dt*v(j,i);这里说错误，该怎么解决啊
   end                                                                                                                                                                    
end
figure;
plot(t,x);
xlabel('Time (s)');
ylabel('Position (m)');
title('FVD Model Simulatio')