二阶ODE初值问题调整：车辆跟驰场景下距离与减速参数引入方法

前置注意：单位统一修正 #

你当前参数存在单位不兼容的问题：

• 减速度da_lc = -3 mi/s²不符合物理常识，推测应为-3 mph/s
• 长度单位ft需要转换为英里（和速度单位mph匹配），1 ft = 1/5280 英里

核心修改点 #

1. 补充初始位置定义、单位转换逻辑 #

根据你给定的初始车头间距、车长参数，初始化前后车的初始位置，同时统一所有参数的单位：

import numpy as np
from scipy.integrate import solve_ivp

# 单位转换常量
FT_TO_MI = 1 / 5280  # 英尺转英里
# 车辆初始速度 单位：mph
v_lc0 = 30
v_fc0 = 30
# 前导车减速参数 单位：mph/s（修正原不合理的mi/s²）
da_lc = -3 
v_lc_stop = 0
# 灵敏度常数k 单位：1/s（原参数单位标注错误，k的量纲应为时间倒数）
k = 17
# 初始车头间距 单位：ft转英里
d0 = 500 * FT_TO_MI
# 车辆平均长度 单位：ft转英里（计算净间距时使用，车头间距计算可按需调用）
car_len = 15 * FT_TO_MI

# 初始位置定义：设后车初始位置为0，前导车初始位置 = 后车初始位置 + 初始车头间距
x_fc0 = 0
x_lc0 = x_fc0 + d0

2. 在ODE函数中动态计算前导车的实时速度、位置 #

前导车的运动状态是随时间变化的，需要在每个求解时刻先计算前导车的v_lc(t)和x_lc(t)，再代入二阶ODE的计算：

def dxfc_dt(t, y):
    x_fc, v_fc = y
    # 计算当前时刻前导车的速度：减速到0后保持静止
    t_stop = (v_lc_stop - v_lc0) / da_lc  # 前导车刹停所需时间
    if t < t_stop:
        v_lc = v_lc0 + da_lc * t
        # 匀减速位移公式，除以3600是把秒转小时，匹配mph的距离单位
        x_lc = x_lc0 + v_lc0 * t + 0.5 * da_lc * t * t / 3600
    else:
        v_lc = 0
        # 刹停后位置不变，先算刹停时的位移
        x_lc_stop = x_lc0 + v_lc0 * t_stop + 0.5 * da_lc * t_stop * t_stop / 3600
        x_lc = x_lc_stop
    
    # 代入你给出的二阶ODE推导后的一阶方程组
    d2xfc_dt2 = k * (v_lc - v_fc) / (x_lc - x_fc)
    return [v_fc, d2xfc_dt2]

3. 修正solve_ivp的调用参数 #

你原代码的solve_ivp参数顺序错误，需要正确传入时间范围、初始状态、求解时间点：

# 求解时间范围：0到10秒，可按需调整
t_span = [0, 10]
# 初始状态：[后车初始位置，后车初始速度]
y0 = [x_fc0, v_fc0]
# 求解时间点
time = np.linspace(t_span[0], t_span[1], 100)
# 调用求解器
z2 = solve_ivp(dxfc_dt, t_span, y0, t_eval=time)

逻辑说明 #

• 两车间距参数已经通过初始位置赋值、ODE中实时计算x_lc - x_fc的方式纳入求解逻辑
• 前导车减速参数已经通过实时计算前导车运动状态的逻辑纳入求解，自动处理刹停前/后的状态切换

内容的提问来源于stack exchange，提问作者qhumphrey