单级倒立摆系统的极点配置与状态观测器设计

单级倒立摆系统的极点配置与状态观测器设计

14122156 杨郁佳

（1）倒立摆的运动方程并将其线性化

选取小车的位移z ，及其速度z 、摆的角位置θ及其角速度θ

作为状态变量，即T

x z z θθ??=???

? 则系统的状态空间模型为 01000100000010()1000mg M M x u M m g Ml Ml x ????????????-????=+????????+-???????????? []1000y x

= 设M=2kg ，m=0.2kg ，g=9.81m/2

s ，则单级倒立摆系统的状态方程为

(101001010)

01020.500013030

011040.54x x x x u x x x x ??????????????????-????????=+????????????????-????????????

[]12100034x x y x x ????

??=??????

（2）状态反馈系统的极点配置。

首先，使用MATLAB ，判断系统的能控性矩阵是否为满秩。

MATLAB 程序如下：

A=[0 1 0 0; 0 0 -1 0; 0 0 0 1; 0 0 11 0];

B=[0; 0.5; 0; -0.5];

C=[1 0 0 0];

D=0;

rct=rank(ctrb(A,B))

[z,p,k]=ss2zp(A,B,C,D)

MATLAB 程序执行结果如下：

系统能控，系统的极点为

1=0

λ

2=0

λ

3=3.3166

λ

4=-3.3166

λ

可以通过状态反馈来任意配置极点，将极点配置在

1=-3

λ*

2=-4

λ*

3=-5

λ*

4=-6

λ*

MATLAB程序如下：

A=[0 1 0 0; 0 0 -1 0; 0 0 0 1; 0 0 11 0]; B=[0; 0.5; 0; -0.5];

P=[-3 -4 -5 -6];

K=place(A,B,P)

MATLAB程序执行结果如下：

因此，求出状态反馈矩阵为

K=[-72.0 -68.4 -332.0 -104.4]

采用MATLAB/Simulink构造单级倒立摆状态反馈控制系统的仿真模型。

首先，在MATLAB的Command Window中输入各个矩阵的值，并且在模型中的积分器中设置非零初值：2。运行仿真程序，显示仿真曲线，如下。

处。

仿真结果表明倒立摆的杆子与数值方向的偏角从初值2，经过控制稳定在=0

（3）状态观测器实现状态反馈极点配置。

MATLAB程序如下：

A=[0 1 0 0; 0 0 -1 0; 0 0 0 1; 0 0 11 0];

B=[0; 0.5; 0; -0.5];

C=[1 0 0 0];

rob=rank(obsv(A,C))

MATLAB程序执行结果如下：

Rob=4说明系统能观，可以设计状态观测器。取状态观测器的特征值为-3，-4，-5，-6 MATLAB 程序如下：

A=[0 1 0 0; 0 0 -1 0; 0 0 0 1; 0 0 11 0];

A1=A ’;

C=[1 0 0 0];

C1=C ’;

P=[-3 -4 -5 -6];

H1=place(A1,C1,P);

H=H1’

MATLAB 程序执行结果如下：

状态观测器矩阵[]18 130 540 1790T

H =--

采用MATLAB/Simulink 构造具有状态观测器的单级倒立摆状态反馈控制系统的仿真模型。

首先，在MATLAB的Command Window中输入各个矩阵的值，并且在模型中的积分器中设置非零初值：2。运行仿真程序，显示仿真曲线，如下。

对比两个仿真图，可以发现加上状态观测器对单级倒立摆的控制效果基本上无影响。