基于模型预测mpc汽车轨迹跟踪算法。 matlab和carsim联合仿真。 支持matlab2018和carsim2019。 理论pdf+源码
想让汽车在弯道里丝滑过弯?传统PID控制就像拿着遥控玩具车——方向盘打大了容易画龙,打小了又跟不上。这时候就得搬出模型预测控制(MPC)这个老司机,它能提前预判未来几步的走位,边开边调整方向盘。
!MPC轨迹跟踪示意图
先搞明白汽车怎么"思考"
在Matlab里搭建个二自由度单车模型,这就是MPC的"大脑"。来看看这个简化版动力学模型的核心代码:
function dx = vehicle_model(x, u) % 参数定义 m = 1573; % 整车质量kg Iz = 2873; % 转动惯量 lf = 1.1; % 前轴到质心距离 lr = 1.58; % 后轴到质心距离 Caf = 80000; % 前轮侧偏刚度 Car = 80000; % 后轮侧偏刚度 beta = atan((lr*tan(u(1)))/(lf+lr)); % 等效转向角 alpha_f = beta - x(4) - lf*x(2)/x(1); % 前轮侧偏角 alpha_r = beta - x(4) + lr*x(2)/x(1); % 后轮侧偏角 % 动力学方程 dx = zeros(4,1); dx(1) = x(2); % X方向速度 dx(2) = (Caf*alpha_f + Car*alpha_r)/m - x(3)*x(1); % Y方向加速度 dx(3) = (lf*Caf*alpha_f - lr*Car*alpha_r)/Iz; % 横摆角加速度 dx(4) = x(3); % 横摆角 end这段代码把汽车变成了微分方程——输入方向盘转角和当前状态,输出下一刻的车辆姿态。注意这里用的小角度近似在实际仿真中要当心,车速太低或者转向太猛会翻车!
CarSim和Matlab的梦幻联动
配置联合仿真时最容易卡在版本兼容性上。记住这两个关键操作:
- 在CarSim的Simulink模板里找到carsimss.slx
- 把S-Function模块拖到你的控制模型里,像这样连接信号:
% 加载CarSim数据文件 cs_loadfile('MPC_Tracking.par'); % 设置仿真时长 csport = csroot; csport.paramvalue = 10; % 10秒仿真 % 生成S-Function csbuild('MPC_Interface');注意检查Matlab的工作目录是否包含CarSim生成的datalink.sdf文件,不然会报"找不到路径"的错误。遇到过有人因为中文路径名折腾一整天,血的教训啊!
MPC控制器的灵魂代码
核心是求解这个滚动优化的QP问题:
function [u_opt, cost] = mpc_solver(x0, ref_traj) % 定义预测时域 N = 10; % 构造代价函数矩阵 Q = diag([10, 5, 2, 1]); % 状态权重 R = 0.1; % 控制量权重 % 生成约束矩阵 A_ineq = [1 0; -1 0]; % 方向盘转角限制 b_ineq = [30*pi/180; 30*pi/180]; % 调用quadprog求解 options = optimoptions('quadprog','Display','off'); [u_opt, cost] = quadprog(H, f, A_ineq, b_ineq, [], [], [], [], [], options); end这里有个调参小技巧:增大Q矩阵中横向误差的权重会让车辆更"强迫症"地贴紧参考线,但可能引发方向盘抽搐。实践中可以给方向盘转角加个低通滤波,就像老司机不会猛打方向一样。
仿真翻车现场实录
第一次跑联合仿真时出现的蛇形走位:
!仿真异常轨迹
问题出在预测时域太短——N=10相当于只预瞄1秒(0.1s步长)。改成N=20后车辆立即老实了,这说明做人啊,眼光还是得放长远点!
最后说个冷知识:CarSim的轮胎模型比我们用的简化版复杂100倍,所以实际调试时要给MPC加个转向速率限制,否则真车可能会把方向盘电机烧了。仿真世界里的完美控制,放到现实里可能分分钟教你做人。