一种充分光滑分数阶投影算子的构造方法
A Construction Method of Sufficiently Smooth Fractional-Order Projection Operator
DOI: 10.12677/AAM.2023.1211462, PDF, HTML, XML, 下载: 113  浏览: 160  科研立项经费支持
作者: 孙 珂, 马志垚*:辽宁工业大学理学院,辽宁 锦州
关键词: 分数阶投影算子自适应律Fractional-Order Projection Operator Adaptive Laws
摘要: 本文首先介绍了Caputo型分数阶导数的定义及其满足的相关引理。其次为了构造自适应律,提出了一种充分光滑的分数阶参数投影算子构造方法。最后给出了所提投影算子满足的一些性质,并提供了严密的证明。
Abstract: This paper first introduces the definition of Caputo type fractional-order derivative and the correla-tion lemma about it. Secondly, a construction method of sufficiently smooth fractional-order projec-tion operator is proposed for adaptive laws. Finally, some properties satisfied by the proposed pro-jection operator are given, and a strict proof is provided.
文章引用:孙珂, 马志垚. 一种充分光滑分数阶投影算子的构造方法[J]. 应用数学进展, 2023, 12(11): 4697-4701. https://doi.org/10.12677/AAM.2023.1211462

1. 引言

投影算子是一种常用的自适应控制算法,它的主要作用是将参数投影到一个预定义的空间或集合中,从而限制参数的范围,以满足特定的约束条件。一般来说,基于投影的参数自适应律是不连续的,这违反了微分方程经典解存在的Lipschitz条件。因此,从现实的角度来看,这种不连续性是不可取的。为了消除这个限制,Pomet等人在 [1] 中通过在凸集周围引入边界层,从而得到Lipschitz连续投影算子。在此基础上,Cai等人在 [2] 中利用更强的任意多次连续可微性取代了Lipschitz连续性,同时对其他投影性质进行了微调。进而上述投影算子的结论在理论研究和实际工程领域中也得到越来越广泛的应用 [3] [4] [5] 。然而,上述投影算子只适用于整数阶系统的参数自适应控制中,分数阶系统是对传统控制理论中的整数阶系统的扩展,可以更准确地描述具有非刚性的动态系统,如热力学系统、柔性系统等。进而在车辆悬架系统、分形与混沌、电池管理系统、机器人等领域都有广泛的应用 [6] [7] [8] 。因此,分数阶系统控制设计问题目前已经成为学者们研究的热点问题。但目前还没有分数阶投影算子用于抑制参数的范围,因此本文基于分数阶相关定义和引理,给出了一种充分光滑的分数阶参数投影算子构造方法。总结了所提投影算子满足的四条性质,并予以证明。

2. Caputo型分数阶导数定义及相关引理

定义1 [6] :假设 P : [ t 0 , + ) R 是连续可微函数,它的 α 阶Caputo型分数阶导数可定义:

D 0 C t α P ( t ) = 1 Γ ( ω α ) 0 t P ( ω ) ( τ ) ( t τ ) α + 1 ω d τ (1)

其中 ω 1 α ω ω 为正整数。 Γ ( ) = 0 + τ 1 e τ d τ Gamma 函数,且满足 Γ ( 1 ) = 1

引理1 [9] :若 ( t ) 是一个光滑函数,则下面不等式成立

1 2 D 0 C t α ( T ( t ) ( t ) ) T ( t ) D 0 C t α ( t ) (2)

引理2 [10] :假设 V = 1 2 z 2 ,若 z D 0 C t α z < 0 ,则有

3. 分数阶投影算子

θ Ω : { θ | θ θ 0 } 是一个未知的常参数,其中 θ 0 > 0 。令 θ ˜ = θ θ ^ 为参数估计误差,我们给出一

种充分光滑的分数阶投影算子构造方法,具体表示如下:

D 0 C t α θ ^ = Proj ( μ , θ ^ ) = μ η 1 η 2 2 ( ε 2 + 2 ε θ 0 ) n + 1 θ 0 2 p ( θ ^ ) (3)

其中

p ( θ ^ ) = θ ^ T θ ^ θ 0 2 (4)

η 1 = { p n + 1 ( θ ^ ) , if p ( θ ^ ) > 0 0 , else , η 2 = 1 2 p ( θ ^ ) T μ + ( 1 2 p ( θ ^ ) T μ ) 2 + δ 2 (5)

其中 ε > 0 δ > 0 ,符号表示梯度。

定理1:若 θ ^ ( 0 ) Ω ,所提分数阶投影算子(3)~(5)满足如下性质:

1. θ ^ ( t ) θ 0 + ε , t 0 .

2. θ ˜ Proj ( μ , θ ^ ) θ ˜ T μ .

3. Proj ( μ , θ ^ ) μ ( 1 + ( θ 0 + ε θ 0 ) 2 ) + θ 0 + ε θ 0 2 δ .

4. Proj ( μ , θ ^ ) is C n .

证明:1. 当 p ( θ ^ ) 0 ,由(4)可得 θ ^ θ 0 < θ 0 + ε 。取 V θ ^ = 1 2 θ ^ T θ ^ ,根据引理1,可得

D 0 C t α V θ ^ θ ^ T D 0 C t α θ ^ = θ ^ T μ η 1 η 2 2 ( ε 2 + 2 ε θ 0 ) n + 1 θ 0 2 θ ^ T θ ^ (6)

p ( θ ^ ) > 0 p ( θ ^ ) T μ = 2 θ ^ T μ 0 时,可得下面不等式

D 0 C t α V θ ^ η 1 η 2 2 ( ε 2 + 2 ε θ 0 ) n + 1 θ 0 2 θ ^ T θ ^ < 0 (7)

由于 η 1 , η 2 > 0 θ ^ 收敛到了 Ω 中,且 θ ^ ( 0 ) Ω 。当 p ( θ ^ ) T μ = 2 θ ^ T μ > 0 ,可推导出

D 0 C t α V θ ^ θ ^ T μ [ ( θ ^ T θ ^ θ 0 2 ) n + 1 ( ε 2 + 2 ε θ 0 ) n + 1 θ 0 2 θ ^ T θ ^ ] θ ^ T μ ( θ ^ T θ ^ θ 0 2 ) n + 1 ( ε 2 + 2 ε θ 0 ) n + 1 θ 0 2 θ ^ T θ ^ ( θ ^ T μ ) 2 + δ 2 [ 1 ( θ ^ T θ ^ θ 0 2 ) n + 1 ( ε 2 + 2 ε θ 0 ) n + 1 θ 0 2 θ ^ T θ ^ ] θ ^ T μ

其中 p ( θ ^ ) = θ ^ T θ ^ θ 0 2 > 0 θ ^ 2 > θ 0 2 。假设 θ ^ 2 = θ 0 + ε ,其中 ε 为一任意正常数,则有

D 0 C t α V θ ^ [ 1 ( ε 2 + 2 ε θ 0 ) n + 1 ( ε 2 + 2 ε θ 0 ) n + 1 θ 0 2 θ ^ T θ ^ ] θ ^ T μ = ( 1 θ ^ T θ ^ θ 0 2 ) θ ^ T μ (8)

其中 θ ^ T μ > 0 θ ˜ T θ ^ > θ 0 2 。因此 D 0 C t α V θ ^ < 0 θ ^ 被吸引到 Ω 中。

2. 当 p ( θ ^ ) 0 ,可得 η 1 = 0 ,则有 θ ˜ T Proj ( μ , θ ^ ) = θ ˜ T μ 。当 p ( θ ^ ) > 0 ,有

θ ˜ T Proj ( μ , θ ^ ) θ ˜ T μ θ ˜ T ( Proj ( μ , θ ^ ) μ ) θ ˜ T θ ^ [ η 1 η 2 ( ε 2 + 2 ε θ 0 ) n + 1 θ 0 2 ] 0 (9)

因为 η 1 η 2 ( ε 2 + 2 ε θ 0 ) n + 1 θ 0 2 > 0 ,则要求 θ ˜ T θ ^ 0 。显然,下面的不等式是成立的。

θ ˜ T θ ^ = ( θ θ ^ ) T θ ^ = 1 2 ( θ ^ T θ ^ θ T θ ) 1 2 ( θ ^ θ ) T ( θ ^ θ ) 1 2 ( θ ^ 2 θ 2 ) < 0 (10)

3. 当 p ( θ ^ ) 0 Proj ( μ , θ ^ ) = μ 。当 p ( θ ^ ) > 0 θ 0 2 < θ ^ T θ ^ ( θ 0 + ε ) 2 ,根据性质1,可得

θ ^ T θ ^ θ 0 2 ε 2 + 2 ε θ 0 1 (11)

计算 Proj ( μ , θ ^ ) 得到

Proj ( μ , θ ^ ) μ + ( θ ^ T θ ^ θ 0 2 ε 2 + 2 ε + θ 0 ) n + 1 η 2 θ ^ θ 0 2 μ + η 2 θ ^ θ 0 2 (12)

因为 η 2 = θ ^ T μ + ( θ ^ T μ ) 2 + δ 2 ,(12) 满足

Proj ( μ , θ ^ ) μ + θ 0 + ε θ 0 2 ( 2 θ ^ μ + δ ) μ ( 1 + 2 ( θ 0 + ε ) 2 θ 0 2 ) + θ 0 + ε θ 0 2 δ (13)

4. 由于 C n 函数的乘积依然是 C n ,假设投影算子(3)中的函数是 C n ,显然 p ( θ ^ ) η 2 ( θ ^ , μ ) C n 。.对于函数 f ( x ) = { x n + 1 , if x > 0 0 , else ,有

f ( i ) ( x ) = { ( n + 1 ) ! ( n + 1 i ) ! x n + 1 i , if x > 0 0 , else (14)

因为

lim x 0 + f ( i 1 ) ( x ) f ( i 1 ) ( 0 ) x 0 = lim x 0 + ( n + 1 ) ! ( n + 2 i ) ! x n + 2 i x = 0

lim x 0 f ( i 1 ) ( x ) f ( i 1 ) ( 0 ) x 0 = lim x 0 0 0 x = 0

则一定可得 η 1 ( θ ^ ) C n 。至此,定理1证毕。

4. 结论

本文基于Caputo型分数阶导数相关定义和引理,提出了一种充分光滑的分数阶参数投影算子构造方法。进而总结出了所提投影算子满足的四条性质,并给出了详细的证明过程。所提分数阶投影算子一方面可用于分数阶系统中对于参数的抑制问题,还由于其充分光滑的特性适用于对参数反复求导的情况,更利于实际应用。本文的不足之处在于只是给出了数学推导,缺少实际仿真算例的验证。在未来的工作中,我们会对所提投影算子进行实际仿真验证,并推广到高阶参数系统的控制设计中。

基金项目

辽宁工业大学博士科研启动基金(XB2021015)。

参考文献

NOTES

*通讯作者。

参考文献

[1] Pomet, J.B. and Praly, L. (1992) Adaptive Nonlinear Regulation: Estimation from the Lyapunov Equation. IEEE Trans-actions on Automatic Control, 37, 729-740.
https://doi.org/10.1109/9.256328
[2] Cai, Z., de Queiroz, M.S. and Dawson, D.M. (2006) A Sufficiently Smooth Projection Operator. IEEE Transactions on Automatic Control, 51, 135-139.
https://doi.org/10.1109/TAC.2005.861704
[3] Tong, S., Min, X. and Li, Y. (2020) Observer-Based Adaptive Fuzzy Tracking Control for Strict-Feedback Nonlinear Systems with Unknown Control Gain Functions. IEEE Transactions on Cybernetics, 50, 3903-3913.
https://doi.org/10.1109/TCYB.2020.2977175
[4] 颜文旭, 黄杰, 许德智. 考虑边端效应的直线感应电机投影自适应指令滤波反推控制[J]. 控制与决策, 2018, 33(8): 1399-1406.
[5] 李静, 王晨, 张家旭. 基于自适应快速终端滑模的车轮滑移率跟踪控制[J]. 浙江大学学报(工学版), 2021, 55(1): 169-176.
[6] Podlubny, I. (1998) Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications. Elsevier, Amsterdam.
[7] 纪玉德. 关于分数阶系统的稳定性与反馈控制研究[D]: [博士学位论文]. 石家庄: 河北师范大学, 2016.
[8] Ben Makhlouf, A. (2022) On the Stability of Caputo Fractional-Order Systems: A Survey. In: Naifar, O. and Ben Ma-khlouf, A. Eds., Fractional Order Systems—Control Theory and Applications. Studies in Systems, Decision and Control, Vol. 364, Springer, Cham, 1-8.
https://doi.org/10.1007/978-3-030-71446-8_1
[9] Aguila-Camacho, N., Duarte-Mermoud, M.A. and Gallegos, J.A. (2014) Lyapunov Functions for Fractional Order Systems. Communications in Nonlinear Science and Numerical Simulation, 19, 2951-2957.
https://doi.org/10.1016/j.cnsns.2014.01.022
[10] Baleanu, D., Machado, J.A.T. and Luo, A.C. (2011) Fractional Dynamics and Control. Springer Science & Business Media, Berlin.
https://doi.org/10.1007/978-1-4614-0457-6