基于.一阶倒立摆的matlab仿真实验 下载本文

.

[A,B,C,D]=tf2ss(num,den); Q=ctrb(A,B)

rank(Q) %系统能控性分析

Q =

1.0000 -0.0882 27.9247 0 1.0000 -0.0882 0 0 1.0000

ans = 3

由得到的rank(Q)的值可知,原系统的能控性矩阵为3,所以我们可知原系统是不能控的。 m=obsv(A,C) rank(m) m =

. .

.

2.3565 0 0 -0.2078 65.7876 5.4423 65.8059 -0.3580 -0.4798 ans = 3

由得到的rank(m)的值可知,原系统的能观性矩阵为3,所以我们可知原系统是能观的。 eig(A) ans =

5.2810 -5.2864 -0.0827

由eig(A)的值可知系统是不稳定的。

. .

.

4 系统设计

4.1极点配置与控制器的设计

采用极点配置法设计多输出的倒立摆系统的控制方案。可以用完全状态反馈来解决,控制摆杆和小车的位置。

设计状态反馈阵时,要使系统的极点设计成两个主导极点和两个非主导极点,用二阶系统的分析方法确定参数。

根据系统性能要求: 最大超调量10%,调节时间为 1s

-??运用超调量计算公式:?%??w?5.66 得到??0.707n

1-?2 t?3

s??n得到两个主导极点为:P1=-4+4.33j P2=-4-4.33j 选取两个非主导极点:

P3=-20 P4=-20

根据MATLAB求取状态增益矩阵,程序如下: a=[0 1 0 0

0 -0.0883167 0.629317 0 0 0 0 1

0 -0.23655 27.8285 0]; b=[0

0.883167 0

2.35655];

. .

.

p=[-4+4.33j -6-4.33j -20 -20]; k=acker(a,b,p) k =

-740.4267 -247.3685 646.9576 113.8866

4.2系统仿真:

根据状态空间表达式建立一阶倒立摆SimuLink仿真图,如下:

. .