过程控制系统仿真试验指导

1、过程控制系统 Matlab/Simulink仿真实验实验一过程控制系统建模i实验二PID控制2实验三串级控制6实验四比值控制13实验五解耦控制系统19实验一 过程控制系统建模指导内容：(略)作业题目一：常见的工业过程动态特性的类型有哪几种？通常的模型都有哪些？在 Simulink中建立 相应模型，并求单位阶跃响应曲线。作业题目二：W 2某二阶系统的模型为 G(s)= n2 ,二阶系统的性能主要取决于匚，B nS 2 nS n两个参数。试利用 Simulink仿真两个参数的变化对二阶系统输出响应的影响，加深对二阶 系统的理解，分别进行下列仿真：(1) Bn =2不变时，；分别为0.1,0.8,

2、1.0, 2.0时的单位阶跃响应曲线；(2) 5 =0.8不变时，Bn分别为2, 5, 8, 10时的单位阶跃响应曲线。31实验二PID控制指导内容：PID控制器的参数整定是控制系统设计的核心内容，它根据被控过程的特征确定PID控制器的比例系数、积分时间和微分时间。PID控制器参数整定的方法很多，概括起来有两大类：(1) 理论计算整定法主要依据系统的数学模型，经过理论计算确定控制器参数。这种方法所得到的计算数据 未必可以直接使用，还必须通过工程实际进行调整和修改。(2) 工程整定方法主要有Ziegler-Nichols整定法、临界比例度法、衰减曲线法。这三种方法各有特点，其 共同点都是通过实验

3、，然后按照工程实验公式对控制器参数进行整定。但无论采用哪一种方法所得到的控制器参数，都需要在实际运行中进行最后调整与完善。工程整定法的基本特点是：不需要事先知道过程的数学模型，直接在过程控制系统中进行现场整定；方法简单，计算简便，易于掌握。a. Ziegler-Nichols 整定法Ziegler-Nichols整定法是一种基于频域设计PID控制器的方法。基于频域的参数整定是需要考虑模型的，首先需要辨识出一个能较好反映被控对象频域特性的二阶模型。根据这样的模型，结合给定的性能指标可推导出公式，而后用于PID参数的整定。基于频域的设计方法在一定程度上回避了精确的系统建模，而且有较为明确的物理意

4、义，比常规的 PID控制可适应的场合更多。目前已经有一些基于频域设计PID控制器的方法，如Ziegler-Nichols法、Cohen-Coon法等。Ziegler-Nichols法是最常用的整定 PID参数的 方法。如果单位阶跃响应曲线看起来是一条S形的曲线，则可用此法，否则不能用。S形曲线用延时时间L和时间常数T来描述，则对象的传递函数可以近似为：C(s) KeLsR(s) Ts 1利用延时时间L、放大系数K和时间常数T,根据表一中的公式确定Kp , T和T的值。表一 Ziegler-Nichols 整定法控制器类型比例度5 / %积分时间t微分时间TPK* LTO00PIK* L1.1T

5、L0.30PIDK* L 0.85T2.2L0.5Lb. 临界比例度法临界比例度法适用于已知对象传递函数的场合，在闭环的控制系统里，将调节器置于纯比例作用下，从大到小逐渐改变调节器的比例度，得到等幅振荡的过渡过程。此时的比例度称为临界比例度5k ,两个相邻波峰间的时间间隔，称为临界振荡周期Tk。采用临界比例度法时，系统产生临界振荡的条件是系统的阶数是3阶或3阶以上。临界比例度法的步骤如下：(1) 将调节器的积分时间 Ti置于最大(T =* )，微分时间置零，比例度 8适当，平衡操作一段时间，把系统投入自动运行；(2) 将比例度&逐渐减小，得到等幅振荡过程，记下临界比例度8k和临界震荡周

6、期Tk的值；(3) 根据孩和Tk值，采用表二的经验公式，计算出调节器的各个参数，即5、I和e的值。表二临界比例度法整定控制器参数控制器类型比例度6 / %积分时间t微分时间TP2&k00PI2.2祝0.88览0PID17k0.50Tk0.125Tk按“先P后I最后D”的操作程序将调节器整定参数调到计算数值上。若还不够满意， 可再作一步调整。临界比例度法注意事项：(1) 有的过程控制系统，临界比例度很小，调节阀不是全关就是全开，对工业生产不 利；(2) 有的过程控制系统，当调节器比例度8调到最小刻度值时，系统仍然不产生等幅振荡，对此，将最小刻度的比例度作为临界比例度dk进行调节器参数整定

7、。c. 曲线衰减法衰减曲线法根据衰减频率特性整定控制器参数。先把控制系统中调节器参数置成纯比例作用(Ti =°°,7=0 )，使系统投入运行，再把比例度从大到小逐渐调小，直到出现4:1衰减过程曲线。此时比例度为4:1，衰减比例度为Ss,上升时间为tr,两个相邻波峰间的时间间隔为Ts,称为4:1衰减振荡周期。根据§ , L, Ts ,使用表三的经验公式可以计算出调节器的各个整定参数值。表三临界比例度法整定控制器参数控制器类型比例度5 / %积分时间T微分时间TP们CO0PI12s2tr 或 0.5Ts0PID0.861.2tr 或 0.3Ts0.4tr 或 0.1T

8、s按“先P后I最后D”的操作程序将调节器整定参数调到计算数值上。若还不够满意,可再作一步调整。衰减曲线法的注意事项：(1) 对于反应较快的系统，要认定 4:1衰减曲线和读出Ts比较困难，此时，可以认为记录指针来回摆动两次就达到稳定是4:1衰减过程。(2) 在生产过程中，负荷变化会影响过程特性。当负荷变化较大时，必须重新整定调 节器参数。(3) 若认为4:1衰减太慢，可采用10:1衰减过程。对于10:1衰减曲线整定调节器参数 的步骤与上述完全相同，仅仅是计算公式不同。作业题目：建立如下所示Simulink仿真系统图。利用Simulink仿真软件进行如下实验:1. 建立如图所示的实验 Simuli

9、nk原理图。2. 双击原理图中的 PID模块，出现参数设定对话框，将PID控制器的积分增益和微分增益改为0,使其具有比例调节功能，对系统进行纯比例控制。3. 进行仿真，观测系统的响应曲线，分析系统性能；然后调整比例增益，观察响应曲 线的变化，分析系统性能的变化。4. 重复（步骤2,3）,将控制器的功能改为比例微分控制，观测系统的响应曲线，分析 比例微分的作用。5. 重复（步骤2,3）,将控制器的功能改为比例积分控制，观测系统的响应曲线，分析 比例积分的作用。6. 重复（步骤2,3）,将控制器的功能改为比例积分微分控制，观测系统的响应曲线， 分析比例积分微分的作用。7. 将PID控制器的积分微分

10、增益改为0,对系统进行纯比例控制。不断修改比例增益，使系统输出的过度过程曲线的衰减比n=4，记下此时的比例增益值。8. 修改比例增益，使系统输出的过度过程曲线的衰减比n=2，记下此时的比例增益值。9. 修改比例增益，使系统输出呈现临界振荡波形，记下此时的比例增益。10. 将PID控制器的比例、积分增益进行修改，对系统进行比例积分控制。不断修改比例、积分增益，使系统输出的过渡过程曲线的衰减比n=2,4,10,记下此时比例和积分增益。11. 将PID控制器的比例、积分、微分增益进行修改，对系统进行比例积分控制。不断修改比例、积分、微分增益，使系统输出的过度过程曲线的衰减比n=2,4,10,记下此时

11、比例、积分、微分增益。实验三串级控制指导内容；例一：串级与单回路控制对比仿真某隧道窑炉系统，主副对象的传递函数 G°1,考虑烧成带温度为主变量,Go2分别为：1燃烧室温度为副变量的串级控制系统，其Goi (s)(30s 1)(3s 1)Go2 (S)202(10s 1)(s 1)2主副控制器的传递函数Gc1 , Gc2 分另1J为:1Gc1(s) =Kc1(1+ ) , Gc2(s) = Kc2 试分别T1s°采用单回路控制和串级控制设计主副PID控制器的参数，并给出整定后系统的阶跃响应的特性响应曲线和阶跃扰动的响应曲线，并说明不同控制方案对系统的影响。解：串级控制设计是一

12、个反复调整测试的过程，使用 Simulink能大大简化这一过程。根据题意，首先建立如图的Simulink模型。图中采用单回路控制的 Simulink图，其中，q1为一次扰动，取阶跃信号； q2为二次扰 动，取阶跃信号；G°2为副对象；G°1为主对象；r为系统输入，取阶跃信号； c为系统输出， 它连接到示波器上，可以方便地观测输出。图中的PID C1为单回路PID控制器，它是按照 PID原理建立的Simulink中的子模块， 其内部结构如下图所示:PID控制器模块子系统的参数设置如下图:经过不断的试验，当输入比例系数为3.7，积分系数为38,微分系数为0时，系统阶跃响应达到比

13、较满意的效果，系统阶跃响应如下图：采用这套PID参数时，二次扰动作用下，系统的输出响应如下图:采用这套PID参数时，一次扰动作用下，系统输出响应如下图:综合以上各图可以看出采用单回路控制，系统的阶跃响应达到要求时，系统对一次,次扰动的抑制效果不是很好。下面考虑采用串级控制时的情况，下图为串级控制时的Simulink模型图：图中，q1为一次扰动，取阶跃信号；q2为二次扰动，取阶跃信号；PID C1为主控制器， 采用PID控制，PID C2为副控制器，采用 PID控制；Go2为副对象；Goi为主对象；r为系 统输入，取阶跃信号；c为系统输出，它连接到示波器上，可以方便地观测输出。经过不断试验，当P

14、ID C1为主控制器输入比例系数为 8.4,积分系数为12.8,微分系数 为0时；当PID C2为主控制器输入比例系数为 10,积分系数为0,微分系数为0时；系统 阶跃响应达到比较满意的效果，系统阶跃响应如下图所示：采用这套PID参数时，二次扰动作用下，系统的输出响应如下图:采用这套PID参数时，一次扰动作用下，系统的输出响应如下图:综合以上各图可以看出， 采用串级控制，系统的阶跃响应达到要求时， 系统对一次扰动,二次扰动的抑制也能达到很好的效果。综合单回路控制和串级控制的情况，系统的控制性能对比如下表所示。表一系统采用单回路控制和串级控制的对比控制品质指标单回路控制Kc1=3.7, Tc1=

15、38串级控制Kc1=8.4, Tc1=12.8, Kc2=10衰减率0.750.75调节时间18060残偏差00二次阶跃扰动下的系统短期 最大偏差0.270.013一次阶跃扰动下的系统短期 最大偏差0.340.13从表中可以看出系统的动态过程改善更为明显，可见对二次扰动的最大动态偏差可以减小约20倍，对一次扰动的最大动态偏差也可以减小约2.6倍，系统的调节时间提高了3倍。PID参数的整定例二：串级控制的参数整定仿真某隧道窑炉系统，考虑烧成带温度为主变量，统中，其主副对象的传递函数Goi, Go2分别为:1燃烧室温度为副变量所构成的串级控制系Goi(s)=，(30s 1)(3s 1)G°

16、;2(s) = 2(10s 1)(s 1)2试采用逐次逼近法整定 PID控制器的参数，并给出整定后系统的阶跃响应曲线。解：按照逐次逼近法的步骤1先主回路开环，按单回路方法整定副控制器，建立的Simulink框图如下图:不断地头验，当 K c2=10时，副回路阶跃响应如下图:从图中可以大致地看出，此衰减比约为4:1,因此此时可进入步骤2,即主回路闭环,取Kc2=10整定主控制器，此时的 Simulink框图如下：不断地实验，调节主调节器，当Kci=12 ,主回路的阶跃响应如下图所示:由图可知此时的Tsi=16，按步骤3,将主回路闭环的条件下，重新整定副控制器参数, 根据衰减曲线法，取 Kc1=1

17、0, Ki=8,此时系统主回路阶跃响应曲线如下图所示：由图可知系统的阶跃响应效果不理想，超调太大，需要进行精调。按照步骤4反复实验当Kci=8.4 , Ki =12.8, Kc2=10时，系统阶跃响应如下图所示：由图可知系统的阶跃响应效果比较理想，此时整定的主副回路的参数比较合理。当然，还可以调整不同的参数组合，以取得满意的输出响应曲线。这是一个反复的实验过程，利用Simulink这一便捷的工具具有很大的优势。作业题目：串级控制系统仿真。已知某串级控制系统的主副对象的传递函数Goi , Go2分别为:Goi(s)=1一，G°2=，副回路干扰通道的传递函数为：Gd2(s)= 1。100

18、s 110s 1s 20s 1(1) 画出串级控制系统的方框图及相同控制对象下的单回路控制系统方框图。(2) 用Simulink画出上述两个系统的仿真框图(3) 选用PID调节器，整定主副控制器的参数，使该串级控制系统性能良好，并绘制相应的单位阶跃响应曲线。(4) 比较单回路控制系统及串级控制系统在相同的副扰动下的单位阶跃响应曲线，并说明原因。实验四比值控制指导内容：单闭环比例控制系统仿真综合实例。3例一：假设系统从动量传递函数为G(s) =e项，设计该从动对象的单闭环比值15s - 1控制系统。解：本题的基本步骤如下： 分析从动量无调节器的开环系统稳定性。由控制理论知，开环稳定性分析是系统校

19、正的前提。系统稳定性的分析可利用Bode图进行，编制 MA TLAB Bode图绘制程序(M-dile )如下：clear allclose allT=15;K0=3;tao=5;num=K0;den=T,1;G=tf(num,den, 'inputdelay' ,tao);margin(G)执行该程序得系统的Bode图如图所示，可见系统是稳定的。稳定裕量为5.05dB,对应增益为1.8。Bode DiagramGm = 5.05 dB (at 0.352 rad/sec) , P m = 55.5 deg (at 0.189 rad/sec) 10.XJBacAdMrnaa

20、M0-10-20-30-40 0mearesanp-720-1440-2160-288010-310-210-1100101Frequency (rad/sec) 选择从动量控制器形式及整定其参数。根据工程整定的论述，选择 PI形式的控制器，即G(s) = Kp +垃 。s本处采用稳定边界法整定系统。先让KI =0,调整Kp使系统等幅振荡(由稳定性分析图知在Kp=1.8附近时系统震荡），即使系统处于临界稳定状态。 pK_此时的震荡周期为Tcr =19s ,比例系数为Kpcr出1.88 ,则 10. 82, p 2.2K =M=0.05s。0.88Tcr系统Simulink框图如下所示其中的PI

21、D控制器结构如下蛀写5 翼蛀4 Link; imT i t led/PTDEil« Edit Zi ew S.imulati m F£.rn» - - lools Kelp Q21 I |io O|»LV I 用因囹Readylow整定后从动闭环系统的单位阶跃响应如图所示:可见系统有约25%30%的超调量，在比值控制中应进一步调整使之处于振荡与不振荡 的边界。调节Kp =0.3,K| =0.02时，系统响应图如下所示，基本达到了振荡临界要求。 系统过程仿真。单闭环比值控制过程相当于从动量变化的随动控制过程。假定主动量由一常值10加幅度为0.3的随机扰动构

22、成，从动量受均值为0、方差为1的随机干扰。主动量和从动量的比值根据工艺要求及测量仪表假定为3.系统的控制过程 Simulink仿真框图如图所示。其中控制常量及随机扰动采用封装形式。主动控制量的封装结构如下:运行结果如下所示（图中曲线从上往下分别为从动量跟踪结果、主动量给定值和随机干扰）:3530 ,-25 ,-20 ,15 ,-10-5 -0 -. r 二 /-5 11111050100150200250300可见除初始时间延时外，从动量较好地跟随主动量变化而变化，并且基本维持比值3,有效地克服了主动量和从动量的扰动。例二：（续例一）单闭环比值控制系统鲁棒性。由控制理论知，在例一中延时是影响系

23、统稳定性的关键环节。假设延时时间常数变化±10%，即延时时间为 4.55.5,分析系统鲁棒性。解：系统仿真框图如图 a所示，图b为延时选择模块 Subsystem的展开图，图c为延时 常数设定方法，改变数为4.5,4.6,4.7, 4.8 - 5.5共11个值。经过运行后在工作空间绘图（使用语句：plot (tout, simout); hold on; grid on)即可见到图 d的仿真结果。分析图d仿真结果并与例一的图进行比较可见，随着延时环节的变化， 从动量跟随主动量的规律有较小变化，但并未改变系统稳定性及精度，说明系统在延时发生 ±10%变化时仍 能正常工作，系统

24、的鲁棒性较强。图a系统仿真框图图b延时选择模块统封装结构Olock Faxaseters: Tiaitsport DelaysTmnsgrt Deli?App? FP<tclf3d d«l«r lo h« incut slEn«, wccqscr if achieved vdren the de Lar is 芒更thax the s-irfe-ulat ion step si zcS'a. EELnet crsTina daI nr：Ini Iiftl input:F9m < itl butf4¥ sisti1024Fad

25、e crier (for 1 inuriza.tion):F Blitfft feedihrouElt Qf 1皿过上 durlnE llihealzaUon邸： I C an calK«lpAjply图c延时参数的设置图d仿真结果作业题目：在例一中如系统传递函数为G(s)二伟、es ,其他参数不变，试对其进行单闭环比值控制系统仿真分析，并讨论 G(s)=房£/分母中“ 15”变化±10%时控制系统的鲁 棒性。实验五解耦控制系统指导内容:对角阵解耦控制仿真综合实例纯原料量与含水量是影响混凝土快干性和强度的两个因素。系统输入控制量为纯原料量和含水量，系统输出量为混凝

26、土的快干性和强度，采用对角阵解耦方法对该系统进行控制仿真。解：设某双输入双输出系统，初步选择输入xi、X2分别对应输出y、V2。系统输入、输出之间的传递关系为：110.5Y（s）_ 7s+1 3s + 1 Xi（s）1 k（s）| 30.3!x2（s）_：.11s + 1 5s + 1 j求系统相对增益以及系统耦合分析（1）由式（1）得系统静态放大系数矩阵为灯 k12 _ 11 0.5虹原口-3 0.3_即系统的第一放大系数矩阵为:-P11P12kn k1211 0.5P = I=|=|p21p22 -_k21 k22 -3 0.3_系统的相对增益矩阵为:0.69 0.31 A =110.31

27、 0.69由相对增益矩阵可以看出，控制系统输入、输出的配对选择是正确的； 通道间存在较强的相互耦合，应对系统进行解耦分析。系统的输入、输出结构如下图所示。确定解耦调节器根据解耦数学公式求解对角矩阵，即Gpii(s) Gpi2(s)'Gp2i(s) Gp22(s)1&L4矿 +64岛+ 4,皿9厂+4心f +W.2j + 4.8lg.9f + 仙一21+3一3IM.lc +fr4 24 x161,4/4 6*1 五j+ 4 81Gpii(S)Gp22(S)-Gpi2(S)Gp22(S)Gii(s)Gp22(s) Gpi2(s)Gp2i(s) |-Gpii(s)Gp2i(s) Gpii(s)Gp22(s)1i08.9s2 46.2s 3.3 -ii.55s2 - 2.7s -0.i5222i6i.4s +64.2s+4.8495s2 +264s + 33i08.9s2 十46.2s + 3.3采用对角矩阵解耦后，系统的结构如下图所示。解耦前后系统的Simulink阶跃仿真框图及结果如下:D不存在耦合时的仿真框图和结果图a不存在耦合时的仿真框图（上）和结果（下）