过程控制系统建模方法.pptx

《过程控制系统建模方法.pptx》由会员分享，可在线阅读，更多相关《过程控制系统建模方法.pptx（117页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、v 建模概念建模概念v过程控制系统建模的两个基本方法过程控制系统建模的两个基本方法42.1过程控制系统建模概念过程控制系统建模概念第1页/共117页建模的概念建模的概念 建立被控对象的数学模型建立被控对象的数学模型建立被控对象的数学模型建立被控对象的数学模型 例如，典型的例如，典型的例如，典型的例如，典型的RLCRLC电路电路电路电路u uc cu u用方程列写的数学模型用方程列写的数学模型用方程列写的数学模型用方程列写的数学模型用图形表示的数学模型用图形表示的数学模型用图形表示的数学模型用图形表示的数学模型输出输出输出输出u uc c输入输入输入输入 u u第2页/共117页建模需要三类主要

2、的信息源建模需要三类主要的信息源1、要确定明确的输入量与输出量、要确定明确的输入量与输出量通通常常选选一一个个可可控控性性良良好好，对对输输出出量量影影响响最最大大的的一一个个输输入入信信号号作作为为输输入入量量，其其余余的的输输入入信信号号则则为为干干扰量。扰量。建模的概念建模的概念第3页/共117页2、要有先验知识、要有先验知识在在建建模模中中，被被控控对对象象内内部部所所进进行行的的物物理理、化化学学过过程符合已经发现的许多定理、原理及模型。程符合已经发现的许多定理、原理及模型。在建模中必须掌握建模对象所要用到的先验知识。在建模中必须掌握建模对象所要用到的先验知识。3、试验数据、试验数据

3、过程的信息也能通过对对象的试验与测量而获得。过程的信息也能通过对对象的试验与测量而获得。合适的实验数据是验证模型和建模的重要依据。合适的实验数据是验证模型和建模的重要依据。建模需要三类主要的信息源（续）建模需要三类主要的信息源（续）第4页/共117页建模的概念建模的概念图图2.1 2.1 数学建模的信息源数学建模的信息源第5页/共117页用于控制的数学模型要求准确可靠，但并非越准确越好，这是因为用于控制的数学模型要求准确可靠，但并非越准确越好，这是因为：准确准确复杂复杂实时性差实时性差影响在线运用影响在线运用闭闭环环控控制制本本身身具具有有一一定定的的鲁鲁棒棒性性，模模型型的的误误差差可可以以

4、视视为为干干扰扰，闭闭环环控控制制在某种程度上具有自动消除干扰影响的能力。在某种程度上具有自动消除干扰影响的能力。所以所以适用性适用性在在建建立立数数学学模模型型时时，要要抓抓住住主主要要因因素素，忽忽略略次次要要因因素素，需需要要做做很很多多近近似似处理，如：线性化、分布参数系统和模型降阶处理等。处理，如：线性化、分布参数系统和模型降阶处理等。建模的概念建模的概念关于过程系统建模的说明关于过程系统建模的说明模型的准确性模型的准确性第6页/共117页42.1过程控制系统建模概念过程控制系统建模概念 建模概念建模概念 过程控制系统建模的两个基本方法过程控制系统建模的两个基本方法第7页/共117页

5、机理法建模机理法建模机理法建模机理法建模 过程系统过程系统过程系统过程系统建模方法建模方法建模方法建模方法测试法建模测试法建模 过程控制系统建模的两个基本方法过程控制系统建模的两个基本方法第8页/共117页1、机理法建模、机理法建模根根据据生生产产过过程程中中实实际际发发生生的的变变化化机机理理，写写出出各各种有关的平衡方程：种有关的平衡方程：物物质质平平衡衡方方程程；能能量量平平衡衡方方程程；动动量量平平衡衡方方程程以以及及反反映映流流体体流流动动、传传热热、传传质质、化化学学反反应应等等基基本本规规律律的的运运动动方方程程，物物性性参参数数方方程程和和某某些些设设备备的的特特性性方方程程等

6、等。如如RLC电路，遵循基尔霍夫定理电路，遵循基尔霍夫定理如单容水槽，遵循物料平衡关系如单容水槽，遵循物料平衡关系过程控制系统建模的两个基本方法过程控制系统建模的两个基本方法第9页/共117页1、机理法建模（续）、机理法建模（续）机理法建模的首要条件机理法建模的首要条件生产过程的机理为人们充分掌握，并且可以比较确切地加以数学描述生产过程的机理为人们充分掌握，并且可以比较确切地加以数学描述机理法建模的问题机理法建模的问题烦琐烦琐模型中某些参数难以确定模型中某些参数难以确定机理建模的准则机理建模的准则简单、适用、满足合理的精度和实时性要求简单、适用、满足合理的精度和实时性要求过程控制系统建模的两个

7、基本方法过程控制系统建模的两个基本方法第10页/共117页2、测试法建模、测试法建模根据工业过程的输入和输出的实测数据进行数根据工业过程的输入和输出的实测数据进行数学处理后得到的模型。学处理后得到的模型。特点是把被研究的工业过程视为一个黑匣子，特点是把被研究的工业过程视为一个黑匣子，完全从外特性上测试和描述它的完全从外特性上测试和描述它的动态动态性质，不性质，不需要深入掌握其内部机理。需要深入掌握其内部机理。过程控制系统建模的两个基本方法过程控制系统建模的两个基本方法输入输入输入输入输出输出输出输出二二二二 阶阶阶阶 系系系系 统统统统第11页/共117页测试建模中施加扰动的必要性测试建模中施

8、加扰动的必要性为了获得动态特性，必须使被研究的过程处于被激励的状态，如施加一个为了获得动态特性，必须使被研究的过程处于被激励的状态，如施加一个阶跃扰动或脉冲扰动等。阶跃扰动或脉冲扰动等。测试建模中对过程机理了解的必要性测试建模中对过程机理了解的必要性测试建模较机理建模简单测试建模较机理建模简单测试建模的分类测试建模的分类经典辨识法和现代辨识法经典辨识法和现代辨识法过程控制系统建模的两个基本方法过程控制系统建模的两个基本方法2、测试法建模（续）、测试法建模（续）第12页/共117页2.2机理建模方法机理建模方法v单容对象的传递函数单容对象的传递函数v具有纯迟延的单容对象特性具有纯迟延的单容对象特

9、性v无自平衡能力的单容对象特性无自平衡能力的单容对象特性v多容对象的动态特性多容对象的动态特性第13页/共117页1、单容水槽、单容水槽涉及变量涉及变量流入量流入量Qi调节阀开度调节阀开度u流出量流出量Qo水位水位h研究目标研究目标分分析析水水位位在在调调节节阀阀开开度扰动下的动态特性度扰动下的动态特性u h单容对象的传递函数单容对象的传递函数图图2.2 单容水槽单容水槽1控制阀门控制阀门 2水槽水槽 3负载阀负载阀(液阻液阻R)第14页/共117页各量定义如下：各量定义如下：Qi输入水流量的稳态值输入水流量的稳态值(m3s)Qi输入水流量的增量输入水流量的增量(m3s)Qo输出水流量的稳态值

10、输出水流量的稳态值(m3s)Qo输出水流量的增量输出水流量的增量(m3s)h0液位的稳态值液位的稳态值(m)h液位的增量液位的增量(m)u调节阀的开度调节阀的开度(m2)A液槽横截面积液槽横截面积(m2)R流出侧负载阀门的液阻流出侧负载阀门的液阻(sm2)单容对象的传递函数单容对象的传递函数单容对象的传递函数单容对象的传递函数第15页/共117页根据物料平衡关系，有：根据物料平衡关系，有：初始时刻，水槽处于平衡状态：初始时刻，水槽处于平衡状态：nQo=Qi，hh0进水阀开度发生阶跃变化进水阀开度发生阶跃变化u时：时：QiQi+Qihh+hQoQo+Qo于是有于是有下一步设法消去下一步设法消去Q

11、o，并将，并将Qi折算到折算到u单容对象的传递函数单容对象的传递函数单容对象的传递函数单容对象的传递函数第16页/共117页v将将Qi折算到折算到uQi 是是由由控控制制阀阀开开度度变变化化u引引起起的的，当当阀阀前后压差不变时，前后压差不变时，Qi与与u成正比关系：成正比关系：其中，其中，Ku为阀门流量系数为阀门流量系数(ms)单容对象的传递函数单容对象的传递函数单容对象的传递函数单容对象的传递函数第17页/共117页设法消去设法消去Qo流出量与液位高度的关系为：流出量与液位高度的关系为：在液位稳态值在液位稳态值(Qo,h0)附近线性化，得：附近线性化，得：单容对象的传递函数单容对象的传递函

12、数单容对象的传递函数单容对象的传递函数第18页/共117页单容水槽的传递函数单容水槽的传递函数单容水槽的传递函数单容水槽的传递函数单容对象的传递函数单容对象的传递函数单容对象的传递函数单容对象的传递函数第19页/共117页被控参数为炉内温度被控参数为炉内温度T,控控制量为电热丝两端电压制量为电热丝两端电压u由热力学知识，有：由热力学知识，有：其中，其中，M为加热丝质量为加热丝质量C为比热为比热H为传热系数为传热系数A为传热面积为传热面积单容对象的传递函数单容对象的传递函数2、电加热炉、电加热炉一阶惯性环一阶惯性环一阶惯性环一阶惯性环节节节节第20页/共117页 单容对象特性分析单容对象特性分析

13、单容对象特性分析单容对象特性分析 单容对象的动态特性都是一阶惯性环节单容对象的动态特性都是一阶惯性环节单容对象的动态特性都是一阶惯性环节单容对象的动态特性都是一阶惯性环节 解微分方程解微分方程解微分方程解微分方程 得得得得单容对象的传递函数单容对象的传递函数单容对象的传递函数单容对象的传递函数第21页/共117页 单容对象特性分析（续）单容对象特性分析（续）单容对象特性分析（续）单容对象特性分析（续）对时间常数对时间常数对时间常数对时间常数T T的说明的说明的说明的说明 T T反映对象响应速度的快慢反映对象响应速度的快慢反映对象响应速度的快慢反映对象响应速度的快慢 对单容水槽，对单容水槽，对单

14、容水槽，对单容水槽，T=RAT=RA 对放大系数对放大系数对放大系数对放大系数K K的说明的说明的说明的说明 K K是系统的稳态指标是系统的稳态指标是系统的稳态指标是系统的稳态指标 K K大，系统的灵敏度高大，系统的灵敏度高大，系统的灵敏度高大，系统的灵敏度高单容对象的传递函数单容对象的传递函数单容对象的传递函数单容对象的传递函数第22页/共117页2.2机理建模方法机理建模方法v单容对象的传递函数单容对象的传递函数v具有纯迟延的单容对象特性具有纯迟延的单容对象特性v无自平衡能力的单容对象特性无自平衡能力的单容对象特性v 多容对象的动态特性多容对象的动态特性第23页/共117页有有一一储储水水

15、槽槽调调节节阀阀1距距水水槽槽有有一一段段较较长长的的距距离离。调调节节阀阀1开开度度变变化化所所引引起起的的流流入入量量变变化化Qi，需需要要经经过过一一段段传传输输时时间间T0，才才能能对对水水槽槽液液位位产产生生影响，影响，T0是纯延迟时间。是纯延迟时间。具有纯延迟的单容对象特性具有纯延迟的单容对象特性纯延迟现象产生的原因是由于纯延迟现象产生的原因是由于纯延迟现象产生的原因是由于纯延迟现象产生的原因是由于扰动发生的地点与测定被控参扰动发生的地点与测定被控参扰动发生的地点与测定被控参扰动发生的地点与测定被控参数位置有一定距离。数位置有一定距离。数位置有一定距离。数位置有一定距离。第24页/

16、共117页无延迟的单容对象的微分方程无延迟的单容对象的微分方程无延迟的单容对象的微分方程无延迟的单容对象的微分方程有有有有纯纯纯纯延延延延迟迟迟迟的的的的单单单单容容容容对对对对象象象象的的的的微微微微分分分分方方方方程程程程纯延迟的单容对象的传递函数纯延迟的单容对象的传递函数纯延迟的单容对象的传递函数纯延迟的单容对象的传递函数具有纯延迟的单容对象特性具有纯延迟的单容对象特性第25页/共117页2.2机理建模方法机理建模方法v单容对象的传递函数单容对象的传递函数v具有纯迟延的单容对象特性具有纯迟延的单容对象特性v无自平衡能力的单容对象特性无自平衡能力的单容对象特性v多容对象的动态特性多容对象的

17、动态特性第26页/共117页自平衡过程自平衡过程自平衡过程自平衡过程受扰后被调量能够自动地稳定在新的平衡点上的过受扰后被调量能够自动地稳定在新的平衡点上的过程程如，用惯性环节描述的单容对象如，用惯性环节描述的单容对象自平衡过程是一种稳定的过程自平衡过程是一种稳定的过程无自平衡过程无自平衡过程受扰后，无法自动恢复平衡的过程受扰后，无法自动恢复平衡的过程如，用积分环节描述的单容对象如，用积分环节描述的单容对象无自平衡能力的单容对象特性无自平衡能力的单容对象特性第27页/共117页无自平衡能力的单容水槽无自平衡能力的单容水槽无自平衡能力的单容对象特性无自平衡能力的单容对象特性无自平衡能力的单容对象无

18、自平衡能力的单容对象 流出端采用容积式计量泵流出端采用容积式计量泵流出端采用容积式计量泵流出端采用容积式计量泵 排出恒定的流量排出恒定的流量排出恒定的流量排出恒定的流量QoQo 输入流量受扰后，水位或一直上升或一直下降输入流量受扰后，水位或一直上升或一直下降输入流量受扰后，水位或一直上升或一直下降输入流量受扰后，水位或一直上升或一直下降 无法通过控制使其平衡无法通过控制使其平衡无法通过控制使其平衡无法通过控制使其平衡第28页/共117页无自平衡能力的单容对象其动态方程为无自平衡能力的单容对象其动态方程为将上式改写为将上式改写为无自平衡能力的单容对象的特性分析无自平衡能力的单容对象的特性分析无自

19、平衡能力的单容对象特性无自平衡能力的单容对象特性A A 液槽截面积液槽截面积响应速度响应速度Ta响应时间响应时间第29页/共117页图图2.9 无自平衡能力单容对象阶跃响应曲线无自平衡能力单容对象阶跃响应曲线 无自平衡能力的单容对象特性无自平衡能力的单容对象特性无自平衡能力的单容对象的特性分析（续）无自平衡能力的单容对象的特性分析（续）无自平衡能力的单容对象的传递函数无自平衡能力的单容对象的传递函数这是一个积分环节这是一个积分环节第30页/共117页2.2机理建模方法机理建模方法v单容对象的传递函数单容对象的传递函数v具有纯迟延的单容对象特性具有纯迟延的单容对象特性v无自平衡能力的单容对象特性

20、无自平衡能力的单容对象特性v多容对象的动态特性多容对象的动态特性第31页/共117页具有两个水槽具有两个水槽双容对象双容对象1）有自衡能力的双容对象）有自衡能力的双容对象2）有自衡能力的多容对象）有自衡能力的多容对象3）无自衡能力的双容对象）无自衡能力的双容对象4）相互作用的双容对象）相互作用的双容对象多容对象的动态特性多容对象的动态特性第32页/共117页1）具有自平衡能力的双容对象）具有自平衡能力的双容对象两个串联对象的模型两个串联对象的模型流入：阀门开度的微小扰动流入：阀门开度的微小扰动u被控参数：下水槽的水位变化被控参数：下水槽的水位变化h2u与与h2间的动态方程间的动态方程多容对象的

21、动态特性多容对象的动态特性第33页/共117页其中：其中：C C1 1、C C2 2两液槽的容量系数；两液槽的容量系数；R R1 1、R R2 2两液槽的出水端的阻力；两液槽的出水端的阻力；T T1 1=R=R1 1C C1 1第一个容器的时间常数；第一个容器的时间常数；T T2 2=R=R2 2C C2 2第二个容器的时间常数；第二个容器的时间常数；K=KK=Ku uR R2 2双容对象的放大系数。双容对象的放大系数。多容对象的动态特性多容对象的动态特性v具有自平衡能力的双容对象的动态特性方程具有自平衡能力的双容对象的动态特性方程第34页/共117页多容对象的动态特性多容对象的动态特性 具有

22、自平衡能力的双容对象的传递函数具有自平衡能力的双容对象的传递函数具有自平衡能力的双容对象的传递函数具有自平衡能力的双容对象的传递函数 有纯延迟时有纯延迟时有纯延迟时有纯延迟时第35页/共117页图图图图2.112.11具有自平衡能力的双容对象的阶跃响应具有自平衡能力的双容对象的阶跃响应具有自平衡能力的双容对象的阶跃响应具有自平衡能力的双容对象的阶跃响应多容对象的动态特性多容对象的动态特性第36页/共117页2）具有自平衡能力的多容对象）具有自平衡能力的多容对象有有n个个相相互互独独立立的的多多容容对对象象的的时时间间常常数数为为T1、T2.Tn，总放大系数为，总放大系数为K，则传递函数为：，则

23、传递函数为：多容对象的动态特性多容对象的动态特性第37页/共117页若若T1=T2=Tn=T则 G(s)=若还有纯延迟，则若还有纯延迟，则 G(s)=2）具有自平衡能力的多容对象）具有自平衡能力的多容对象多容对象的动态特性多容对象的动态特性第38页/共117页3）无自平衡能力的双容对象）无自平衡能力的双容对象一个有自平衡能力的单容对象和一个无自平一个有自平衡能力的单容对象和一个无自平衡能力的单容对象的串联衡能力的单容对象的串联多容对象的动态特性多容对象的动态特性一阶惯性环节一阶惯性环节积分环节积分环节由由于于多多了了一一个个中中间间液液糟糟,作作为为被被 控控 参参 数数 的的h2，并并不不能

24、能立立即即以以最最大大速速度变化度变化第39页/共117页 R1C1 +=多容对象的动态特性多容对象的动态特性多容对象的动态特性方程：多容对象的动态特性方程：第40页/共117页令令T=R1C1,Ta=C2/Ku，则得则得T +=其对应的传递函数为其对应的传递函数为G(s)=多容对象的动态特性多容对象的动态特性第41页/共117页有纯延迟的情况则有纯延迟的情况则 G(s)=多容对象的动态特性多容对象的动态特性第42页/共117页图图2.13无自平衡能力双容对象的阶跃响应曲无自平衡能力双容对象的阶跃响应曲线线多容对象的动态特性多容对象的动态特性第43页/共117页4）相互作用的双容对象）相互作用

25、的双容对象两个水槽中，一水槽液位的高低会影响另一两个水槽中，一水槽液位的高低会影响另一水槽液位变化，两者之间有相互作用，结果水槽液位变化，两者之间有相互作用，结果会改变水槽的等效时间常数。会改变水槽的等效时间常数。多容对象的动态特性多容对象的动态特性图图2.14 具有相互作用的双容模型具有相互作用的双容模型第44页/共117页设被控参数为设被控参数为输入扰动为输入扰动为原来平衡时原来平衡时Q0=Q1=Qi,H10=h20当输入有扰动当输入有扰动后后多容对象的动态特性多容对象的动态特性4）相互作用的双容对象）相互作用的双容对象第45页/共117页多容对象的动态特性多容对象的动态特性（1）（2）（

26、4）（3）第46页/共117页可得对应的传递函数为可得对应的传递函数为 多容对象的动态特性多容对象的动态特性4）相互作用的双容对象）相互作用的双容对象v若以若以h2为被控参数，则为被控参数，则第47页/共117页v2.1 过程控制系统建模概念过程控制系统建模概念v2.2 机理建模方法机理建模方法v2.3 测试建模方法测试建模方法第二章第二章过程控制系统建模方过程控制系统建模方法法第48页/共117页当生产过程机理不明、模型参数难以确定时，当生产过程机理不明、模型参数难以确定时，需要用过程辩识方法把数学模型估计出来。需要用过程辩识方法把数学模型估计出来。复杂的工业过程对象常由高阶非线性微分方复杂

27、的工业过程对象常由高阶非线性微分方程描述，求解困难。程描述，求解困难。机理建模得到的近似数学模型也需用试验测机理建模得到的近似数学模型也需用试验测量加以验证量加以验证2.3测试建模方法测试建模方法采用试验测定的原因采用试验测定的原因第49页/共117页2.3测试建模方法测试建模方法v对象特性的试验测定方法对象特性的试验测定方法v测定动态特性的时域法测定动态特性的时域法v测定动态特性的频域法测定动态特性的频域法v测定动态特性的统计相关法测定动态特性的统计相关法v最小二乘法最小二乘法第50页/共117页对象特性的实验测定方法对象特性的实验测定方法试验测定的方法试验测定的方法施加激励施加激励试验测定

28、方法的分类试验测定方法的分类（根据加入的激励信号和结果的分析方法分）（根据加入的激励信号和结果的分析方法分）n经典辨识法（非参数模型辨识）经典辨识法（非参数模型辨识）n时域法时域法n频域法频域法n相关法相关法n现代辨识法现代辨识法n最小二乘最小二乘n梯度校正梯度校正n极大似然法极大似然法由以时间或频率为自变量的试验曲线得到：由以时间或频率为自变量的试验曲线得到：1 1）数学模型结构）数学模型结构 2 2）模型参数）模型参数首先假定模型的结构；首先假定模型的结构；然后极小化模型与过程间的误差准则函数以然后极小化模型与过程间的误差准则函数以确定参数。确定参数。第51页/共117页（1）时域法）时域

29、法对被控对象施加对被控对象施加阶跃输入阶跃输入，测绘出对象输出，测绘出对象输出量随时间变化的响应曲线，或施加量随时间变化的响应曲线，或施加脉冲输入脉冲输入测绘出输出的脉冲响应曲线。测绘出输出的脉冲响应曲线。由响应曲线的结果分析，确定出被控对象的由响应曲线的结果分析，确定出被控对象的传递函数。传递函数。该方法测试设备简单，测试工作量小、应用该方法测试设备简单，测试工作量小、应用广泛，缺点是测试精度不高。广泛，缺点是测试精度不高。对象特性的实验测定方法对象特性的实验测定方法第52页/共117页（2）测定动态特性的）测定动态特性的频域方法频域方法对被控对象施加不同频率的正弦波，测出输对被控对象施加不

30、同频率的正弦波，测出输入量与输出量的幅值比和相位差，获得对象入量与输出量的幅值比和相位差，获得对象的频率特性，来确定被控对象的传递函数。的频率特性，来确定被控对象的传递函数。对象特性的实验测定方法对象特性的实验测定方法第53页/共117页（3）测定动态特性的）测定动态特性的统计相关法统计相关法对对被被控控对对象象施施加加某某种种随随机机信信号号或或直直接接利利用用对对象象输输入入端端本本身身存存在在的的随随机机噪噪音音进进行行观观察察和和记记录，录，可可以以在在生生产产过过程程正正常常运运行行状状态态下下进进行行，在在线线辨识，精度也较高。辨识，精度也较高。统统计计相相关关法法要要求求积积累累

31、大大量量数数据据，并并要要用用相相关关仪和计算机对这些数据进行计算和处理。仪和计算机对这些数据进行计算和处理。对象特性的实验测定方法对象特性的实验测定方法第54页/共117页2.3测试建模方法测试建模方法v对象特性的试验测定方法对象特性的试验测定方法v测定动态特性的时域法测定动态特性的时域法v测定动态特性的频域法测定动态特性的频域法v测定动态特性的统计相关法测定动态特性的统计相关法v最小二乘法最小二乘法第55页/共117页在在被被控控对对象象上上，人人为为地地加加非非周周期期信信号号后后，测测定定被被控控对对象象的的响响应应曲曲线线，然然后后再再根根据据响响应应曲曲线，求出被控对象的传递函数。

32、线，求出被控对象的传递函数。测定动态特性的时域法测定动态特性的时域法 图图2.15 测试响应曲线的框图测试响应曲线的框图1被测对象被测对象 2变换器变换器 3记录仪记录仪第56页/共117页1、输入信号选择及实验注意事项、输入信号选择及实验注意事项阶跃输入信号阶跃输入信号是时域法首选的输入信号是时域法首选的输入信号 对象的阶跃响应曲线对象的阶跃响应曲线直观地反映对象的动态特性直观地反映对象的动态特性直接来自原始的记录曲线，无需转换直接来自原始的记录曲线，无需转换实验也比较简单实验也比较简单从响应曲线中也易于直接求出其对应的传递函数从响应曲线中也易于直接求出其对应的传递函数测定动态特性的时域法测

33、定动态特性的时域法第57页/共117页1、输入信号选择及实验注意事项（续）、输入信号选择及实验注意事项（续）有时生产现场运行条件受到限制，不允许被有时生产现场运行条件受到限制，不允许被控对象的被控参数有较大幅度变化，无法测控对象的被控参数有较大幅度变化，无法测出一条完整的阶跃响应曲线，可改用出一条完整的阶跃响应曲线，可改用矩形脉矩形脉冲作为输入信号冲作为输入信号，得到脉冲响应后，得到脉冲响应后，再将其再将其换成一条阶跃响应曲线。换成一条阶跃响应曲线。测定动态特性的时域法测定动态特性的时域法第58页/共117页测定动态特性的时域法测定动态特性的时域法图图2.16 2.16 由矩形脉冲响应确定阶跃

34、响应由矩形脉冲响应确定阶跃响应 设设：u u1 1(t).(t).u u2 2(t)(t)作作用用下下的的阶阶跃跃响响应应曲曲线线为为y y1 1(t)(t)和和y y2 2(t),(t),则脉冲响应曲线为则脉冲响应曲线为:而阶跃响应曲线为而阶跃响应曲线为:第59页/共117页将时间轴按将时间轴按 t分成分成n等份等份在在区间内区间内在在区区间间内内依此类推，最后得到完整依此类推，最后得到完整的阶跃响应曲线。的阶跃响应曲线。测定动态特性的时域法测定动态特性的时域法由矩形脉冲响应确定阶跃响应方法（续）由矩形脉冲响应确定阶跃响应方法（续）第60页/共117页2、实验结果的数据处理、实验结果的数据处

35、理如如何何将将实实验验所所获获得得的的各各种种不不同同响响应应曲曲线线进进行行处处理理，以以便便用用一一些些简简单单的的典典型型微微分分方方程程或或传传递递函函数数来来近近似似表表达达，既既适适合合工工程程应应用用，又又有有足够的精度，这就是数据处理要解决的问题。足够的精度，这就是数据处理要解决的问题。测定动态特性的时域法测定动态特性的时域法微分方程微分方程传递函数传递函数（含结构和参数）（含结构和参数）第61页/共117页典型的工业过程的传递函数（自平衡过程）：典型的工业过程的传递函数（自平衡过程）：测定动态特性的时域法测定动态特性的时域法一阶惯性环节加纯延迟一阶惯性环节加纯延迟二阶惯性环节

36、加纯延迟二阶惯性环节加纯延迟n n个相同极点的个相同极点的n n阶惯性环节加纯延迟阶惯性环节加纯延迟 第62页/共117页v典型的工业过程的传递函数（无自平衡过程）：典型的工业过程的传递函数（无自平衡过程）：其传递函数中应含有一个积分环节其传递函数中应含有一个积分环节测定动态特性的时域法测定动态特性的时域法第63页/共117页选择何种传递函数？选择何种传递函数？测试者对被控对象的验前知识测试者对被控对象的验前知识本人的经验本人的经验与标准的一阶、二阶阶跃响应曲线比较与标准的一阶、二阶阶跃响应曲线比较选择传递函数阶次的原则：选择传递函数阶次的原则：低阶低阶数据处理简单，计算量也小，但准确程度较低

37、数据处理简单，计算量也小，但准确程度较低高阶高阶数据处理麻烦，计算量大，但拟合精度也较高数据处理麻烦，计算量大，但拟合精度也较高测定动态特性的时域法测定动态特性的时域法第64页/共117页以下讨论两种传递函数的参数确定方法以下讨论两种传递函数的参数确定方法带纯延迟的一阶惯性环节带纯延迟的一阶惯性环节带纯延迟的二阶惯性环节带纯延迟的二阶惯性环节测定动态特性的时域法测定动态特性的时域法第65页/共117页传递函数参数的确定：传递函数参数的确定：1）如果传函为带纯迟延的一阶惯性环节：）如果传函为带纯迟延的一阶惯性环节：K的确定的确定用参数用参数T、的确定的确定作图法作图法两点法两点法测定动态特性的时

38、域法测定动态特性的时域法第66页/共117页对于传函对于传函作图法确定参数作图法确定参数T、v简单，但误差大简单，但误差大v适用于适用于PID参数的工程整定参数的工程整定测定动态特性的时域法测定动态特性的时域法第67页/共117页对于传函对于传函两点法确定参数两点法确定参数T、y(t)归一化后的归一化后的y*(t)有有测定动态特性的时域法测定动态特性的时域法 根据阶跃响应曲线上两点数据根据阶跃响应曲线上两点数据(t1,y1),(t2,y2)确定参确定参数数T、（两个未知数两个方程）（两个未知数两个方程）第68页/共117页两点法结论：两点法结论：两点法确定参数两点法确定参数T、（续）（续）测定

39、动态特性的时域法测定动态特性的时域法第69页/共117页为了计算方便，可取，为了计算方便，可取，则可得则可得（2-47）测定动态特性的时域法测定动态特性的时域法第70页/共117页校验校验用用另另外外两两组组数数据据检检验验计计算算出出T、准准确确与与否否，即即 测定动态特性的时域法测定动态特性的时域法第71页/共117页传递函数参数的确定：传递函数参数的确定：2）如果传函为带纯延迟的二阶惯性环节）如果传函为带纯延迟的二阶惯性环节增益增益K值仍可按式值仍可按式2-43计算计算时间纯延迟时间纯延迟可根据阶跃响应曲线来确定；可根据阶跃响应曲线来确定；在考虑截去纯延迟部分后，阶跃响应在考虑截去纯延迟

40、部分后，阶跃响应y(t)为为无量纲形式无量纲形式y*(t);根据阶跃响应曲线上两个点的数据确定根据阶跃响应曲线上两个点的数据确定T1和和T2测定动态特性的时域法测定动态特性的时域法第72页/共117页测定动态特性的时域法测定动态特性的时域法第73页/共117页式式（2-39）截截去去纯纯延延迟迟并并化化为为无无量量纲纲形形式式后后，所对应的传递函数形式为所对应的传递函数形式为G(s)=(2-49)测定动态特性的时域法测定动态特性的时域法第74页/共117页与上式对应的阶跃响应为与上式对应的阶跃响应为或或 (2-50)同样取两点，联立方程组，求出同样取两点，联立方程组，求出T1T1，T2T2测定

41、动态特性的时域法测定动态特性的时域法第75页/共117页测定动态特性的时域法测定动态特性的时域法第76页/共117页求近似解：求近似解：(2-52)测定动态特性的时域法测定动态特性的时域法第77页/共117页表表2.1列列出出高高阶阶惯惯性性对对象象1/(Ts+1)n中中阶阶数数n与与t1/t2比值的关系比值的关系测定动态特性的时域法测定动态特性的时域法nt1/t2nt1/t210.3280.65820.46930.53100.7140.581150.62120.73560.651370.67140.75第78页/共117页表表中中t1和和t2应应是是y*(t)取取0.4和和0.8所所对对应应

42、的的t1、t2；中中的的时时间间常常数数T可可由由下式求得下式求得 测定动态特性的时域法测定动态特性的时域法第79页/共117页两点法求两点法求n n阶惯性环节的时间常数的步骤：阶惯性环节的时间常数的步骤：1 1）取阶跃响应曲线上两点）取阶跃响应曲线上两点 2 2）根据）根据t2/t1t2/t1的值确定系统的阶次和时间常数的值确定系统的阶次和时间常数 若若 一阶一阶 若若 二阶二阶 若若 根据根据t2/t1t2/t1确定阶次确定阶次n n测定动态特性的时域法测定动态特性的时域法第80页/共117页测定动态特性的时域法测定动态特性的时域法拟合非自平衡过程阶跃响应曲线拟合非自平衡过程阶跃响应曲线用

43、带延迟的积分环节来近似该响应曲线：用带延迟的积分环节来近似该响应曲线：其中：其中：=t2误差大误差大u阶跃输阶跃输入幅值入幅值第81页/共117页（2）用用式式（2-42）来来近近似似图图2.19响响应应曲曲线线、公式重写如下：公式重写如下：由响应曲线有由响应曲线有：延迟时间延迟时间=t1t1至至A之间以惯性环节作用为主之间以惯性环节作用为主故故T=t2-t1曲线稳态之后是积分环节作用为主曲线稳态之后是积分环节作用为主故故测定动态特性的时域法测定动态特性的时域法第82页/共117页测定动态特性的频域法测定动态特性的频域法对象特性的试验测定方法对象特性的试验测定方法测定动态特性的时域法测定动态特

44、性的时域法测定动态特性的频域法测定动态特性的频域法测定动态特性的统计相关法测定动态特性的统计相关法最小二乘法最小二乘法第83页/共117页被控对象的动态特性可用频率特性来描述被控对象的动态特性可用频率特性来描述与与传传递递函函数数及及微微分分方方程程一一样样，也也表表征征系系统统的的运运动规律动规律。在所研究对象的输入端加入某个频率的正弦波在所研究对象的输入端加入某个频率的正弦波信号，同时记录输入和输出的稳定振荡波形，信号，同时记录输入和输出的稳定振荡波形，可测得该被控对象的频率特性。可测得该被控对象的频率特性。测定动态特性的频域法测定动态特性的频域法第84页/共117页2.3.4测定动态特性

45、的统计相关法对象特性的试验测定方法对象特性的试验测定方法测定动态特性的时域法测定动态特性的时域法测定动态特性的频域法测定动态特性的频域法测定动态特性的统计相关法测定动态特性的统计相关法最小二乘法最小二乘法第85页/共117页原原理理：但但输输入入为为白白噪噪声声（一一种种随随机机信信号号）时时，输输出出与与输输入入的的互互相相关关函函数数与与脉脉冲响应函数成正比。冲响应函数成正比。2.3.4 测定动态特性的统计相关法测定动态特性的统计相关法 相关分析法求对象脉冲响应函数的方框图相关分析法求对象脉冲响应函数的方框图第86页/共117页2.3.5最小二乘法对象特性的试验测定方法对象特性的试验测定方

46、法测定动态特性的时域法测定动态特性的时域法测定动态特性的频域法测定动态特性的频域法测定动态特性的统计相关法测定动态特性的统计相关法最小二乘法最小二乘法第87页/共117页最小二乘的基本思想最小二乘的基本思想“未知量的最可能未知量的最可能的值，是使各次实际观测值和计算值之间的的值，是使各次实际观测值和计算值之间的差值的平方和为最小的值差值的平方和为最小的值”。2.3.5 最小二乘法最小二乘法第88页/共117页设设y是是一一根根金金属属轴轴的的长长度度，T是是金金属属轴轴的的温度，求轴长温度，求轴长y和温度和温度T之间的关系。之间的关系。第第一一步步:确确定定y和和T之之间间关关系系的的数数学学

47、模模型型的的类型和结构。类型和结构。第第二二步步:确确定定变变量量之之间间函函数数关关系系中中的的未未知知参数参数。2.3.5 最小二乘法最小二乘法第89页/共117页y和和T之间是线性关系。之间是线性关系。y=y0(1+T)=a+bTy00时金属轴的长度；时金属轴的长度；-膨胀系数。膨胀系数。2.3.5 最小二乘法最小二乘法第90页/共117页在加噪声干扰下，写成在加噪声干扰下，写成Zi=Yi（真值）（真值）+ni(随机误差随机误差)即即Zi=a+bTi+nii=1，2，,N2.3.5 最小二乘法最小二乘法第91页/共117页根根据据N次次（2）观观测测数数据据（TiZi）来来估估计出未知参

48、数计出未知参数a和和b的值。的值。2.3.5 最小二乘法最小二乘法第92页/共117页a和和b的的值值确确定定能能使使观观测测值值和和模模型型的的计计算算值之间误差为最小。每次观测误差为值之间误差为最小。每次观测误差为ni=zi-yi=zi-(a+b)相加起来的误差为相加起来的误差为 2.3.5 最小二乘法最小二乘法第93页/共117页常用误差平方和作为总误差。即常用误差平方和作为总误差。即 2.3.5 最小二乘法最小二乘法第94页/共117页按照求极值的原理，要使按照求极值的原理，要使J最小，只要将最小，只要将J分别对分别对a和和b求偏导数，令其等于零，求偏导数，令其等于零，a和和b的估计值

49、的估计值和和满足下面的条件满足下面的条件:2.3.5 最小二乘法最小二乘法第95页/共117页即即和和由下列方程组所确定由下列方程组所确定 2.3.5 最小二乘法最小二乘法第96页/共117页由式（由式（2-76）可解出）可解出 2.3.5 最小二乘法最小二乘法第97页/共117页假假 定定 一一 个个 变变 量量 y与与 一一 组组 n维维 变变 量量X=（x1,x2-xn）有线性关系。即）有线性关系。即y=2.3.5 最小二乘法最小二乘法第98页/共117页 N个参数的线性系统个参数的线性系统2.3.5 最小二乘法最小二乘法第99页/共117页假假设设用用 y(i)和和x1(i)、x2(i

50、)-xn(i),i=1,2,.,m表表示示实实测测数数据据。可可以以通通过过m个个线线性性方方程的方程组表示数据之间的关系。程的方程组表示数据之间的关系。i=1,2,.,m回归函数，回归函数，是回归系数。是回归系数。2.3.5 最小二乘法最小二乘法第100页/共117页 可以用矩阵形式表示如下可以用矩阵形式表示如下2.3.5 最小二乘法最小二乘法第101页/共117页2.3.5 最小二乘法最小二乘法第102页/共117页若若m=n,根据方程通过下式求解根据方程通过下式求解 只要只要即方阵即方阵X的逆存在，则能够唯一地求解的逆存在，则能够唯一地求解，表示，表示的估计值。的估计值。2.3.5 最小