西安电子科技大学，陕西 西安

收稿日期：2022年1月1日；录用日期：2022年1月14日；发布日期：2022年1月28日

摘要

本文将针对分数阶( $\alpha >1$ )系统，同时考虑最常见的两种不确定类型，即多胞不确定和无穷范数。首先，给出了不确定分数阶系统无状态反馈时鲁棒稳定性条件。其次，采用线性矩阵不等式的形式给出了不确定分数阶系统状态反馈控制器的设计方法。最后，通过两个数值算例验证了理论的正确性。

关键词

分数阶系统，多胞不确定，范数有界不确定，鲁棒稳定性

Study on the Robust Stability of Fractional-Order System

Xiaozhong Guo

Xidian University, Xi’an Shaanxi

Received: Jan. 1^st, 2022; accepted: Jan. 14^th, 2022; published: Jan. 28^th, 2022

ABSTRACT

In this paper, we will consider the two most common types of uncertainty for fractional-order systems, which are subjected to poly-topic and H-infinity bounded uncertainties with $\alpha >1$. Firstly, the robust stability conditions of uncertain fractional-order systems without state feedback are given. Secondly, the design method of state feedback controller for the fractional system with uncertainty is given in the form of linear matrix inequality. Finally, two numerical examples are given to verify the correctness of the theory.

Keywords:Fractional-Order System, Poly-Topic Uncertainty, Norm Bounded Uncertainty, Robust Stability

Copyright © 2022 by author(s) and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

1. 引言

不确定性研究有多种形式。在 [1] [2] 中，分数阶微积分作为整数阶微积分的推广和延伸，由于可以更好地刻画函数的历史发展过程，有更强的记忆性和全局相关性，得到了国内外学者的广泛关注，利用多胞和范数有界描述不确定性，给出了一个基于线性矩阵不等式稳定性分析和控制器设计的充分条件。在 [3] 中，考虑了区间不确定性下的分数阶系统的鲁棒稳定性，分数阶与其他模型参数之间存在确定的线性耦合关系。在 [4] 中，分数阶的研究在控制领域中主要集中于物理对象的建模。 [5] [6] 中，基于线性矩阵不等式分别对 $1<\alpha <2$ 和 $0<\alpha <1$ 的情况，给出了具有区间不确定性的分数阶系统稳定的充要条件。在 [7] 中，同时考虑了多胞和二范数有界不确定性。在 [8] 中，研究了在多胞和范数有界不确定条件下的分数阶系统的稳定性和镇定问题。在 [9] 中，研究了不确定分数阶的鲁棒稳定界问题，并且系统具有范数有界不确定性。文献 [10] [11] [12] [13] 研究不确定性矩阵谱范数的上界，考虑的稳定区域要么是复平面的左半平面，要么是单一的关于实轴对称的凸域。这些上界都是基于Lyapunov定理。文献 [14] 考虑了当分数阶系统具有多胞型不确定时的鲁棒镇定控制器的设计方法。

在 [15] 中，作者从控制输入能效的角度研究了分数阶控制系统，给出了控制输入能量函数的发散条件，对作用于对象上的控制器能量优化起到了作用。

基于上述研究，本文将研究两种不确定性与稳定性之间的关系，推导出了分数阶( $\alpha >1$ )系统在多胞和范数有界不确定情况下，系统的稳定性条件。设计状态反馈控制器使分数阶系统变得稳定，并通过数值算例，验证理论的正确性。

2. 预备知识

在本文的后续中，如果没有特别强调或指出，一般符号的含义如下所示。

$X^{\text{T}}$ 表示矩阵X的转置， $X^{\text{*}}$ 表示矩阵X的共轭转置， $X>0$ 表示X是正定矩阵。

$sym\left\{X\right\}$ 等于 $X+X^{*}$，矩阵中的 $•$ 表示这是对称矩阵中的对称部分。 $\otimes$ 表示矩阵的Kronecker积。 $D^{\alpha }$ 是 $\alpha$ 阶Caputo微分算子。

${\Vert X\Vert }_{\infty }$ 是X的无穷大范数， $I_n$ 是n阶单位矩阵。

引理1. 当 $1<\alpha <2$ 时，系统稳定的充分必要条件是存在正定矩阵 $X=X^{\text{T}}>0$，$X\in R^{n\times n}$ 使得如下线性矩阵不等式成立：

$sym\left\{\Lambda \otimes \left(AX\right)\right\}<0$ (1)

其中 $\otimes$ 表示矩阵的Kronecker积：

$\Lambda =\left[\begin{array}{cc}\sin \theta & -\cos \theta \\ \cos \theta & \sin \theta \end{array}\right]$，$\theta =\pi -\alpha \frac{\pi }{2}$ (2)

引理2. 对于具有适当维数的正定矩阵 $\chi >0$ 和矩阵 $Q,P$，有以下不等式成立：

$Q^{\text{T}}P+P^{\text{T}}Q\le Q^{\text{T}}\chi Q+P^{\text{T}}{\chi }^{-1}P$ (3)

引理3. 对于给定的对称矩阵：

$\Lambda =\left[\begin{array}{cc}A& B\\ B^{\text{T}}& C\end{array}\right]$ (4)

$A,B,C$ 都具有恰当的维数，则以下三个条件是等价的：

1. $\Lambda <0$ ；

2. $A<0$，$C-B^{\text{T}}{A}^{-1}B<0$ ；

3. $C<0$，$A-B{C}^{-1}{B}^{\text{T}}<0$ ；

引理4. 当 $0<\alpha <2$ 时，分数阶系统渐近稳定的充分必要条件是：

$\left|Arg\left(eig\left(A\right)\right)\right|>\alpha \pi /2$ (5)

其中 $Arg\left(\right)$ 是辐角函数， $eig\left(A\right)$ 是矩阵A的所有特征值。

3. 结果与分析

3.1. 鲁棒稳定性分析

定理3.1 当 $u\left(t\right)=0$ 时，多胞和范数有界不确定分数阶系统具有鲁棒渐近稳定性的充分条件是：存在矩阵 $X_j=X_j^{\text{T}}>0$，$X_j\in R^{n\times n}$，$X_{ij}\in Z^{+}$，$\left(i,j=1,2,\cdots ,N\right)$ 满足约束：

$\left[\begin{array}{cc}sym\left(\Lambda \otimes {\hat{A}}_i{X}_j\right)+{\displaystyle {\sum }_{k_1=1}^m\left\{{\chi }_{ii{k}_1}\left[I_2\otimes \left(F_{i{k}_1}{F}_{i{k}_1}^{\text{T}}\right)\right]\right\}}& .\\ I_2\otimes X_j& -\frac{{\left({\displaystyle {\sum }_{k_1=1}^m{\chi }_{ii{k}_1}^{-1}}\right)}^{-1}}{{\delta }_{1\max }^2}\left(I_2\otimes I_n\right)\end{array}\right]<0$ (6)

证明. 根据引理1，如果存在正定矩阵 $X_j$ 及 $X\left(\gamma \right)={\displaystyle {\sum }_{j=1}^N{\gamma }_j{X}_j}$，使得分数阶系统渐近稳定，那么则有：

$sym\left\{\Lambda \otimes \left(AX\right)\right\}<0$

$\iff {\displaystyle {\sum }_{i=1}^N{\displaystyle {\sum }_{j=1}^N{\gamma }_i{\gamma }_{j}sym\left(\Lambda \otimes \left(A_i{X}_j\right)\right)}}<0$

$\begin{array}{l}\iff {\displaystyle {\sum }_{i=1}^N{\gamma }_i^2}sym\left\{\Lambda \otimes \left[{\hat{A}}_i+{\displaystyle {\sum }_{k_1=1}^m{\delta }_{i{k}_1}{F}_{i{k}_1}}\right]X_i\right\}\\ +{\displaystyle {\sum }_{i=1}^{N-1}{\displaystyle {\sum }_{j=i+1}^N{\gamma }_i{\gamma }_{j}sym\left\{\Lambda \otimes \left[\left({\hat{A}}_i+{\displaystyle {\sum }_{k_1}^m{\delta }_{i{k}_1}{F}_{i{k}_1}}\right)X_j\right]+\left({\hat{A}}_j+{\displaystyle {\sum }_{k_1}^m{\delta }_{j{k}_1}{F}_{j{k}_1}}\right)X_i\right\}}}<0\end{array}$ (7)

根据引理2，存在矩阵 ${\chi }_{ii{k}_1}>0\left(i=1,2,\cdots ,N\right)$，

${\chi }_{ij{k}_1}>0$，${\chi }_{ji{k}_1}>0$，$\left(i=1,2,\cdots ,N-1\right)$，$\left(j=i+1,\cdots ,N\right)$，

且有 $\Vert {\delta }_{i{k}_1}\Vert \le {\delta }_{1\max }\left(i=1,2,\cdots ,N;k_1=1,2,\cdots ,m\right)$，可以得到：

$\begin{array}{l}{\displaystyle {\sum }_{i=1}^N{\gamma }_i^2}\left\{sym\left(\Lambda \otimes {\hat{A}}_i{X}_i\right)+{\displaystyle {\sum }_{k_1=1}^m\left\{{\chi }_{ii{k}_1}\left[I_2\otimes \left(F_{i{k}_1}{F}_{i{k}_1}^{\text{T}}\right)\right]+{\chi }_{ii{k}_1}^{-1}{\delta }_{1\max }^2\left[I_2\otimes X_i^{\text{T}}{X}_i\right]\right\}}\right\}\\ +{\displaystyle {\sum }_{i=1}^N{\displaystyle {\sum }_{j=i+1}^N{\gamma }_i}{\gamma }_j}\{sym\left(\Lambda \otimes {\hat{A}}_i{X}_j\right)+{\displaystyle {\sum }_{k_1=1}^m\left\{{\chi }_{ij{k}_1}\left[I_2\otimes \left(F_{i{k}_1}{F}_{i{k}_1}^{\text{T}}\right)\right]+{\chi }_{ij{k}_1}^{-1}{\delta }_{1\max }^2\left[I_2\otimes X_j^{\text{T}}{X}_j\right]\right\}}\\ +sym\left(\Lambda \otimes {\hat{A}}_j{X}_i\right)+{\displaystyle {\sum }_{k_1=1}^m\left\{{\chi }_{ji{k}_1}\left[I_2\otimes \left(F_{j{k}_1}{F}_{j{k}_1}^{\text{T}}\right)\right]+{\chi }_{ji{k}_1}^{-1}{\delta }_{1\max }^2\left[I_2\otimes X_i^{\text{T}}{X}_i\right]\right\}}\}<0\end{array}$ (8)

使用引理3，可以将式(8)写为：

$\left[\begin{array}{cc}sym\left(\Lambda \otimes {\hat{A}}_i{X}_j\right)+{\displaystyle {\sum }_{k_1=1}^m\left\{{\chi }_{ii{k}_1}\left[I_2\otimes \left(F_{i{k}_1}{F}_{i{k}_1}^{\text{T}}\right)\right]\right\}}& .\\ I_2\otimes X_j& -\frac{{\left({\displaystyle {\sum }_{k_1=1}^m{\chi }_{ii{k}_1}^{-1}}\right)}^{-1}}{{\delta }_{1\max }^2}\left(I_2\otimes I_n\right)\end{array}\right]<0$，$i,j=1,2,\cdots ,N$ (9)

根据引理1，系统的鲁棒稳定性得证。

3.2. 状态反馈控制器设计

定理3.2 当 $u\left(t\right)=Kx\left(t\right)$ 时，多胞和范数有界不确定分数阶系统具有鲁棒渐近稳定性的充分条件是：存在矩阵 $X_j=X_j^{\text{T}}>0$，$X_j\in R^{n\times n}$，$R_j\in R^{m\times n}$，${\chi }_{1ij},{\chi }_{2ij}\in Z^{+}$，$\left(i,j=1,2,\cdots ,N\right)$ 满足约束：

$\left[\begin{array}{ccc}\begin{array}{c}sym\left[\Lambda \otimes \left({\hat{A}}_i{X}_j+{\hat{B}}_i{R}_j\right)\right]\\ +{\displaystyle {\sum }_{k_1=1}^m{\chi }_{1ij{k}_1}\left[I_2\otimes \left(F_{i{k}_1}{F}_{i{k}_1}^{\text{T}}\right)\right]}\\ +{\displaystyle {\sum }_{k_2=m+1}^I{\chi }_{2ij{k}_2}\left[I_2\otimes \left(F_{i{k}_2}{F}_{i{k}_2}^T\right)\right]}\end{array}& .& .\\ I_2\otimes X_j& -\frac{{\left({\displaystyle {\sum }_{k_1=1}^m{\chi }_{1ij{k}_1}^{-1}}\right)}^{-1}}{{\delta }_{1\max }^2}\left(I_2\otimes I_n\right)& .\\ I_2\otimes R_j& 0& -\frac{{\left({\displaystyle {\sum }_{k_2=m\text{+}1}^I{\chi }_{2ij{k}_2}^{-1}}\right)}^{-1}}{{\delta }_{2\max }^2}\left(I_2\otimes I_n\right)\end{array}\right]<0$ (10)

式中状态反馈器可描述为：

$u\left(t\right)=Kx\left(t\right)$，$K\left(\gamma \right)={\displaystyle {\sum }_{i=1}^N{\gamma }_i{K}_i}={\displaystyle {\sum }_{i=1}^N{\gamma }_i{R}_i{\displaystyle {\sum }_{i=1}^N{X}_i^{-1}}}$

证明. 根据引理1，如果存在正定矩阵 $X_j$，$X_i$ 及 $X\left(\gamma \right)={\displaystyle {\sum }_{j=1}^N{\gamma }_j{X}_j}$，$K\left(\gamma \right)={\displaystyle {\sum }_{i=1}^N{\gamma }_i{K}_i}$ 使得分数阶系统渐近稳定，那么则有：

$sym\left\{\Lambda \otimes \left(A+BK\right)X\right\}<0$

$\begin{array}{l}\iff {\displaystyle {\sum }_{i=1}^N{\gamma }_i^2}sym\left\{\Lambda \otimes \left[\left({\hat{A}}_i+{\displaystyle {\sum }_{k_1=1}^m{\delta }_{i{k}_1}{F}_{i{k}_1}}\right)X_i\right]+\Lambda \otimes \left[\left({\hat{B}}_i+{\displaystyle {\sum }_{k_2=m+1}^I{\delta }_{i{k}_2}{F}_{i{k}_2}}\right)R_i\right]\right\}\\ +{\displaystyle {\sum }_{i=1}^{N-1}{\displaystyle {\sum }_{j=i+1}^N{\gamma }_i{\gamma }_{j}sym}}\{\Lambda \otimes \left[\left({\hat{A}}_i+{\displaystyle {\sum }_{k_1}^m{\delta }_{i{k}_1}{F}_{i{k}_1}}\right)X_j+\left({\hat{A}}_j+{\displaystyle {\sum }_{k_2=m+1}^I{\delta }_{j{k}_1}{F}_{j{k}_1}}\right)X_i\right]\\ +\Lambda \otimes \left[\left({\hat{B}}_i+{\displaystyle {\sum }_{k_2=m+1}^I{\delta }_{i{k}_2}{F}_{i{k}_2}}\right)R_j+\left({\hat{B}}_j+{\displaystyle {\sum }_{k_2=m+1}^I{\delta }_{j{k}_2}{F}_{j{k}_2}}\right)R_i\right]\}<0\end{array}$ (11)

根据引理2，存在矩阵 ${\chi }_{1ii{k}_1},{\chi }_{2ii{k}_2}>0$，${\chi }_{1ij{k}_1},{\chi }_{2ij{k}_2}>0$，${\chi }_{1ji{k}_1},{\chi }_{2ji{k}_2}>0$，$\Vert {\delta }_{i{k}_1}\Vert \le {\delta }_{1\max }$，$\left(i=1,2,\cdots ,N-1;j=i+1,\cdots ,N\right)$，且有 $\Vert {\delta }_{i{k}_2}\Vert \le {\delta }_{2\max }$，$\left(i=1,2,\cdots ,N\right)$，$\left(k_1=1,2,\cdots ,m\right)$ 可以推导出：

$\begin{array}{l}{\displaystyle {\sum }_{i=1}^N{\gamma }_i^2}\{sym\left[\Lambda \otimes \left({\hat{A}}_i{X}_i+{\hat{B}}_i{R}_i\right)\right]+{\displaystyle {\sum }_{k_1=1}^m\left\{{\chi }_{1ii{k}_1}\left[I_2\otimes \left(F_{i{k}_1}{F}_{i{k}_1}^{\text{T}}\right)\right]+{\chi }_{1ii{k}_1}^{-1}{\delta }_{1\max }^2\left[I_2\otimes X_i^{\text{T}}{X}_i\right]\right\}}\\ +{\displaystyle {\sum }_{k_2=m+1}^I\left\{{\chi }_{2ii{k}_2}\left[I_2\otimes \left(F_{i{k}_2}{F}_{i{k}_2}^{\text{T}}\right)\right]+{\chi }_{2ii{k}_2}^{-1}{\delta }_{2\max }^2\left[I_2\otimes R_i^{\text{T}}{R}_i\right]\right\}}\}\\ +{\displaystyle {\sum }_{i=1}^{N-1}{\displaystyle {\sum }_{j=i+1}^N{\gamma }_i{\gamma }_j\{sym\left[\Lambda \otimes \left({\hat{A}}_i{X}_j+{\hat{B}}_i{R}_j\right)\right]}}+{\displaystyle {\sum }_{k_1=1}^m\{{\chi }_{1ij{k}_1}\left[I_2\otimes \left(F_{i{k}_1}{F}_{i{k}_1}^{\text{T}}\right)\right]}\\ +{\chi }_{1ij{k}_1}^{-1}{\delta }_{1\max }^2\left[I_2\otimes X_j^{\text{T}}{X}_j\right]\}+{\displaystyle {\sum }_{k_2=m+1}^I\left\{{\chi }_{2ij{k}_2}\left[I_2\otimes \left(F_{i{k}_2}{F}_{i{k}_2}^{\text{T}}\right)\right]+{\chi }_{2ij{k}_2}^{-1}{\delta }_{2\max }^2\left[I_2\otimes R_j^{\text{T}}{R}_j\right]\right\}}\\ +sym\left[\Lambda \otimes \left({\hat{A}}_j{X}_i+{\hat{B}}_j{R}_i\right)\right]+{\displaystyle {\sum }_{k_1=1}^m\left\{{\chi }_{1ji{k}_1}\left[I_2\otimes \left(F_{j{k}_1}{F}_{j{k}_1}^{\text{T}}\right)\right]+{\chi }_{1ji{k}_1}^{-1}{\delta }_{1\max }^2\left[I_2\otimes X_i^{\text{T}}{X}_i\right]\right\}}\\ +{\displaystyle {\sum }_{k_2=m+1}^I\left\{{\chi }_{2ji{k}_2}\left[I_2\otimes \left(F_{j{k}_2}{F}_{j{k}_2}^{\text{T}}\right)\right]+{\chi }_{2ii{k}_2}^{-1}{\delta }_{2\max }^2\left[I_2\otimes R_i^{\text{T}}{R}_i\right]\right\}}\}<0\end{array}$ (12)

使用引理3，可以将式(12)写为：

$\left[\begin{array}{ccc}\begin{array}{c}sym\left[\Lambda \otimes \left({\hat{A}}_i{X}_j+{\hat{B}}_i{R}_j\right)\right]\\ +{\displaystyle {\sum }_{k_1=1}^m{\chi }_{1ij{k}_1}\left[I_2\otimes \left(F_{i{k}_1}{F}_{i{k}_1}^{\text{T}}\right)\right]}\\ +{\displaystyle {\sum }_{k_2=m+1}^I{\chi }_{2ij{k}_2}\left[I_2\otimes \left(F_{i{k}_2}{F}_{i{k}_2}^{\text{T}}\right)\right]}\end{array}& .& .\\ I_2\otimes X_j& -\frac{{\left({\displaystyle {\sum }_{k_1=1}^m{\chi }_{1ij{k}_1}^{-1}}\right)}^{-1}}{{\delta }_{1\max }^2}\left(I_2\otimes I_n\right)& .\\ I_2\otimes R_j& 0& -\frac{{\left({\displaystyle {\sum }_{k_2=m\text{+}1}^I{\chi }_{2ij{k}_2}^{-1}}\right)}^{-1}}{{\delta }_{2\max }^2}\left(I_2\otimes I_n\right)\end{array}\right]<0$，$i,j=1,2,\cdots ,N$

(13)

根据引理1，系统的鲁棒稳定性得证。

4. 数值仿真

在本节中，通过两个数值示例来证明上述理论的正确性，由于分数阶微分不能直接仿真，因此使用 [11] 中提出的改进的Oustaloup方法进行了以下仿真。在第一个示例中研究了不确定分数阶系统的稳定性问题，在第二个实例中，设计状态反馈器，使得系统稳定。

4.1. 示例1

考虑具有多胞和范数有界不确定性系统的鲁棒稳定性问题，其中系统参数给出如下：

$N=2$，$\alpha =1.35$

$\text{SS-I}:{\hat{A}}_1=\left[\begin{array}{cc}-1.8& 1.3\\ 0.44& -0.75\end{array}\right],F_{11}=\left[\begin{array}{cc}0.35& 0.41\\ 0.61& 0.13\end{array}\right],F_{12}=\left[\begin{array}{cc}0.33& 0.56\\ 0.47& 0.3\end{array}\right]$

$\text{SS-II}:{\hat{A}}_2=\left[\begin{array}{cc}-1.5& 0.45\\ 0.87& -0.98\end{array}\right],F_{21}=\left[\begin{array}{cc}0.23& 0.31\\ 0.34& 0.51\end{array}\right],F_{22}=\left[\begin{array}{cc}0.12& 0.43\\ 0.55& 0.49\end{array}\right]$

根据定理3.1，使用MATLB提供的LMI工具箱，得到LMI式(6)的解：

$X_1=\left[\begin{array}{cc}0.5565& -0.1182\\ -0.1182& 0.0669\end{array}\right]$，$X_2=\left[\begin{array}{cc}36.6824& -33.2225\\ -33.2225& 37.9283\end{array}\right]$

当 ${\gamma }_1=0.1$ ； ${\gamma }_2=0.9$ 时，给系统加入不确定参数，随机选取不确定参数的值，由图1可以看出随着最大不确定性界 ${\delta }_{1\max }$ 的增加，不确定性因素的影响也会增加，系统的输出响应也将由收敛变的发散。这与定理3.1相一致。

4.2. 示例2

考虑具有多胞和范数有界不确定性系统的鲁棒稳定性问题，其中系统参数给出如下：

$N=2$，$\alpha =1.35$

$\text{SS-I}:{\hat{A}}_1=\left[\begin{array}{ccc}-3& 5& -1\\ -3& 5& -1\\ 2& 2& -5\end{array}\right],F_{11}=\left[\begin{array}{ccc}1& 1& 0\\ 0& 0& 0\\ 0& 0& 0\end{array}\right],F_{12}=\left[\begin{array}{ccc}0& 0& 0\\ 0& 1& 1\\ 0& 0& 0\end{array}\right],B_1=\left[\begin{array}{c}1\\ 2\\ 1\end{array}\right],F_{13}=\left[\begin{array}{c}1\\ 0\\ 0\end{array}\right],F_{14}=\left[\begin{array}{c}0\\ 1\\ 0\end{array}\right]$

图1. 系统在不同不确定界值时的输出响应

$\text{SS-II}:{\hat{A}}_2=\left[\begin{array}{ccc}-3& 4& 0\\ -3& 5& -1\\ 1& 2& -5\end{array}\right],F_{21}=\left[\begin{array}{ccc}1& 1& 0\\ 0& 0& 1\\ 0& 0& 0\end{array}\right],F_{22}=\left[\begin{array}{ccc}0& 0& 0\\ 0& 1& 1\\ 0& 0& 1\end{array}\right],B_2=\left[\begin{array}{c}2\\ 1\\ 1\end{array}\right],F_{23}=\left[\begin{array}{c}0\\ 1\\ 0\end{array}\right],F_{24}=\left[\begin{array}{c}1\\ 0\\ 0\end{array}\right]$

根据定理3.2，使用MATLB提供的LMI工具箱，得到LMI式(10)的解：

$X_1=\left[\begin{array}{ccc}0.0572& -0.0548& 0.0103\\ -0.0548& -0.0055& -0.0598\\ 0.0103& -0.0598& 0.1526\end{array}\right]$，$X_2=\left[\begin{array}{ccc}0.9001& 3.2029& -2.7717\\ 3.2029& -8.7520& 0.9859\\ -2.7717& 0.9859& 5.9512\end{array}\right]$，

$R_1={\left[\begin{array}{c}0.1080\\ -0.3451\\ 0.2591\end{array}\right]}^{\text{T}}$，$R_2={\left[\begin{array}{c}-21.6417\\ 19.5770\\ -6.3547\end{array}\right]}^{\text{T}}$

该系统中，当系统的控制输入 $u\left(t\right)=0$，系统不是渐近稳定的，由图2系统的时间响应可以得出，原系统的状态是发散的。

根据定理3.2，为了保证系统的稳定性，加入状态反馈控制器 $u\left(t\right)=Kx\left(t\right)$ 使得系统渐近稳定，随机选取系统的加权系数 ${\gamma }_1=0.6$ ； ${\gamma }_2=0.4$，系统不确定参数的无穷大范数 ${\delta }_{1\max }={\delta }_{2\max }=0.2$，求解出控制器的增益 $K=\left[\begin{array}{ccc}-3.6615& -1.9968& -1.1921\end{array}\right]$。由图3系统的时间响应可以看出，系统时间响应由图2的发散状态最终变了收敛状态从而说明了定理3.2的正确性。

图2. 系统的特征值和时间响应

图3. 系统的特征值和时间响应

5. 结论

本文推导了在 $\alpha >1$ 时，多胞和无穷范数有界不确定分数阶系统的鲁棒稳定性的充分条件。在不满足鲁棒稳定条件时，设计了相应的线性状态反馈稳定控制器。最后，通过两个数值例子验证了理论的正确性。所有的结果都是根据MATLAB提供的LMI计算得出。仿真结果表明随着不确定性无穷范数的增大，其对系统的稳定性影响逐步增大，系统逐渐失稳。但是对于具体的系统，其无穷范数的上界到底如何确定，这也是后续需要进行进一步研究的问题。同时基于该方法设计的状态反馈控制器可以有效镇定不确定分数阶系统，且控制代价小，适合于工程实现。

文章引用

郭校中. 分数阶系统的鲁棒稳定性研究
Study on the Robust Stability of Fractional-Order System[J]. 动力系统与控制, 2022, 11(01): 19-26. https://doi.org/10.12677/DSC.2022.111003

参考文献