过程控制系统仿真实验指导资料

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

2、.0 时的单位阶跃响应曲线；2) 0.8不变时，n分别为 2, 5, 8, 10 时的单位阶跃响应曲线。实验二 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 来描述，则对象的传递函数可以近似为：LsC(s)Ke LsR(s)Ts 1利用延时时间L 、放大系数K 和时间常数T，根据表一中的公式确定K ， Ti 和的值。p表一 Ziegler-Nichols 整定法控制器类型比例度/%积分时间Ti微分时间PK*LT0PIK*L

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

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

7、数整定。c曲线衰减法衰减曲线法根据衰减频率特性整定控制器参数。先把控制系统中调节器参数置成纯比例作用（Ti,0 ） ，使系统投入运行，再把比例度从大到小逐渐调小，直到出现4:1衰减过程曲线。此时比例度为4:1 ，衰减比例度为s，上升时间为t r，两个相邻波峰间的时Ts，称为 4:1 衰减振荡周期。根据 s， tr ， Ts，使用表三的经验公式可以计算出调节器的各个整定参数值。表三 临界比例度法整定控制器参数控制器类型比例度/%积分时间Ti微分时间Ps0PI1.2 s2tr 或 0.5Ts0PID0.8 s1.2tr或 0.3Ts0.4tr或 0.1Ts按“先 P 后 I 最后 D ”的操作程序

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

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

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

11、例、积分、微分增益。实验三 串级控制指导内容；例一：串级与单回路控制对比仿真某隧道窑炉系统，考虑烧成带温度为主变量，主副对象的传递函数Go1， Go2分别为：1Go1 (s)，o1 (30s 1)(3s 1)燃烧室温度为副变量的串级控制系统，Go2 (s)1(10s 1)(s 1)21主副控制器的传递函数Gc1， Gc2分别为：Gc1(s) Kc1(1 T11s)，Gc2(s) Kc2。 试分别采用单回路控制和串级控制设计主副PID 控制器的参数，并给出整定后系统的阶跃响应的特性响应曲线和阶跃扰动的响应曲线，并说明不同控制方案对系统的影响。解：串级控制设计是一个反复调整测试的过程，使用Simu

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

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

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

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

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

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

18、出串级控制系统的方框图及相同控制对象下的单回路控制系统方框图。（2） 用Simulink画出上述两个系统的仿真框图（3） 选用PID调节器，整定主副控制器的参数，使该串级控制系统性能良好，并绘制相应的单位阶跃响应曲线。（4） 比较单回路控制系统及串级控制系统在相同的副扰动下的单位阶跃响应曲线，并说明原因。例一：假设系统从动量传递函数为单闭环比例控制系统仿真综合实例。实验四 比值控制G(s)315s 1e 5t由控制理论知，开环稳定性分析是系统校正的前提。系统稳定性的分析可利用Bode 图MA TLAB Bode 图绘制程序(M-dile )如下：clear allclose allT=15;K

19、0=3;tao=5;num=K0;den=T,1;G=tf(num,den, inputdelay ,tao);margin(G)执行该程序得系统的Bode 图如图所示，可见系统是稳定的。稳定裕量为5.05dB，对应1.8。Bode Diagram)Bd(edutingaMGm = 5.05 dB (at 0.352 rad/sec) , P m = 55.5 deg (at 0.189 rad/sec) 10 0- 10- 20- 30- 40 0)ged(esah-720-1440-2160-288010-310-210-110 010 1Frequency (rad/sec)根据工程整定

20、的论述，选择PI 形式的控制器，即G(s) K KI 。ps本处采用稳定边界法整定系统。先让KI =0，调整Kp使系统等幅振荡(由稳定性分析K p =1.8 附近时系统震荡），即使系统处于临界稳定状态。K此 时 的 震 荡 周 期 为Tcr1 9s ， 比 例 系 数 为 Kpcr 1.88 ， 则 K pcr 0. 82，p 2. 2KI0.820.88Tcr0.05s。系统 Simulink 框图如下所示其中的 PID 控制器结构如下整定后从动闭环系统的单位阶跃响应如图所示：可见系统有约25%30%的超调量，在比值控制中应进一步调整使之处于振荡与不振荡调节 Kp 0.3, KI 0.02

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

22、在例一中延时是影响系统稳定性的关键环节。假设延时时间常数变化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% 变化时仍能正常工作，系

23、统的鲁棒性较强。图 a 系统仿真框图图 b 延时选择模块统封装结构作业题目：在例一中如系统传递函数为图 c 延时参数的设置图 d 仿真结果G(s) 3 e 4s，其他参数不变，试对其进行单闭环比15s 1值控制系统仿真分析，并讨论G(s)3 e 4s 分母中“ 15”变化15s 110% 时控制系统的鲁棒性。实验五 解耦控制系统指导内容：对角阵解耦控制仿真综合实例纯原料量与含水量是影响混凝土快干性和强度的两个因素。系统输入控制量为纯原料量和含水量，系统输出量为混凝土的快干性和强度，采用对角阵解耦方法对该系统进行控制仿真。解：设某双输入双输出系统，初步选择输入x1、 x2分别对应输出y1、 y2

24、。系统输入、输出之间的传递关系为：0.53s 1 X1 (s)0.3 X2(s)5s 11）11Y1(s)7s 1Y2 (s)311s 1 求系统相对增益以及系统耦合分析k11 k1211 0.5由式（ 1 ）得系统静态放大系数矩阵为k21 k223 0.3p11 p12k11 k1211 0.5即系统的第一放大系数矩阵为：Pp21 p22k21 k223 0.30.69 0.31系统的相对增益矩阵为：。0.31 0.69由相对增益矩阵可以看出，控制系统输入、输出的配对选择是正确的；通道间存在较强的相互耦合，应对系统进行解耦分析。系统的输入、输出结构如下图所示。 确定解耦调节器 根据解耦数学公式求解对角矩阵，即Gp11(s) Gp12(s)1Gp11(s)Gp22(s)Gp12(s)Gp22(s)Gp21(s) Gp22(s)G11(s)Gp22(s) Gp12(s)Gp21(s) Gp11(s)Gp21(s) Gp11(s)Gp22(s)1108.9s2 46.2s 3.311.55s2 2.7s 0.15222161.4s2 64.2s 4.8 495s2 264s 33108.9s2 46.2s 3.3采用对角矩阵解耦后，系统的结构如下图所示。解耦前后系统的Simulink 阶跃仿真框图及结果如下：1）不存在耦合时的仿真框图和结果图 a 不存在耦合时