过程控制第4章被控对象数学模型

1、过程控制过程控制第四章第四章 被控对象的数学模型被控对象的数学模型目目 录录4.1 被控对象的数学模型被控对象的数学模型 4.2 被控对象数学模型的建立被控对象数学模型的建立4.3 机理法建立被控对象的数学模型机理法建立被控对象的数学模型4.4 实验法建立被控对象的数学模型实验法建立被控对象的数学模型定值控制系统、程序控制系统、随动系统（伺服控制系统） 线性系统线性系统和非线性系统 连续系统与离散系统 单输入单输出单输入单输出系统与多输入多输出系统 4.1 被控对象的数学模型被控对象的数学模型 要实现过程控制，首先要了解和掌握被控对象的过程特性，而用数学语言对被控对象的数学模型特性进行描述。数

2、学模型的分类:数学模型（对象特性）：对象输入量与输出量之间的关系(数学表达式)输入量：控制变量各种各样的干扰变量通道：对象的输入变量至输出变量的信号联系控制通道：控制变量至被控变量的信号联系通道干扰通道：干扰至被控变量的信号联系通道 被控对象控制通道干扰通道干扰变量控制变量被控变量对象输出为控制通道输出与各干扰通道输出之和被控对象：需要控制的设备、机械或生产过程。输出量：被控变量对象的数学模型：稳态数学模型和动态数学模型。 稳态数学模型描述对象的输入量与输出量之间的稳态关系 动态数学模型描述对象在输入量改变以后输出量的变化情况。稳态数学模型是动态数学模型在对象达到平衡时的特例。数学模型的表示方

3、法：参量模型：通过数学方程式表示常用的描述形式：微分方程、传递函数、差分方程等参量模型的微分方程的一般表达式：( )(1)()11010( )( )( )( )( )( )( )nnmnmytayta y ta y tb xtb x tb x ty(t)表示输出量，x(t)表示输入量，通常输出量的阶次不低于输入量的阶次(nm) 通常n=1，称该对象为一阶对象模型一阶对象模型；n=2，称二阶对象模型二阶对象模型。 非参量模型：采用曲线、表格等形式表示。例如“阶跃响应曲线、脉冲响应曲线等。优点：形象、清晰、定性。缺点：直接利用来分析系统较困难（必要时须进行数学处理获得参量模型）。设计过程控制方案（

4、被控变量及检测点选择，控制变量的确定，控制结构形式都与对象特性有关）整定控制器参数（控制规律的选择）指导设计生产工艺设备 进行仿真试验研究 培训系统运行操作人员 建模目的：建模的方法：机理建模、实验建模、混合建模机理建模机理建模：内部机理平衡方程（物料平衡、能量平衡、化学反应等基本方程） 数学模型优点：明确的物理意义缺点：只能用于简单、线性过程，需试验验证。实验建模实验建模在所要研究的对象上，人为的施加一个输入作用(输入量) ，然后用仪表记录表征对象特性的物理量(输出量)随时间变化的规律，得到一系列实验数据或曲线。这些数据或曲线就可以用来表示对象特性。 这种应用对象输入输出的实测数据来决定其模

5、型的方法，通常称为系统辨识系统辨识。系统辨识分为：过程辨识过程辨识和参数辨识参数辨识。有时，为进一步分析对象特性，可对这些数据或曲线进行处理，使其转化为描述对象特性的解析表达式。 混合建模混合建模将机理建模与实验建模结合起来机理分析（数学模型的结构）实验的方法（某些未知的或不确定的参数，即参数估计）。 4.24.2被控对象数学模型的建立被控对象数学模型的建立4.3 机理法建立被控对象的数学模型问题：处于平衡状态的对象加入干扰以后，不经控制系统能否自行达到新的平衡状态？iq0q0qiq图：假设初始为平衡状态qi=qo，水箱水位保持不变。当发生变化时(qiqo)，此时水箱的水位开始升高根据流体力学

6、原理，水箱出口流量与H是存在一定的对应关系的：0/qH R因此，qi H qo，直至qi=qo可见该系统受到干扰以后，即使不加控制，最终自身是会回到新的平衡状态，这种特性称为“自衡特性自衡特性”。 右图：如果水箱出口由泵打出，其不同之处在于：qi当发生变化时，qo不发生变化。如果qiqo ，水位H将不断上升，直至溢出，可见该系统是“无自衡能力无自衡能力”。绝大多数对象都有自衡能力，一般而言有自衡能力的系统比无自衡能力的系统容易控制。 R：出口阀门的阻力系数、液阻（与阀门开度有关）1、自平衡单容过程（一阶线性对象）问题：求右图所示的对象模型（输入输出模型）解：该对象的输入量为qi 输出量（被控变

7、量）为液位h根据物料平衡方程：单位时间流入水槽的物料 单位时间流出水槽的物料 水槽物料储藏量的变化率dVqqiodtV ChdhqqCiodt由于出口流量可以近似地表示为：ohqR CihdhqRdt RCidhhR qdtiq0qChR：出口阀门的阻力系数、液阻（与阀门开度有关）消去qo：TRCKR令、 idhThK qdt已知：R idhThK qdt对上式作拉氏变换：( )( )( )iTsH sH sK Q s一阶对象的传递函数：( )( )1iH sKQ sTs该对象的阶跃响应：如果qi为幅值为a的阶跃响应，则( )iaQ ss( )( )1(1)iKKaH sQ sTss Ts11

8、( )( )(1)Kah tLH sLs Ts11KaKaTLsTs11*()1TKaLsTs(1)tTKae一阶线性对象（总结）idhThK qdt( )( )1iH sKQ sTs( )(1)tTh tKaeh()h(t)T0.632h()qita01(0)(1)0( )(1)( )(1)0.632 ( )TThKaehKaeKah TKaeh K放大系数放大系数，在阶跃输入作用下，对象输出达到新的稳定值时，输出变化量与输入变化量之比，表征对象的稳态特性稳态特性。K越大，表示输入量对输出量的影响越大。T时间常数时间常数，在阶跃输入作用下，对象输出达到最终稳态变化量的63.2所需要的时间，时

9、间常数T是反映响应变化快慢或响应滞后的重要参数，表征对象的动态特性。用T表示的响应滞后称阻容滞后（容量滞后）。 T大，响应速度慢，控制不及时； T小，响应速度快，控制及时，T太小易引起振荡。2 2无自平衡单容过程无自平衡单容过程 受扰过程的平衡状态被破坏后，在没有操作人员或仪表等干预下，依靠被控过程自身能力不能重新回到平衡状态。dtdhCQQ211QdthdC因为Q2=0，则可得：CssQsHsG1)()()(1令：T=C（积分时间常数）Ts1无自平衡单容过程阶跃响应曲线如图所示。TssQsHsG1)()()(1当过程具有纯滞后时：seTsKsQsHsG1)()()(1：过程的纯滞后时间 自平

10、衡过程： 无自平衡过程：seTssQsHsG1)()()(1带纯滞后的自平衡过程的阶跃响应：seTsKsG1)(3 3、多容过程（二阶线性对象）、多容过程（二阶线性对象）问题：求右图所示的对象模型（输入输出模型）该对象的输入量为qi 输出量为液位h2（同样利用物料平衡方程）槽1：111111 ( )idhhqqCqdtR2212 ()o2odhhqqCqdtR槽2：22211221122222()id hdhR C R CR CR ChRqdtdt联立方程求解：222121222() id hdhTTTThK qdtdt传递函数：221 21212( )( )()1(1)(1)iHsKKQ s

11、TT sTT sTsT s1Riq1C2C2Rh1h2q1qo1112222 TRCTR C KR令：已知：221 21212( )( )()1(1)(1)iHsKKQ sTT sTT sTsT s1211221222112211212122121( )( )(1)(1)111()111ttTTKah tLHsLs T sT sTTLKasTT T sTT T sTTKaeeTTTT( )iaQ ss二阶线性对象（总结）二阶线性对象（总结）典型的微分方程：典型的传递函数：典型的阶跃响应函数：qita222121222()id hdhT TTThK qdtdt221 21212( )( )()1

12、(1)(1)iHsKKQ sTT sTT sTsT s121222121( )1ttTTTTh tKaeeTTTT2( )h t例题：试分析下图所示串联液体贮槽对象的动态特性，写出以q1的变化量为输入、液位h2的变化量为输出的微分方程式及传递函数表达式。图中，C1、C2分别为贮槽1及贮槽2的截面积，阀门的阻力系数分别为R1、R2。 （提示： ）1221hhqR解：1121dhqqCdt1221 hhqR2232dhqqCdt232 hqR22211221122212212()d hdhR C R CR CR CR ChRqdtdt得：22222hdhqCRdt代入：2211212()hdhdh

13、qCCRdtdt代入得：1222212hhhdhCRRdt22211122112221( )( )()1HsRQ sRC R C sRCR CR C sq在工业生产过程中，对象特性往往比较复杂，很难得到对象特性的数学描述；即使得到也是高阶微分方程或偏微分方程，难以求解 ，因此常用实验方法来分析对象特性。 4.4 实验法建立被控对象的数学模型实验法建立被控对象的数学模型 实验方法主有以下几种： 1、加专门信号 时域方法 例如：阶跃反应曲线法、矩形脉冲法 频域方法 例如：正弦波、梯形波 统计相关法 （施加随机信号） 例如：白噪声2、不加专门信号：正常操作时记录信号q所谓实验测取对象特性，就是在所研

14、究的对象上施加人为输入作用(输入量)，用仪表测取并记录表征对象特性的物理量(输出量)随时间变化的规律，得到一系列实验数据或曲线(非参量模型)，有时再加以必要的数据处理（得到参量模型）。通常称为系统辨识系统辨识。1、阶跃反应曲线法、阶跃反应曲线法Q1h Q2Q1th tt0t0水槽对象及其阶跃响应曲线p测量对象在阶跃输入信号作用下，输出量随时间的变化规律。 优点：方法简单。不需要特殊的信号发生器；一般也不需要增加特殊的记录仪表。测量工作量不大，数据处理也较方便。 缺点：精度较差。对象在阶跃信号作用下，从不稳定到稳定一般所需时间较长，受其它扰动因素影响，精度受到限制；若加大阶跃信号幅值，又会影响正

15、常的生产。阶跃信号一般取额定值的5%15%。加测试信号前，要求系统尽可能保持稳定稳定状态，否则会影响测试结果在相同条件下重复测试多次重复测试多次以抽取其共性减少随机干扰因素的影响测试过程注意：3、记录输入量开始作阶跃信号变化时至反应曲线的起始点的时间 （滞后时间）5、分别做阶跃输入信号为正、反方向两种变化情况的试验对比，以反映非线性对被控过程的影响。2、合理选择阶跃信号值（一般阶跃变化在正常输入信号最大幅值的 515之间，大多取10）2、矩形脉冲法、矩形脉冲法x ty tt0t1矩形脉冲特性曲线矩形脉冲波信号x tt0t1t2正弦信号x tt0 除此之外，还可采用矩矩形脉冲波形脉冲波和正弦正弦

16、信号，分别称为矩形脉冲波法矩形脉冲波法和频率特性法频率特性法。q对象在t0时刻突然加一阶跃干扰，到t1时刻突然除去阶跃干扰，输出量随时间的变化规律。 q由于加在对象上的干扰经过一段时间后被除去，干扰幅值可以取的比较大，以提高实验精度，而对象输出又不致于长期偏离给定值，对正常生产影响较小。 测试结果具有较高的精度，但数据处理较为复杂，需要进行相应的转换。方法：根据阶跃响应曲线的形状选定模型传递函数的形式确定具体参数。1)(TsKsG) 1)(1()(21sTsTKsGnTsKsG) 1()(seTsKsG1)(sesTsTKsG) 1)(1()(21sneTsKsG) 1()(或或或3 3、由阶

17、跃响应曲线确定过程传递函数、由阶跃响应曲线确定过程传递函数自平衡能力过程传递函数近似形式：无自平衡特性的被控对象传递函数近似形式：seTssG1)(sesTsTKsG) 1()(21snesTsTKsG) 1()(21TssG1)() 1()(21sTsTKsGnsTsTKsG) 1()(21或或或两方面考虑：u根据被控过程的先验知识选用合适的传递函数形式；u根据建立数学模型的目的及对模型的准确性要求，选用合适的传递函数形式。 在满足精度要求的情况下，尽量选用低阶传递函数低阶传递函数的形式。大量的实际工业过程一般都采用一、二阶传递函数的形式来描述。 确定了传递函数形式之后，由阶跃响应曲线来求取

18、被控对象动态特性的特征参数（即放大系数K、时间常数T、滞后时间等），被控过程的数学模型（传递函数）就可确定。4、传递函数的选用、传递函数的选用由阶跃响应曲线确定一阶惯性环节的特性参数由阶跃响应曲线确定一阶惯性环节的特性参数阶跃响应曲线如下图时可用一阶惯性环节来近似。要确定的参数只有T和K。 K= y()/x05 5、确定传递函数的方法、确定传递函数的方法( )( )( )1y sKG sx sTs在阶跃响应曲线的起点处做切线，该切线与y()的交点所对应的时间即为T。若阶跃变化量为X0）：0( )(1)tTy tKXe00( )(1)TyKXeKX 0ty t/001()t TtKXeT0KXT

19、tan根据测试数据直接计算T：令则2/Tt )(%39)2/(yTyTt )(%63)(yTy)(%5 .86)2(yTyTt2 在阶跃响应曲线上找到上述几个数据所对应的时间t1、t2、t3，则可计算出T。如果由t1、t2 和t3分别取的数值T有差异，可以用求平均值的方法对T加以修正。 0( )(1)( )(1)ttTTy tKXeye若阶跃变化量为X0）： ( )( )( )1y sKG sx sTs若误差较大，则选用其它形式近似。2. 由阶跃响应曲线确定一阶惯性加滞后环节的特性参数由阶跃响应曲线确定一阶惯性加滞后环节的特性参数阶跃响应曲线是S形单调曲线时，选用有纯滞后的一阶惯性环节作为该过

20、程的传递函数。seTsKsG1)(放大系数：K= y()/x0YtY()确定被控过程时间常数T与滞后时间方法：作图法和两点计算法。 作图法 曲线斜率最大（A点）处作一条切线，该切线与时间轴交于B点，与y(t)的稳态值y()交于C点，C点在时间轴上的投影为D点，则 T: BD: OBseTsKsG1)( 计算法 首先将y(t)转换成无量纲形式y0(t)。)()()()(00ytyKxtyty其阶跃响应无量纲的形式如下：tettyTt10)(0 在图中选取二个不同时刻t1和t2，以及对应的y0(t1)和y0(t2)，其中t2 t1 ，通过计算可确定T和。 TtTtetyety211)(1)(201

21、0)(1ln)(1ln201012tytyttT)(1ln)(1ln)(1ln)(1ln2010201102tytytyttyt 为了方便计算，可选y0(t1)=0.39、y0(t2)=0.632，代入上两式可得 21122)(2ttttT 计算出T和后，还应该把计算的结果与实测曲线进行比较，以检验所得模型的准确性。 若计算结果与实测值的差距可以接受，表明所求得的一阶惯性加滞后环节传递函数满足要求。否则，表明用一阶惯性加滞后环节近似被控过程的传递函数不合适，应选用高阶传递函数。 )(1ln)(1ln201012tytyttT)(1ln)(1ln)(1ln)(1ln2010201102tytyt

22、yttyt例：40.8:tT04( )0.55y t3. 由阶跃响应曲线确定二阶环节的特性参数 对无滞后的二阶环节，只需确定参数K、T1和T2。其相应的阶跃响应曲线如图所示。) 1)(1()(21sTsTKsG 化为无量纲形式的阶跃响应y0(t)后，传递函数如下表示： ) 1)(1(1)(21sTsTsG其相应的单位阶跃响应为：2112221101)(TtTteTTTeTTTty 根据上式，可以利用阶跃响应曲线上两个点的数据求出T1和T2。 若选取y0(t1)=0.4、y0(t2)=0.8两点，再从曲线上确定对应的t1和t2，即可得到方程组2 . 06 . 022122111212211212

23、211TtTtTtTteTTTeTTTeTTTeTTT上式的近似解为 )55. 074. 1 ()(16. 212121212121ttTTTTttTT注意：采用上式确定T1和T2时，应满足0.32t1/t20.46：被控过程数学模型应用高于二阶的环节近似，即G(s)=K/(Ts+1)n ；T (t1+t2)/2.16 n。式中的n可根据t1/t2由下表查出。 n n1 12 23 34 45 56 68 8101012121414t t1 1/t/t2 20.32 0.32 0.46 0.46 0.53 0.53 0.58 0.58 0.62 0.62 0.65 0.65 0.6850.68

24、50.71 0.71 0.7350.735 0.75 0.75 思考题与习题48 实验测得某液位控制对象的阶跃响应数据如下：t/(s)01020406080100140180250300400500600h/(cm)000.20.82.03.65.48.811.814.416.618.419.219.6（1）画出液位的阶跃响应曲线；(近似画一下，确定大致形状)（2）用一阶惯性环节加纯滞后近似描述该过程的动态特性，确定K， ，T。已知阶跃扰动为u=20%。（只要思路和求解方法，不要具体过程）seTsKsG1)( K= y()/u 首先将输出数据转换为无量纲形式，在图中选取二个不同时刻t1和t2，以及对应的y0(t1)和y0(t2)，通过计算可确定T和。)(1ln)(1ln201012tytyttT)(1ln)(1ln)(1ln)(1ln2010201102tytytyttyt一般可选y0(t1)=0.39