一阶倒立摆含观测器的状态反馈控制系统综合与设计专业特制

1、广 西 大 学 实 验 报 告实验内容：一阶倒立摆含观测器的状态反馈控制系统综合与设计 张凯强 (0902100202) 、毛世博（0902100110）、李季（0902100206） 一、 实验时间：2013.3.8二、 实验地点：综合楼702三、 实验目的1、 理解并掌握线性状态反馈控制的原理和方法；2、 理解并掌握线性观测器的设计方法；3、 练习控制性能比较与评估的方法。四、实验设备与软件：1、倒立摆试验台2、matlab软件五、实验原理1、被控对象模型及其线性化根据牛顿定律建立系统垂直和水平方向的动力学方程，计及u=f，得 (1) (2)保留低阶项，项，忽略微小的高次项，在竖直位置处进

2、行线性化。由(1)(2)得 (3) (4)令，输入为，则状态方程为 （5）代入参数，忽略摩擦得 （6）该状态方程输入是加速度，输出是小车位置和摆杆角度。2、 时不变线性连续系统的状态反馈控制与观测器对时不变线性连续系统 以系统状态为反馈变量产生控制 这种控制方式称为状态反馈控制，但状态作为系统内部变量，一般很难直接测出，为此引入状态观测器。全维状态观测器的动态方程为 若输出矩阵c为满秩时，可设计较简单的降维状态观测器，其最小维数为n-m（n代表状态个数，m代表输出个数）。6、 实验内容1、 状态反馈及极点配置（1） 能控性检查：输入代码：clear;a = 0 1 0 0; 0 0 0 0;

3、0 0 0 1; 0 0 29.4 0;b = 0 1 0 3'c = 1 0 0 0; 0 1 0 0;d = 0 0'uc = ctrb(a,b);rank(uc)输出：ans = 4系统能控性矩阵满秩，即系统状态完全能控。（2） 系统极点配置选取系统主导极点：，闭环非主导极点距虚轴的距离为主导极点的5倍以上，则取：，输入代码：clear;a = 0 1 0 0; 0 0 0 0; 0 0 0 1; 0 0 29.4 0;b = 0 1 0 3' p = -10-0.0001*j,-10+0.0001*j,-2-2*sqrt(3)*j,-2+2*sqrt(3)*j;

4、k = place(a,b,p)输出：k = -54.4218 -24.4898 93.2739 16.1633（3） 极点配置系统仿真根据系统空间表达式，搭建模型。 仿真波形如图从仿真结果可以看出，小车最终稳定，小车速度，摆杆角度，角速度最终都稳定在0位置，小车位置超调5%，调整时间2s，基本符合控制要求。2、 采用状态观测器的状态反馈系统设计（1） 闭环观测器极点配置 <1> 判断可观性输入代码： a = 0 1 0 0;0 0 0 0;0 0 0 1;0 0 29.4 0;b = 0;1;0;3;c = 1 0 0 0;0 0 1 0;d = 0;sys = ss(a,b,c

5、,d);observe_matrix = obsv(a,c);rank_of_obsv = rank(observe_matrix)输出： rank_of_obsv = 4系统完全可观。输出矩阵c的秩为2，所以降维观测器的最小维数为4-2=2。 <2> 设定降维观测器的期望极点观测器特征值的选取一般是状态反馈配置极点2-3倍，所以选取状态观测器为-5，-5。输入代码：r = 0 1 0 0;0 0 0 1;p = c;r;invp = inv(p);p = -5;-5;<3> 求取等价系统的模型输入代码：aa = p*a*invpa11 = aa(1:2,1:2);a1

6、2 = aa(1:2,3:4);a21 = aa(3:4,1:2);a22 = aa(3:4,3:4);bb = p*bb1 = bb(1:2);b2 = bb(3:4);cc = c*invp输出：aa = 0 0 1.0000 0 0 0 0 1.0000 0 0 0 0 0 29.4000 0 0bb = 0 0 1 3cc = 1 0 0 0 0 1 0 0<4> 求取矩阵l输入代码：syms ssystem_eq = expand(s-p(1)*(s-p(2) syms l_1 l_2 l_3 l_4syms sl = l_1 0;0 l_4;eq = collect(d

7、et(s*eye(2)-(a22-l*a12),s) 输出：system_eq = s2 + 10*s + 25 eq = s2 + (l_1 + l_4)*s + l_1*l_4选取l= ll = 5 0;0 5;<5>求取降维观测器的动态方程 输入代码：aw = (a22 - ll*a12)bu = (b2 - ll*b1)by = (a21 - ll*a11) + (a22 - ll*a12)*llcw = invp(1:4,3:4)dy = invp(1:4,1:2)+invp(1:4,3:4)*ll输出：aw = -5 0 0 -5bu = 1 3by = -25.0000 0 0 4.4000cw = 0 0 1 0 0 0 0 1dy = 1 0 5 0 0 1 0 5（2） 系统仿真 仿真波形如图与不带观测器的