自適應反步控制
文章目錄
- 自適應反步控制
- 1. 基本思想
- A. 第一步
- B. 第二步
1. 基本思想
基于傳統反步法,考慮了系統方程中以線性形式出現的未知參數。核心思想包括參數估計率和控制率。
考慮二階系統:
{x˙1=x2+φ1T(x1)θx˙2=u+φ2T(x1,x2)θ(1)\begin{cases} \dot{x}_1 = x_2 + \varphi_1^{T}(x_1)\theta \\ \dot{x}_2 = u + \varphi_2^{T}(x_1, x_2)\theta \end{cases} \tag{1} {x˙1?=x2?+φ1T?(x1?)θx˙2?=u+φ2T?(x1?,x2?)θ?(1)
其中 θ\thetaθ 是未知的常數參數(列向量),φ1(x1)\varphi_1(x_1)φ1?(x1?) 和 φ2(x1,x2)\varphi_2(x_1,x_2)φ2?(x1?,x2?) 為已知的非線性函數。uuu 是控制輸入,xrx_rxr? 是 x1x_1x1? 期望跟蹤到的軌跡。
設計目標有二:確保閉環系統的所有狀態保持有界;讓 x1x_1x1? 漸進跟蹤到 xrx_rxr?。
假設:參考軌跡 xrx_rxr? 及其導數有界并可預測;φ1(x1)\varphi_1(x_1)φ1?(x1?) 和 φ2(x1,x2)\varphi_2(x_1,x_2)φ2?(x1?,x2?) 是已知且連續可微的;參數不確定性是線性形式出現的。這些假設的目的是保證在設計自適應律時,Lyapunov 導函數可以被有效約束,并且更新率是可實現的。
如果參數 θ\thetaθ 已知,可以使用靜態積分分反步法設計虛擬控制率 α1\alpha_1α1?,但此處的參數 θ\thetaθ 并非已知,所以需要設計合適的 θ\thetaθ 更新規則以保證閉環穩定。
A. 第一步
定義跟蹤誤差
z1=x1?xr(2)z_1 = x_1 - x_r \tag{2} z1?=x1??xr?(2)
z2=x2?α1?x˙r(3)z_2 = x_2 - \alpha_1 - \dot{x}_r \tag{3} z2?=x2??α1??x˙r?(3)
z2z_2z2? 表示“實際 x2x_2x2? 與理想目標值 α1+x˙r\alpha_1 + \dot{x}_rα1?+x˙r? 之間的偏差”。
對(2)求導數,并將(1)、(3)帶入:
z˙1=x˙1?x˙r=x2+φ1T(x1)θ?x2+z2+α1=z2+α1+φ1Tθ(4)\dot{z}_1 = \dot{x}_1 - \dot{x}_r = x_2 + \varphi_1^{T}(x_1)\theta - x_2 + z_2 + \alpha_1 = z_2 + \alpha_1 + \varphi_1^{T}\theta \tag{4} z˙1?=x˙1??x˙r?=x2?+φ1T?(x1?)θ?x2?+z2?+α1?=z2?+α1?+φ1T?θ(4)
其中 θ\thetaθ 是未知的,不能直接使用已知參數的反步法。設計 θ^1\hat{\theta}_1θ^1? 用以近似 θ\thetaθ;同時設計一個參數更新率來在線調整 θ^1\hat{\theta}_1θ^1?。
設計 Lyapunov 候選函數:
V1=12z12+12θ~1TΓ?1θ~1(5)V_1 = \frac{1}{2} z_1^2 + \frac{1}{2} \tilde{\theta}_1^{T} \Gamma^{-1} \tilde{\theta}_1 \tag{5} V1?=21?z12?+21?θ~1T?Γ?1θ~1?(5)
其中 Γ>0\Gamma > 0Γ>0 是參數更新率的自適應增益矩陣,θ~1=θ?θ^1\tilde{\theta}_1 = \theta - \hat{\theta}_1θ~1?=θ?θ^1? 表示參數估計誤差。
對其求導:
V˙1=z1z˙1+θ~1TΓ?1θ~˙1(6)\dot{V}_1 = z_1 \dot{z}_1 + \tilde{\theta}_1^{T} \Gamma^{-1} \dot{\tilde{\theta}}_1 \tag{6} V˙1?=z1?z˙1?+θ~1T?Γ?1θ~˙1?(6)
將(4)帶入:
V˙1=z1(z2+α1+φ1Tθ)+θ~1TΓ?1θ~˙1(7)\dot{V}_1 = z_1(z_2 + \alpha_1 + \varphi_1^{T}\theta) + \tilde{\theta}_1^{T} \Gamma^{-1} \dot{\tilde{\theta}}_1 \tag{7} V˙1?=z1?(z2?+α1?+φ1T?θ)+θ~1T?Γ?1θ~˙1?(7)
展開:
V˙1=z1z2+z1α1+z1φ1Tθ^1+z1φ1Tθ~1?θ~1TΓ?1θ^˙1(8)\dot{V}_1 = z_1 z_2 + z_1 \alpha_1 + z_1 \varphi_1^{T} \hat{\theta}_1 + z_1 \varphi_1^{T} \tilde{\theta}_1 - \tilde{\theta}_1^{T} \Gamma^{-1} \dot{\hat{\theta}}_1 \tag{8} V˙1?=z1?z2?+z1?α1?+z1?φ1T?θ^1?+z1?φ1T?θ~1??θ~1T?Γ?1θ^˙1?(8)
由于 φ1Tθ~=θ~1Tφ1\varphi_1^{T}\tilde{\theta} = \tilde{\theta}_1^{T}\varphi_1φ1T?θ~=θ~1T?φ1?,可得:
V˙1=z1z2+z1α1+z1φ1Tθ^1+θ~1T(φ1z1?Γ?1θ^˙1)(9)\dot{V}_1 = z_1 z_2 + z_1 \alpha_1 + z_1 \varphi_1^{T} \hat{\theta}_1 + \tilde{\theta}_1^{T}(\varphi_1 z_1 - \Gamma^{-1} \dot{\hat{\theta}}_1) \tag{9} V˙1?=z1?z2?+z1?α1?+z1?φ1T?θ^1?+θ~1T?(φ1?z1??Γ?1θ^˙1?)(9)
定義虛擬控制 α1\alpha_1α1? 和參數更新率:
α1=?c1z1?φ1Tθ^1\alpha_1 = -c_1 z_1 - \varphi_1^{T} \hat{\theta}_1 α1?=?c1?z1??φ1T?θ^1?
θ^˙1=Γφ1z1(10)\dot{\hat{\theta}}_1 = \Gamma \varphi_1 z_1 \tag{10} θ^˙1?=Γφ1?z1?(10)
帶入可得:
V˙1=?c1z12+z1z2(11)\dot{V}_1 = -c_1 z_1^2 + z_1 z_2 \tag{11} V˙1?=?c1?z12?+z1?z2?(11)
B. 第二步
設計實際控制輸入 uuu,使 z2→0z_2 \to 0z2?→0。
由(3)可得:
z˙2=x˙2?α˙1?x¨r(12)\dot{z}_2 = \dot{x}_2 - \dot{\alpha}_1 - \ddot{x}_r \tag{12} z˙2?=x˙2??α˙1??x¨r?(12)
注意 α1\alpha_1α1? 實際上是 x1x_1x1?、θ^1\hat{\theta}_1θ^1?、xrx_rxr? 的函數,故采用復合函數求導方式。
選取新的 Lyapunov 函數:
V2=V1+12z22(14)V_2 = V_1 + \frac{1}{2} z_2^2 \tag{14} V2?=V1?+21?z22?(14)
求導:
V˙2=?c1z12+z2(z1+z˙2)(15)\dot{V}_2 = -c_1 z_1^2 + z_2(z_1 + \dot{z}_2) \tag{15} V˙2?=?c1?z12?+z2?(z1?+z˙2?)(15)
為消除未知參數項 θ\thetaθ,設計:
u=?z1?c2z2+?α1?x1x2?(φ2??α1?x1φ1)Tθ^2+?α1?θ^1Γφ1z1+?α1?xrx˙r+x¨r(18)u = -z_1 - c_2 z_2 + \frac{\partial \alpha_1}{\partial x_1} x_2 - \left( \varphi_2 - \frac{\partial \alpha_1}{\partial x_1} \varphi_1 \right)^T \hat{\theta}_2+ \frac{\partial \alpha_1}{\partial \hat{\theta}_1} \Gamma \varphi_1 z_1 + \frac{\partial \alpha_1}{\partial x_r} \dot{x}_r + \ddot{x}_r \tag{18} u=?z1??c2?z2?+?x1??α1??x2??(φ2???x1??α1??φ1?)Tθ^2?+?θ^1??α1??Γφ1?z1?+?xr??α1??x˙r?+x¨r?(18)
此時(13)變為:
z˙2=?z1?c2z2+(φ2??α1?x1φ1)T(θ?θ^2)(19)\dot{z}_2 = -z_1 - c_2 z_2 + \left( \varphi_2 - \frac{\partial \alpha_1}{\partial x_1} \varphi_1 \right)^T (\theta - \hat{\theta}_2) \tag{19} z˙2?=?z1??c2?z2?+(φ2???x1??α1??φ1?)T(θ?θ^2?)(19)
定義新的 Lyapunov 函數:
V2=V1+12z22+12θ~2TΓ?1θ~2(20)V_2 = V_1 + \frac{1}{2} z_2^2 + \frac{1}{2} \tilde{\theta}_2^{T} \Gamma^{-1} \tilde{\theta}_2 \tag{20} V2?=V1?+21?z22?+21?θ~2T?Γ?1θ~2?(20)
求導:
V˙2=?c1z12?c2z22?θ~2TΓ?1θ^˙2+θ~2T(φ2??α1?x1φ1)z2(21)\dot{V}_2 = -c_1 z_1^2 - c_2 z_2^2 - \tilde{\theta}_2^{T} \Gamma^{-1} \dot{\hat{\theta}}_2 + \tilde{\theta}_2^{T} \left( \varphi_2 - \frac{\partial \alpha_1}{\partial x_1} \varphi_1 \right) z_2 \tag{21} V˙2?=?c1?z12??c2?z22??θ~2T?Γ?1θ^˙2?+θ~2T?(φ2???x1??α1??φ1?)z2?(21)
令:
θ^˙2=Γ(φ2??α1?x1φ1)z2(22)\dot{\hat{\theta}}_2 = \Gamma \left( \varphi_2 - \frac{\partial \alpha_1}{\partial x_1} \varphi_1 \right) z_2 \tag{22} θ^˙2?=Γ(φ2???x1??α1??φ1?)z2?(22)
最終得到:
V˙2=?c1z12?c2z22(23)\dot{V}_2 = -c_1 z_1^2 - c_2 z_2^2 \tag{23} V˙2?=?c1?z12??c2?z22?(23)
可保證系統全局漸近穩定。