Продольная динамика автомобиля

# --- Извлечение сигналов из результата модели m1_longitudinal_vehicle ---
v          = data["Longitudinal Vehicle.v"]
accel      = data["Longitudinal Vehicle.a"]
v_rel      = data["Longitudinal Vehicle.v_headwind"]

speed_kmh  = data["Скорость автомобиля, км/ч"]
wind_ms    = data["Скорость ветра, м/с"]
slope_rad  = data["Уклон дороги, рад"]

f_total    = data["Longitudinal Vehicle.F_total"]
f_drive    = data["Longitudinal Vehicle.F_drive"]
f_resist   = data["Longitudinal Vehicle.F_resistance"]
f_aero     = data["Longitudinal Vehicle.C_d_force"]

clutch_A   = data["Муфта A.fundamental_friction_clutch.T"]
clutch_B   = data["Муфта B.fundamental_friction_clutch.T"]
clutch_C   = data["Муфта C.fundamental_friction_clutch.T"]
clutch_D   = data["Муфта D.fundamental_friction_clutch.T"]
clutch_E   = data["Муфта E.fundamental_friction_clutch.T"]

k_rpm      = 60/(2*pi)/1000
w_engine   = data["Гидротрансформатор.impeller_flange.w"]
w_ring     = data["Ravigneaux Gear.ring_flange.w"]
w_carrier  = data["Ravigneaux Gear.carrier_flange.w"]
w_wheel    = data["Longitudinal Vehicle.w"]

gear       = data["Номер передачи"]
gratio     = data["Передаточное число коробки передач"]
k_speed    = data["Гидротрансформатор.converter.speed_ratio"]
k_torque   = data["Гидротрансформатор.converter.torque_ratio"]
fuel       = data["Расход топлива, л/100км"]

gr()

p1 = plot(speed_kmh.time,
    [3.6 .* wind_ms.value (180/pi) .* slope_rad.value],
    xlabel="Время, с", ylabel="Скорость, км/ч  ·  Уклон, град",
    label=["Скорость ветра" "Уклон дороги"],
    title="Условия движения",
    legend=:bottomright,
    left_margin = 10Plots.Measures.mm)

p2 = plot(gear.time, gear.value,
    xlabel="Время, с", ylabel="Номер передачи",
    label="Номер передачи",
    title="Переключение передач",
    left_margin = 10Plots.Measures.mm)

p3 = plot(v.time, [v.value accel.value],
    xlabel="Время, с", ylabel="Скорость, м/с  ·  Ускорение, м/с²",
    label=["Скорость автомобиля" "Ускорение"],
    title="Кинематика автомобиля",
    left_margin = 10Plots.Measures.mm)

p4 = plot(fuel.time, fuel.value,
    xlabel="Время, с", ylabel="Расход топлива, л/100 км",
    label="Расход топлива", ylimits=(0,20),
    title="Топливная экономичность",
    left_margin = 10Plots.Measures.mm)

plot(
    p1, p2,
    p3, p4,
    layout=(2,2), size=(750,700), lw=2,
    titlefontsize=10, guidefontsize=9, legendfontsize=8
)