经典：计算机仿真技术基础5

1、电力电子与电力传动实验室电力电子与电力传动实验室2021/8/61第第2章章 系统建模的基本方法与模型处理技术系统建模的基本方法与模型处理技术 3. 数学模型数学模型离散时间系统是同连续时间系统相对应的，它的输入、输出均是离散时间信号。同连续时间系统一样，它的数学模型也分为四种形式。 a)差分方程差分方程)() 1() 1()()() 1() 1()(11011kubkubmkbmkubkyakyankyankymmnn)(nm 设系统的输入序列u(k)，输出序列y(k)，系统输入、输出满足)()() 1(kukyky前向差分或) 1() 1()(kukyky后向差分一般的或)() 1()()

2、() 1()(101mkubkbkubnkyakyakymn(1)电力电子与电力传动实验室电力电子与电力传动实验室2021/8/62b)传递函数传递函数)()()()1 ()1(1110)1(111zUzbzbzbbzYzazazammmmnnnnnm nnnnmmmmzazazazbzbzbbzUzYzH)1(111)1(11101)()()(对式(1)两边取Z变换，若系统的初始条件为零，y(k)=0，u(k)=0，k0，则可得其中 Y(z)是序列y(k)的z变换；U(k)是序列u(k)的z变换。系统的离散传递函数为(2)电力电子与电力传动实验室电力电子与电力传动实验室2021/8/63c)

3、状态方程状态方程)()() 1()()() 1(kukkykukkDCXGFXX前面列出的模型只描述了系统的输入序列和输出序列之间的关系，为了进行仿真，通常要采用系统的内部模型，即离散状态空间模型。通常引进状态变量序列x(k)，构造系统的状态空间模型。一般的形式为(3)d)结构图结构图 离散系统的结构图表示和连续系统的相似，只要将每一个方框图内的连续系统的传递函数s换成离散系统的z函数即可。 电力电子与电力传动实验室电力电子与电力传动实验室2021/8/642.2.4 采样系统的数学模型采样系统的数学模型 随着计算机科学与技术的发展，人们不仅采用数字计算机而且利用计算机进行控制系统的分析与设计

4、，形成数字控制系统（或计算机控制系统），其控制器是由数字计算机组成的。它的输入变量和控制变量只是在采样点（时刻）取值的间断的脉冲序列信号，描述它的数学模型是离散的差分方程或离散状态方程；而被控对现象是连续的，其数学模型是连续时间模型，所以整个系统实际是一个连续离散混合系统。它主要有连续的控制对象、离散的控制器、采样器和保持器等几个环节组成，这就是采样系统的典型形式。 描述采样系统的模型就是连续离散混合模型，采样系统的框图如图1所示。采样控制系统里，采样开关和保持器是作为物理实体存在的。电力电子与电力传动实验室电力电子与电力传动实验室2021/8/65数 字 控制器保持器 控制对象r(t)e(t

5、)Te(kT)Tu(kT)u(t) y(t) 数字控制器把系统的模拟信号e(t)经过采样器及A/D转换器变成计算机可以接受的数字信号，经过计算机处理以后以数字量输出，再经过D/A转换器变成模拟量输入到被控对象。一般地，D/A转换器要将计算机第k次的输出保持一段时间，直到计算机第k+1次计算结果给它以后，其值才改变，因此通常把D/A转换器看成零阶保持器。电力电子与电力传动实验室电力电子与电力传动实验室2021/8/66严格地来讲，A/D转换器、计算机处理器、D/A转换器这三者并不是同步并行地进行工作，而是一种串行流水的工作方式，通常三者完成各自的任务所花费的时间并不严格相等，但如果三个时间的总和

6、与采样周期T相比可以忽略不计时，一般就认为数字控制器对控制信号的处理是瞬时完成的，采样开关是同步进行的，如果要考虑完成任务的时间的话，可以在系统中增加一个纯滞后环节。显然D/A转换的作用相当于一个零阶的信号重构器。 我们如果将连续的控制对象同保持器一起进行离散化，那么采样系统就简化为离散系统，采样系统的数学模型可以采用离散系统的四种数学模型表示 。 详细的内容会在采样控制系统仿真一节中介绍。电力电子与电力传动实验室电力电子与电力传动实验室2021/8/672.3非线性模型的线性化处理非线性模型的线性化处理 在实际工作中，纯粹的线性系统几乎是不存在的，说一个物理系统是线性的，实际上是看它的某些主

7、要物理性能可以充分精确地用一个线性模型加以描述而已。所谓“充分精确”是指实际系统与理想的线性系统之间的差别，相对于所研究的问题而言，已经小到忽略不计的程度。由于非线性模型的性质一般比线性模型的性质复杂的多，所以工程上常常用线性的关系近似地代替非线性关系。 在进行模型线性化处理时，一般把受控量与输入量之间的函数关系分成两类：一类函数的函数值与各阶导数值都是连续的，至少在工作范围内是连续的，称这类函数是光滑的，光滑函数的非线性是不严重的非线性，或可以在一个小的范围内用线性函数来近似；另一类函数是不光滑函数，不光滑函数是严重的非线性，一般来说不能用线性函数来近似，而只能视系统的物理性质来采取特定的线

8、性化方法。 电力电子与电力传动实验室电力电子与电力传动实验室2021/8/682.3.1微偏线性化方法1. 线性化的基本概念线性化的基本概念 所谓非线性数学模型的线性化就是对一个非线性系统的数学模型找出其稳定的平衡点，如果在工作过程中，代表系统属性的各物理量只在该平衡点附近产生微小的变化，非线性系统模型就能够以此平衡点为基础，表示成一个线性模型，关于线性系统的控制理论都能适用于该模型。这便是自动控制理论里关于小偏差线性化方法或称增量线性化方法的概念。2. 非线性数学模型的线性化的基本方法非线性数学模型的线性化的基本方法 对于非线性系统，当系统变量偏离工作点的偏差值很小时，由级数理论可知，若变量

9、在给定的工作区间内其各阶导数存在，便可在给定工作点的邻域内将非线性特性展开为泰勒级数，当偏差的范围很小时，可以忽略级数中偏差的高次项，得到只包含偏差的一次项的线性方程。电力电子与电力传动实验室电力电子与电力传动实验室2021/8/69 20000! 2)()()()(xxfxxfxfxxf0 xxxxxfxfxxf)()()(000 xxfxf)()(0对于光滑函数f(x)，常用微偏线性化的方法进行线性化处理，把f(x)在x0处展开成taylor 级数，如式(1)所示：当充分小的时候，式(1)可以写成或(1)(2)式(2)是一个线性模型，在工程上常常把)(),(xfxf函数近似地当作光滑函数，

10、在工作点处求函数的微分，完成线性化。由于是在工作点附近小范围内进行的，因而这样的线性化也被称为增量线性化或微偏线性化。 都连续的函数y=f(x)在工作点(x0，y0)处的微分是xxfyd)(d把dx、dy改写成x、y就可以得到式(2)，完成函数y=f(x)的微偏线性化。 只要电力电子与电力传动实验室电力电子与电力传动实验室2021/8/610nnnnxxfxxfxxfxxxfxxxf00022110201021),(),(xxx),(21nxxxf),(020100nxxxx对多变量非线性模型可以做类似的处理。对多变量非线性系统函数，把它在工作点成taylor级数可以得到处展开就可以把函数进行

11、微偏线性化。3. 求线性化微分方程的步骤求线性化微分方程的步骤(1)按物理和化学定律，列出系统的原始方程式，确定平衡点处各变量的数值。(2)找出原始方程式中间变量与其它因素的关系，若为非线性函数，在原平衡点邻域内，各阶导数存在并且是唯一的，则可进行线性化处理。(3)将非线性特性展开为泰勒级数，忽略偏差量的高次项，留下一次项，求出它的系数值。(4)消去中间变量，在原始方程式中，将各变量用平衡点的值加偏差量来表示。 电力电子与电力传动实验室电力电子与电力传动实验室2021/8/6114. 微偏线性化的几何意义)(xfy xxfy)(在希望的工作点用切线方程来近似代替曲线方程，是微偏线性化方法的几何

12、含义电力电子与电力传动实验室电力电子与电力传动实验室2021/8/612uxxaxxx21221sin0,例1 将下面的非线性系统在平衡点附近线性化。其中，输入u为常数。021 xx0arcsin2010 xuax，1.在线性化之前，首先需要确定平衡点的值。因为是，得到解：平衡点，所以2.将非线性特性展开为泰勒级数,忽略偏差量的高次项,留下一次项，求出它的系数值。2110222112222211111cos0000 xxxxxfxxfxxxxfxxfxxxxx电力电子与电力传动实验室电力电子与电力传动实验室2021/8/6133.将平衡点带入增量方程中21221arcsincosxxuxxxx

13、uaaxcos(arcsin104.写成标准的状态方程其中21xxx21xxx电力电子与电力传动实验室电力电子与电力传动实验室2021/8/6142.4 高阶模型的降阶处理高阶模型的降阶处理 在系统分析、设计和仿真中，常常会遇到一些复杂系统，这些系统的状态变量很多，阶次很高。例如电力系统，研究电磁暂态时的发电机系统等。 对于高阶系统进行仿真或设计是很麻烦的，从仿真的角度来讲，高阶系统的仿真要占用较多的内存和机时。从系统设计的角度来看，高阶系统的控制器往往十分复杂，有的甚至是不可能实现的。因此需要对高阶系统进行简化降阶，使其变得比较易于计算机和工程上实现，同时又要能在一定的精度范围内表现原系统的

14、特性。 所谓的模型简化，就是说为高阶复杂的系统准备一个低阶的近似模型，它在计算上、分析上都比原来的高阶系统模型简单，而且还可以提供关于原系统的足够多的信息。衡量一个模型的简化方法的可行性通常有四条标准：准确性、稳定性、简便型和灵活性。电力电子与电力传动实验室电力电子与电力传动实验室2021/8/615准确性：要求简化的模型与原型的主要特征一致，如主导极点一致，静态增益一致，频率响应与时间响应基本一致等。稳定性：要求简化的模型的稳定性与原型一致，而且具有相近的稳定裕量。简便性：要求从原型获得简化模型的过程简单,计算量小。灵活性：要求根据实际情况方便地进行调整，并得出所有侧重的简化模型。 通常以上

15、几个要求是难以同时满足的，有的方法准确性好，但计算量大；有的计算方便，但未必能保证稳定性等等，在实际中往往要综合考虑才行。 一般模型的简化技术分为两大类：一类是在状态空间模型上进行简化，称为时域模型简化法，另一类是在传递函数模型上进行简化，称为频域简化法。具体的简化方法在近几年不断涌现，在工程实际中都有应用。 电力电子与电力传动实验室电力电子与电力传动实验室2021/8/6162.5连续系统模型的离散化处理连续系统模型的离散化处理 在前面介绍采样系统的数学模型时，已经介绍了数字控制系统的模型。在仿真时可以先将连续的数学模型进行离散化得到离散的数学模型，然后再进行仿真。进行的每一步计算都是以这个

16、离散的模型为基础的，原来的连续系统的模型不再参与计算，这些结构上比较简单的离散系统，比较便于在计算机上求解。 系 统仿真模型建立仿真实验建模模 型计算机电力电子与电力传动实验室电力电子与电力传动实验室2021/8/617连续系统数学模型离散系统数学模型S域传递函数Z域传递函数微分方程差分方程连续系统状态方程离散系统状态方程连续系统结构图离散系统结构图电力电子与电力传动实验室电力电子与电力传动实验室2021/8/618基本思想基本思想:传递函数是控制工程领域中应用最为广泛的模型描述形式。对于给定的函数G(s)，设法找到s域到z域的某种映射关系，它将s域的变量s映射到z平面上，由此得到与连续系统传

17、递函数G(s)相对应的离散传递函数G(z)。2.5.1 替换法)(zfs )(sG)(| )()(zfssGzG即首先求出，然后直接把它带入就可以得出进而根据G(z)由Z反变换求得系统的时域离散模型差分方程，据此便可以进行快速求解。 电力电子与电力传动实验室电力电子与电力传动实验室2021/8/619S域传递函数Z域传递函数S=f(Z)(| )()(zGsGzfs根据Z变换定义，可得S域到Z域的最基本的映射关系是zTsezTsln1或 其中T为采样周期。电力电子与电力传动实验室电力电子与电力传动实验室2021/8/6201.简单替换法简单替换法)(zfs )(sGNkkkMrrrsasbsG1

18、01)()(sGNiiisAsG1)( 为了求出的函数关系，让我们先从传递函数开始。将展开成部分分式，设无重极点，则因为上式右侧各项具有相同的形式，所以只研究其中一项的转换对G(s)转换为G(z)就具有普遍的意义。（1）（2）电力电子与电力传动实验室电力电子与电力传动实验室2021/8/621设sAsG)()()()(tAutydttdy这个传递函数意味着模拟输入u(t)和输出y(t)之间有如下关系。即dtdy(t)Tnyny )() 1(用)()(t y ny)()(t u nu便可写出与上式相联系的差分方程 )()()() 1(1nAunynynyT（3）（4）（5）电力电子与电力传动实验

19、室电力电子与电力传动实验室2021/8/622ty(t)nn+1y(n)y(n+1)T电力电子与电力传动实验室电力电子与电力传动实验室2021/8/623)()()1()(1zAUzYTzzYT方程两边取Z变换，并进行整理得进一步得比较G(s) 与G(z)sAsG)() 1(1)()()(zTAzUzYzGTzs1TzssGzG1| )()(如果G(s)中的则简单替换法（Euler法）! 2/)( )1 (2TsTsezTs（6）（7）（8）电力电子与电力传动实验室电力电子与电力传动实验室2021/8/624)()()() 1(1nAunynynyT由(5)式可得)()(T)() 1(nynA

20、unyny（9） )()()(tytAudttdy由(4)式dttytAunTydttynTyTnyTnnTTnnT)1()1()()()()()() 1(（10） 可得对比(9)式、(10)式可见简单替换法的几何意义为几何意义几何意义电力电子与电力传动实验室电力电子与电力传动实验室2021/8/625tnn+1T)(ty)(nydttytAuTnnT)1()()()(tydttyTnnT)1()(电力电子与电力传动实验室电力电子与电力传动实验室2021/8/626例例1. 已知二阶系统的传递函数为231)()()(2sssUsYsG用简单替换法确定它的z传递函数和差分方程。 解：Tzs/1)

21、21)(1 ()23(21311)()()(222TTzTzTTzTzzUzYzG：做变换得进行z反变换得)()()21)(1 () 1()32()2(2nuTnyTTnyTny电力电子与电力传动实验室电力电子与电力传动实验室2021/8/627sAsG)(2.双线性替换法双线性替换法同上对于传递函数最基本的形式这个传递函数意味着模拟输入u(t)和输出y(t)之间有如下关系。即)()()(tAutydttdydttdyTnyny)( )() 1()( )() 1(21tynyny)( )() 1(21tununu利用（1） （2） 电力电子与电力传动实验室电力电子与电力传动实验室2021/8/

22、628ty(t)nn+1y(n)y(n+1)Ty(n+1)+y(n)/2电力电子与电力传动实验室电力电子与电力传动实验室2021/8/629便可写出与上式相联系的差分方程)() 1(2)() 1(2)() 1(1nunuAnynynynyT方程两边取Z变换，并进行整理得)() 1(2)()12()()21(zUzAzYTzzYT进一步得112)()()(zzTAzUzYzG比较G(s) 与G(z)可知sAsG)(双线性替换法（Tustin变换 ）112| )()(zzTssGzG112zzTs如果则（3） （4） （5） （6） 电力电子与电力传动实验室电力电子与电力传动实验室2021/8/6

23、30几何意义几何意义由(3)式得（7） )()()(tytAudttdy由(2)式dttytAunTydttynTyTnyTnnTTnnT)1()1()()()()()() 1(（8） 可得对比(7)式、(8)式可见简单替换法的几何意义为)() 1(2)() 1(2)() 1(1nunuAnynynynyT )()() 1() 1(2)() 1(nynAunynAuTnyny电力电子与电力传动实验室电力电子与电力传动实验室2021/8/631tnn+1T)(ty)(nydttytAuTnnT)1()()()(tydttyTnnT)1()(电力电子与电力传动实验室电力电子与电力传动实验室2021

24、/8/632这一变换也可以由下面的关系直接推得2/12/12/2/TsTseeezTsTsTs例例2. 利用Tustin变化求上例中的二阶系统的差分模型。 原系统为231)()()(2sssUsYsG则)2)(1(2)2(4)2)(1(2) 12(2)11)(2(3)11)(2(1)()()(22222TTzTzTTzzTzzTzzTzUzYzG进行反变换得到相应的差分模型)() 1(2)2()()2)(1(2) 1()2(4)2)(1(21)2(2222nuTnuTnuTnyTTnyTTTny电力电子与电力传动实验室电力电子与电力传动实验室2021/8/6333. 根匹配法根匹配法由控制理论

25、可知，连续系统的动态特性是由其传递函数的增益及零点、极点的分布情况所决定的，根匹配的基本思想就是构造一个相应于系统传递函数的离散传递函数，使两者的零点、极点相匹配，并且两者具有相同的终态响值，具体的做法是 ,假定线性系统的传递函数为)(,)()()()()(2121nmpspspsqsqsqsKsGnm构造离散系统的传递函数为)(,)()()()()(2121nmpzpzpzzqzqzqzKzGnmnmz(1)(2)电力电子与电力传动实验室电力电子与电力传动实验室2021/8/634), 2 , 1(nippii与), 2 , 1(miqqii与Tsez mnz其中，满足某种匹配关系将s平面上

26、的零、极点映射到z平面上）是因为当nm时，根据经典的根轨迹法可知，此时G(s)存在n-m个位于无穷远点的零点，可以认为n-m个无穷远零点位于负实轴上的无穷远处，则根据映射关系可知在z平面上必须再配n-m个位于原点的零点；Kz是根据G(z)与G(s)的终值相同的条件而确定的增益。 显然，按照上述思想构造的离散模型不仅具有良好的算法稳定性，而且保持了原系统的动态特性。（利用关系式如上所述，实现动态匹配往往从以下几个方面加以考虑：G(z)与G(s)具有相同数目的极点和零点；G(z)具有与G(s)的极点、零点相匹配的极点和零点；G(z)具有与G(s)的终值相匹配的终值1. 调节相位，使G(z)与G(s)的动态响应达到最佳匹配。 电力电子与电力传动实验室电力电子与电力传动实验室2021/8/635一般来说，在做系统的离散化时首先要假设原连续系统满足以下条件：系统是线性的；系统必须是可以进行拉普拉斯变换的；系统必须是渐进稳定的，并满足终值定理；终值必须是非零的。然后，通过以下的步骤来进行系统的离散化： 确定连续系统的传递函数G(s)；将G(s) 写成式(1)形式，以确定零、极点，即pi (i=1,2,m)，qj (j=1,2,n)；利用关系式Tsez 利用上一步求得的z平面上的零、极点写