《自动控制设计》由会员分享，可在线阅读，更多相关《自动控制设计（17页珍藏版）》请在金锄头文库上搜索。
1、自动控制原理课程设计本课程设计的目的着重于自动控制基本原理与设计方法的综合实际应用。主 要内容包括：古典自动控制理论(PID)设计、现代控制理论状态观测器的设计、 自动控制 MATLAB 仿真。通过本课程设计的实践，掌握自动控制理论工程设计的 基本方法和工具。I 内容某生产过程设备如图 1 所示，由液容为 C1 和 C2 的两个液箱组成，图中 Q为稳态液体流量(m3 /s)，AQ为液箱A输入水流量对稳态值的微小变化(m3 /s)，iAQ为液箱A到液箱B流量对稳态值的微小变化(m3 /s)，AQ为液箱B输出水流12量对稳态值的微小变化(m3 /s)，h为液箱A的液位稳态值(m)，Ah为液箱A液I
2、I面高度对其稳态值的微小变化(m)，h为液箱B的液位稳态值(m)，Ah为液箱B22液面高度对其稳态值的微小变化 (m)， R ,R 分别为 A， B 两液槽的出水管液阻12(m/(m3/s)。设u为调节阀开度(m2)。已知液箱A液位不可直接测量但可观，液箱B液位可直接测量。要求1. 建立上述系统的数学模型；2. 对模型特性进行分析，时域指标计算，绘出bode,乃示图，阶跃反应曲线3. 对B容器的液位分别设计：P，PI, PD, PID控制器进行控制；4. 对原系统进行极点配置，将极点配置在一1 + j和一1 j;（极点可以不一样）5. 设计一观测器，对液箱 A 的液位进行观测（此处可以不带极点
3、配置）；6. 如果要实现液位 h2 的控制，可采用什么方法，怎么更加有效？试之。用 MATLAB 对上述设计分别进行仿真。（提示：流量Q=液位h/液阻R,液箱的液容为液箱的横断面积，液阻R=液面差变化Ah/流量变化AQ。2 双容液位对象的数学模型的建立及 MATLAB 仿真过程一、对系统数学建模如图一所示，被控参数Ah的动态方程可由下面几个关系式导出:2液箱 A： AQ - AQ = C dAhii 11 dt液箱 B： AQ - AQ = C dAh2122 dtR =Ah / AQ1 1 1R = Ah / AQ2 2 2AQ 二 K Auiu 消去中间变量，可得TT d2Ah2 + (T
4、 + T )dAh + Ah = KAu1 2 dt 212 dt2式中， C ,C 两液槽的容量系数12R , R 两液槽的出水端阻力12T =RC第一个容积的时间常数1 1 1T =R C 第二个容积的时间常数2 2 2K二K R 双容对象的放大系数u2其传递函数为：G( s ) 2-AU(S) 一 TT S2 + (T + T )S +11 2 1 2二对模型特性进行分析，绘出bode，奈氏图，阶跃反应曲线 当输入为阶跃响应时的 Matlab 仿真：令 T1=T2=6 ； K=1G(S)AH (S)AU (S)136S 2 + 12S +11(6S +1)2单位阶跃响应的MATLAB程序
5、：num1=1;den1=36 12 1;G1=tf(num1，den1);figure(1);step(G1);xlabel(时间(sec)；ylabel(输出响应)；title (二阶系统单位阶跃响应);step(G1，100);运行结果如下：阶跃反应曲线：图1c)=l； c(t)=l； t= 45.5s；ppt =10s； dt =45.5s； s最大超调量：(5(t )=pc(t )- pC(8)/C(8)*100%=0%稳态误差分析：(S)_开环传递函数 G (S)-AUs-36 s 2 +12 s +1稳态误差e 1 ；ss用MATLAB绘制的奈氏图如下图2所示，其程序如下:nyq
6、uist(1,conv(6 1,6 1)r.-|*iFigure 1r xFile Edit View Insert Tools Desktop Window HelpD Q S | E| h Nyquist Diagram.8_.b.40.0.o.2o.24. _.Ko.- - -总 Aeu-EilriE-LI-0.8-0.6-0.4-0.200.20.40.60.81Real Axis图2在工程实践中，一般希望正相角裕度r为45o 60o，增益裕度K 10dB,即 gK 3。当系统为单位负反馈时的Bode图：用MATLAB绘制的奈氏图如下图3g所示，其程序如下：sys=tf(1,conv(
7、6 1,6 1);margin(sys);figureL Figure 1File Edit View Insert Tools Desktop Window HelpBod已 DiagramGm = Inf dB (at Inf rad/sec), Pm = -180 deg (at 0 rad/sec)10 10 101Frequency (rad/sec)图3三：对B容器的液位分别设计：P, PI, PD, PID控制器进行控制PID 控制的原理和特点(1)P 控制：取 P=9;I=0;D=0;(2) PI 控制：P=6, 1=0.4, D=0；3) PD 控制： P=9， I=0， D
8、=5；(4) PID 控制：P=5, 1=0.3, D=4；四系统极点配置在T+j; T-j根据传递函数G (S) =2=AU(S) TTS2 + (T + T )S +11 2 1 2 d 2 AhdAh得微分万程 TT2 + (T + T )2 + Ah = KAu1 2 dt 21 2 dt 2令 Ah = x , Ah = x , Ah = x2 1 2 2 2 2得状态万程12即：1xT T 112T + T一 T 2 XTT 212+ KuX1 X2121T + T一 2TT12X-0 _1XK2输出： y= X1极点配置：令 K=1； T1=T2=2；X = X12X = - X
9、 - 2 x + u2120即：X1X2一 0.25一1Xi +X2用MATLAB确定状态反馈矩阵K,使得系统闭环极点配置在（T+j,-l-j）,程序如 下：A=0 1；-0.25 -1；B=0；1；P=-1+j；-1-j；K=place(A,B,P)运行结果为K =1.75001.0000仿真：五设计一观测器，对液箱A的液位进行观测 建立状态观测器：根据传递函数G (S) =2=AU(S) TTS2 + (T + T )S +11 2 1 2d 2 AhdAh得微分万程 TT - 亠+ (T + T )一a + Ah = KAu1 2 dt 21 2 dt 2令 Ah = x , Ah =
10、x , Ah = x2 1 2 2 2 2得状态万程x = x12+ KuT + TTT-12 xT T 2即：2x1 x2011212x-0 _1L x K21T + T2TT12u输出：y=x1 全维观测器的建立： 令G= g1 ,得Lg 2A GC =TT12T + T1TT12g1 01=2g2TT121 g11T +TT 2 gTT 212det 九 I -(A GC )_= detTT12g1 11T + T+ tTT12y T + T =入 2 + (-TT12g )九 + 1 g121TT12期望特征式：（a为设定值）2f * （九）=（九 + a）（九 + a）二九2 + 2
11、a + a 2对比 1 式和 2 式，得T + T12 + g = 2aTT 212TT121 - g 201TT121 - g1- +-T * 2 g- 212g1 f g2本实验中，需观测的状态为水箱 A 溶液的液位 h1 建立数学模型R1=R2=1; c1=c2=1;Au -Ah =丝1 dt12 dtx1 = Ah ; x2 = Ah ;12X1-10 _x1+丁x.21-1x 20x1y = 10 cx 2x1y 二 g 2x 2令状态观测器的极点为(-6-j,-6+j)设计此给定系统状态观测器的 MATLAB 程序如下A=-1 0;1 -1;B=1 0;C=1 1;A1=A;B1=
12、C;C1=B;P=-6-j -6+j;K=acker(A1,B1,P);G=K运行结果为G =26-16仿真：oG a i n313ain220.51020304050607080901000.50010203040506070809010003 a i rr10ain4Sc:pe2仿真图:Ga inGainBIntegrator-0.5 kIntegratorl六、如果要实现液位h2的控制，可采用什么方法，怎么更加有效？试之 前馈反馈控制方法这种调节系统中要直接测量干扰量的变化，液位h2作为反馈量，流量Q作为前 馈量，可以克服流量Q干扰量的偏差，同时可以加快控制的速度，使调节更加及 时有效。Welcome !欢迎您的下载, 资料仅供参考!TT12g =1-TTa2 2T + T g，= 2a - _22G=所以全维状态观测器得方程是X = （A - GC） X + G （y） + bu
