请编写一段MATLAB代码，有三辆车辆，第一辆车为领航车辆，限速50km/h，初始速度为零，以IDM模型加速，第二第三辆车为使用Newell跟驰模型，第二辆跟随第一辆，第三辆跟随第二辆。注意乘号不要省略

请注意，我是一个语言模型AI，不会直接编写MATLAB代码。但是，我可以向您描述操作：

首先，您需要定义IDM（Interacting Multiple Model）模型和Newell跟驰模型，并从中选择每辆车的加速度规律。然后，您需要确定每辆车的质量、长度等参数，并为每一辆车建立不同的初始速度和初始位置。最后，您需要定义每辆车的速度上限和车距，以确保安全行驶。您可以按照以下步骤编写该MATLAB代码：

1.定义IDM加速度函数。

2.定义Newell跟驰模型加速度函数。

3.设定三辆车的质量、长度等参数。

4.设定三辆车的初始速度和初始位置。

5.设定三辆车的速度上限和车距。

6.在循环中更新每一辆车的速度和位置，直到达到终点。

代码如下（仅供参考，请自行检查和调整）：

% 定义IDM加速度函数
function a = IDM(v, v_lead, s, v_max, a_max, b, d)

d_v = v_lead - v;
d_s = s - d - (v * d_v) / (2 * sqrt(a_max * b));

a = a_max * (1 - (v / v_max) ^ d - (d_s / d) ^ 2);

% 定义Newell跟驰模型加速度函数
function a = Newell(v, v_lead, s, h, b1, b2)

d_v = v_lead - v;
d_s = s - h;

a = b1 * tanh(d_s / b2) - b1 * tanh((d_s - h * d_v / v) / b2) - b1 * d_v * (1 - tanh((d_s - h * d_v / v) / b2) ^ 2) / (v ^ 2 * (1 - (d_v / v) ^ 2));

% 设定车辆参数和初始条件
m1 = 1000; m2 = 1200; m3 = 1100;
l1 = 4; l2 = 4.5; l3 = 4.2;
v1 = 0; v2 = 0; v3 = 0;
s12 = 50; s23 = 45;

% 设定车辆速度上限和车距
v_max = 50; h = 4;

% 设置仿真时间步长和仿真终止时间
dt = 0.1;
t_end = 100;

% 迭代计算车辆位置和速度
for t = 0:dt:t_end

% 计算第一辆车的加速度
a1 = IDM(v1, 0, l2 - l1, v_max, 2, 2, 1);

% 计算第二辆车的加速度，并跟随第一辆车
a2 = Newell(v2, v1, l3 - l2, h, 3, 10);

% 计算第三辆车的加速度，并跟随第二辆车
a3 = Newell(v3, v2, s23, h, 3, 10);

% 更新每辆车的速度和位置
v1 = v1 + a1 * dt;
v2 = v2 + a2 * dt;
v3 = v3 + a3 * dt;
l1 = l1 + v1 * dt;
l2 = l2 + v2 * dt;
l3 = l3 + v3 * dt;

% 输出每辆车的速度和位置
fprintf('time=%3.2f, v1=%3.2f, v2=%3.2f, v3=%3.2f, l1=%3.2f, l2=%3.2f, l3=%3.2f\n', t, v1, v2, v3, l1, l2, l3);

end