《过程控制与自动化仪表》第4章被控过程的数学模型课件

1、 第4章 被控过程的数学模型第1页，共56页。本章要点1）掌握被控过程机理建模的方法与步骤； 2）熟悉被控过程的自衡和非自衡特性； 3）熟悉单容过程和多容过程的阶跃响应曲线及解析表达式； 4）重点掌握被控过程基于阶跃响应的建模步骤、作图方法和数据处理； 5）熟悉被控过程的一次完成最小二乘建模方法，学会用MATLAB语言编写算法程序。 6）熟悉被控过程的递推最小二乘建模方法，学会用MATLAB语言编写算法程序。 第2页，共56页。第一节过程建模的基本概念一、 被控过程的数学模型及其作用 被控过程的数学模型是指过程的输入变量与输出变量之间定量关系的描述其中：过程的输入变量至输出变量的信号联系称为通

2、道控制作用至输出变量的信号联系称为控制通道干扰作用至输出变量的信号联系称为干扰通道过程的输出为控制通道与干扰通道的输出之和 第3页，共56页。过程的数学模型静态数学模型动态数学模型 过程输出实际上是多个输入的函数，每个输入对被控量的影响都是不一样的，通常选用一个可控性良好的输入作为控制量，而其他输入被统称为外部扰动。第4页，共56页。二、被控过程的特性依据过程特性的不同分为自衡特性与无自衡特性、单容特性与多容特性、振荡与非振荡特性等 1有自衡特性和无自衡特性 当原来处于平衡状态的过程出现干扰时，其输出量在无人或无控制装置的干预下，能够自动恢复到原来或新的平衡状态，则称该过程具有自衡特性，否则，

3、该过程则被认为无自衡特性。 第5页，共56页。无自衡过程及其阶跃响应曲线 自平衡特性其传递函数的典型形式有：一阶惯性环节 二阶惯性环节 二阶惯性+纯滞后环节 一阶惯性+纯滞后环节 具有自衡特性的过程及其响应曲线 无平衡特性其传递函数的典型形式有：一阶环节 二阶环节 二阶+纯滞后环节 一阶+纯滞后环节 第6页，共56页。3振荡与非振荡过程的特性在阶跃输入作用下，输出会出现多种形式。图中，a)、b)和c)为振荡过程，d)和e)为非振荡过程。 4具有反向特性的过程 对过程施加一阶跃输入信号，若在开始一段时间内，过程的输出先降后升或先升后降，即出现相反的变化方向，则称其为具有反向特性的被控过程。第7页

4、，共56页。4具有反向特性的过程锅炉汽包水位的变化过程即为典型的具有反向特性的过程。 冷水量 蒸发率水位进水量大 汽泡溃灭 水位第8页，共56页。三、过程建模方法1机理演绎法 （解析法） 根据被控过程的内部机理，运用已知的静态或动态平衡关系，用数学解析的方法求取被控过程的数学模型。2试验辨识法 先给被控过程人为地施加一个输入作用，然后记录过程的输出变化量，得到一系列试验数据或曲线，最后再根据输入输出试验数据确定其模型的结构（包括模型形式、阶次与纯滞后时间等）与模型的参数。 3. 混合法 机理演绎法与试验辩识法的相互交替使用的一种方法第9页，共56页。第二节 解析法建立过程的数学模型一、解析法建

5、模的一般步骤1）明确过程的输出变量、输入变量和其他中间变量；2）以物料或能量平衡方程为基本依据,列写基本平衡方程；3）以基本平衡方程减去起始平衡点的平衡方程（初始 条件），得到增量方程；4） 整理得到以增量为变量的微分方程标准形式。二、单容过程的解析法建模例1：某单容液位过程，如右图。贮罐中液位高度h为被控参数,流入贮罐的体积流量q1为过程的输入量并可通过阀门1的开度来改变；流出贮罐的体积流量q2为过程的干扰，其大小可以通过阀门2的开度来改变。试确定q1与h之间的数学关系?第10页，共56页。解 根据动态物料平衡关系，即在单位时间内贮罐的液体流入量与单位时间内贮罐的液体流出量之差应等于贮罐中液

6、体贮存量的变化率则有：写为增量形式为其中分别为偏离某平衡状态的增量。A为贮罐的截面积假定近似成正比而与阀门2的液阻成反比与则有 ,带入增量式中可得单容液位过程的微分方程增量式 进行拉普拉斯变换，并写成传递函数形式 其中：为被控过程的时间常数 为被控过程的放大系数 为被控过程的容量系数，或称 过程容量，这里第11页，共56页。一阶对象的特性参数都具有明显的物理意义：ttx0y( )yx放大倍数K的物理意义K表明了稳态时，输出对输入的放大倍数 。求法： K = y( ) / x0K 越大，表示对象的输入对输出的影响越大。阶跃响应函数：第12页，共56页。时间常数T的物理意义对象受到阶跃输入后，输出

7、达到新的稳态值的63.2所需的时间，就是时间常数T。或对象受到阶跃输入后，输出若保持初始速度变化到新的稳态值所需时间就是时间常数。tT0.632y)求法：第13页，共56页。在相同的阶跃输入作用下，对象的时间常数不同时，被控变量的响应曲线如图所示 。 T反映了对象输出对输入的响应速度T越大，响应越慢。如水槽对象中 T=AR ，说明水槽面积越大，水位变化越慢。第14页，共56页。T也反映了过渡过程时间被控变量变化到新的稳态值所需要的时间理论上需要无限长。当t时，才有yKx0 ，但是当t =3T 时，便有： 即：经过3T时间，输出已经变化了满幅值的95。这时，可以近似地认为动态过程基本结束。tTy

8、)3T第15页，共56页。例2积分时间常数恒量泵第16页，共56页。输入输出特性： 我们把这种平衡被打破没有能力自动达到新平衡位置的系统称为： “无自衡能力过程”。第17页，共56页。例3 求T1 和Q1的关系R热阻； C热容第18页，共56页。第19页，共56页。例4第20页，共56页。单容过程模型总结 总结以上几个例子可以发现每个被控过程中都有一个储存“能量”的环节。液容 液位控制系统中水箱的储水量 系数为水箱截面积A热容 温度控制系统中热水所含的热量 系数为热水质量与比热的乘积Gcp电容 电压控制系统系统中的电容 系数为电容量C 所谓单容过程指的是只有一个储能环节的过程。 有自衡特性的单

9、容过程是一阶惯性环节，无自衡特性的是一个积分环节。第21页，共56页。三、多容过程的解析法建模以自衡特性的双容过程为例，如图设q1为过程输入量，第二个液位槽的液位h2为过程输出量，若不计第一个与第二个液位槽之间液体输送管道所形成的时间延迟，试求q1与h2之间的数学关系。 解根据动态平衡关系，列出以下增量方程 进行拉普拉斯变换，整理得到传递函数、数学模型 为槽1的时间常数为槽2的时间常数 其中第22页，共56页。四、多容过程的解析法建模二阶自平衡能力过程 与自平衡单容过程的阶跃响应（如曲线1）相比，双容过程的阶跃响应（如曲线2）一开始变化较慢，其原因是槽与槽之间存在液体流通阻力而延缓了被控量的变

10、化。显然，若串接容器越多，则过程容量越大，时间延缓越长。 第23页，共56页。二阶无自平衡能力过程第24页，共56页。例6输入变量输出变量第25页，共56页。思路：提高阶次并消除中间变量公式一相加得：将第26页，共56页。公式一公式二第27页，共56页。令：第28页，共56页。如果过程为多个容器依次连接则传递函数为：第29页，共56页。过程模型中的纯滞后第30页，共56页。第三节 实验法建立过程的数学模型试验辨识法可分为经典辨识法与现代辨识法两大类。在经典辨识法中，最常用的有基于响应曲线的辨识方法；在现代辨识法中，又以最小二乘辨识法最为常用。 一、响应曲线法 响应曲线法是指通过操作调节阀，使被

11、控过程的控制输入产生一阶跃变化或方波变化，得到被控量随时间变化的响应曲线或输出数据，再根据输入输出数据，求取过程的输入输出之间的数学关系。响应曲线法又分为阶跃响应曲线法和方波响应曲线法 第31页，共56页。1、阶跃响应曲线法1）试验测试前，被控过程应处于相对稳定的工作状态; 一、注意事项2）在相同条件下应重复多做几次试验 ，减少随机干扰的影响;3）对正、反方向的阶跃输入信号进行试验，以衡量过程的非线性程度;4）一次试验后，应将被控过程恢复到原来的工况并稳定一段时间 再做第二次试验; 5）输入的阶跃幅度不能过大，以免对生产的正常进行产生不利影响, 但也不能过小，以防其它干扰影响的比重相对较大而影

12、响试验结果。 一般取阶跃变化在正常输入信号最大幅值的5%10%之间，多取10%。 第32页，共56页。控制器被控过程测量变送+-ryuym调节阀信号发生器记录仪实验信号的测取框图第33页，共56页。二、模型结构的确定在完成阶跃响应试验后，应根据试验所得的响应曲线确定模型的结构 对于大多数过程，数学模型和传递函数分别为:一阶惯性一阶惯性+纯滞后 二阶惯性+纯滞后 二阶惯性 对于某些无自衡特性过程， 其对应的传递函数为：注意： 对于更高阶或其它较复杂的系统，应在保证辨识精度的前提下，数学模型结构应尽可能简单 第34页，共56页。三、模型参数的确定（1）确定一阶环节的参数 该响应曲线可近似为无时延的

13、一阶环节则其输入与输出的关系为：为过程的放大系数，为时间常数。 其中上式中，当以上式为斜率在t=0处作切线，切线方程为当则有：和时第35页，共56页。由以上分析可知，图解法：根据再在阶跃响应曲线的起点t=0处作切线，该切线与的交点所对应的时间（上图中阶跃响应曲线上的OB段）即为 先确定的确定还可以使用计算法：令t分别为时，则有以及在阶跃响应曲线上求得三个状态下的时间t1、t2、t3，不难计算出一般求多个值的平均第36页，共56页。T第37页，共56页。（2）确定一阶时延环节的参数 一阶惯性+时延环节的传递函数 需要确定转化为标幺值和的确定步骤是：先将阶跃响应即：相应的阶跃响应表达式为 选取两个

14、不同时刻t1，t2（两点法），代入两边取自然对数，求解化简可得：求出和第38页，共56页。补充例题1利用两点法，取A(20,2.8)和B(50,5.9)两点进行计算：解：第39页，共56页。补充例题1可知其参数为：第40页，共56页。（3）确定二阶环节的参数 传递函数为：三个需要确定的参数的确定与一阶环节确定方法相同 的确定采用两点法。二阶无时延环节的输入、输出关系为 含两个未知数的函数只要任意求两点就可以解出T1和T2，经常取稳态值的40%和80 % 两点来计算。求解可得第41页，共56页。注意：用这种方法确定T1和T2时，应满足的条件 因为，说明有一个时间常数远大于另一个，系统可以用一个一

15、阶环节近似。其中其中时，应为二阶以上环节。 当当时当时这时这个函数取得极大值的条件是：T1=T2 ，则用一个二阶环节近似：第42页，共56页。对于n阶环节传递函数可以按近似计算的大小由下表确定其中n可以根据n12345678101214t1/t20.320.460.530.580.620.650.670.6850.710.7350.75高阶过程的n与的关系h ()Othn=1 n=2 n=3 n=4 n=5 容量个数越多（阶数n越多），阶跃响应曲线上升越慢。第43页，共56页。2 非自衡对象高阶非自衡对象的近似作图法作图求解第44页，共56页。第45页，共56页。（4）确定二阶时延环节的参数

16、二阶时延环节阶跃响应曲线如右图： 传递函数为：需确定参数4个通过拐点F作切线 得纯滞后时间 ，容量滞后时间 以及、而总的纯滞后时间 可以证明：与的关系为其中在的约束条件下，可以解得和这个方程为超越方程，求解比较复杂，通常采用图解法 自学图解法第46页，共56页。有些工艺对象不允许长时间施加较大幅度的扰动，改为施加脉宽为t的方波脉冲，得到的响应曲线称为“方波响应”。2、方波响应曲线法方波响应曲线法是在正常输入的基础上，施加一方波输入，并测取相应输出的变化曲线，据此估计过程参数。第47页，共56页。 一个是在t = 0时加入的正阶跃信号x1（t） 另 一个是在 t =t 时加入的负阶跃信号x2（t） x（t）= x1（t）+ x2（t）其中， x2（t）= - x1（t -t）原理：方波信号是两个阶跃信号的代数和。tttttxxxx0 x0 x0第48页，共56页。 根据此式，方波响应可逐点拆分为阶跃飞升曲线y1（t）和 y2（t）。 对应的响应也为两个阶跃响应之和： y（t）= y1（t）+ y2（t） = y1（t）- y1（t-t）ty2(t)OtxOO tttx1(t)x2(t)= x1(t-t)txy1(t)y(t)y(t)第49页，共56页。如图输出响应由两个时间相差t0、极性相反、形状完全相同的