carfollowingmodel 在车联网研究之中

在车联网研究之中，对于先进的车辆跟驰模型的仿真，以及绘制时空图

%% Parameters
y_neutral=25; p_y=23.3; a=2; T=0.1; y_min=7.02; v_max=33.6; r=v_max/p_y; p=5;
%k_a=（a*T*（1-r*T）-1）/a*T*（1-2*r*T）; k_b=（（a*T*（1-r*T）-1）^3）/（（a^2）*（T^2）*（（1-2*r*T）^2））; 
%k_c=1; k_d=（（a*T*（1-r*T）-1）^2）/（a*T*（1-2*r*T））;   %加反馈
%% 
v_sta=20; y_sta=v_sta/（v_max/p_y）-p_y/2+y_neutral;
x=zeros（110051）; y=x; v=x; x_re=x; %51辆车，1100个T（110s）程序中i车实际上是编号为51-i的车
v（:51）=20;
v（100:12951）=0; % 头车在100-102s停
v（1:）=20;  %初始时刻速度都为20m/s
x（1:）=2*y_sta:y_sta:52*y_sta;   %初始时刻位置，末车设为2*y_sta，即以末车后2*y_sta的点为绝对坐标零点
y（11:50）=y_sta*ones（150）; y（:51）=0;   %初始时刻车间距为稳态值y_sta
x_re（1:）=x（151）*ones（151）-x（1:）;   %与头车的距离

for t=2:1100
    x（t:）=v（t-1:）*T+x（t-1:）;    %t时刻每辆车的绝对位移
    x_re（t:）=x（t51）*ones（151）-x（t:）;  %t时刻每辆车与头车的距离
  
    for i=1:50
        y（ti）=x（ti+1）-x（

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1512  2020-10-01 05:13  CM_kkh_50.m
     文件         260  2020-10-01 05:13  Vop.m