状态观测器的设计与实践
发布日期:2021-07-01 04:03:05 浏览次数:2 分类:技术文章

本文共 1845 字,大约阅读时间需要 6 分钟。

文章目录

声明:本文内容整理自在B站的课程内容,版权归原作者所有,博主仅作学习笔记记录。

1. 系统建模

在这里插入图片描述

假设只有位置可测、速度不可测。

z 1 = x z 2 = x ˙ , 该 量 不 可 测 u = F z_1=x\\ z_2=\dot x,该量不可测\\ u=F z1=xz2=x˙u=F

建立系统状态方程为:

在这里插入图片描述

2. 观测器设计

观测器(observer): 根据系统的输入和输出来估计系统的状态。

下面来介绍观测器(本文介绍的即大名鼎鼎的龙伯格观测器)的设计步骤:

x ˙ = A x + B u ① y = C x + D u ② \dot x = Ax + Bu\qquad\qquad\qquad①\\ y=Cx+Du\qquad\qquad\qquad② x˙=Ax+Buy=Cx+Du

x ^ \hat{x} x^为估计值, y ^ \hat{y} y^为估计的输出,则

x ^ ˙ = A x ^ + B u + L ( y − y ^ ) ③ y ^ = C x ^ + D u ④ \dot{ \hat x} = A\hat x + Bu + L(y - \hat y)\qquad③\\ \hat y=C \hat x+Du\qquad\qquad\qquad④ x^˙=Ax^+Bu+L(yy^)y^=Cx^+Du

将④代入③,得到:

x ^ ˙ = A x ^ + B u + L y − L ( C x ^ + D u ) \dot {\hat x} = A\hat x + Bu + Ly - L(C\hat x + Du) x^˙=Ax^+Bu+LyL(Cx^+Du)

即:

x ^ ˙ = ( A − L C ) x ^ + ( B − L D ) u + L y ⑤ \dot {\hat x} = (A-LC)\hat x + (B-LD)u + Ly\qquad⑤ x^˙=(ALC)x^+(BLD)u+Ly
上式即为观测器核心。

令①-⑤得,

x ˙ − x ^ ˙ = [ A x + B u ] − [ ( A − L C ) x ^ + ( B − L D ) u + L y ] \dot x - \dot {\hat x} = [Ax + Bu] - [(A - LC)\hat x + (B - LD)u + Ly] x˙x^˙=[Ax+Bu][(ALC)x^+(BLD)u+Ly]
将②代入上式,得:
x ˙ − x ^ ˙ = A x + B u − ( A − L C ) x ^ − ( B − L D ) u − L C x − L D u \dot x - \dot {\hat x} = Ax + Bu - (A - LC)\hat x - (B - LD)u - LCx - LDu x˙x^˙=Ax+Bu(ALC)x^(BLD)uLCxLDu
继续整理,得:
x ˙ − x ^ ˙ = ( A − L C ) ( x − x ^ ) \dot x - \dot {\hat x} = (A - LC)(x - \hat x) x˙x^˙=(ALC)(xx^)

即误差系统方程为:

e ˙ = ( A − L C ) e \dot e = (A - LC)e e˙=(ALC)e
其中, e = x − x ^ e=x - \hat x e=xx^

我们的目标是要保证 e → 0 e→0 e0,也就是说 ( A − L C ) (A-LC) (ALC)的特征值要小于0才行!即:

R e [ e i g ( A − L C ) ] < 0 Re[eig(A-LC)]<0 Re[eig(ALC)]<0

由此可得该观测器的参数矩阵L,观测器到此即设计完毕。

对于第1部分中的系统,观测器设计如下:


在这里插入图片描述


3. 实践仿真

由以上模型,在Simulink中搭建模型如下:

在这里插入图片描述

输出结果如下:

在这里插入图片描述在这里插入图片描述
从上图可见,观测值跟系统几乎重合,观测器有效!

给定观测器 z ^ 1 \hat z_1 z^1一个初值,运行结果如下:

在这里插入图片描述

在这里插入图片描述

致谢

感谢所发布的所有视频,让博主重新拾起了扔到脑后多少年的东西,严谨的态度,手把手的带你手撕公式,这才是真正的doctor,向他看齐!与君共勉!!!


参考文献:

  • https://zhuanlan.zhihu.com/p/108466027

转载地址:https://miracle.blog.csdn.net/article/details/111039778 如侵犯您的版权,请留言回复原文章的地址,我们会给您删除此文章,给您带来不便请您谅解!

上一篇:各路大神对于观测器的文章总结【持续更新】
下一篇:MATLAB - plotFreqSpectrum函数的使用

发表评论

最新留言

第一次来,支持一个
[***.219.124.196]2024年04月11日 05时09分10秒