滑模控制(Sliding mode control)快速入门

news2024/11/24 16:38:21

0. 简介

最近作者受到邀请,让我帮忙给刚入门的学弟讲讲滑模控制。可是作者也不知道怎么向未入门的学弟讲解这些基础知识,所以作者翻了翻近几年写的很好的文章以及视频。综合起来,来总结出一套比较基础,且适用于初学者的文章吧。这里我们先贴一下王崇卫同学的笔记。
在这里插入图片描述
对应的视频连接在下面:

【Advanced控制理论】17

1. 滑模控制目的

对于滑模控制而言,我觉得我们先要明白其目的再来学习。一开始我们对滑动控制的定义是:滑动模式是先使用受控系统产生两个以上的子系统,然后再刻意加入一些切换条件产生滑动模式,以达成控制目标的一种技术。

滑模控制(sliding mode control, SMC)也叫变结构控制,其本质上是一类特殊的非线性控制,且非线性表现为控制的不连续性。这种控制策略与其他控制的不同之处在于系统的“结构”并不固定,而是可以在动态过程中,根据系统当前的状态(如偏差及其各阶导数等)有目的地不断变化,迫使系统按照预定“滑动模态”的状态轨迹运动。

例如滑动模式控制中存在滑动曲面 s = 0 s=0 s=0,一开始时,系统会在有限时间内到达滑动曲面,之后就会沿着滑动曲面移动。在滑动模式的理论叙述中,系统会约束在滑动曲面上,因此只需将系统视为在滑动曲面上滑动。不过实际系统的实现是用高频切换来让系统近似在滑动曲面上滑动,高频切换的控制信号让系统在很邻近滑动曲面的范围内切跳(chatter),而且其频率是不固定的。虽然整体系统是非线性的,不过下图中,当系统到达滑动曲面后,理想(没有切跳)系统会限制在 s = 0 s=0 s=0的滑动曲面上,滑动曲面是线性时不变系统,在原点处指数稳定。
在这里插入图片描述

2. 滑模控制优缺点

2.1 滑模控制的优点:

滑动模态可以设计且与对象参数和扰动无关,具有快速响应、对参数变化和扰动不灵敏( 鲁棒性)、无须系统在线辨识、物理实现简单

2.2 滑模控制的缺点:

当状态轨迹到达滑动模态面后,难以严格沿着滑动模态面向平衡点滑动,而是在其两侧来回穿越地趋近平衡点,从而产生抖振——滑模控制实际应用中的主要障碍。国内外主要通过改进滑模趋近律达到减弱抖振的目的。

3. 滑模控制需要条件

上文讲到滑模变结构控制器设计也包括两部分,一是能从状态空间的任何位置有限时间到达滑模面 s = 0 s = 0 s=0,二是在滑模面上可以收敛到原点(平衡点)。这也就代表我们要存在有一个稳定的滑模面,且该滑模面是可达的。为此有以下四个条件:

  • 稳定性条件:在s=0的滑模面上,状态是收敛的,即滑动模态存在;
  • 可达性条件:在切换面s=0以外的运动点将于有限时间内到达切换面;
  • 保证滑模运动的稳定性;
  • 达到控制系统运动品质要求。

下面将按照四个条件来叙述如何设计滑模控制的控制器,这里的部分内容借鉴了文章滑动模型控制(Sliding Mode Control),并结合作者的理解进行写作。

3.1 被控系统的滑模面生成

首先第一步就是我们需要明白,我们需要找到一个滑模面来让被控系统在滑模面上维持稳定。
例如假设存在一个被控系统:
x ˙ 1 = x 2 x ˙ 2 = u \begin{aligned} \dot{x}_1 &= x_2 \\ \dot{x}_2 &= u \end{aligned} x˙1x˙2=x2=u

这个时候我们就需要根据被控系统设计一个滑模面,滑模面一般可以设计为如下的形式
s ( x ) = ∑ i = 1 n − 1 c i x i + x n s(x) = \sum_{i=1}^{n-1} c_i x_i + x_n s(x)=i=1n1cixi+xn
因为在滑模控制中,要保证多项式 p n − 1 + c n p n − 2 + ⋯ + c 2 p + c 1 p^{n − 1} + c_n p^{n − 2} + \cdots + c_2 p + c_1 pn1+cnpn2++c2p+c1为Hurwitz (简单来说这条条件是为了满足状态在 s = 0 s=0 s=0的滑模面上可以收敛)。

什么是Hurwitz,即上述多项式的特征值的实数部分在左半平面,即为负。

我们可以看到上述的被控系统是存在有两个变量,所以需要取 n = 2 n=2 n=2,即 s ( x ) = c 1 x 1 + x 2 s ( x ) = c_1 x_1 + x_2 s(x)=c1x1+x2,为了保证多项式 p + c 1 p+c_1 p+c1为Hurwitz,需要多项式 p + c 1 = 0 p+c_1=0 p+c1=0的特征值实数部分为负,即 c 1 > 0 c_1>0 c1>0

我们知道滑模控制需要使得状态 x 1 x_1 x1 x 2 x_2 x2的导数均达到零,我们令 s = 0 s=0 s=0,分析一下结果有
{ c x 1 + x 2 = 0 x ˙ 1 = x 2    ⇒    c x 1 + x ˙ 1 = 0    ⇒    { x 1 ( t ) = e − c t x 1 ( 0 ) x 2 ( t ) = x ˙ 1 ( t ) = − c x 1 ( 0 ) e − c t \left\{\begin{aligned} &cx_1 + x_2 = 0 \\ &\dot{x}_1 = x_2 \end{aligned}\right. ~~ \Rightarrow ~~ c x_1 + \dot{x}_1 = 0 ~~ \Rightarrow ~~ \left\{\begin{aligned} &x_1(t) = \text{e}^{-ct} x_1(0) \\ &x_2(t) = \dot{x}_1(t) = -c x_1(0) \text{e}^{-ct} \end{aligned}\right. {cx1+x2=0x˙1=x2    cx1+x˙1=0    {x1(t)=ectx1(0)x2(t)=x˙1(t)=cx1(0)ect
通过上式可以看到状态 x 1 x_1 x1 x 2 x_2 x2 最终都是趋向于零的,而且速度是以指数速率趋紧的。指数速率意味着当 t = 1 / c t=1/c t=1/c时,趋零过程完成 63.2 % 63.2\% 63.2%,当 t = 3 / c t=3/c t=3/c时,趋零过程完成 95.021 % 95.021\% 95.021%。那么我们通过调节参数 c c c的大小即可实现对趋零速度的调节, c c c 越大,速度越快。

因此如果满足了 s = c x 1 + x 2 = 0 s=cx_1 + x_2=0 s=cx1+x2=0,那么系统的状态 x 1 x_1 x1 x 2 x_2 x2也将沿着滑模面趋近于零 ( s = 0 s=0 s=0称之为滑模面)。

3.2 可达性控制器设计

在拿到滑模面后则证明被控系统的稳定性条件成立,下面一步就是可达性条件,即状态 x x x 从状态空间中任意一点出发,可以在有限时间到达 s = 0 s=0 s=0 的滑模面上,此时我们可以采用李雅普诺夫间接法来分析,从前面可知,切换函数 s s s 是状态变量 x x x 的函数,取以下的李雅普诺夫函数

V = 1 2 s 2 V = \frac{1}{2} s^2 V=21s2

对时间求导可得
V ˙ = s s ˙ = s ( − sgn ( s ) − s ) = − sgn ( s ) s − s 2 = − ( ∣ s ∣ + s 2 ) < 0 \begin{aligned} \dot{V} &= s \dot{s} \\ &= s (-\text{sgn}(s) - s) \\ &= -\text{sgn}(s) s - s^2 \\ &= -(|s| + s^2) < 0 \end{aligned} V˙=ss˙=s(sgn(s)s)=sgn(s)ss2=(s+s2)<0

为了使系统稳定,我们需要使 V ˙ < 0 \dot{V}<0 V˙<0,即 s s ˙ < 0 s \dot{s}<0 ss˙<0。此时系统对于 s s s而言是渐进稳定,不能保证其有限时间到 s = 0 s=0 s=0 的滑模面上(渐进稳定是当 t t t趋于无穷时,状态变量 x x x 趋于 0 0 0,即无限时间到达),因此需要 s s ˙ < − σ s \dot{s}<-\sigma ss˙<σ σ \sigma σ是一个极小的正数。 以上就是可达性条件成立的必要依据 \color{red}{以上就是可达性条件成立的必要依据} 以上就是可达性条件成立的必要依据

但是实际上每次设计总不能都用李雅普诺夫函数判断,于是人们就提出了趋近律这一概念,常用的趋近律有如下几种,其中 sgn ( s ) \text{sgn}(s) sgn(s) 是符号函数, s > 0 , sgn ( s ) = 1 ; s < 0 , sgn ( s ) = − 1 ; s = 0 , sgn ( s ) = 0 s>0,\text{sgn}(s)=1; s<0, \text{sgn}(s)=-1; s=0, \text{sgn}(s)=0 s>0,sgn(s)=1;s<0,sgn(s)=1;s=0,sgn(s)=0:

  1. 等速趋近律: s ˙ = − ϵ  sgn ( s ) ,      ϵ > 0 \dot{s} = -\epsilon ~\text{sgn}(s), ~~~~\epsilon > 0 s˙=ϵ sgn(s),    ϵ>0

  2. 指数趋近律: s ˙ = − ϵ  sgn ( s ) − k s ,      ϵ > 0 , k > 0 \dot{s} = -\epsilon ~\text{sgn}(s) - k s, ~~~~\epsilon > 0, k>0 s˙=ϵ sgn(s)ks,    ϵ>0,k>0

  3. 幂次趋近律: s ˙ = − k ∣ s ∣ α  sgn ( s ) − k s ,      k > 0 , 0 < α < 1 \dot{s} = -k |s|^\alpha ~\text{sgn}(s) - k s, ~~~~k>0, 0<\alpha<1 s˙=ksα sgn(s)ks,    k>0,0<α<1

一般在使用时候我们需要完成这些参数的调整,一般我们使用的是指数趋近率,并将 ϵ \epsilon ϵ k k k的值均设为1,简化为:

s ˙ =  sgn ( s ) − s \dot{s} = ~\text{sgn}(s) - s s˙= sgn(s)s

然后我们可知 s ( x ) = c 1 x 1 + x 2 s ( x ) = c_1 x_1 + x_2 s(x)=c1x1+x2,则 s ˙ =  sgn ( s ) − s = c 1 x 1 ˙ + x 2 ˙ = c 1 x 2 + u \dot{s} = ~\text{sgn}(s) - s = c_1 \dot{x_1} + \dot{x_2} = c_1x_2+u s˙= sgn(s)s=c1x1˙+x2˙=c1x2+u。则我们可以得到控制器 u u u为:

u =  sgn ( s ) − s − c 1 x 2 u = ~\text{sgn}(s) - s - c_1x_2 u= sgn(s)sc1x2

这就得到了我们必要的两个条件即,存在滑模面 s s s以及可达性控制器 u u u.

4. 滑模控制Python代码

下面是最简单的python代码

…详情请参照古月居

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/392064.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

Node.js的简单学习一-----未完待续

文章目录前言学习目标一、初识Node.js1.1 回顾与思考1.1.1 需要掌握那些技术1.1.2 浏览器中的JavaScript的组成部分1.2 Node.js简介1 什么是Node.js2 Node.js中的JavaScript运行环境3 Node.js 可以做什么1.3 Node.js环境的安装1.4 在Node.js环境中执行JavaScript 代码终端中的快…

掌握好Framework 才是王道~

现在面试对Android开发者的要求越来越高了&#xff01;从最开始的阿里、头条、腾讯等大厂&#xff0c;到现在的互联网车企&#xff0c;面试总喜欢问道 Framework底层原理的相关问题 Android Framework的三大核心功能&#xff1a; 1、View.java:View工作原理&#xff0c;实现包…

海尔三翼鸟:生态聚拢的密度,决定场景落地速度

最近学到一个新词&#xff0c;叫做涌现能力。 怎么理解呢&#xff1f;我们以当下最火的ChatGPT为例&#xff0c;GPT1模型是1.17亿参数&#xff0c;GPT2有15亿参数&#xff0c;GPT3有1750亿个参数。研究人员在放大模型规模的进程中发现一个惊人的现象&#xff0c;模型参数达到一…

Vue基础学习 -- Vue数据与方法 v-指令补充 vue生命周期

Vue数据与方法 实例化Vue对象时&#xff0c;需要在data中定义数据&#xff0c;data选项是一个函数&#xff0c;它应该返回一个对象&#xff0c;然后Vue通过响应性系统将其包裹起来&#xff0c;并以$data形式储存在实例中。为了方便起见&#xff0c;该归降的任何顶级property可…

【蓝桥杯刷题】盗版Huybery系列之手抓饼赛马

【蓝桥杯刷题】—— 盗版Huybery系列之手抓饼赛马&#x1f60e;&#x1f60e;&#x1f60e; 目录 &#x1f4a1;前言&#x1f31e;&#xff1a; &#x1f49b;盗版Huybery系列之手抓饼赛马题目&#x1f49b; &#x1f4aa; 解题思路的分享&#x1f4aa; &#x1f60a;题…

计算机SCI期刊投稿,除了投稿信,还要做什么准备? - 易智编译EaseEditing

投稿信的准备 期刊的编辑往往需要一些有关作者及其论文的信息。 而作者也希望给编辑提供一些有助于其全文送审及决策的信息。 这些信息都应该包括在投稿信中。 投稿信应包括以下几方面的内容&#xff1a; 文题和所有作者的姓名;稿件适宜的栏目; 为什么此论文适合于在该刊而…

大数据框架之Hive:第2章 Hive安装

第2章 Hive安装 2.1 Hive安装地址 1&#xff09;Hive官网地址 http://hive.apache.org/ 2&#xff09;文档查看地址 https://cwiki.apache.org/confluence/display/Hive/GettingStarted 3&#xff09;下载地址 http://archive.apache.org/dist/hive/ 4&#xff09;gith…

OpenCV入门(四)快速学会OpenCV3画基本图形

OpenCV入门&#xff08;四&#xff09;快速学会OpenCV3画基本图形 1.画点 在OpenCV中&#xff0c;点分为2D平面中的点和3D平面中的点&#xff0c;区别就是3D中点多了一个z坐标。我们首先介绍2D中的点&#xff0c;坐标为整数的点可以直接用(x, y)代替&#xff0c;其中x是横坐标…

快速幂----快速求解底数的n次幂

目录 一.快速幂 1.问题的引入 2.快速幂的介绍 3.核心思想 4.代码实现 2.猴子碰撞的方法数 1.题目描述 2.问题分析 3.代码实现 一.快速幂 1.问题的引入 问题:求解num的n次幂,结果需要求余7 对于这个问题我们可能就是直接调用函数pow(a,b)来直接求解a的b次幂问题,但是如果…

灵巧手抓持<分类><仿真>

获取灵巧手抓取物体时的抓持类型&#xff0c;需要考虑&#xff1a;手本身的结构、被抓取物体的形状尺寸、抓持操作任务的条件。 研究方法&#xff1a;基于模型的方法、基于数据驱动的方法 基于模型的方法&#xff1a;建立灵巧手抓持相关的运动学和动力学模型建立目标函数求解…

【操作系统】进程和线程的区别

文章目录1. 概述2. 进程3. 线程4. 协程5. 进程与线程区别1. 概述 进程和线程这两个名词天天听&#xff0c;但是对于它们的含义和关系其实还有点懵的&#xff0c;其实除了进程和线程&#xff0c;还存在一个协程&#xff0c;它们的关系如下&#xff1a; 首先&#xff0c;我们需要…

C++基础——Ubuntu下编写C++环境配置总结(C++基本简介、Ubuntu环境配置、编写简单C++例程)

【系列专栏】&#xff1a;博主结合工作实践输出的&#xff0c;解决实际问题的专栏&#xff0c;朋友们看过来&#xff01; 《QT开发实战》 《嵌入式通用开发实战》 《从0到1学习嵌入式Linux开发》 《Android开发实战》 《实用硬件方案设计》 长期持续带来更多案例与技术文章分享…

优思学院|质量管理简史:从质检到质控,到六西格玛管理

在现代的工业化生产中&#xff0c;质量管理已成为企业生产运营的不可或缺的一部分。经过长时间的探索和实践&#xff0c;质量管理的理论和实践已经逐步成熟&#xff0c;并形成了完整的体系。下面&#xff0c;我将简要介绍一下质量管理学的发展历程。 质量管理的萌芽期&#xff…

socket聊天室--socket的建立

socket聊天室–socket实现 文章目录 socket聊天室--socket实现socket()bind()listen()accept()connect()发送接收read()函数recv()函数write()函数send()函数close()关闭套接字IP 地址格式转换函数socket() #include <sys/types

web餐饮开源程序

简介 一款专门针对餐饮行业而开发桌面应用程序 技术 借助Panuon.UI.Silver控件库&#xff0c;开发的一款餐饮软件。 运行环境&#xff1a;.NETFramework,Versionv4.8。 运行数据库&#xff1a;MySql。 ORM框架&#xff1a;SqlSugar。 第三方插件&#xff1a;Panuon.UI.Silv…

OpenKruise 开发者不容错过的带薪实习机会!马上加入 LFX Mentorship 计划

LFX Mentorship 计划由 Linux Foundation 组织发起&#xff0c;为像 OpenKruise 这样的 CNCF 托管项目提供了激励开源贡献、扶植社区发展的优秀土壤。参与其中的开发者不仅有机会在经验丰富的社区 Mentor 指导下贡献开源项目、为职业生涯加分&#xff0c;完成工作后还能获得 $3…

【JAVA】Java中集合的遍历

&#x1f3c6;今日学习目标&#xff1a;集合遍历 &#x1f603;创作者&#xff1a;颜颜yan_ ✨个人主页&#xff1a;颜颜yan_的个人主页 ⏰本期期数&#xff1a;第五期 &#x1f389;专栏系列&#xff1a;JAVA 集合遍历前言一、Iterator接口二、foreach循环总结前言 在实际开发…

JavaScript基础二、变量

零、文章目录 文章地址 个人博客-CSDN地址&#xff1a;https://blog.csdn.net/liyou123456789个人博客-GiteePages&#xff1a;https://bluecusliyou.gitee.io/techlearn 代码仓库地址 Gitee&#xff1a;https://gitee.com/bluecusliyou/TechLearnGithub&#xff1a;https:…

UE4 材质学习 (焚烧材质)

效果步骤随便从网上下载一张图片&#xff08;地址&#xff1a;链接/链接&#xff09;&#xff0c;导入UE中新建一个材质函数这里命名为“E_Function”双击打开该材质函数&#xff0c;由于需要输出变发光和变透明两种效果&#xff0c;因此这里需要两个输出节点&#xff1a;分别命…

MIPI协议

MIPI调试指南Rev.0.1 June 18, 2019 © 2018 Horizon Robotics. All rights reserved.Revision HistoryThissection tracks the significant documentation changes that occur fromrelease-to-release. The following table lists the technical content changes foreach …