模型参考自适应控制系统结构

模型参考自适应控制系统的典型结构如下图所示,系统由参考模型、控制器和自适应率组成。

  • 控制器包括被控对象的前馈控制器和反馈控制器,可以根据自适应率进行调整;
  • 参考模型实际上是一种理想控制系统,其输出代表了期望的性能,对调节系统的特性要求,如超调量、阻尼时间等等;
  • 自适应率用来消除被控对象输出和参考模型期望输出的误差,改变控制器参数或者生成辅助输入。

img

基于 L y a p u n o v Lyapunov Lyapunov稳定性原理的二阶 S I S O SISO SISO系统自适应控制器设计

考虑如下二阶 S I S O SISO SISO系统模型:

x ¨ + θ T Φ ( x , x ˙ ) = u \ddot{x}+\theta ^T \Phi(x,\dot{x})=u x¨+θTΦ(x,x˙)=u

其中 θ = [ θ 1 θ 2 θ 3 ] T ∈ R 3 \theta = [\theta _1\quad \theta _2\quad \theta _3]^T \in \mathbb{R}^3 θ=[θ1θ2θ3]TR3是未知常数向量, Φ ( x , x ˙ ) = [ Φ 1 ( x , x ˙ ) Φ 2 ( x , x ˙ ) Φ 3 ( x , x ˙ ) ] T ∈ R 3 \Phi(x,\dot{x}) = [\Phi _1(x,\dot{x}) \quad \Phi _2(x,\dot{x}) \quad \Phi _3(x,\dot{x})]^T \in \mathbb{R}^3 Φ(x,x˙)=[Φ1(x,x˙)Φ2(x,x˙)Φ3(x,x˙)]TR3为已知的有节基函数, θ T Φ ( x , x ˙ ) \theta ^T \Phi(x,\dot{x}) θTΦ(x,x˙)代表系统的结构匹配不确定性。

选取状态变量 x 1 = x , x 2 = x ˙ x_1=x,x_2=\dot{x} x1=x,x2=x˙,则有:

{ x ˙ 1 = x 2 x ˙ 2 = − θ T Φ ( x , x ˙ ) + u \begin{cases} \dot{x}_1=x_2 \\ \dot{x}_2=-\theta ^T \Phi(x,\dot{x})+u \end{cases} {x˙1=x2x˙2=θTΦ(x,x˙)+u

设计如下控制器:

u = − ( λ + 1 ) x ˙ − λ x + θ ^ T Φ ( x , x ˙ ) u=- (\lambda+1) \dot{x}-\lambda x +\hat{\theta}^T \Phi(x,\dot{x}) u=(λ+1)x˙λx+θ^TΦ(x,x˙)

= − ( λ + 1 ) x 2 − λ x 1 + θ ^ T Φ ( x 1 , x 2 ) =- (\lambda+1) x_2-\lambda x_1 +\hat{\theta}^T \Phi(x_1,x_2) =(λ+1)x2λx1+θ^TΦ(x1,x2)

其中 θ ^ \hat{\theta} θ^是参数 θ \theta θ的估计值, λ \lambda λ为控制器输入参数且应该大于等于0。

接下来设计估计参数 θ ^ \hat{\theta} θ^的更新率,这里结合 L y a p u n o v Lyapunov Lyapunov稳定性进行设计。假设参数估计的误差 θ ~ = θ ^ − θ \tilde{\theta}=\hat{\theta}-\theta θ~=θ^θ,因为 θ \theta θ为常数,所以 θ ~ ˙ = θ ^ ˙ \dot{ \tilde{\theta} } = \dot{\hat{\theta}} θ~˙=θ^˙,定义 L y a p u n o v Lyapunov Lyapunov 函数:

V ( x 1 , x 2 , θ ~ ) = 1 2 ( λ x 1 + x 2 ) 2 + 1 2 θ ~ T θ ~ V(x_1,x_2,\tilde{\theta})=\frac{1}{2} (\lambda x_1+x_2 )^2+\frac{1}{2}\tilde{\theta}^T \tilde{\theta} V(x1,x2,θ~)=21(λx1+x2)2+21θ~Tθ~

求导:

V ˙ ( x 1 , x 2 , θ ~ ) = ( λ x 1 + x 2 ) ( λ x ˙ 1 + x ˙ 2 ) + θ ~ T θ ^ ˙ \dot{V}(x_1,x_2,\tilde{\theta})=(\lambda x_1 +x_2)(\lambda \dot{x}_1 +\dot{x}_2)+\tilde{\theta}^T\dot{\hat{\theta}} V˙(x1,x2,θ~)=(λx1+x2)(λx˙1+x˙2)+θ~Tθ^˙

= ( λ x 1 + x 2 ) ( λ x 2 − θ T Φ ( x 1 , x 2 ) + u ) + θ ~ T θ ^ ˙ =(\lambda x_1 +x_2)(\lambda x_2 -\theta ^T \Phi(x_1,x_2)+u)+\tilde{\theta}^T\dot{\hat{\theta}} =(λx1+x2)(λx2θTΦ(x1,x2)+u)+θ~Tθ^˙

= ( λ x 1 + x 2 ) ( λ x 2 − θ T Φ ( x 1 , x 2 ) − ( λ + 1 ) x 2 − λ x 1 + θ ^ T Φ ( x 1 , x 2 ) ) + θ ~ T θ ^ ˙ =(\lambda x_1 + x_2 )(\lambda x_2 - \theta ^T\Phi(x_1,x_2)- (\lambda+1) x_2-\lambda x_1 +\hat{\theta}^T \Phi(x_1,x_2))+\tilde{\theta}^T\dot{\hat{\theta}} =(λx1+x2)(λx2θTΦ(x1,x2)(λ+1)x2λx1+θ^TΦ(x1,x2))+θ~Tθ^˙

= − ( λ x 1 + x 2 ) 2 + ( λ x 1 + x 2 ) θ ~ T Φ ( x 1 , x 2 ) + θ ~ T θ ^ ˙ = -(\lambda x_1 +x_2)^2 + (\lambda x_1 +x_2) \tilde{\theta}^T \Phi(x_1,x_2)+\tilde{\theta}^T\dot{\hat{\theta}} =(λx1+x2)2+(λx1+x2)θ~TΦ(x1,x2)+θ~Tθ^˙

= − ( λ x 1 + x 2 ) 2 + θ ~ T ( ( λ x 1 + x 2 ) Φ ( x 1 , x 2 ) + θ ^ ˙ ) =-(\lambda x_1 +x_2)^2 +\tilde{\theta}^T ((\lambda x_1 +x_2)\Phi (x_1,x_2)+\dot{\hat{\theta}} ) =(λx1+x2)2+θ~T((λx1+x2)Φ(x1,x2)+θ^˙)

L y a p u n o v Lyapunov Lyapunov理论: V ˙ ≤ − ρ V + C \dot{V}\leq -\rho V +C V˙ρV+C V → C 2 ρ V \rightarrow \sqrt{\frac{C}{2 \rho}} V2ρC ,为了达到系统的稳定,设计 θ ^ \hat{\theta } θ^的自适应更新率为:

θ ^ ˙ = − ( λ x 1 + x 2 ) Φ ( x 1 , x 2 ) − k θ ^ \dot{\hat{\theta}} =-(\lambda x_1 +x_2)\Phi (x_1,x_2) - k\hat{\theta} θ^˙=(λx1+x2)Φ(x1,x2)kθ^

其中 k k k为自适应率输入参数,将 θ ^ \hat{\theta } θ^的更新率带入则有:

V ˙ ( x 1 , x 2 , θ ~ ) = − ( λ x 1 + x 2 ) 2 − k θ ~ T θ ^ \dot{V}(x_1,x_2,\tilde{\theta})=-(\lambda x_1 +x_2)^2-k\tilde{\theta}^T\hat{\theta} V˙(x1,x2,θ~)=(λx1+x2)2kθ~Tθ^

= − ( λ x 1 + x 2 ) 2 − k θ ~ T ( θ ~ + θ ) =-(\lambda x_1 +x_2)^2-k\tilde{\theta}^T(\tilde{\theta}+\theta) =(λx1+x2)2kθ~T(θ~+θ)

= − ( λ x 1 + x 2 ) 2 − k θ ~ T θ ~ − k θ ~ T θ = -(\lambda x_1 + x_2)^2 - k\tilde{\theta}^T \tilde{\theta} - k\tilde{\theta} ^T \theta =(λx1+x2)2kθ~Tθ~kθ~Tθ

≤ − ( λ x 1 + x 2 ) 2 − k θ ~ T θ ~ + k 2 θ ~ T θ ~ + k 2 θ T θ \leq -(\lambda x_1 + x_2)^2 - k\tilde{\theta}^T \tilde{\theta} +\frac{k}{2}\tilde{\theta}^T\tilde{\theta}+\frac{k}{2}\theta ^T \theta (λx1+x2)2kθ~Tθ~+2kθ~Tθ~+2kθTθ

≤ − ( λ x 1 + x 2 ) 2 − k 2 θ ~ T θ ~ + k 2 θ T θ \leq -(\lambda x_1 + x_2)^2 - \frac{k}{2}\tilde{\theta}^T \tilde{\theta} + \frac{k}{2}\theta ^T \theta (λx1+x2)22kθ~Tθ~+2kθTθ

由于 θ \theta θ为常数,所以按照以上方式设计的控制器 u u u和自适应率 θ ^ \hat \theta θ^满足 L y a p u n o v Lyapunov Lyapunov稳定性。

Logo

为开发者提供学习成长、分享交流、生态实践、资源工具等服务,帮助开发者快速成长。

更多推荐