马宾《过程控制与自动化仪表》第4章 被控过程的数学模型.ppt

第4章 被控过程的数学模型,本章要点,1）掌握被控过程机理建模的方法与步骤； 2）熟悉被控过程的自衡和非自衡特性； 3）熟悉单容过程和多容过程的阶跃响应曲线及解析表达式； 4）重点掌握被控过程基于阶跃响应的建模步骤、作图方法和数据处理； 5）熟悉被控过程的一次完成最小二乘建模方法，学会用MATLAB语言编写算法程序。 6）熟悉被控过程的递推最小二乘建模方法，学会用MATLAB语言编写算法程序。,第一节 过程建模的基本概念,一、 被控过程的数学模型及其作用,被控过程的数学模型是指过程的输入变量与输出变量之间定量关系的描述 其中： 过程的输入变量至输出变量的信号联系称为通道 控制作用至输出变量的信号联系称为控制通道 干扰作用至输出变量的信号联系称为干扰通道 过程的输出为控制通道与干扰通道的输出之和,过程的数学模型,静态数学模型,动态数学模型,过程输出实际上是多个输入的函数，每个输入对被控量的影响都是不一样的，通常选用一个可控性良好的输入作为控制量，而其他输入被统称为外部扰动。,二、被控过程的特性,依据过程特性的不同分为自衡特性与无自衡特性、单容特性与多容特性、振荡与非振荡特性等,1有自衡特性和无自衡特性,当原来处于平衡状态的过程出现干扰时，其输出量在无人或无控制装置的干预下，能够自动恢复到原来或新的平衡状态，则称该过程具有自衡特性，否则，该过程则被认为无自衡特性。,无自衡过程及其阶跃响应曲线,自平衡特性其传递函数的典型形式有：,一阶惯性环节,二阶惯性环节,二阶惯性+ 纯滞后环节,一阶惯性+ 纯滞后环节,具有自衡特性的过程及其响应曲线,无平衡特性其传递函数的典型形式有：,一阶环节,二阶环节,二阶+纯滞后环节,一阶+纯滞后环节,3振荡与非振荡过程的特性,在阶跃输入作用下，输出会 出现多种形式。图中，a)、 b)和c)为振荡过程，d)和e) 为非振荡过程。,4具有反向特性的过程,对过程施加一阶跃输入信号，若在开始一段时间内，过程的输出先降后升或先升后降，即出现相反的变化方向，则称其为具有反向特性的被控过程。,4具有反向特性的过程,锅炉汽包水位的变化过程即为典型的具有反向特性的过程。,汽泡溃灭,三、过程建模方法,1机理演绎法 （解析法）,根据被控过程的内部机理，运用已知的静态或动态平衡关系，用数学解析的方法求取被控过程的数学模型。,2试验辨识法,先给被控过程人为地施加一个输入作用，然后记录过程的输出变化量，得到一系列试验数据或曲线，最后再根据输入输出试验数据确定其模型的结构（包括模型形式、阶次与纯滞后时间等）与模型的参数。,3. 混合法,机理演绎法与试验辩识法的相互交替使用的一种方法,第二节 解析法建立过程的数学模型,一、解析法建模的一般步骤,1）明确过程的输出变量、输入变量和其他中间变量； 2）以物料或能量平衡方程为基本依据,列写基本平衡方程； 3）以基本平衡方程减去起始平衡点的平衡方程（初始 条件），得到增量方程； 4） 整理得到以增量为变量的微分方程标准形式。,二、单容过程的解析法建模,例1：某单容液位过程，如右图。贮 罐中液位高度h为被控参数,流入贮罐 的体积流量q1为过程的输入量并可通 过阀门1的开度来改变；流出贮罐的 体积流量q2为过程的干扰，其大小可 以通过阀门2的开度来改变。试确定q1 与h之间的数学关系?,解,根据动态物料平衡关系，即在单位时间内贮罐的液体流入量与单位 时间内贮罐的液体流出量之差应等于贮罐中液体贮存量的变化率,则有：,写为增量形式为,其中,分别为偏离某平衡状态的增量。A为贮罐的截面积,假定,近似成正比而与阀门2的液阻,成反比,与,则有,带入增量式中可得单容液位过程的微分方程增量式,进行拉普拉斯变换， 并写成传递函数形式,其中：,为被控过程的时间常数,为被控过程的放大系数,为被控过程的容量系数，或称,过程容量，这里,一阶对象的特性参数都具有明显的物理意义：,放大倍数K的物理意义 K表明了稳态时，输出对输入的放大倍数 。求法： K = y( ) / x0,K 越大，表示对象的输入对输出的影响越大。,阶跃响应函数：,时间常数T的物理意义 对象受到阶跃输入后，输出达到新的稳态值的63.2所需的时间，就是时间常数T。 或对象受到阶跃输入后，输出若保持初始速度变化到新的稳态值所需时间就是时间常数。,求法：,在相同的阶跃输入作用下，对象的时间常数不同时，被控变量的响应曲线如图所示 。,T反映了对象输出对输入的响应速度 T越大，响应越慢。如水槽对象中 T=AR ，说明水槽面积越大，水位变化越慢。,T也反映了过渡过程时间 被控变量变化到新的稳态值所需要的时间理论上需要无限长。,当t时，才有yKx0 ，但是当t =3T 时，便有：,即：经过3T时间，输出已经变化了满幅值的95。这时，可以近似地认为动态过程基本结束。,例2,积分时间常数,恒量泵,输入输出特性：,我们把这种平衡被打破没有能力自动达到新平衡位置的系统称为： “无自衡能力过程”。,例3 求T1 和Q1的关系,R热阻； C热容,例4,单容过程模型总结,总结以上几个例子可以发现每个被控过程中都有一个储存“能量”的环节。,液容 液位控制系统中水箱的储水量 系数为水箱截面积A,热容 温度控制系统中热水所含的热量 系数为热水质量与比热的乘积Gcp,电容 电压控制系统系统中的电容 系数为电容量C,所谓单容过程指的是只有一个储能环节的过程。 有自衡特性的单容过程是一阶惯性环节，无自衡特性的是一个积分环节。,三、多容过程的解析法建模,以自衡特性的双容过程为例，如图设q1为过程输入量，第二个液位槽的液位h2为过程输出量，若不计第一个与第二个液位槽之间液体输送管道所形成的时间延迟，试求q1与h2之间的数学关系。,解,根据动态平衡关系， 列出以下增量方程,进行拉普拉斯变换，整理 得到传递函数、数学模型,为槽1的时间常数,为槽2的时间常数,其中,四、多容过程的解析法建模,二阶自平衡能力过程 与自平衡单容过程的阶跃响应（如曲线1）相比，双容过程的阶跃响应（如曲线2）一开始变化较慢，其原因是槽与槽之间存在液体流通阻力而延缓了被控量的变化。显然，若串接容器越多，则过程容量越大，时间延缓越长。,二阶无自平衡能力过程,例6,输入变量,输出变量,思路：提高阶次并消除中间变量,公式一,公式一,公式二,令：,如果过程为多个容器依次连接则传递函数为：,过程模型中的纯滞后,第三节 实验法建立过程的数学模型,试验辨识法可分为经典辨识法与现代辨识法两大类。 在经典辨识法中，最常用的有基于响应曲线的辨识方法； 在现代辨识法中，又以最小二乘辨识法最为常用。,一、响应曲线法,响应曲线法是指通过操作调节阀，使被控过程的控制输入产生一阶跃变化或方波变化，得到被控量随时间变化的响应曲线或输出数据，再根据输入输出数据，求取过程的输入输出之间的数学关系。响应曲线法又分为阶跃响应曲线法和方波响应曲线法,1、阶跃响应曲线法,1）试验测试前，被控过程应处于相对稳定的工作状态;,一、注意事项,2）在相同条件下应重复多做几次试验 ，减少随机干扰的影响;,3）对正、反方向的阶跃输入信号进行试验，以衡量过程的非线性程度;,4）一次试验后，应将被控过程恢复到原来的工况并稳定一段时间 再做第二次试验;,5）输入的阶跃幅度不能过大，以免对生产的正常进行产生不利影响, 但也不能过小，以防其它干扰影响的比重相对较大而影响试验结果。 一般取阶跃变化在正常输入信号最大幅值的5%10%之间，多取10%。,控制器,被控过程,测量变送,+,-,r,y,u,调节阀,信号发生器,记录仪,实验信号的测取框图,二、模型结构的确定,在完成阶跃响应试验后，应根据试验所得的响应曲线确定模型的结构,对于大多数过程，数学模型和传递函数分别为:,一阶惯性,一阶惯性+纯滞后,二阶惯性+纯滞后,二阶惯性,对于某些无自衡特性过程， 其对应的传递函数为：,注意： 对于更高阶或其它较复杂的系统，应在保证辨识精度的前提下， 数学模型结构应尽可能简单,三、模型参数的确定,（1）确定一阶环节的参数,该响应曲线可近似为无时延的一阶环节 则其输入与输出的关系为：,为过程的放大系数，,为时间常数。,其中,上式中，当,以上式为斜率在t=0处作切线，切线方程为,当,则有：,和,时,由以上分析可知，图解法：,根据,再在阶跃响应曲线的起点t=0处作切线，该切线与,的交点所对应的时间（上图中阶跃响应曲线上的OB段）即为,先确定,的确定还可以使用计算法：,令t分别为,时，则有,以及,在阶跃响应曲线上求得,三个状态下的时间t1、t2、t3，不难计算出,一般求多个值的平均,（2）确定一阶时延环节的参数,一阶惯性+时延环节的传递函数,需要确定,选取两个不同时刻t1， t2（两点法），代入,两边取自然对数， 求解化简可得：,补充例题1,利用两点法，取A(20,2.8)和B(50,5.9)两点进行计算：,解：,补充例题1,可知其参数为：,（3）确定二阶环节的参数,传递函数为：,三个需要确定的参数,的确定与一阶环节确定方法相同,的确定采用两点法。,二阶无时延环节的输入、输出关系为,含两个未知数的函数只要任意求两点就可以解出T1和T2，经常取稳态值的40%和80 % 两点来计算。,求解可得,注意：用这种方法确定T1和T2时，应满足,的条件,因为，,说明有一个时间常数远大于另一个，系统可以用一个一阶环节近似。,其中,其中,时，应为二阶以上环节。,当,这时这个函数取得极大值的条件是：T1=T2 ，则用一个二阶环节近似：,对于n阶环节传递函数,可以按,近似计算,的大小由下表确定,其中n可以根据,容量个数越多（阶数n越多），阶跃响应曲线上升越慢。,2 非自衡对象,高阶非自衡对象的近似,作图法,作图求解,（4）确定二阶时延环节的参数,二阶时延环节阶跃响应曲线如右图：,传递函数为：,需确定参数4个,通过拐点F作切线,得纯滞后时间,，容量滞后时间,以及,、,而总的纯滞后时间,在,的约束条件下，可以解得,和,这个方程为超越方程，求解比较复杂，通常采用图解法,自学图解法,有些工艺对象不允许长时间施加较大幅度的扰动，改为施加脉宽为t的方波脉冲，得到的响应曲线称为“方波响应”。,2、方波响应曲线法,方波响应曲线法是在正常输入的基础上，施加一方波输入，并测取相应输出的变化曲线，据此估计过程参数。,一个是在t = 0时加入的正阶跃信号x1（t） 另 一个是在 t =t 时加入的负阶跃信号x2（t） x（t）= x1（t）+ x2（t）,其中， x2（t）= - x1（t -t）,原理：方波信号是两个阶跃信号的代数和。,根据此式，方波响应可逐点拆分为阶跃飞升曲线y1（t）和 y2（t）。,对应的响应也为两个阶跃响应之和： y（t）= y1（t）+ y2（t） = y1（t）- y1（t-t）,如图输出响应由两个时