自动控制设计(自动控制原理课程设计)

自动控制原理课程设计本课程设计的目的着重于自动控制基本原理与设计方法的综合实际应用。主要内容包括：古典自动控制理论（PID）设计、现代控制理论状态观测器的设计、自动控制MATLAB仿真。通过本课程设计的实践，掌握自动控制理论工程设计的基本方法和工具。1内容某生产过程设备如图1所示，由液容为C1和C2的两个液箱组成，图中Q为稳态液体流量)/(3sm，iQ为液箱A输入水流量对稳态值的微小变化)/(3sm，1Q为液箱A到液箱B流量对稳态值的微小变化)/(3sm，2Q为液箱B输出水流量对稳态值的微小变化)/(3sm，1h为液箱A的液位稳态值)(m，1h为液箱A液面高度对其稳态值的微小变化)(m，2h为液箱B的液位稳态值)(m，2h为液箱B液面高度对其稳态值的微小变化)(m，21,RR分别为A，B两液槽的出水管液阻))//((3smm。设u为调节阀开度)(2m。已知液箱A液位不可直接测量但可观，液箱B液位可直接测量。h20+△h2Q10+△Q1Qi+△Qih10+△h1Q20+△Q2AB要求图1某生产过程示意图1.建立上述系统的数学模型；2.对模型特性进行分析，时域指标计算，绘出bode,乃示图，阶跃反应曲线3.对B容器的液位分别设计：P，PI，PD，PID控制器进行控制；4.对原系统进行极点配置，将极点配置在－1＋j和－1－j；（极点可以不一样）5.设计一观测器，对液箱A的液位进行观测（此处可以不带极点配置）；6.如果要实现液位h2的控制，可采用什么方法，怎么更加有效？试之。用MATLAB对上述设计分别进行仿真。（提示：流量Q=液位h/液阻R，液箱的液容为液箱的横断面积，液阻R=液面差变化h/流量变化Q。）2双容液位对象的数学模型的建立及MATLAB仿真过程一、对系统数学建模如图一所示，被控参数2h的动态方程可由下面几个关系式导出：液箱A：dthdCQQi111液箱B：dthdCQQ2221111/QhR222/QhRuKQui消去中间变量，可得：uKhdthdTTdthdTT222122221)(式中，21,CC——两液槽的容量系数21,RR——两液槽的出水端阻力111CRT——第一个容积的时间常数222CRT——第二个容积的时间常数2RKKu_双容对象的放大系数其传递函数为：1)()()()(212212STTSTTKSUSHSG二．对模型特性进行分析，绘出bode,奈氏图，阶跃反应曲线当输入为阶跃响应时的Matlab仿真：令T1=T2=6；K=1112361)()()(22SSSUSHSG2)16(1S单位阶跃响应的MATLAB程序：num1=[1];den1=[36121];G1=tf(num1,den1);figure(1);step(G1);xlabel('时间(sec)');ylabel('输出响应');title('二阶系统单位阶跃响应');step(G1,100);运行结果如下：阶跃反应曲线：图1c(∞)=1;c(tp)=1;tp=45.5s;td=10s;ts=45.5s;最大超调量：δ(tp)=[c(tp)-c(∞)]/c(∞)*100%=0%稳态误差分析：开环传递函数112361)()()(22SSSUSHSG，稳态误差1sse；用MATLAB绘制的奈氏图如下图2所示，其程序如下：nyquist([1],conv([61],[61]))图2在工程实践中，一般希望正相角裕度r为45o~60o,增益裕度Kg10dB,即Kg3。当系统为单位负反馈时的Ｂode图：用MATLAB绘制的奈氏图如下图3所示，其程序如下：sys=tf([1],conv([61],[61]));margin(sys);figure图3三：对B容器的液位分别设计：P，PI，PD，PID控制器进行控制PID控制的原理和特点(1)P控制：取P=9;I=0;D=0;（2）PI控制：P=6，I=0.4，D=0；（3）PD控制：P=9，I=0，D=5；（4）PID控制：P=5，I=0.3，D=4；四．系统极点配置在-1+j;-1-j根据传递函数1)()()()(212212STTSTTKSUSHSG得微分方程uKhdthdTTdthdTT222122221)(令222212,,xhxhxh得状态方程uKxxTTTTTTxxKuxTTTTxTTxxx01101212121212122121121221＋即：输出：1xy极点配置：令K=1;T1=T2=2;uxxxxuxxxxx10125.0102212121221＋即：用MATLAB确定状态反馈矩阵K，使得系统闭环极点配置在（-1+j,-1-j）,程序如下：A=[01;-0.25-1];B=[0;1];P=[-1+j;-1-j];K=place(A,B,P)运行结果为K=1.75001.0000仿真：仿真图五．设计一观测器，对液箱A的液位进行观测建立状态观测器：根据传递函数1)()()()(212212STTSTTKSUSHSG得微分方程uKhdthdTTdthdTT222122221)(令222212,,xhxhxh得状态方程uKxxTTTTTTxxKuxTTTTxTTxxx01101212121212122121121221＋即：输出：1xy全维观测器的建立：令21ggG，得1112122212121212121212122121212120101011111detdet()11ggAGCTTTTggTTTTTTTTgTTgIAGCgTTgTTTTTTTT期望特征式：2(a2))(()(22为设定值）aaaaf对比1式和2式，得22112212112aTTgagTTTT得21212221121TTTTagaTTg21ggG所以全维状态观测器得方程是1112221212ˆˆ()()010ˆ()1xAGCxGybuggxyuTTggKTTTT本实验中，需观测的状态为水箱A溶液的液位1h，建立数学模型R1=R2=1;c1=c2=1;;2;12122111hxhxdthdhhdthdhu21]01[0121110121xxyuxxxx21]11[xxy令状态观测器的极点为（-6-j,-6+j）设计此给定系统状态观测器的MATLAB程序如下A=[-10;1-1];B=[10];C=[11];A1=A';B1=C';C1=B';P=[-6-j-6+j];K=acker(A1,B1,P);G=K'运行结果为G=26-16仿真：仿真图：六、如果要实现液位h2的控制，可采用什么方法，怎么更加有效？试之前馈反馈控制方法这种调节系统中要直接测量干扰量的变化，液位h2作为反馈量，流量Q作为前馈量，可以克服流量Q干扰量的偏差，同时可以加快控制的速度，使调节更加及时有效。