基于反演设计和容错控制研究Conformable分数阶单机无穷大电力系统的混沌同步问题，在设计过程中引入虚拟的误差函数，利用李雅普诺夫稳定性理论确保带有执行器故障的被控系统的鲁棒同步，最后，模拟结果表明所提控制策略的有效性。

单机无穷大电力系统，容错同步，Conformable分数阶导数，反演设计

Qunli Zhang, Xuerong Fu, Ying Wei

College of Mathematics and Statistics, Heze University, Heze Shandong

Received: Dec. 26^th, 2021; accepted: Jan. 16^th, 2022; published: Jan. 29^th, 2022

A synchronization strategy based on backstepping design and tolerant control is studied for conformable fractional order unipolar infinite power system. The virtual control error functions are introduced in the design, and the Lyapunov stability theory is used to ensure the robust synchronization of the controlled systems with the actuator faults. Finally, the simulation results show the effectiveness of the proposed control strategy.

Keywords:Unipolar Infinite Power System, Tolerant Synchronization, Conformable Fractional Order Derivative, Backstepping Design

Copyright © 2022 by author(s) and beplay安卓登录

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

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

分数阶微积分是当前非线性研究领域的热点，能够较准确简洁地阐述具有历史记忆性和空间全域相关性等力学与物理过程，在量子力学、电磁振动、材料力学等领域有着十分广泛的应用 [ 1 ] [ 2 ] [ 3 ]。Khalil等 [ 4 ] 提出一种表现形式简洁易懂的Conformable分数阶微分定义，与整数阶微分有着良好分析、计算等统一 [ 5 ]，受到学者们的极大关注。Atangana等 [ 6 ] 继续发展了Conformable分数阶微分理论，为后继学者的研究奠定了基础。Feng和Meng [ 7 ] 根据Conformable分数阶微积分性质，研究了一类Conformable分数阶动力学方程的振荡和渐近性质。当前国内外关于Conformable分数阶微分的研究虽然有一些代表性成果 [ 7 ] [ 8 ] [ 9 ] [ 10 ] [ 11 ]，但尚处于初期阶段。

近年来，研究成果为数不多的容错控制已经成为一种有效的增加系统安全和可靠性的方法 [ 12 ] [ 13 ] [ 14 ] [ 15 ]，其设计方法一般有：被动方法、主动方法，即使故障发生，被控系统仍能保持渐近稳定。反演法 [ 16 ] [ 17 ] [ 18 ] [ 19 ] [ 20 ] 可将复杂的系统分解成不超过原系统阶数的子系统，通过设计中间虚拟控制量、修正虚拟控制律确保系统的稳定性，是一种性能优良的控制方法。Zhang等 [ 16 ] 设计反演控制器，对Morse吸引子做了混沌同步分析。Deng等 [ 17 ] 通过引入辅助信号，引入基于期望轨迹的前馈补偿概念，提出了一种实用的输出反馈反演控制器，用于液压执行器的运动控制。Zhang等 [ 19 ] 针对磁浮列车系统的稳定悬浮控制问题，设计了基于反步法的非线性控制器。Su等 [ 20 ] 在控制设计中引入虚拟控制函数，使用改进的自适应定律估计外部扰动的未知边界，利用未知边界的更新值，为不确定的非线性系统合成一类连续自适应鲁棒状态反馈控制器。

有效对电力系统进行稳定分析、评估是保障电力系统安全稳定运行的前提，Zhang等 [ 21 ] 基于轨迹断面特征根研究了电力系统暂态稳定定量评估方法，Zhao等 [ 22 ] 考虑二次系统具有时滞，用新的Lyapunov-Krasovskii泛函构造函数和解析方法得到了稳定判据，Wang等 [ 23 ] 首先利用分岔图及相图对系统的基本动力学行为进行分析，然后借助递归型滑模控制方法为混沌电力系统设计励磁控制器并给出具体表达式。如今单机无穷大系统模式在现实中的应用较为关注，文献 [ 24 ] [ 25 ] [ 26 ] [ 27 ] [ 28 ] 采用多种控制方法对系统性能指标进行稳定性分析。Yan [ 29 ] 采用Lyapunov指数谱、分岔图和吸引子相图分析了Conformable分数阶单机无穷大电力系统的分岔和混沌现象。整数阶电力系统模型已经无法满足研究需要，分析Conformable分数阶单机无穷大电力系统混沌振荡产生的机理有着重要的理论和实际意义 [ 29 ] [ 30 ] [ 31 ] [ 32 ]。

本文基于容错控制理论，设计中间虚拟控制量，通过修正虚拟控制律，基于反演方法设计带有执行器故障的控制器，研讨Conformable分数阶单机无穷大电力系统混沌同步问题。

定义：函数 x ( t ) : [ a , + ∞ ) → R 的conformable分数阶导数为

d q x ( t ) d t q = lim ε → 0 x ( t + ε t 1 − q ) − x ( t ) ε ,     t > 0 ,   0 < q ≤ 1 ,

d q x ( 0 ) d t q = lim t → 0 + d q x ( t ) d t q .

注：函数 x ( t ) 可导，则有 d q x ( t ) d t q = x ′ ( t ) t 1 − q 。

定理1：若 d q x ( t ) d t q = f ( x ( t ) ) ，则有 x ( t ) = x ( 0 ) + ∫ 0 t f ( x ( s ) ) s q − 1 d s 。

证明：

d q x ( t ) d t q = lim ε → 0 x ( t + ε t 1 − q ) − x ( t ) ε = lim ε → 0 ∫ 0 t + ε t 1 − q f ( x ( s ) ) s q − 1 d s − ∫ 0 t f ( x ( s ) ) s q − 1 d s ε = lim ε → 0 ∫ t t + ε t 1 − q f ( x ( s ) ) s q − 1 d s ε = lim ε → 0 f ( x ( t + ε t 1 − q ) ) ( t + ε t 1 − q ) q − 1 t 1 − q = f ( x ( t ) ) .

所以有 x ( t ) = x ( 0 ) + ∫ 0 t f ( x ( s ) ) s q − 1 d s 成立。

定理2：设 d q x ( t ) d t q = f ( x ( t ) ) ，则 ∀ r > 0 ，有

d q d t q ‖ x ( t ) − y ( t ) ‖ ≤ ( x ( t ) − y ( t ) ) T ( f ( x ( t ) ) − f ( y ( t ) ) ) ‖ x ( t ) − y ( t ) ‖ 2 ‖ x ( t ) − y ( t ) ‖ ,

其中函数 y ( t ) 满足 d q y ( t ) d t q = f ( y ( t ) ) 。

证明： ∀ r > 0 ， d q d t q ( e r t q q x ( t ) ) = r e r t q q x ( t ) + e r t q q f ( x ( t ) ) 。

于是有

d q d t q ( e r t q q ( x ( t ) − y ( t ) ) ) = r e r t q q ( x ( t ) − y ( t ) ) + e r t q q ( f ( x ( t ) ) − f ( y ( t ) ) ) .

由定理1，得

e r t q q ( x ( t ) − y ( t ) ) = ( x ( 0 ) − y ( 0 ) ) + ∫ 0 t e r s q q ( r ( x ( s ) − y ( s ) ) + ( f ( x ( s ) ) − f ( y ( s ) ) ) ) s q − 1 d s .

e r t q q ‖ x ( t ) − y ( t ) ‖ ≤ ‖ x ( 0 ) − y ( 0 ) ‖ + ∫ 0 t e r s q q s q − 1 ‖ r ( x ( s ) − y ( s ) ) + ( f ( x ( s ) ) − f ( y ( s ) ) ) ‖ d s .

r e r t q q ‖ x ( t ) − y ( t ) ‖ + e r t q q d q d t q ‖ x ( t ) − y ( t ) ‖ ≤ e r t q q ‖ r ( x ( t ) − y ( t ) ) + ( f ( x ( t ) ) − f ( y ( t ) ) ) ‖ .

d q d t q ‖ x ( t ) − y ( t ) ‖ ≤ ‖ r ( x ( t ) − y ( t ) ) + ( f ( x ( t ) ) − f ( y ( t ) ) ) ‖ − r ‖ x ( t ) − y ( t ) ‖ = 2 r ( x ( t ) − y ( t ) ) T ( f ( x ( t ) ) − f ( y ( t ) ) ) + ‖ f ( x ( t ) ) − f ( y ( t ) ) ‖ ‖ r ( x ( t ) − y ( t ) ) + ( f ( x ( t ) ) − f ( y ( t ) ) ) ‖ + r ‖ x ( t ) − y ( t ) ‖ .

由于 2 r ( x ( t ) − y ( t ) ) T ( f ( x ( t ) ) − f ( y ( t ) ) ) + ‖ f ( x ( t ) ) − f ( y ( t ) ) ‖ ‖ r ( x ( t ) − y ( t ) ) + ( f ( x ( t ) ) − f ( y ( t ) ) ) ‖ + r ‖ x ( t ) − y ( t ) ‖ 关于 r > 0 单调递减，所以

d q d t q ‖ x ( t ) − y ( t ) ‖ ≤ lim r → + ∞ 2 r ( x ( t ) − y ( t ) ) T ( f ( x ( t ) ) − f ( y ( t ) ) ) + ‖ f ( x ( t ) ) − f ( y ( t ) ) ‖ ‖ r ( x ( t ) − y ( t ) ) + ( f ( x ( t ) ) − f ( y ( t ) ) ) ‖ + r ‖ x ( t ) − y ( t ) ‖ = ( x ( t ) − y ( t ) ) T ( f ( x ( t ) ) − f ( y ( t ) ) ) ‖ x ( t ) − y ( t ) ‖ = ( x ( t ) − y ( t ) ) T ( f ( x ( t ) ) − f ( y ( t ) ) ) ‖ x ( t ) − y ( t ) ‖ 2 ‖ x ( t ) − y ( t ) ‖ .

带有执行器故障函数得Conformable分数阶单机无穷大电力系统混沌微分方程 [ 29 ] [ 30 ] 为驱动系统：

{ d q x 1 ( t ) d t q = x 2 ( t ) , d q x 2 ( t ) d t q = − 0.5 x 2 ( t ) + 2.4 x 3 ( t ) + 0.2 − sin x 1 ( t ) , d q x 3 ( t ) d t q = x 4 ( t ) , d q x 4 ( t ) d t q = − x 3 ( t ) + m ( t ) , (1)

其中 q = 0.6 ， m ( t ) 是执行器故障函数。

响应系统为：

{ d q y 1 ( t ) d t q = y 2 ( t ) , d q y 2 ( t ) d t q = − 0.5 y 2 ( t ) + 2.4 y 3 ( t ) + 0.2 − sin y 1 ( t ) , d q y 3 ( t ) d t q = y 4 ( t ) , d q y 4 ( t ) d t q = − y 3 ( t ) + m ⌢ ( t ) − p ( e 4 ( t ) + k ( e 3 ( t ) + e 2 ( t ) + e 1 ( t ) ) , (2)

其中 m ⌢ ( t ) 是 m ( t ) 的估计。

设误差 e i ( t ) = y i ( t ) − x i ( t ) , i = 1 , 2 , 3 , 4 ， m ˜ ( t ) = m ⌢ ( t ) − m ( t ) ，则有如下结论：

定理3：若执行器故障误差 m ˜ ( t ) 是关于t的可导函数，满足

d q m ˜ ( t ) d t q = − ( e 4 ( t ) + k ( e 3 ( t ) + e 2 ( t ) + e 1 ( t ) ) ) − r m ˜ ( t ) , (3)

当存在常数 k > 0 , p > 0 , r > 0 ，满足

W = ( − 1 3 4 3 4 3 4 − 1.9 1 4 3 4 1 4 − ( k − 2.4 ) ) < 0 ， Q = ( − 1 3 4 3 4 3 4 k 0 3 4 − 1.9 1 4 1 − 1.9 k 2 0 3 4 1 4 2.4 − k 2.4 k − k 2 2 0 3 4 k 1 − 1.9 k 2 2.4 k − k 2 2 k − p 0 0 0 0 0 − r ) < 0 (4)

时，驱动系统(1)和响应系统(2)混沌同步。

证明：由(1)和(2)得误差系统：

{ d q e 1 ( t ) d t q = e 2 ( t ) , d q e 2 ( t ) d t q = − 0.5 e 2 ( t ) + 2.4 e 3 ( t ) − 2 sin e 1 ( t ) 2 cos x 1 ( t ) + y 1 ( t ) 2 , d q e 3 ( t ) d t q = e 4 ( t ) , d q e 4 ( t ) d t q = − e 3 ( t ) + m ˜ ( t ) − p ( e 4 ( t ) + k ( e 3 ( t ) + e 2 ( t ) + e 1 ( t ) ) ) , (5)

由定理2和误差系统(5)中的第一个方程得：

d q ‖ e 1 ( t ) ‖ d t q ≤ e 1 ( t ) e 2 ( t ) ‖ e 1 ( t ) ‖ 2 ‖ e 1 ( t ) ‖ .

令 e 2 ( t ) = − e 1 ( t ) ，则 d q ‖ e 1 ( t ) ‖ d t q ≤ − ‖ e 1 ( t ) ‖ ，

‖ e 1 ( t ) ‖ ≤ ‖ e 1 ( 0 ) ‖ e − t q q → 0 ,       t → + ∞ .

设 α 2 ( t ) = e 2 ( t ) + e 1 ( t ) ，则结合(5)中第一、二个方程有

{ d q e 1 ( t ) d t q = α 2 ( t ) − e 1 ( t ) , d q α 2 ( t ) d t q = 0.5 α 2 ( t ) − 0.5 e 1 ( t ) + 2.4 e 3 ( t ) − 2 sin e 1 ( t ) 2 cos x 1 ( t ) + y 1 ( t ) 2 .

再次用定理2得：

d q ‖ ( e 1 ( t ) , α 2 ( t ) ) ‖ d t q ≤ − e 1 2 ( t ) + 1.5 | e 1 ( t ) | | α 2 ( t ) | + 0.5 α 2 2 ( t ) + 2.4 α 2 ( t ) e 3 ( t ) ‖ ( e 1 ( t ) , α 2 ( t ) ) ‖ 2 ‖ ( e 1 ( t ) , α 2 ( t ) ) ‖ .

令 e 3 ( t ) = − α 2 ( t ) ，则有

d q ‖ ( e 1 ( t ) , α 2 ( t ) ) ‖ d t q ≤ − e 1 2 ( t ) + 1.5 | e 1 ( t ) | | α 2 ( t ) | − 1.9 α 2 2 ( t ) ‖ ( e 1 ( t ) , α 2 ( t ) ) ‖ 2 ‖ ( e 1 ( t ) , α 2 ( t ) ) ‖ ≤ − 0.25 ‖ ( e 1 ( t ) , α 2 ( t ) ) ‖ .

于是有 ‖ ( e 1 ( t ) , α 2 ( t ) ) ‖ ≤ ‖ ( e 1 ( 0 ) , α 2 ( 0 ) ) ‖ e − 0.25 t q q → 0 ， t → + ∞ 。

设 α 3 ( t ) = e 3 ( t ) + α 2 ( t ) ，则结合(5)中第一、二、三个方程有

{ d q e 1 ( t ) d t q = α 2 ( t ) − e 1 ( t ) , d q α 2 ( t ) d t q = − 1.9 α 2 ( t ) − 0.5 e 1 ( t ) + 2.4 α 3 ( t ) − 2 sin e 1 ( t ) 2 cos x 1 ( t ) + y 1 ( t ) 2 , d q α 3 ( t ) d t q = 2.4 α 3 ( t ) − 1.9 α 2 ( t ) − 0.5 e 1 ( t ) − 2 sin e 1 ( t ) 2 cos x 1 ( t ) + y 1 ( t ) 2 + e 4 ( t ) .

再次利用定理2得：

d q ‖ ( e 1 ( t ) , α 2 ( t ) , α 3 ( t ) ) ‖ d t q ≤ [ − e 1 2 ( t ) + 1.5 | e 1 ( t ) | | α 2 ( t ) | − 1.9 α 2 2 ( t ) + 0.5 α 2 ( t ) α 3 ( t ) ‖ ( e 1 ( t ) , α 2 ( t ) , α 3 ( t ) ) ‖ 2     + 2.4 α 3 2 ( t ) + 1.5 | e 1 ( t ) | | α 3 ( t ) | + e 4 ( t ) α 3 ( t ) ‖ ( e 1 ( t ) , α 2 ( t ) , α 3 ( t ) ) ‖ ] ‖ ( e 1 ( t ) , α 2 ( t ) , α 3 ( t ) ) ‖ .

令 e 4 ( t ) = − k α 3 ( t ) ，则有

d q ‖ ( e 1 ( t ) , α 2 ( t ) , α 3 ( t ) ) ‖ d t q ≤ ( | e 1 ( t ) | | α 2 ( t ) | | α 3 ( t ) | ) T W ( | e 1 ( t ) | | α 2 ( t ) | | α 3 ( t ) | ) ‖ ( e 1 ( t ) , α 2 ( t ) , α 3 ( t ) ) ‖ 2 ‖ ( e 1 ( t ) , α 2 ( t ) , α 3 ( t ) ) ‖ .

由条件(4)中 W < 0 知， ‖ ( e 1 ( t ) , α 2 ( t ) , α 3 ( t ) ) ‖ → 0 ， t → + ∞ 。

设 α 4 ( t ) = e 4 ( t ) + k α 3 ( t ) ，则结合(5)中第一、二、三、四个方程有

{ d q e 1 ( t ) d t q = α 2 ( t ) − e 1 ( t ) , d q α 2 ( t ) d t q = − 1.9 α 2 ( t ) − 0.5 e 1 ( t ) + 2.4 α 3 ( t ) − 2 sin e 1 ( t ) 2 cos x 1 ( t ) + y 1 ( t ) 2 , d q α 3 ( t ) d t q = ( 2.4 − k ) α 3 ( t ) − 1.9 α 2 ( t ) − 0.5 e 1 ( t ) − 2 sin e 1 ( t ) 2 cos x 1 ( t ) + y 1 ( t ) 2 + α 4 ( t ) , d q α 4 ( t ) d t q = ( k − p ) α 4 ( t ) + ( 2.4 k − k 2 − 1 ) α 3 ( t ) + ( 1 − 1.9 k ) α 2 ( t ) − 0.5 k e 1 ( t )                               − 2 k sin e 1 ( t ) 2 cos x 1 ( t ) + y 1 ( t ) 2 + m ˜ ( t ) .

再次利用定理2和(3)得：

d q ‖ ( e 1 ( t ) , α 2 ( t ) , α 3 ( t ) , α 4 ( t ) , m ˜ ( t ) ) ‖ d t q ≤ Z T ( t ) Q Z ( t ) ‖ ( e 1 ( t ) , α 2 ( t ) , α 3 ( t ) , α 4 ( t ) , m ˜ ( t ) ) ‖ 2 ‖ ( e 1 ( t ) , α 2 ( t ) , α 3 ( t ) , α 4 ( t ) , m ˜ ( t ) ) ‖ ,

其中 Z ( t ) = ( | e 1 ( t ) | , | α 2 ( t ) | , | α 3 ( t ) | , | α 4 ( t ) | , | m ˜ ( t ) | ) T 。由条件(4)知，

‖ ( e 1 ( t ) , α 2 ( t ) , α 3 ( t ) , α 4 ( t ) , m ˜ ( t ) ) ‖ → 0 , t → + ∞ .

我们取 k = 4 , r = 3 , p = 200 , m ( t ) = 0.6 + sin ( 2 x 4 ( t ) ) ，满足条件(4)，得到conformable分数阶单机无穷大电力系统(1)和(2)的混沌同步情况(见图1~4)。

图1. 状态分量 x 1 ( t ) , y 1 ( t ) 同步

图2. 状态分量 x 2 ( t ) , y 2 ( t ) 同步

图3. 状态分量 x 3 ( t ) , y 3 ( t ) 同步

图4. 状态分量 x 4 ( t ) , y 4 ( t ) 同步

注：当系统(1)、(2)中的第一二三个方程也存在执行器故障函数时，可模仿定理的推导得到类似结论。

本文基于容错控制、反演设计研究了Conformable分数阶单机无穷大电力系统的混沌同步问题。在故障存在的情况下，合理假设中间误差函数用反演法设计出控制器，仍能使系统具有渐近稳定性。

本文十分感谢山东省自然科学基金(项目编号：ZR2020QA002，ZR2017MA029)的资助。

张群力,付雪荣,魏 莹. Conformable分数阶单机无穷大电力系统的混沌容错同步Tolerant Synchronization Analysis for Conformable Fractional Order Unipolar Infinite Power System[J]. 应用数学进展, 2022, 11(01): 590-599. https://doi.org/10.12677/AAM.2022.111066