过程控制系统仿真实验指导.doc

过程控制系统Matlab/Simulink仿真实验实验一 过程控制系统建模1实验二 PID控制2实验三 串级控制6实验四 比值控制13实验五 解耦控制系统19实验一 过程控制系统建模指导内容：（略）作业题目一：常见的工业过程动态特性的类型有哪几种？通常的模型都有哪些？在Simulink中建立相应模型，并求单位阶跃响应曲线。作业题目二：某二阶系统的模型为，二阶系统的性能主要取决于，两个参数。试利用Simulink仿真两个参数的变化对二阶系统输出响应的影响，加深对二阶系统的理解，分别进行下列仿真：（1）不变时，分别为0.1, 0.8, 1.0, 2.0时的单位阶跃响应曲线；（2）不变时，分别为2, 5, 8, 10时的单位阶跃响应曲线。实验二 PID控制指导内容：PID控制器的参数整定是控制系统设计的核心内容，它根据被控过程的特征确定PID控制器的比例系数、积分时间和微分时间。PID控制器参数整定的方法很多，概括起来有两大类：（1） 理论计算整定法主要依据系统的数学模型，经过理论计算确定控制器参数。这种方法所得到的计算数据未必可以直接使用，还必须通过工程实际进行调整和修改。（2） 工程整定方法主要有Ziegler-Nichols整定法、临界比例度法、衰减曲线法。这三种方法各有特点，其共同点都是通过实验，然后按照工程实验公式对控制器参数进行整定。但无论采用哪一种方法所得到的控制器参数，都需要在实际运行中进行最后调整与完善。工程整定法的基本特点是：不需要事先知道过程的数学模型，直接在过程控制系统中进行现场整定；方法简单，计算简便，易于掌握。a Ziegler-Nichols整定法Ziegler-Nichols整定法是一种基于频域设计PID控制器的方法。基于频域的参数整定是需要考虑模型的，首先需要辨识出一个能较好反映被控对象频域特性的二阶模型。根据这样的模型，结合给定的性能指标可推导出公式，而后用于PID参数的整定。基于频域的设计方法在一定程度上回避了精确的系统建模，而且有较为明确的物理意义，比常规的PID控制可适应的场合更多。目前已经有一些基于频域设计PID控制器的方法，如Ziegler-Nichols法、Cohen-Coon法等。Ziegler-Nichols法是最常用的整定PID参数的方法。如果单位阶跃响应曲线看起来是一条S形的曲线，则可用此法，否则不能用。S形曲线用延时时间L和时间常数T来描述，则对象的传递函数可以近似为：利用延时时间L、放大系数K和时间常数T，根据表一中的公式确定，和的值。表一 Ziegler-Nichols整定法控制器类型比例度积分时间微分时间P0PI0PID b 临界比例度法临界比例度法适用于已知对象传递函数的场合，在闭环的控制系统里，将调节器置于纯比例作用下，从大到小逐渐改变调节器的比例度，得到等幅振荡的过渡过程。此时的比例度称为临界比例度，两个相邻波峰间的时间间隔，称为临界振荡周期。采用临界比例度法时，系统产生临界振荡的条件是系统的阶数是3阶或3阶以上。临界比例度法的步骤如下：（1）将调节器的积分时间置于最大（），微分时间置零，比例度适当，平衡操作一段时间，把系统投入自动运行；（2）将比例度逐渐减小，得到等幅振荡过程，记下临界比例度和临界震荡周期的值；（3）根据和值，采用表二的经验公式，计算出调节器的各个参数，即、和的值。表二 临界比例度法整定控制器参数控制器类型比例度积分时间微分时间P0PI0PID按“先P后I最后D”的操作程序将调节器整定参数调到计算数值上。若还不够满意，可再作一步调整。临界比例度法注意事项：（1）有的过程控制系统，临界比例度很小，调节阀不是全关就是全开，对工业生产不利；（2）有的过程控制系统，当调节器比例度调到最小刻度值时，系统仍然不产生等幅振荡，对此，将最小刻度的比例度作为临界比例度进行调节器参数整定。c 曲线衰减法衰减曲线法根据衰减频率特性整定控制器参数。先把控制系统中调节器参数置成纯比例作用（），使系统投入运行，再把比例度从大到小逐渐调小，直到出现4:1衰减过程曲线。此时比例度为4:1，衰减比例度为，上升时间为，两个相邻波峰间的时间间隔为，称为4:1衰减振荡周期。根据，使用表三的经验公式可以计算出调节器的各个整定参数值。表三 临界比例度法整定控制器参数控制器类型比例度积分时间微分时间P0PI或0PID或或按“先P后I最后D”的操作程序将调节器整定参数调到计算数值上。若还不够满意，可再作一步调整。衰减曲线法的注意事项：（1）对于反应较快的系统，要认定4:1衰减曲线和读出比较困难，此时，可以认为记录指针来回摆动两次就达到稳定是4:1衰减过程。（2）在生产过程中，负荷变化会影响过程特性。当负荷变化较大时，必须重新整定调节器参数。（3）若认为4:1衰减太慢，可采用10:1衰减过程。对于10:1衰减曲线整定调节器参数的步骤与上述完全相同，仅仅是计算公式不同。作业题目：建立如下所示Simulink仿真系统图。利用Simulink仿真软件进行如下实验：1. 建立如图所示的实验Simulink原理图。2. 双击原理图中的PID模块，出现参数设定对话框，将PID控制器的积分增益和微分增益改为0，使其具有比例调节功能，对系统进行纯比例控制。3. 进行仿真，观测系统的响应曲线，分析系统性能；然后调整比例增益，观察响应曲线的变化，分析系统性能的变化。4. 重复（步骤2,3），将控制器的功能改为比例微分控制，观测系统的响应曲线，分析比例微分的作用。5. 重复（步骤2,3），将控制器的功能改为比例积分控制，观测系统的响应曲线，分析比例积分的作用。6. 重复（步骤2,3），将控制器的功能改为比例积分微分控制，观测系统的响应曲线，分析比例积分微分的作用。7. 将PID控制器的积分微分增益改为0，对系统进行纯比例控制。不断修改比例增益，使系统输出的过度过程曲线的衰减比n=4，记下此时的比例增益值。8. 修改比例增益，使系统输出的过度过程曲线的衰减比n=2，记下此时的比例增益值。9. 修改比例增益，使系统输出呈现临界振荡波形，记下此时的比例增益。10. 将PID控制器的比例、积分增益进行修改，对系统进行比例积分控制。不断修改比例、积分增益，使系统输出的过渡过程曲线的衰减比n=2,4,10，记下此时比例和积分增益。11. 将PID控制器的比例、积分、微分增益进行修改，对系统进行比例积分控制。不断修改比例、积分、微分增益，使系统输出的过度过程曲线的衰减比n=2,4,10，记下此时比例、积分、微分增益。实验三 串级控制指导内容；例一：串级与单回路控制对比仿真某隧道窑炉系统，考虑烧成带温度为主变量，燃烧室温度为副变量的串级控制系统，其主副对象的传递函数Go1，Go2分别为：， 主副控制器的传递函数Gc1，Gc2分别为：，。试分别采用单回路控制和串级控制设计主副PID控制器的参数，并给出整定后系统的阶跃响应的特性响应曲线和阶跃扰动的响应曲线，并说明不同控制方案对系统的影响。解：串级控制设计是一个反复调整测试的过程，使用Simulink能大大简化这一过程。根据题意，首先建立如图的Simulink模型。图中采用单回路控制的Simulink图，其中，q1为一次扰动，取阶跃信号；q2为二次扰动，取阶跃信号；Go2为副对象；Go1为主对象；r为系统输入，取阶跃信号；c为系统输出，它连接到示波器上，可以方便地观测输出。图中的PID C1为单回路PID控制器，它是按照PID原理建立的Simulink中的子模块，其内部结构如下图所示： PID控制器模块子系统的参数设置如下图：经过不断的试验，当输入比例系数为3.7，积分系数为38，微分系数为0时，系统阶跃响应达到比较满意的效果，系统阶跃响应如下图：采用这套PID参数时，二次扰动作用下，系统的输出响应如下图：采用这套PID参数时，一次扰动作用下，系统输出响应如下图：综合以上各图可以看出采用单回路控制，系统的阶跃响应达到要求时，系统对一次，二次扰动的抑制效果不是很好。下面考虑采用串级控制时的情况，下图为串级控制时的Simulink模型图：图中，q1为一次扰动，取阶跃信号；q2为二次扰动，取阶跃信号；PID C1为主控制器，采用PID控制，PID C2为副控制器，采用PID控制；Go2为副对象；Go1为主对象；r为系统输入，取阶跃信号；c为系统输出，它连接到示波器上，可以方便地观测输出。经过不断试验，当PID C1为主控制器输入比例系数为8.4，积分系数为12.8，微分系数为0时；当PID C2为主控制器输入比例系数为10，积分系数为0，微分系数为0时；系统阶跃响应达到比较满意的效果，系统阶跃响应如下图所示：采用这套PID参数时，二次扰动作用下，系统的输出响应如下图：采用这套PID参数时，一次扰动作用下，系统的输出响应如下图：综合以上各图可以看出，采用串级控制，系统的阶跃响应达到要求时，系统对一次扰动，二次扰动的抑制也能达到很好的效果。综合单回路控制和串级控制的情况，系统的控制性能对比如下表所示。表一 系统采用单回路控制和串级控制的对比控制品质指标单回路控制Kc1=3.7，Tc1=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参数的整定例二：串级控制的参数整定仿真某隧道窑炉系统，考虑烧成带温度为主变量，燃烧室温度为副变量所构成的串级控制系统中，其主副对象的传递函数Go1，Go2分别为： ，试采用逐次逼近法整定PID控制器的参数，并给出整定后系统的阶跃响应曲线。解：按照逐次逼近法的步骤1先主回路开环，按单回路方法整定副控制器，建立的Simulink框图如下图：不断地实验，当Kc2=10时，副回路阶跃响应如下图：从图中可以大致地看出，此衰减比约为4:1，因此此时可进入步骤2，即主回路闭环，取Kc2=10整定主控制器，此时的Simulink框图如下：不断地实验，调节主调节器，当Kc1=12，主回路的阶跃响应如下图所示：由图可知此时的Ts1=16，按步骤3，将主回路闭环的条件下，重新整定副控制器参数，根据衰减曲线法，取Kc1=10，Ki=8，此时系统主回路阶跃响应曲线如下图所示：由图可知系统的阶跃响应效果不理想，超调太大，需要进行精调。按照步骤4反复实验当Kc1=8.4，Ki =12.8，Kc2=10时，系统阶跃响应如下图所示：由图可知系统的阶跃响应效果比较理想，此时整定的主副回路的参数比较合理。当然，还可以调整不同的参数组合，以取得满意的输出响应曲线。这是一个反复的实验过程，利用Simulink这一便捷的工具具有很大的优势。作业题目：串级控制系统仿真。已知某串级控制系统的主副对象的传递函数Go1，Go2分别为：，副回路干扰通道的传递函数为：。（1） 画出串级控制系统的方框图及相同控制对象下的单回路控制系统方框图。（2） 用Simulink画出上述两个系统的仿真框图（3） 选用PID调节器，整定主副控制器的参数，使该串级控制系统性能良好，并绘制相应的单位阶跃响应曲线。（4） 比较单回路控制系统及串级控制系统在相同的副扰动下的单位阶跃响应曲线，并说明原因。实验四 比值控制指导内容：单闭环比例控制系统仿真综合实例。例一：假设系统从动量传递函数为，设计该从动对象的单闭环比值控制系统。解：本题的基本步骤如下：分析从动量无调节器的开环系统稳定性。由控制理论知，开环稳定性分析是系统校正的前提。系统稳定性的分析可利用Bode图进行，编制MATLAB 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。选择从动量控制器形式及整定其参数。根据工程整定的论述，选择PI形式的控制器，即。本处采用稳定边界法整定系统。先让=0，调整使系统等幅振荡（由稳定性分析图知在=1.8附近时系统震荡），即使系统处于临界稳定状态。此时的震荡周期为，比例系数为，则，。系统Simulink框图如下所示其中的PID控制器结构如下整定后从动闭环系统的单位阶跃响应如图所示：可见系统有约25%30%的超调量，在比值控制中应进一步调整使之处于振荡与不振荡的边界。调节时，系统响应图如下所示，基本达到了振荡临界要求。系统过程仿真。单闭环比值控制过程相当于从动量变化的随动控制过程。假定主动量由一常值10加幅度为0.3的随机扰动构成，从动量受均值为0、方差为1的随机干扰。主动量和从动量的比值根据工艺要求及测量仪表假定为3. 系统的控制过程Simulink仿真框图如图所示。其中控制常量及随机扰动采用封装形式。主动控制量的封装结构如下：运行结果如下所示（图中曲线从上往下分别为从动量跟踪结果、主动量给定值和随机干扰）：可见除初始时间延时外，从动量较好地跟随主动量变化而变化，并且基本维持比值3，有效地克服了主动量和从动量的扰动。例二：（续例一）单闭环比值控制系统鲁棒性。由控制理论知，在例一中延时是影响系统稳定性的关键环节。假设延时时间常数变化，即延时时间为4.55.5，分析系统鲁棒性。解：系统仿真框图如图a所示，图b为延时选择模块Subsystem的展开图，图c为延时常数设定方法，改变数为4.5,4.6,4.7，4.85.5共11个值。经过运行后在工作空间绘图（使用语句：plot（tout，simout）；hold on；grid on）即可见到图d的仿真结果。 分析图d仿真结果并与例一的图进行比较可见，随着延时环节的变化，从动量跟随主动量的规律有较小变化，但并未改变系统稳定性及精度，说明系统在延时发生变化时仍能正常工作，系统的鲁棒性较强。图a 系统仿真框图图b 延时选择模块统封装结构图c 延时参数的设置图d 仿真结果作业题目：在例一中如系统传递函数为，其他参数不变，试对其进行单闭环比值控制系统仿真分析，并讨论分母中“15”变化时控制系统的鲁棒性。实验五 解耦控制系统指导内容：对角阵解耦控制仿真综合实例纯原料量与含水量是影响混凝土快干性和强度的两个因素。系统输入控制量为纯原料量和含水量，系统输出量为混凝土的快干性和强度，采用对角阵解耦方法对该系统进行控制仿真。解：设某双输入双输出系统，初步选择输入x1、x2分别对应输出y1、y2。系统输入、输出之间的传递关系为： （1） 求系统相对增益以及系统耦合分析由式（1）得系统静态放大系数矩阵为即系统的第一放大系数矩阵为：系统的相对增益矩阵为：。由相对增益矩阵可以看出，控制系统输入、输出的配对选择是正确的；通道间存在较强的相互耦合，应对系统进行解耦分析。系统的输入、输出结构如下图所示。 确定解耦调节器根据解耦数学公式求解对角矩阵，即采用对角矩阵解耦后，系统的结构如下图所示。解耦前后系统的Simulink阶跃仿真框图及结果如下： 1）不存在耦合时的仿真框图和结果图a 不存在耦合时的仿真框图（上）和结果（下）2）系统耦合Simulink仿真框图和结果图b 系统耦合Simulink仿真框图（上）和结果（下）3）对角矩阵解耦后的仿真框图和结果图c对角矩阵解耦后的仿真框图（上）和结果（下）对比图a和图b可知，本系统的耦合影响主要体现在幅值变化和响应速度上，但影响不显著。其实不进行解耦通过闭环控制仍有可能获得要求品质。对比图a和图c可知，采用前馈解耦器后系统的响应和不存在耦合结果一样，采用前馈实现了系统解耦。解耦后系统可按两个独立的系统进行分析和控制。 控制器形式选择与参数整定通过解耦，原系统已可看成两个独立的单输入输出系统。考虑到PID应用的广泛性和系统无静差要求，控制器形式采用PI形式。PI参数整定通过解耦的两个单输入输出系统进行，整定采取试误法进行。当x1y1通道Kp=15，Ki=2时系统的阶跃响应如图：当x2y2通道Kp=35，Ki=5时系统阶跃响应如图： 系统仿真采用对角矩阵解耦时，控制系统如图所示。为了比较解耦和不解耦两种情况，分别列出两种情况的Simulink框图和仿真结果。解耦时系统的Simulink仿真框图及结果（第二幅图中的响应曲线从上往下依次是通道x2y2的输入波形和响应波形、通道x1y1的输入波形和响应波形以及随机扰动波形）：不解耦时系统的Simulink仿真框图及结果（第二幅图中的响应曲线在t=1s处从上往下依次是通道x2y2的输入波形和响应波形、通道x1y1的输入波形和响应波形以及随机扰动波形）：由图对比结果可知，系统解耦后系统的动态响应有一定改善，但改善不大，这是由于耦合较弱所致。因此当要求不高时，系统可以不采取解耦措施。作业题目：在例题中若输入输出之间传递关系改为，其他参数不变，试利用对角阵解耦方法实现系统的过程控制。