滑模控制算法——基本原理（附MATLAB程序）

滑模控制算法（Sliding Mode Control, SMC）是一种非线性控制策略，广泛用于处理系统的非线性、扰动以及建模不确定性。它的核心思想是通过设计一个滑模面，使得系统状态在滑模面上滑动，从而实现系统的稳定性和期望的动态性能。

一、基本概念

滑模面（Sliding Surface）：选择一个合适的滑模面，使得系统在该面上运动时能实现期望的动态特性。滑模面通常定义为系统状态变量的线性组合，例如 s=Cx+D，其中 C 和 D 是设计参数。

滑模控制律（Sliding Mode Control Law）：在滑模面上，控制律的设计目的是使得系统状态在滑模面上滑动，并且保持在该面上。常见的控制律形式包括：u=−ksign(s)，其中 k是控制增益，sign(s) 是滑模面 s的符号函数。

二、公式推导

滑模控制算法的推导涉及几个关键步骤：定义滑模面、设计滑模控制律、以及验证系统稳定性。以下是滑模控制算法的主要推导过程及公式。

2.1 系统模型

考虑一个动态系统，其状态空间方程为：

$\dot{x}=Ax+Bu$

其中：

是系统状态向量。

是控制输入。

和

是系统矩阵。

2.2 定义滑模面

选择一个滑模面 $s\left ( x \right )$ ：

$s\left ( x \right )=Cx+D$

其中：

和

是设计参数。 $s\left ( x \right )$ 是滑模面函数，通常我们希望系统状态在 $s\left ( x \right )=0$ 的面上滑动。

2.3 滑模控制律设计

控制律的设计目的是使得系统状态在滑模面上滑动。选择控制律如下：

$u=u_{eq}-Ksign\left ( s\left ( x \right ) \right )$

其中：

$u_{eq}$ 是等效控制项。

是控制增益。 $sign\left ( s\left ( x \right ) \right )$ 是符号函数，它为滑模面 $s\left ( x \right )$ 提供正负方向的信息。

2.4 等效控制

为了使系统状态在滑模面上稳定滑动，需要计算等效控制 $u_{eq}$ 。在滑模面上，系统应满足 $\dot{s}=0$ 。计算 $\dot{s}$ ：

$\dot{s}=C\dot{x}$

由于 $\dot{x}=Ax+Bu$ ，代入得到

$\dot{s}=C\left ( Ax+Bu \right )=CAx+CBu$

为了使系统状态在滑模面上滑动，即 $\dot{s}=0$ ，我们需要：

CAx+CBu=0

所以，等效控制 $u_{eq}$ 为：

$u_{eq}=-\left ( CB \right )^{-1}CAx$

2.5 滑模面稳定性分析

系统的控制目标是使得状态在滑模面 $s\left ( x \right )=0$ 上稳定。考虑控制律：

$u=u_{eq}-Ksign\left ( s\left ( x \right ) \right )$

代入系统方程：

$\dot{x}=Ax+B\left [ u_{eq}-Ksign\left ( s\left ( x \right ) \right ) \right ]$

将等效控制 $u_{eq}$ 代入，得到：

$\dot{x}=Ax+B\left [ -\left ( CB \right )^{-1}CAx-Ksign\left ( s\left ( x \right ) \right ) \right ]$

因此，系统的状态方程变为：

$\dot{x}=\left ( A-BK\left ( CB \right )^{-1}CA \right )x-BKsign\left ( s\left ( x \right ) \right )$

在滑模面上，状态的动态行为由以下方程描述：

$\dot{s}=C\dot{x}=C\left ( A-BK\left ( CB \right ) ^{-1}CA\right )x-CBKsign\left ( s\left ( x \right ) \right )$

简化为：

$\dot{s}=-CBKsign\left ( s\left ( x \right ) \right )$

在滑模面上，为了确保系统的稳定性，选择合适的使得滑模面 $s\left ( x \right )=0$ 的稳定性条件满足。通常，我们要求足够大，使得滑模面上的运动是稳定的，并且状态变量能够快速趋向于滑模面。

三、优势及挑战

3.1优势

鲁棒性：滑模控制对系统的不确定性和外部扰动具有较强的鲁棒性。即使系统模型存在误差，滑模控制依然能保证系统的稳定性。

适应性：滑模控制可以处理系统的非线性特性，通过设计适当的滑模面和控制律，能够实现对多种非线性动态的控制。

3.2 挑战

抖振现象：滑模控制可能会引起系统的抖振现象（chattering），这主要是由于控制律的高频切换引起的。为了解决这个问题，常常需要使用平滑的滑模控制方法，如高阶滑模控制或自适应滑模控制。

设计复杂性：对于复杂系统，滑模控制的设计和实现可能比较复杂，需要进行详细的数学分析和调试。

3.3 方法改进

由于滑模控制中的符号函数可能导致高频切换（抖振现象），可以使用平滑的替代方法，例如：

平滑控制：使用连续函数替代符号函数，如饱和函数或符号函数的平滑近似。

高阶滑模控制：设计高阶滑模面以减小抖振。

四、MATLAB仿真程序

下面是一个基本的 MATLAB 程序示例，演示如何使用滑模控制算法对一个简单的线性系统进行仿真。这个例子将以一个二阶系统为基础，并实现一个简单的滑模控制器。程序中包含了系统建模、滑模控制器设计、仿真及结果展示的步骤。

4.1 定义系统参数

我们以一个简单的二阶系统作为示例，其状态空间方程为：

$\dot{x}=Ax+Bu$

y=Cx

其中：

是系统矩阵

是输入矩阵

是输出矩阵

假设我们的系统如下：

$A=\begin{bmatrix} 0 &1 \\ -2& -2 \end{bmatrix}$ $B=\begin{bmatrix} 0\\ 1 \end{bmatrix}$ $C=\begin{bmatrix} 1 & 0 \end{bmatrix}$

4.2 编写MATLAB程序

% 清除工作区和命令窗口clear; clc;% 系统参数A = [0 1; -2 -2]; % 系统矩阵 AB = [0; 1];       % 输入矩阵 BC = [1 0];        % 输出矩阵 CD = 0;            % 系统 D 矩阵（此处为零）% 滑模控制设计lambda = 2;       % 选择滑模面参数k = 10;           % 控制增益% 初始化x0 = [1; 0];      % 初始状态tspan = [0 10];   % 仿真时间% 定义状态空间方程及滑模控制律% 采用匿名函数system = @(t, x) A*x + B*(-k*sign(C*x)); % 运行仿真[t, X] = ode45(system, tspan, x0);% 计算系统输出Y = C * X';% 绘制结果figure;subplot(2,1,1);plot(t, X(:,1), 'r', 'LineWidth', 1.5);hold on;plot(t, X(:,2), 'b', 'LineWidth', 1.5);xlabel('Time (s)');ylabel('States');legend('x1', 'x2');title('System States');subplot(2,1,2);plot(t, Y, 'k', 'LineWidth', 1.5);xlabel('Time (s)');ylabel('Output');title('System Output');% 滑模面 s(x) 和控制律 u(t)s = @(x) C*x;u = @(x) -k*sign(s(x));% 打印控制输入disp('Control input u(t):');for i = 1:length(t)    fprintf('Time: %.2f s, Control Input: %.2f\n', t(i), u(X(i,:)'));end

代码解释：

系统参数：定义系统的状态矩阵、输入矩阵和输出矩阵。这里选择了一个简单的二阶系统。