线性系统理论大作业

  • 格式:doc
  • 大小:406.50 KB
  • 文档页数:12

线性系统理论大作业(总11页)

-CAL-FENGHAI.-(YICAI)-Company One1

-CAL-本页仅作为文档封面,使用请直接删除 目录

题目一 ............................................. 错误!未指定书签。

(一)状态反馈加积分器校正的输出反馈系统设计 ....... 错误!未指定书签。

(1)建立被控对象的状态空间模型,并判断系统性质 ...... 错误!未指定书签。

(2)状态反馈增益矩阵和积分增益常数的设计 ............ 错误!未指定书签。

(3)全维观测器设计 .................................. 错误!未指定书签。

(4)如何在闭环调速系统中增加限流环节 ................ 错误!未指定书签。

(二)二次型最优全状态反馈控制和按负载扰动前馈补偿的复合控制系统设计错误!未指定书签。

(1)线性二次型最优全状态反馈设计 .................... 错误!未指定书签。

(2)降维观测器设计 .................................. 错误!未指定书签。

题目二 ............................................. 错误!未指定书签。

(1)判断系统是否存在最优控制律 ...................... 错误!未指定书签。

(2)非零给定点的最优控制设计和仿真分析 .............................. 错误!未指定书签。

(3)权矩阵的各权值对动态性能影响分析 .................................. 错误!未指定书签。

题目一

(一)状态反馈加积分器校正的输出反馈系统设计

(1)建立被控对象的状态空间模型,并判断系统性质

1)画出与题目对应的模拟结构图,如图1所示:

图1 原始系统结构图

取状态变量为1x=n,2x=dI,3x=du,控制输入u=cu

将已知参数代人并设输出y=n=1x,得被控对象的状态空间表达式为

其中,237500039.768011=-3.696-17.85727.05600-588.235100TelalalasCGDCARTTRTT, 000=023529.41ssBKT,2375-30.4880=000GDE,100C

2)检查被控系统的结构性质

判断系统能控性、能观性、稳定性

程序如下:

A=[0 39.768 0;-3.696 -17.857 27.056;0 0 -588.235];

B=[0;0;23529.41];C=[1 0 0];

Qc=ctrb(A,B);

Qo=obsv(A,C);

L=length(A);

if rank(Qc)==L

disp('系统是状态完全能控');

else

disp('系统是状态不完全能控');

end

if rank(Qo)==L

disp('系统是状态完全能观');

else

disp('系统是状态不完全能观');

end

disp(eig(A))%利用A的特征值判断系统稳定性

运行结果:

系统是状态完全能控

系统是状态完全能观

1.0e+02 *

-0.0893 + 0.0820i

-0.0893 - 0.0820i

-5.8823 + 0.0000i

由于矩阵A全部特征值均具有负实部,因此系统渐近稳定。

原系统

设负载转矩为0,输入为阶跃信号,系统simulink仿真如下:

图2 原始开环系统结构框图

图3 原始开环系统仿真

分析:由系统仿真图可以看出,调节时间大于0.5s,不满足性能指标。

(2)状态反馈增益矩阵和积分增益常数的设计

由于原系统能控,可以使用状态反馈。为满足设计指标,采用状态反馈加积分器校正的输出反馈系统。因增广系统能控,故可采用线性状态反馈控制律12uKxKw。

将闭环系统极点配置到复平面左半开平面的任意期望位置且可消除阶跃扰动及阶跃参考输入作用下的稳态误差。

式中:1111213=KKKK,123=Txxxx,wvy,v为系统参考输入。 由经典控制理论,闭环极点为21,21nnj的欠阻尼二阶线性定常系统的超调量及调节时间为21%=100%e,3.5snt。

系统需满足%%≤8,0.5sst≤,计算可得0.6266≥,n≥11.1714,取=0.7,=n12,设计指标的期望闭环主导极点对为*1,28.48.57j。选择2个期望的闭环非主导极点离虚轴为主导极点5倍以上,取*3,450,据期望闭环极点,采用MATLAB极点配置函数可求出增广系统状态反馈增益阵,程序如下:

A=[0 39.768 0;-3.696 -17.857 27.056;0 0 -588.235];

B=[0;0;23529.41];C=[1 0 0];

Az=[A [0;0;0];-C 0];Bz=[B;0];Cz=[C 0];

P=[-8.4+j*8.57;-8.4-j*8.57;-50;-50];

Km=acker(Az,Bz,P);

K=[Km(1,1),Km(1,2),Km(1,3),-Km(1,4)]

运行程序可得:

系统simulink仿真如下:

图4 状态反馈加积分器校正系统结构框图

图5 状态反馈加积分器校正系统仿真

由图可知,超调量%(1.0420-1)/14.2%8%,调节时间0.5sst,满足要求。

图6 加负载扰动后系统状态反馈加积分器校正系统仿真

0时刻扰动,最终系统稳定在1,因此系统稳态误差为0。

(3)全维观测器设计

由于系统能观,可以使用状态观测器。

237500039.768011=-3.696-17.85727.05600-588.235100TelalalasCGDCARTTRTT,

000=023529.41ssBKT,100C,

新系统的特征根为:-61,27.807212.1936j,基于通常选择观测器的响应速度比所考虑的状态反馈闭环系统快2-5倍这一经验规则,取观测器期望极点为:-150,-60,-70。

应用MATLAB极点配置函数求解新系统全维观测器,程序如下:

A=[0 39.768 0;-3.696 -17.857 27.056;-38.8235 -88.9412 -98.8233];

B=[0;0;23529.41];C=[1 0 0];

P=[-150;-60;-70]; Gt=acker(A',C',P);%求对偶系统的状态反馈增益阵G

G=Gt';%求系统的观测器偏差反馈增益矩阵G

运行程序可得:

163.3197G=8.1911129.8454。

带观测器的状态反馈加积分调节系统仿真结构图如下:

图7 带观测器的状态反馈加积分器校正系统结构框图

图8 系统加全维观测器波形图

图9 全维观测器波形图

由仿真图可知,系统的稳态误差为0,动态误差满足超调量%%≤8,调节时间0.5sst≤的要求。

状态估计误差收敛速度与状态观测器极点的配置有关。一般而言状态观测器极点在复平面的左半开平面距离虚轴距离越远,则估计误差收敛速度越快。但是,观测器响应速度过快会产生大量噪声,影响系统的正常工作故不宜取值过大。综合工程实际出发,一般取为比状态反馈闭环系统快2—5倍。

(4)如何在闭环调速系统中增加限流环节

从加快启动电动机的角度来看,闭环调速系统应允许有较大的启动电流,而造成堵转的故障消失后,系统电流应能自动恢复正常。所以常规的熔断器或过流继电器在这里均不能作为限流保护措施。因为它们是通过切断电路来保护设备的,虽然能起到保护作用,但故障消失后,系统无法自动恢复正常。为了充分利用设备的过流能力,又保证设备的安全运行,电流截止负反馈则可以限制电流的大小。

电流截止负反馈的作用是:当电枢电流大于某一截止值时,电流负反馈起作用,限制电流不能过大。当电枢电流小于截止值时,电流反馈被截止,对系统的稳态运行不产生影响。电动机启动时,因为电流截止负反馈作用,从而限制启动电流。正常工作时,电流截止负反馈作用很小。电动机发生堵转时,由于电流截止负反馈的作用,使Ud大大下降,因而使Ia不致过大。允许的堵转电流一般为电动机额定电流的2~2.5倍。系统工作在额定值时,由于电流截止负反馈起作用,从而保证系统设备的安全。电流截止负反馈如图所示:

图10 电流截止负反馈结构图

(二)二次型最优全状态反馈控制和按负载扰动前馈补偿的复合控制系统设计

(1)线性二次型最优全状态反馈设计

1)判断系统是否存在最优控制律

要使系统阶跃响应具有良好的动、静态特性,可按非零给定点的最优控制律设计,即()()gutKxtu,由于输入维数和输出维数相等,所以1()()(0)cwutKxtWy。

由于系统完全能控,因此,最优控制()ut存在。

最优控制性能指标为:01()()()()2TTtJxtQxtutRutdt,其中Q为状态加权系数矩阵,R为控制加权系数矩阵。

2)非零给定点的最优控制设计和仿真分析 由图可知,系统输出响应发散,可引入最优控制1()()(0)cwutKxtWy。选取设112233001000000=01000100qQqq,R=1。

程序如下:

A=[0 39.768 0;-3.696 -17.857 27.056;0 0 -588.235];

B=[0;0;23529.41];

C=[1 0 0];

D=0;

R=1;

Q=[100 0 0;0 1 0;0 0 1];

K=lqr(A,B,Q,R);

ac=A-B*K;

W=inv((-C/(A-B*K))*B);

bc=B*W;

cc=C;

dc=D;

step(ac,bc,cc,dc);

grid

运行结果如下:K=9.86354.87020.9809,110.0009cW。

图11 非零给定点最优控制系统单位阶跃响应

3)权矩阵的各权值对动态性能影响分析

a)固定控制加权系数矩阵R=1,且另22q、33q都为1,11q取不同值时,研究非零给定点的最优控制仿真曲线。程序如下:

a_color=['r','g','b','y','c', 'm','k'];

A=[0 39.768 0;-3.696 -17.857 27.056;0 0 -588.235];

B=[0;0;23529.41];

C=[1 0 0];

D=0;

R=1;

syms Q q11;

N=[1 100,200,500,1000,10000];

syms i K;

for i=1:6

q11=N(i);

Q=[q11 0 0;0 1 0;0 0 1]

K=lqr(A,B,Q,R);

ac=A-B*K;

W=inv(-C/(A-B*K)*B);

bc=B*W;

cc=C;

dc=D;

sys(i)=ss(ac,bc,cc,dc);