【MPC学习笔记】01:MPC简介(Lecture 1_1 Unconstrained MPC)

2024-01-03 21:52

本文主要是介绍【MPC学习笔记】01:MPC简介(Lecture 1_1 Unconstrained MPC),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

本笔记来自北航诸兵老师的课程
课程地址:模型预测控制(2022春)lecture 1-1 Unconstrained MPC

文章目录

  • 0 MPC 简介
    • 0.1 案例引入
    • 0.2 系统模型
    • 0.3 MPC的优点
    • 0.4 MPC的缺点
    • 0.5 MPC的未来
  • 1 详细介绍

0 MPC 简介

0.1 案例引入

MPC(Model Predictive Control)模型预测控制,是预测控制的一种,是基于模型来进行控制的。
老师举了下面这个例子来引入MPC的基本思想:

比方说我们为未来的一段时间制定计划,一天中几点到几点该做什么。但是计划赶不上变化,出现变化,出现拖延,计划就得做相应的调整。过了一段时间,根据计划的实际落实情况,再对接下来的计划进行调整。如此往复。不断地执行计划,也不断地修订计划。

0.2 系统模型

在控制系统中,有惯用表示:输入记作 u u u,状态变量记作 x x x,输出记作 y y y
假设系统是离散的,系统的状态方程为:
x ( k + 1 ) = f ( x ( k ) , u ( k ) ) x(k+1)=f(x(k),u(k)) x(k+1)=f(x(k),u(k))

实际上系统可以是,线性的或非线性的,连续的或离散的或既包含连续又包含离散的,确定的或随机的,只要满足该方程即可

设当前时刻为 k k k,当前状态为 x ( k ) x(k) x(k)
在输入 u ( k ) u(k) u(k) 的作用下,系统的状态将由 x ( k ) x(k) x(k) 变为 x ( k + 1 ) x(k+1) x(k+1)
在输入 u ( k + 1 ) u(k+1) u(k+1) 的作用下,系统的状态将由 x ( k + 1 ) x(k+1) x(k+1) 变为 x ( k + 2 ) x(k+2) x(k+2)
在输入 u ( k + 2 ) u(k+2) u(k+2) 的作用下,系统的状态将由 x ( k + 2 ) x(k+2) x(k+2) 变为 x ( k + 3 ) x(k+3) x(k+3)

由上面的列举,知:输入序列➡️输出序列
但在此时,也就是时刻 k k k ,我们并不知道输入序列 { u ( k ) , u ( k + 1 ) , u ( k + 2 ) , ⋯ } \{u(k),u(k+1),u(k+2),\cdots\} {u(k),u(k+1),u(k+2),} 是多少

自然而然就会想到一个问题——怎么确定输入序列?
答:通过优化的方式,Optimization
状态序列 记为 X ( k ) X(k) X(k)
输入序列 记为 U ( k ) U(k) U(k)
输入序列的求解,可用如下优化问题的公式来描述:
U ∗ ( k ) = a r g m i n ∑ i = k ∞ l ( x ( i ) , u ( i ) ) = { u ∗ ( k ) , u ∗ ( k + 1 ) , … } s . t . x ∈ X , u ∈ U \begin{aligned} U^*(k) &= arg\ min\sum^{\infin}_{i=k}l(x(i),u(i)) \\ &=\{u^*(k),u^*(k+1),\dots\} \\ \\ s.t.\quad &x\in \mathscr{X}, u\in \mathscr {U} \end{aligned} U(k)s.t.=arg mini=kl(x(i),u(i))={u(k),u(k+1),}xX,uU
其中, a r g m i n arg\ min arg min 表示使 **代价函数(目标函数)**取值最小时,输入序列 U ( k ) U(k) U(k) 的取值; ∗ ^* 表示最优解; s . t . s.t. s.t. 表示约束条件; l ( x ( i ) , u ( i ) ) l(x(i),u(i)) l(x(i),u(i)) 称为 “Stage cost”。
u ( k ) = u ∗ ( k ) u(k)=u^*(k) u(k)=u(k) ,舍弃求出的 U ∗ ( k ) U^*(k) U(k) 中后续其他时刻的输入,则由 x ( k + 1 ) = f ( x ( k ) , u ( k ) ) x(k+1) = f(x(k),u(k)) x(k+1)=f(x(k),u(k)) 可以求出时刻 k + 1 k+1 k+1 的状态
接着, k + 1 k+1 k+1 变为当前时刻,重复上述步骤,求出时刻 k + 2 k+2 k+2 的状态 x ( k + 2 ) x(k+2) x(k+2),…
以上就是MPC的基本原理

如果只优化一次,将计算出的 U ( k ) U(k) U(k) 序列依次执行,那么就变成了开环优化;而这里每一时刻优化后都只取 u ∗ ( k ) u^*(k) u(k) 执行( u ∗ ( k ) u^*(k) u(k) x ( k ) x(k) x(k)的函数),并且不断进行优化,构成滚动优化(闭环优化), 因此MPC实际上引入了反馈

0.3 MPC的优点

  • 处理控制输入和系统状态上的约束(Constraints)
    • 约束来源:actuator limits; safety; environmental; economic constraints
    • PID没办法解决约束问题
  • 近似最优控制
    • 与线性系统中的最优控制(LQR, 线性二次型调节器)有区别,在LQR中,我们找到的是最优的增益 k k k(假设,已知系统是线性反馈),MPC找的是 u u u

0.4 MPC的缺点

  • 需要在线优化(online optimization),可能会有较大的计算负载

0.5 MPC的未来

随着计算机算力提升,MPC或替代PID成为工业界控制主流
在这里插入图片描述

1 详细介绍

见【MPC学习笔记】02:MPC详细简介(Lecture 1_1 Unconstrained MPC)

这篇关于【MPC学习笔记】01:MPC简介(Lecture 1_1 Unconstrained MPC)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/567169

相关文章

Java深度学习库DJL实现Python的NumPy方式

《Java深度学习库DJL实现Python的NumPy方式》本文介绍了DJL库的背景和基本功能,包括NDArray的创建、数学运算、数据获取和设置等,同时,还展示了如何使用NDArray进行数据预处理... 目录1 NDArray 的背景介绍1.1 架构2 JavaDJL使用2.1 安装DJL2.2 基本操

Golang的CSP模型简介(最新推荐)

《Golang的CSP模型简介(最新推荐)》Golang采用了CSP(CommunicatingSequentialProcesses,通信顺序进程)并发模型,通过goroutine和channe... 目录前言一、介绍1. 什么是 CSP 模型2. Goroutine3. Channel4. Channe

Java中的Opencv简介与开发环境部署方法

《Java中的Opencv简介与开发环境部署方法》OpenCV是一个开源的计算机视觉和图像处理库,提供了丰富的图像处理算法和工具,它支持多种图像处理和计算机视觉算法,可以用于物体识别与跟踪、图像分割与... 目录1.Opencv简介Opencv的应用2.Java使用OpenCV进行图像操作opencv安装j

HarmonyOS学习(七)——UI(五)常用布局总结

自适应布局 1.1、线性布局(LinearLayout) 通过线性容器Row和Column实现线性布局。Column容器内的子组件按照垂直方向排列,Row组件中的子组件按照水平方向排列。 属性说明space通过space参数设置主轴上子组件的间距,达到各子组件在排列上的等间距效果alignItems设置子组件在交叉轴上的对齐方式,且在各类尺寸屏幕上表现一致,其中交叉轴为垂直时,取值为Vert

Ilya-AI分享的他在OpenAI学习到的15个提示工程技巧

Ilya(不是本人,claude AI)在社交媒体上分享了他在OpenAI学习到的15个Prompt撰写技巧。 以下是详细的内容: 提示精确化:在编写提示时,力求表达清晰准确。清楚地阐述任务需求和概念定义至关重要。例:不用"分析文本",而用"判断这段话的情感倾向:积极、消极还是中性"。 快速迭代:善于快速连续调整提示。熟练的提示工程师能够灵活地进行多轮优化。例:从"总结文章"到"用

【前端学习】AntV G6-08 深入图形与图形分组、自定义节点、节点动画(下)

【课程链接】 AntV G6:深入图形与图形分组、自定义节点、节点动画(下)_哔哩哔哩_bilibili 本章十吾老师讲解了一个复杂的自定义节点中,应该怎样去计算和绘制图形,如何给一个图形制作不间断的动画,以及在鼠标事件之后产生动画。(有点难,需要好好理解) <!DOCTYPE html><html><head><meta charset="UTF-8"><title>06

学习hash总结

2014/1/29/   最近刚开始学hash,名字很陌生,但是hash的思想却很熟悉,以前早就做过此类的题,但是不知道这就是hash思想而已,说白了hash就是一个映射,往往灵活利用数组的下标来实现算法,hash的作用:1、判重;2、统计次数;

零基础学习Redis(10) -- zset类型命令使用

zset是有序集合,内部除了存储元素外,还会存储一个score,存储在zset中的元素会按照score的大小升序排列,不同元素的score可以重复,score相同的元素会按照元素的字典序排列。 1. zset常用命令 1.1 zadd  zadd key [NX | XX] [GT | LT]   [CH] [INCR] score member [score member ...]

hdu 2602 and poj 3624(01背包)

01背包的模板题。 hdu2602代码: #include<stdio.h>#include<string.h>const int MaxN = 1001;int max(int a, int b){return a > b ? a : b;}int w[MaxN];int v[MaxN];int dp[MaxN];int main(){int T;int N, V;s

【机器学习】高斯过程的基本概念和应用领域以及在python中的实例

引言 高斯过程(Gaussian Process,简称GP)是一种概率模型,用于描述一组随机变量的联合概率分布,其中任何一个有限维度的子集都具有高斯分布 文章目录 引言一、高斯过程1.1 基本定义1.1.1 随机过程1.1.2 高斯分布 1.2 高斯过程的特性1.2.1 联合高斯性1.2.2 均值函数1.2.3 协方差函数(或核函数) 1.3 核函数1.4 高斯过程回归(Gauss