FVD跟驰模型Python数值仿真

-- coding: utf-8 --

“”"
Created on Fri Oct 15 10:17:28 2021

@author: Libertas
“”"
import openpyxl
wb=openpyxl.Workbook()
wb.save(‘test.xlsx’)#新建一个Excel用于保存车辆信息
from openpyxl import load_workbook
import numpy as np

wb = load_workbook(‘test.xlsx’)#借助Excel保存车辆位置信息

wb1 = wb.active

car_x=[20.109*(20-i) for i in range(20)]#利用三个列表存放车辆位置、速度、加速度，一一对应
car_v=[20 for i in range(20)]
car_a=[0 for i in range(20)]

for ms in range(1,101):#0.1s为一个步长，前10s
car_v[0]+=0.1car_a[0]
car_x[0]+=0.1car_v[0]
wb1.cell(ms,1,car_x[0])
for i in range(1,20):
car_a[i]=0.4*(14+16np.tanh(0.13(car_x[i-1]-car_x[i]-5)-1.57)-car_v[i])+0.5*(car_v[i-1]-car_v[i])
car_v[i]+=0.1car_a[i]
car_x[i]+=0.1car_v[i]
wb1.cell(ms,i&