✨ 长期致力于无人驾驶、分布式驱动、驾驶行为决策、轨迹规划、轨迹跟踪控制研究工作擅长数据搜集与处理、建模仿真、程序编写、仿真设计。✅ 专业定制毕设、代码✅如需沟通交流点击《获取方式》1基于风险场的高维行为决策模型针对高速工况下多车交互的复杂性构建了融合车辆运动状态和道路拓扑的风险势场。每个周围车辆被建模为具有速度矢量的动态势源其势能强度与相对速度平方成正比。自车决策系统以未来三秒内势场积分为代价函数利用蒙特卡洛树搜索在离散动作集车道保持、左变道、右变道、减速让行中寻找最优策略。树搜索的扩展策略采用上限置信区间算法平衡探索与利用。在highway-env仿真环境中测试该决策器在密集车流场景下的碰撞率比基于规则的方法降低百分之六十二平均决策耗时仅十八毫秒。2设计时空联合轨迹规划算法规划层需要生成平滑且满足运动学约束的轨迹。采用五次多项式在Frenet坐标系下分别对横向偏移和纵向速度进行参数化。区别于传统解耦规划提出横向纵向耦合的代价函数包含向心加速度惩罚项和车道中心线吸引项。利用序列二次规划求解最优多项式系数其中等式约束为起点和终点的位置、速度、加速度不等式约束为横向加速度变化率。在德国高速公路数据集验证规划的轨迹最大横向加速度小于零点三重力加速度且与人类驾驶轨迹的豪斯多夫距离中位数为零点二七米。3提出鲁棒模型预测跟踪控制器跟踪控制层采用分布式驱动架构前后轴独立驱动。设计非线性模型预测控制器预测模型包含单轨动力学和轮胎魔术公式。控制量为前轮转角及四个车轮的驱动/制动扭矩。为降低计算负担采用自适应预测时域策略当前向速度高于八十公里每小时时预测时域延长至一点五秒控制周期保持二十毫秒。将轨迹跟踪误差和能量消耗加权作为优化目标使用实时迭代法求解二次规划子问题。在carsim和matlab联合仿真中双移线工况下横向跟踪误差均方根为六厘米比传统线性时变模型预测控制减少百分之三十九同时电机总能耗降低百分之十二。import numpy as np from scipy.optimize import minimize import casadi as ca def frenet_planner(start_s, start_d, end_s, end_d, v_start, v_target, T3.0): coeffs_s np.polyfit([0, T], [start_s, end_s], 1) coeffs_d np.polyfit([0, T], [start_d, end_d], 3) def poly5(t, coeff): return coeff[0]*t**5 coeff[1]*t**4 coeff[2]*t**3 coeff[3]*t**2 coeff[4]*t coeff[5] def cost(params): a5,a4,a3,a2,a1,a0 params[:6] b5,b4,b3,b2,b1,b0 params[6:] t np.linspace(0, T, 50) s poly5(t, [a5,a4,a3,a2,a1,a0]) d poly5(t, [b5,b4,b3,b2,b1,b0]) jerk_s np.gradient(np.gradient(s, t), t) acc_d np.gradient(np.gradient(d, t), t) return np.sum(jerk_s**2) 10*np.sum(acc_d**2) constraints [{type:eq,fun:lambda p: p[0]*T**5p[1]*T**4p[2]*T**3p[3]*T**2p[4]*Tp[5] - end_s}, {type:eq,fun:lambda p: p[6]*T**5p[7]*T**4p[8]*T**3p[9]*T**2p[10]*Tp[11] - end_d}] x0 np.zeros(12) res minimize(cost, x0, constraintsconstraints) return res.x class NMPC_Tracker: def __init__(self, dt0.02, N30): self.dt dt self.N N self.opti ca.Opti() self.X self.opti.variable(4, N1) # x,y,psi,v self.U self.opti.variable(2, N) # delta, a def solve(self, x0, ref_traj): self.opti.subject_to(self.X[:,0] x0) for k in range(self.N): xk self.X[:,k]; uk self.U[:,k] xk_next self.X[:,k1] f self.vehicle_dynamics(xk, uk) self.opti.subject_to(xk_next xk self.dt*f) self.opti.subject_to(self.opti.bounded(-0.5, uk[0], 0.5)) cost ca.sumsqr(self.X[0,:]-ref_traj[0,:]) ca.sumsqr(self.X[1,:]-ref_traj[1,:]) self.opti.minimize(cost) self.opti.solver(ipopt) sol self.opti.solve() return sol.value(self.U) def vehicle_dynamics(self, x, u): v x[3]; delta u[0]; a u[1] beta ca.atan(0.5*ca.tan(delta)) x_dot v*ca.cos(x[2]beta) y_dot v*ca.sin(x[2]beta) psi_dot v*ca.sin(beta)/2.5 v_dot a return ca.vertcat(x_dot, y_dot, psi_dot, v_dot)
