【模糊控制预测】协同自适应巡航控制MPC控制器设计【含Matlab源码 3818期】

💥💥💥💥💥💥💥💥💞💞💞💞💞💞💞💞💞Matlab领域博客之家💞💞💞💞💞💞💞💞💞💥💥💥💥💥💥💥💥
🚀🚀🚀🚀🚀🚀🚀🚀🚀🚀🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚀🚀🚀🚀🚀🚀🚀🚀🚀🚀

🔊博主简介：985研究生，Matlab领域科研开发者；
🏫个人主页：Matlab领域
🏆代码获取方式：
CSDN Matlab领域—代码获取方式

🚅座右铭：路漫漫其修远兮，吾将上下而求索。
更多Matlab信号处理仿真内容点击👇
①Matlab信号处理（高阶版）
②付费专栏Matlab信号处理（进阶版）
③付费专栏Matlab信号处理（初级版）

⛳️关注CSDN Matlab领域，更多资源等你来！！

⛄一、无人机圆轨迹MPC控制器设计简介

1 MPC简介
模型预测控制（以下简称 MPC）是一种依赖于系统模型进行数学优化的复杂控制器。它利用优化算法计算有限时间范围内一系列的控制输入序列，并优化该序列，但控制器仅执行序列中的第一组控制输入，然后再次重复该循环。MPC 主要分为 3 个关键步骤：模型预测、滚动优化、反馈校正。

2 MPC在轨迹跟踪过程中的应用
在轨迹跟踪应用领域，通常 MPC 建模可根据机器人的控制方式选择基于运动学运动状态方程建模或者基于动力学运动状态方程建模。前者是根据车辆转向的几何学角度关系和速度位置关系来建立描述车辆的运动的预测模型，一般只适用于低速运动场景；后者是对被控对象进行综合受力分析，从受力平衡的角度建模，一般应用在高速运动场景，如汽车无人驾驶。本文基于双轮差速运动学模型，推导实现差速运动的MPC轨迹跟踪控制。 双轮差速运动学线性 MPC 轨迹跟踪实现的基本思路是状态方程–线性化–离散化–预测方程–约束线性化–非线性目标函数转为二次规划–求解最优问题。其中线性化可采用泰勒级数展开忽略高次项的方式，离散化可采用后向差分法，这两个部分都有一定的模型精度丢失。

3 作用机理
MPC作用机理描述为：在每一个采用时刻，根据获得的当前测量信息，在线求解一个有限时间开环优化问题，并将得到的控制序列的第一个元素作用于被控对象。在下一个采样时刻，重复上述过程：用新的测量值作为此时预测系统未来动态的初始条件，刷新优化问题并重新求解 。
即MPC算法包括三个步骤：
（1）预测系统未来动态；
（2）（数值）求解开环优化问题；
（3）将优化解的第一个元素（或者说第一部分）作用于系统
这三步是在每个采样时刻重复进行的，且无论采用什么样的模型，每个采样时刻得到的测量值都作为当前时刻预测系统未来动态的初始条件
在线求解开环优化问题获得开环优化序列是MPC和传统控制方法的主要区别，因为后者通常是离线求解一个反馈控制律，并将得到的反馈控制律一直作用于系统。

⛄二、部分源代码

%求解有约束预测控制问题，函数返回值为系统输出及控制输入
%此函数针对无限预测时域与有限控制时域的情况
%function [Y,U]=cmpcQP()
clear;
%初始化
load dstatespace.mat;%下载辨识出的状态空间模型
%输入输出个数
usize=size(Bd,2);
ysize=size(C,1);
xsize=size(Ad,1);
YW=[10 0 0 0;0 6 0 0;0 0 3 0;0 0 0 0];%输出误差加权阵
% YW=[10 0;0 1];%输出误差加权阵
UW=0;%输入加权阵
m=5;%控制时域
p=10;%预测时域
Ts=0.1;%采样周期
Tend=30;%仿真总时间
Nrow=Tend/Ts+m;
LR1=5*ones(Nrow,1);%L1输出期望值
LR2=zeros(Nrow,1);%L2输出期望值
LR3=zeros(Nrow,1);%L2输出期望值
LR4=zeros(Nrow,1);%L2输出期望值
r=[LR1 LR2 LR3 LR4];
% r=[LR1 LR2];

L=zeros(Nrow,3);
for k=1:1:150
L(k,1)=0;
L(k,2)=3;
L(k,3)=3kTs+2;
end
for k=151:1:160
L(k,1)=-1;
L(k,2)=L(k-1,2)+L(k,1)*Ts;
L(k,3)=L(k-1,3)+L(k-1,2)Ts+0.5L(k,1)TsTs;
end
for k=161:1:Nrow
L(k,1)=0;
L(k,2)=L(k-1,2);
L(k,3)=L(k-1,3)+L(k-1,2)*Ts;
end

v=zeros(1,Nrow);
for i=20:1:24
v(1,i)=0.5;
end
for i=40:1:44
v(1,i)=-1.5;
end

%产生均值为a，标准差为b的白噪声
n=randn(1,Nrow);
n=n/std(n);
n=n-mean(n);
a=0; %均值为0
b=sqrt(0.01); %方差为0.01
n=a+b*n;
%dn=[n;n];
dn=n;
K=[0.1;0.1;0.1;0.1];

%输出约束
YMAX=ones(2m,1);
for i=1:2m
if rem(i,2)==0
YMAX(i,1)=3;
else
YMAX(i,1)=3;
end
end
YMIN=ones(2m,1);
for i=1:2m
if rem(i,2)==0
YMIN(i,1)=-3;
else
YMIN(i,1)=0;
end
end
YMAX=ones(4m,1);
for i=1:4m
if rem(i,4)==0
YMAX(i,1)=40;
elseif rem(i,3)==0
YMAX(i,1)=Inf;
elseif rem(i,2)==0&&rem(i,4)~=0
YMAX(i,1)=30;
else
YMAX(i,1)=10;
end
end
YMIN=ones(4m,1);
for i=1:4m
if rem(i,4)==0
YMIN(i,1)=0;
elseif rem(i,3)==0
YMIN(i,1)=-Inf;
elseif rem(i,2)==0&&rem(i,4)~=0
YMIN(i,1)=-30;
else
YMIN(i,1)=0;
end
end

⛄三、运行结果

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]季冬冬.MPC控制器的FPGA实现及其应用[J].吉林大学学报（信息科学版）. 2013,21(11)

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除

🍅 仿真咨询
1 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化

2 机器学习和深度学习方面
卷积神经网络（CNN）、LSTM、支持向量机（SVM）、最小二乘支持向量机（LSSVM）、极限学习机（ELM）、核极限学习机（KELM）、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断

3 图像处理方面
图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知

4 路径规划方面
旅行商问题（TSP）、车辆路径问题（VRP、MVRP、CVRP、VRPTW等）、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化

5 无人机应用方面
无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配

6 无线传感器定位及布局方面
传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化

7 信号处理方面
信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化

8 电力系统方面
微电网优化、无功优化、配电网重构、储能配置

9 元胞自动机方面
交通流 人群疏散 病毒扩散 晶体生长

10 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合