本文将针对分数阶(α>1)系统,同时考虑最常见的两种不确定类型,即多胞不确定和无穷范数。首先,给出了不确定分数阶系统无状态反馈时鲁棒稳定性条件。其次,采用线性矩阵不等式的形式给出了不确定分数阶系统状态反馈控制器的设计方法。最后,通过两个数值算例验证了理论的正确性。 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 (α>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.
本文将针对分数阶( α > 1 )系统,同时考虑最常见的两种不确定类型,即多胞不确定和无穷范数。首先,给出了不确定分数阶系统无状态反馈时鲁棒稳定性条件。其次,采用线性矩阵不等式的形式给出了不确定分数阶系统状态反馈控制器的设计方法。最后,通过两个数值算例验证了理论的正确性。
分数阶系统,多胞不确定,范数有界不确定,鲁棒稳定性
Xiaozhong Guo
Xidian University, Xi’an Shaanxi
Received: Jan. 1st, 2022; accepted: Jan. 14th, 2022; published: Jan. 28th, 2022
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 α > 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 )系统在多胞和范数有界不确定情况下,系统的稳定性条件。设计状态反馈控制器使分数阶系统变得稳定,并通过数值算例,验证理论的正确性。
在本文的后续中,如果没有特别强调或指出,一般符号的含义如下所示。
X T 表示矩阵X的转置, X * 表示矩阵X的共轭转置, X > 0 表示X是正定矩阵。
s y m { X } 等于 X + X * ,矩阵中的 • 表示这是对称矩阵中的对称部分。 ⊗ 表示矩阵的Kronecker积。 D α 是 α 阶Caputo微分算子。
‖ X ‖ ∞ 是X的无穷大范数, I n 是n阶单位矩阵。
引理1. 当 1 < α < 2 时,系统稳定的充分必要条件是存在正定矩阵 X = X T > 0 , X ∈ R n × n 使得如下线性矩阵不等式成立:
s y m { Λ ⊗ ( A X ) } < 0 (1)
其中 ⊗ 表示矩阵的Kronecker积:
Λ = [ sin θ − cos θ cos θ sin θ ] , θ = π − α π 2 (2)
引理2. 对于具有适当维数的正定矩阵 χ > 0 和矩阵 Q , P ,有以下不等式成立:
Q T P + P T Q ≤ Q T χ Q + P T χ − 1 P (3)
引理3. 对于给定的对称矩阵:
Λ = [ A B B T C ] (4)
A , B , C 都具有恰当的维数,则以下三个条件是等价的:
1. Λ < 0 ;
2. A < 0 , C − B T A − 1 B < 0 ;
3. C < 0 , A − B C − 1 B T < 0 ;
引理4. 当 0 < α < 2 时,分数阶系统渐近稳定的充分必要条件是:
| A r g ( e i g ( A ) ) | > α π / 2 (5)
其中 A r g ( ) 是辐角函数, e i g ( A ) 是矩阵A的所有特征值。
定理3.1 当 u ( t ) = 0 时,多胞和范数有界不确定分数阶系统具有鲁棒渐近稳定性的充分条件是:存在矩阵 X j = X j T > 0 , X j ∈ R n × n , X i j ∈ Z + , ( i , j = 1 , 2 , ⋯ , N ) 满足约束:
[ s y m ( Λ ⊗ A ^ i X j ) + ∑ k 1 = 1 m { χ i i k 1 [ I 2 ⊗ ( F i k 1 F i k 1 T ) ] } . I 2 ⊗ X j − ( ∑ k 1 = 1 m χ i i k 1 − 1 ) − 1 δ 1 max 2 ( I 2 ⊗ I n ) ] < 0 (6)
证明. 根据引理1,如果存在正定矩阵 X j 及 X ( γ ) = ∑ j = 1 N γ j X j ,使得分数阶系统渐近稳定,那么则有:
s y m { Λ ⊗ ( A X ) } < 0
⇔ ∑ i = 1 N ∑ j = 1 N γ i γ j s y m ( Λ ⊗ ( A i X j ) ) < 0
⇔ ∑ i = 1 N γ i 2 s y m { Λ ⊗ [ A ^ i + ∑ k 1 = 1 m δ i k 1 F i k 1 ] X i } + ∑ i = 1 N − 1 ∑ j = i + 1 N γ i γ j s y m { Λ ⊗ [ ( A ^ i + ∑ k 1 m δ i k 1 F i k 1 ) X j ] + ( A ^ j + ∑ k 1 m δ j k 1 F j k 1 ) X i } < 0 (7)
根据引理2,存在矩阵 χ i i k 1 > 0 ( i = 1 , 2 , ⋯ , N ) ,
χ i j k 1 > 0 , χ j i k 1 > 0 , ( i = 1 , 2 , ⋯ , N − 1 ) , ( j = i + 1 , ⋯ , N ) ,
且有 ‖ δ i k 1 ‖ ≤ δ 1 max ( i = 1 , 2 , ⋯ , N ; k 1 = 1 , 2 , ⋯ , m ) ,可以得到:
∑ i = 1 N γ i 2 { s y m ( Λ ⊗ A ^ i X i ) + ∑ k 1 = 1 m { χ i i k 1 [ I 2 ⊗ ( F i k 1 F i k 1 T ) ] + χ i i k 1 − 1 δ 1 max 2 [ I 2 ⊗ X i T X i ] } } + ∑ i = 1 N ∑ j = i + 1 N γ i γ j { s y m ( Λ ⊗ A ^ i X j ) + ∑ k 1 = 1 m { χ i j k 1 [ I 2 ⊗ ( F i k 1 F i k 1 T ) ] + χ i j k 1 − 1 δ 1 max 2 [ I 2 ⊗ X j T X j ] } + s y m ( Λ ⊗ A ^ j X i ) + ∑ k 1 = 1 m { χ j i k 1 [ I 2 ⊗ ( F j k 1 F j k 1 T ) ] + χ j i k 1 − 1 δ 1 max 2 [ I 2 ⊗ X i T X i ] } } < 0 (8)
使用引理3,可以将式(8)写为:
[ s y m ( Λ ⊗ A ^ i X j ) + ∑ k 1 = 1 m { χ i i k 1 [ I 2 ⊗ ( F i k 1 F i k 1 T ) ] } . I 2 ⊗ X j − ( ∑ k 1 = 1 m χ i i k 1 − 1 ) − 1 δ 1 max 2 ( I 2 ⊗ I n ) ] < 0 , i , j = 1 , 2 , ⋯ , N (9)
根据引理1,系统的鲁棒稳定性得证。
定理3.2 当 u ( t ) = K x ( t ) 时,多胞和范数有界不确定分数阶系统具有鲁棒渐近稳定性的充分条件是:存在矩阵 X j = X j T > 0 , X j ∈ R n × n , R j ∈ R m × n , χ 1 i j , χ 2 i j ∈ Z + , ( i , j = 1 , 2 , ⋯ , N ) 满足约束:
[ s y m [ Λ ⊗ ( A ^ i X j + B ^ i R j ) ] + ∑ k 1 = 1 m χ 1 i j k 1 [ I 2 ⊗ ( F i k 1 F i k 1 T ) ] + ∑ k 2 = m + 1 I χ 2 i j k 2 [ I 2 ⊗ ( F i k 2 F i k 2 T ) ] . . I 2 ⊗ X j − ( ∑ k 1 = 1 m χ 1 i j k 1 − 1 ) − 1 δ 1 max 2 ( I 2 ⊗ I n ) . I 2 ⊗ R j 0 − ( ∑ k 2 = m + 1 I χ 2 i j k 2 − 1 ) − 1 δ 2 max 2 ( I 2 ⊗ I n ) ] < 0 (10)
式中状态反馈器可描述为:
u ( t ) = K x ( t ) , K ( γ ) = ∑ i = 1 N γ i K i = ∑ i = 1 N γ i R i ∑ i = 1 N X i − 1
证明. 根据引理1,如果存在正定矩阵 X j , X i 及 X ( γ ) = ∑ j = 1 N γ j X j , K ( γ ) = ∑ i = 1 N γ i K i 使得分数阶系统渐近稳定,那么则有:
s y m { Λ ⊗ ( A + B K ) X } < 0
⇔ ∑ i = 1 N γ i 2 s y m { Λ ⊗ [ ( A ^ i + ∑ k 1 = 1 m δ i k 1 F i k 1 ) X i ] + Λ ⊗ [ ( B ^ i + ∑ k 2 = m + 1 I δ i k 2 F i k 2 ) R i ] } + ∑ i = 1 N − 1 ∑ j = i + 1 N γ i γ j s y m { Λ ⊗ [ ( A ^ i + ∑ k 1 m δ i k 1 F i k 1 ) X j + ( A ^ j + ∑ k 2 = m + 1 I δ j k 1 F j k 1 ) X i ] + Λ ⊗ [ ( B ^ i + ∑ k 2 = m + 1 I δ i k 2 F i k 2 ) R j + ( B ^ j + ∑ k 2 = m + 1 I δ j k 2 F j k 2 ) R i ] } < 0 (11)
根据引理2,存在矩阵 χ 1 i i k 1 , χ 2 i i k 2 > 0 , χ 1 i j k 1 , χ 2 i j k 2 > 0 , χ 1 j i k 1 , χ 2 j i k 2 > 0 , ‖ δ i k 1 ‖ ≤ δ 1 max , ( i = 1 , 2 , ⋯ , N − 1 ; j = i + 1 , ⋯ , N ) ,且有 ‖ δ i k 2 ‖ ≤ δ 2 max , ( i = 1 , 2 , ⋯ , N ) , ( k 1 = 1 , 2 , ⋯ , m ) 可以推导出:
∑ i = 1 N γ i 2 { s y m [ Λ ⊗ ( A ^ i X i + B ^ i R i ) ] + ∑ k 1 = 1 m { χ 1 i i k 1 [ I 2 ⊗ ( F i k 1 F i k 1 T ) ] + χ 1 i i k 1 − 1 δ 1 max 2 [ I 2 ⊗ X i T X i ] } + ∑ k 2 = m + 1 I { χ 2 i i k 2 [ I 2 ⊗ ( F i k 2 F i k 2 T ) ] + χ 2 i i k 2 − 1 δ 2 max 2 [ I 2 ⊗ R i T R i ] } } + ∑ i = 1 N − 1 ∑ j = i + 1 N γ i γ j { s y m [ Λ ⊗ ( A ^ i X j + B ^ i R j ) ] + ∑ k 1 = 1 m { χ 1 i j k 1 [ I 2 ⊗ ( F i k 1 F i k 1 T ) ] + χ 1 i j k 1 − 1 δ 1 max 2 [ I 2 ⊗ X j T X j ] } + ∑ k 2 = m + 1 I { χ 2 i j k 2 [ I 2 ⊗ ( F i k 2 F i k 2 T ) ] + χ 2 i j k 2 − 1 δ 2 max 2 [ I 2 ⊗ R j T R j ] } + s y m [ Λ ⊗ ( A ^ j X i + B ^ j R i ) ] + ∑ k 1 = 1 m { χ 1 j i k 1 [ I 2 ⊗ ( F j k 1 F j k 1 T ) ] + χ 1 j i k 1 − 1 δ 1 max 2 [ I 2 ⊗ X i T X i ] } + ∑ k 2 = m + 1 I { χ 2 j i k 2 [ I 2 ⊗ ( F j k 2 F j k 2 T ) ] + χ 2 i i k 2 − 1 δ 2 max 2 [ I 2 ⊗ R i T R i ] } } < 0 (12)
使用引理3,可以将式(12)写为:
[ s y m [ Λ ⊗ ( A ^ i X j + B ^ i R j ) ] + ∑ k 1 = 1 m χ 1 i j k 1 [ I 2 ⊗ ( F i k 1 F i k 1 T ) ] + ∑ k 2 = m + 1 I χ 2 i j k 2 [ I 2 ⊗ ( F i k 2 F i k 2 T ) ] . . I 2 ⊗ X j − ( ∑ k 1 = 1 m χ 1 i j k 1 − 1 ) − 1 δ 1 max 2 ( I 2 ⊗ I n ) . I 2 ⊗ R j 0 − ( ∑ k 2 = m + 1 I χ 2 i j k 2 − 1 ) − 1 δ 2 max 2 ( I 2 ⊗ I n ) ] < 0 , i , j = 1 , 2 , ⋯ , N
(13)
根据引理1,系统的鲁棒稳定性得证。
在本节中,通过两个数值示例来证明上述理论的正确性,由于分数阶微分不能直接仿真,因此使用 [
考虑具有多胞和范数有界不确定性系统的鲁棒稳定性问题,其中系统参数给出如下:
N = 2 , α = 1.35
SS-I : A ^ 1 = [ − 1.8 1.3 0.44 − 0.75 ] , F 11 = [ 0.35 0.41 0.61 0.13 ] , F 12 = [ 0.33 0.56 0.47 0.3 ]
SS-II : A ^ 2 = [ − 1.5 0.45 0.87 − 0.98 ] , F 21 = [ 0.23 0.31 0.34 0.51 ] , F 22 = [ 0.12 0.43 0.55 0.49 ]
根据定理3.1,使用MATLB提供的LMI工具箱,得到LMI式(6)的解:
X 1 = [ 0.5565 − 0.1182 − 0.1182 0.0669 ] , X 2 = [ 36.6824 − 33.2225 − 33.2225 37.9283 ]
当 γ 1 = 0.1 ; γ 2 = 0.9 时,给系统加入不确定参数,随机选取不确定参数的值,由图1可以看出随着最大不确定性界 δ 1 max 的增加,不确定性因素的影响也会增加,系统的输出响应也将由收敛变的发散。这与定理3.1相一致。
考虑具有多胞和范数有界不确定性系统的鲁棒稳定性问题,其中系统参数给出如下:
N = 2 , α = 1.35
SS-I : A ^ 1 = [ − 3 5 − 1 − 3 5 − 1 2 2 − 5 ] , F 11 = [ 1 1 0 0 0 0 0 0 0 ] , F 12 = [ 0 0 0 0 1 1 0 0 0 ] , B 1 = [ 1 2 1 ] , F 13 = [ 1 0 0 ] , F 14 = [ 0 1 0 ]
图1. 系统在不同不确定界值时的输出响应
SS-II : A ^ 2 = [ − 3 4 0 − 3 5 − 1 1 2 − 5 ] , F 21 = [ 1 1 0 0 0 1 0 0 0 ] , F 22 = [ 0 0 0 0 1 1 0 0 1 ] , B 2 = [ 2 1 1 ] , F 23 = [ 0 1 0 ] , F 24 = [ 1 0 0 ]
根据定理3.2,使用MATLB提供的LMI工具箱,得到LMI式(10)的解:
X 1 = [ 0.0572 − 0.0548 0.0103 − 0.0548 − 0.0055 − 0.0598 0.0103 − 0.0598 0.1526 ] , X 2 = [ 0.9001 3.2029 − 2.7717 3.2029 − 8.7520 0.9859 − 2.7717 0.9859 5.9512 ] ,
R 1 = [ 0.1080 − 0.3451 0.2591 ] T , R 2 = [ − 21.6417 19.5770 − 6.3547 ] T
该系统中,当系统的控制输入 u ( t ) = 0 ,系统不是渐近稳定的,由图2系统的时间响应可以得出,原系统的状态是发散的。
根据定理3.2,为了保证系统的稳定性,加入状态反馈控制器 u ( t ) = K x ( t ) 使得系统渐近稳定,随机选取系统的加权系数 γ 1 = 0.6 ; γ 2 = 0.4 ,系统不确定参数的无穷大范数 δ 1 max = δ 2 max = 0.2 ,求解出控制器的增益 K = [ − 3.6615 − 1.9968 − 1.1921 ] 。由图3系统的时间响应可以看出,系统时间响应由图2的发散状态最终变了收敛状态从而说明了定理3.2的正确性。
图2. 系统的特征值和时间响应
图3. 系统的特征值和时间响应
本文推导了在 α > 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