第五章快速数字仿真法

第五章快速数字仿真法第五章快速数字仿真法本章第一节到第三节介绍了三种常用的方法，以便根据情况选用。这些方法不仅在连续系统仿真时可以使用，也可以用于连续控制器在计算机上的离散实现。本章还在第四节介绍了计算机控制系统的仿真。5.1增广矩阵法5.2替换法5.3零极点匹配法5.4计算机控制系统仿真5.1增广矩阵法一、基本思想假定一个连续系统的状态方程为

(5.1.1)这是一个齐次方程，它的解是

(5.1.2)已知可以证明：如果取前五项，则计算精度与四阶龙格-库塔法相同。这就是说，如果被仿真的系统是一个齐次方程，在选定计算步距为n以后，若只取的前五项，则有5.1增广矩阵法

(5.1.3)由于都可以在仿真前计算出，所以(5.1.3)式所示的递推计算公式中右端的系数项可以事先求出，而仿真计算就变成每次只做一个十分简单的递推推算。但是，实际的物理系统模型大多是一个非齐次方程，即 (5.1.4)其中为系统的控制量，假定它是一个单输入系统。根据控制理论可知，(5.1.4)式的解为t≥0(5.1.5)5.1增广矩阵法显见，求解(5.1.5)式的非齐次方程时，它的解，除了一个自由项之外，还有一个强制项：。由于的任意性，(5.1.5)式一般不容易求解。但是，对于某些特殊的输入函数，如果能将控制量增广到状态变量中去，使(5.1.5)式这样的非齐次方程，变成一个齐次方程(5.1.1)式，就可以避免计算复杂的强制项，而利用类似(5.1.3)式的计算方法。二、典型输入函数时的增广矩阵假定被仿真的系统为

(5.1.6)5.1增广矩阵法其中A为维矩阵，即表示有n个状态变量。(1)阶跃输入时设则定义第n+1个状态变量为故可得增广后的状态方程即输出方程为5.1增广矩阵法(2)斜坡输入设则定义5.1增广矩阵法因此，系统增广后的状态方程为初始条件为5.1增广矩阵法(3)指数输入设定义则有故系统增广后的状态方程为5.2替换法一个连续物理系统最常见的数学表现形式就是s域的传递函数。替换法的基本思想就是，设法找到s域（连续域）的某种对应关系，然后将中的变量s转化成变量z，由此得到与系统传递函数相对应的离散系统脉冲传递函数（脉冲传递函数即采样系统输出脉冲序列的Z变换与输入脉冲序列Z变换之比），进而获得进行数字仿真用的递推算式，以便在计算机上求解计算。5.2替换法5.2.1简单替换法s域与z域的基本关系是，式中T为采样周期，即仿真计算时的计算步长；或者，这是一个超越函数，不可直接用来进行替换。实际上必须要寻找其他近似的表示关系。一种最简单的替换关系可以从一阶差分方程中得到。用传递函数表示系统，在时域内可以用一个微分方程来表示，例如，系统

(5.2.1)5.2替换法的时域表示为

(5.2.2)假若导数计算用下述差分来近似：或间写成 (5.2.3)则微分方程(5.2.2)式即等价为下述差分方程：5.2替换法

或 (5.2.4)当微分方程中导数的计算采用(5.2.3)式差分表达式时，称为向后差分法。当然，导数的近似计算还可以采用向前差分法，此时，导数按下述差分来近似计算：方程(5.2.2)式可等价为下述差分方程：5.2替换法或

(5.2.5)现对差分方程(5.2.4)进行Z变换，则是

(5.2.6)比较(5.2.6)与(5.2.1)式，有

(5.2.7)

或

(5.2.8)5.2替换法关系式(5.2.7)或(5.2.8)就是一种最简单的替换方法。它表明，如果将传递函数中的用(5.2.8)式替换，就可以将s域中的传递函数变换为z域中脉冲传递函数，从而可以得到递推算式(5.2.4)式。这种替换相当于向后差分法。从(5.2.4)式可见，这种向后差分法也就是数值积分法中的超前欧拉法。若对(5.2.5)式做z变换，则可得要下述脉冲传递函数：

(5.2.9)

将(5.2.9)式与(5.2.1)式比较，有5.2替换法 (5.2.10)

或(5.2.11)关系式(5.2.10)或(5.2.11)也是一种最简单的替换方式，这种替换式相当于向前差分法，即数值积分中的欧拉法。上述两种替换法比较简单，但局限性很大，实际工程中很少采用。下边对两种变换作简单的讨论。实际上替换式(5.2.8)及(5.2.11)均可以看做是s平面与z平面之间的相互映射。对(5.2.11)式来说，由于，设，所以5.2替换法对于z平面上的单位圆(及稳定域)，，故即等于 (5.2.12)因此将z域中的单位圆映射到s域平面上，正好是以 为圆心，以为半径的一个圆，如图5.1.1(a)所示。这就是说z平面上的单位圆，按(5.2.11)式变换，它将是s平面 上以为半径的一个圆。5.2替换法反过来说，s平面上只有部分面积通过(5.2.11)式才能映射到z平面的单位圆内。显然，若一个系统是稳定的，其极点分布如图5.1.1(a)所示，通过(5.2.11)变化后有两个极点能映射到在z平面的单位圆内；而极点则映射到z平面的单位圆外。这样，原来为稳定的系统，通过(5.2.11)式的替换，仿真模型变得不稳定了，因此这种仿真模型失真太大。由(5.2.12)式还可 以看出，若要使稳定，就要求增大半径，即减 小计算步长T，从而增加了计算工作量，故采用(5.2.11)式不适合快速数字仿真。5.2替换法 图5.2.1简单替换法的映射关系5.2替换法

由(5.2.8)式的替换关系，可以推得，s平面的左半部将映射到z平面单位圆的局部范围内，如图5.1.1(b)所示。为说明这一点，(5.2.8)可以改写为

所以 (5.2.13)显然，若，，则；若，

则(5.2.13)右端分母大于分子，即；5.2替换法若，则(5.2.13)右端分母小于分子，即。这表明，s平面右半部全部映射在该圆之内。由于这个小圆位于z平面的单位圆内，所以原连续系统是稳定的，利用(5.2.8)式的替换关系得到的数学模型也一定是稳定的，并且与计算步长T无关。但是，由于整个s平面的左半部不是映射到z平面的整个单位圆内，而只是映射到其中一个小圆内，显然与准确z变换的结果相差较大，从而使仿真模型失真较为严重。例5.2.1现给定一个二阶系统的传递函数，使用两种替换，求仿真数学模型。计算步长为T=1s。5.2替换法解：(1)利用前差替换公式(5.2.11)式，得对该式做Z反比变换

或5.2替换法(2)利用后差替换公式(5.2.8)式，得Z反变换后得5.2替换法从上述计算结果可见，同一结果采用不同的计算公式，所得到的数学模型差别是很大的，不仅差分方程系数不同，而且形式也常常是不同的。5.2.2双线性变换上述简单替换法使用起来问题较多，比较适用的替换法是双线性变换法（又称Tustin法）。这种替换关系可以从s到z变量准确的映射关系推得。根据定义，z变量与s变量的关系是

或者(5.2.14)5.2替换法(5.2.14)式可以展成无穷级数

如果取该级数的第一项作为它的近似，则

(5.2.15)

或(5.2.16)数学上称这种变换为双线性变换。利用(5.2.15)就可以把传递函转换为脉冲传递函数。5.2替换法应当指出，这种变换关系相当于数值积分中的梯形法。事实上，积分环节(5.2.1)式对应的微分方程形式，若采用梯形积分法可得

(5.2.17)对(5.2.17)式的差分方程进行Z变换得可见，积分环节从梯形积分法所推得的递推公式与从双线性变换所得的递推公式是相同的。5.2替换法采用双线性变换的一个和重要的优点是，在s域里稳定的传递函数，通过这种变换得到的脉冲传递函数也一定是稳定的。因为这种变换将s域左半平面准确地映射到z域的单位圆内。实际上，将代入关系式(5.2.16)式中，得5.2替换法由该式可知，若，即，则；若，则；若，则。这就是说，采用(5.2.16)式映射，s域左半平面映射到z域的单位圆内，如图所示，所以若原系统稳定，那么利用这种变换得到的脉冲传递函数也一定是稳定的。

图5.2.2双线性变换的映射关系5.2替换法例5.2.2对例5-1传递函数

利用双线性变换求其仿真数学模型，计算步长仍取。解：将替换关系(5.2.15)式代入式中，则得5.2替换法反变换后可得5.2替换法对任意传递函数采用这种变换，它的脉冲传递函数即为

(5.2.18)由(5.2.18)式可见，与是同阶的。这种方法在概念上是很清楚的。实际运算时可能很麻烦，为此对常用的典型环节，与先求出与之间的各系数关系式，列成表格，在使用时查找是很方便的。此外目前许多程序都具有这种变换功能，在实际应用中可供选用。另外，将这部分工作编成计算机程序，交付计算机来完成，也不是一件难事。5.2替换法5.2.3状态方程的双线性变换双线性变换公式不仅可以方便地用于传递函数，同时也可用于系统的状态方程。若系统的状态方程为

(5.2.19)式中分别为维列向量；A为维矩阵；B为维矩阵；C为维矩阵；为维矩阵。对(5.2.19)式左边进行拉氏变换，得5.2替换法该式又可写成

(5.2.20)式中T为计算步长。现将用(5.2.15)式代替，则得进一步整理可得

5.2替换法 (5.2.21)5.2替换法

式中(5.2.22)

(5.2.23)

从(5.2.22)式及(5.2.23)式中消去，可得 即5.2替换法将该式代换(5.2.21)式中的的第二项I，则得

(5.2.24)将(5.2.16)式中第二式做Z变换并将(5.2.24)式代入，则得

(5.2.25)

5.2替换法由(5.2.25)式可见，该式相当于一个离散状态方程的输出方程的Z变换，若令

(5.2.26)

则以(5.2.26)式系数组成的离散状态方程如下：

(5.2.27)5.2替换法

进一步研究证明，方程(5.2.27)式的脉冲传递函数就是(5.2.25)式。而(5.2.27)式即为(5.2.19)式的等价差分方程组，并且各矩阵均为系数状态方程(5.2.19)式的系数矩阵和计算步长的函数，它可在计算机上迭代求解。下面我们对双线性变换的一些性质作简单的讨论。前面已经说明了，若原连续系统是稳定的，采用双线性变换所得离散数学模型也是稳定的。现在再进一步研究其稳定特性。对单位阶跃输入的稳态响应可以由(5.2.24)式求得5.2替换法利用(5.2.22)式及(5.2.23)式，可得5.2替换法原连续系统的稳态值可以从(5.2.20)式令求得，即

显然，连续系统的稳态值与离散数学模型的稳态值是相同的。这说明采用双向性变换并不会改变系统的稳态值。采用双向性变换不仅可以保证系统的稳定性及稳态值不变，同时也具有一定的精度。这精度主要是指离散数学模型的频率特性或输出序列与连续系统的频率特性或输出序列是相近的。为了说明这一点，利用例双线性变换后的系统来进行计算。对该题易得到5.2替换法将代入后即可求得脉冲传递函数的频率特性为

将代入上式得若令，则可以得到的频率特性。连续系统的频率特性可将代入例的中即可求得。比较这两种频率特性可见，在低频段相差较小，当超过0.8rad/s时，两个频率特性相差较大，但频率在高频段，又重复对称出现。5.2替换法 这表明在一定的频率范围内双线性变换是有一定精度的，但在某频段内，双线性变换的频率特性会发生一定的畸变。产生这种畸变的主要原因是，通过(5.2.15)式变换后，s平面上无限长的虚轴被压缩成平面上 的单位圆，等效的说被压缩在之内（，为计算步长），即产生频率翘曲。应当指出，根据(5.2.15)式进行变换，虽也将s平面左半部映射到z平面的单位圆内内，但它与的映射关系并不相同，即按照(5.2.15)式将s平面上的一点映射到z平面上的位置与按映射位置并不相同。5.3零极点匹配法零极点匹配法是获得连续传递函数等价离散数学模型的一种简单有效的方法。时间域函数的采样运算就是直接将s域的函数通过

(5.3.1)关系映射到z平面。零极点匹配法（又称零极点映射法）不仅通过(5.3.1)将极点映射到z平面，而且还将的零点映射到z平面。通常，传递函数的极点多于零点数，此时可以认为该传递函数在无穷远处有一定数目的零点，这些零点通过(5.3.1)式映射关系被映射在z=1处。

5.3零极点匹配法若给定的连续系统传递函数为

(5.3.2)式中:分别为重复零极点的实部与虚部。若，即极点数多于零点数个。按上述意思，通过(5.3.1)式映射可得下述脉冲传递函数：

(5.3.3)

5.3零极点匹配法式中是离散传递函数的增益，可以根据z域频率响应在某些关键频率上能较好地等价连续频率响应的原则来选择。通常对低通滤波器性质的传递函数或一般控制系统来说，它可以按直流()增益相等的原则来选取。归纳起来，零极点匹配法的步骤如下：(1)将的所有极点按映射到z平面，如的实极点数在处，则有极点；倘若有 复数极点，则有极点。(2)用同样的方法将全部有限零点按映射到z平面。(3)将所有在处的零点映射为。

5.3零极点匹配法(4)选择离散数学模型的增益，选择的方法是在给定的关键频率处使连续系统的增益与离散系统的增益相等。若按直流增益相等选择，则有

例现将零极点法用于前述例题，求差分方程的表达式。

解：该传递函数没有限零点，他的两个极点为

5.3零极点匹配法对计算步长T=1s，经过映射后，在z与内的两个极点为

由于没有限零点，故需补充两个在z=-1处的零点，所以

可以通过直流增益相等的原则来求得。因为

所以5.3零极点匹配法故可得=0.209，最后脉冲传递函数为

(5.3.4)该式可以转化为下述差分方程：

应当说明，将无穷远处的零点映射在z=-1处，这样的处 理，理由并不充分，主要的根据是在双线性变换 中就在分子上引入了零点（z+1）。这样做的好处，是有助于消除标准Z变换中产生的频率混叠现象。5.3零极点匹配法当然，如若认为无穷远处的零点是负实数，则由于，故可以将其看作是映射在z=0处的零点。此时可以将(5.3.3)式中用来代替，但经验表明，这种处理方法精度较差。此外，为了使离散模型与原连续系统更为接近，常常还可以将无穷远处的零点映射到（0，-1）之间某处，此时可以将(5.3.3)式中须用来代替，式中。这样在中将有两个参数及，进行实际仿真计算，则是很简单的。5.3零极点匹配法对于某些传递函数，其直流增益可能等于零或无穷大，无法利用直流增益相等的原则来确定传递函数的。此时可以通过加入适当类型的输入信号，利用使其稳态值相等的方法选择值。例如根据上述步骤可得到

由于在s=0处有一零点，稳态增益等于零，加阶跃输入信号时，终值将等于零。为使终值5.3零极点匹配法而由可得

若是稳态值相等，则应选择为

5.3零极点匹配法最后可得离散数学模型为当然，对这种系统还可以不依据稳态增益来选取，而依据某个特殊频率时的增益相等的原则来确定。5.4计算机控制系统仿真计算机控制系统是由离散部分（数字计算机或数字控制器）和连续部分（保持器或数模转换器以及被控对象）两部分组合而成。如图所示。

图5.4.1计算机控制系统组成图5.4计算机控制系统仿真这类系统的模型与第四章所述连续系统离散化模型形式相同，都是差分方程，所以从仿真方法来讲两者都是相同的。但是要注意的是连续系统离散化模型是一种近似模型，因为在建立模型过程中人为的加入了采样开关和保持器，所以它只能是近似描述原来的连续系统。而计算机控制系统仿真的特点是，数字控制器是实际存在的离散量运算，采样开关也是实际存在的不是人为虚拟的。采样周期Ts是确定的，它与连续部分离散化虚拟采样周期T是不同的。下面来说明如何来考虑这两者之间的关系。一、采样周期及计算步距5.4计算机控制系统仿真将图所示的计算机控制系统用传递函数形式表示如图所示。 图5.4.2计算机控制系统图 图中——被控对象的传递函数； ——保持器传递函数； ——数字控制器的脉冲传递函数。5.4计算机控制系统仿真是该计算机控制系统实际存在的采样开关的采样周期，这是计算机控制系统的一个重要参数。通常它是根据被控对象的反应快慢事先选定的。比如，当被控对象的阶跃响应曲线如图所示时，若选择数字控制器为P，I，D形式，那么可选为0.3L左右。对于图像所示的计算机控制系统仿真，一般有两种情况：(1)如果可以求出，而且只要求计算输出量，不要求计算中间变量，那么，可以用同样的采样间隔，分别写出数字控制器的差分方程，然后进行数字仿真。此时，采样周期，也就是计算步距

。5.4计算机控制系统仿真(2)如果不太好求出，或者不仅要求计算输出量，还要求计算出中间的其他状态变量，或者被控对象中有非线性环节，因此在做数字仿真时必须将连续部分分成几部分来分别计算。此时，就要在每个连续部分前分别加上虚拟的采样器和保持器。

图5.4.3被控对象的阶跃响应曲线5.4计算机控制系统仿真为了使计算有一定的精度，这些虚拟采样器的采样间隔T（也就是计算步距）就不能选得过分大，一般来讲。此时，在做仿真计算中，就应对连续部分与离散部分分别计算，而不能用同一步距。为了仿真计算的方便，通常选取T使，其中N为正整数。此时对整个计算机系统做仿真计算时，仿真程序将分为两个循环，一个是对连续部分，一个是对离散部分。由于，要对连续部分做N次计算后，才对离散部分做一次计算，如图所示。 图中， ——连续部分的计算步距了； ——离散部分的采样周期； ——要求的仿真时间。5.4计算机控制系统仿真

图5.4.4计算机控制系统数字仿真程序的流程图5.4计算机控制系统仿真二、计算机控制系统的仿真方法现以图所示的计算机控制系统为例进行说明。要实现计算机控制系统仿真，一般是计算离散部分的输出，其输入为，用作为连续部分的输入，计算其输入，这一部分的仿真方法和连续系统完全相同，只要在确定连接矩阵W的系数时，把去掉，看成是一个阶跃信号加在系统上即可。(1)数字控制器模型

5.4计算机控制系统仿真数字控制可以有各种控制规律，既可以有不同的脉冲传递函数或差分方程，但建立模型的方法基本相同，这里先以常用的PID（比例、积分和微分）控制规律为例来说明建立模型的方法。

(5.4.1)

式中：——比例系数； ——积分时间常数； ——为分时间常数。现若用Z变换方法把变换为，则可得5.4计算机控制系统仿真

(5.4.2)式中

由(5.4.2)式进行反变换，可得差分方程如下：

(5.4.3)

(5.4.3)式及所求的数字控制器的一种模型形式。它与连续部分离散化所得的差分方程形式相同。因此计算机控制系统的仿真就是要研究高阶差分方程的仿真方法。5.4计算机控制系统仿真(2)差分方程仿真设系统闭环脉冲传递函数的一般形式为

(5.4.4)其差分方程为

解上式即可求出采样时刻的输出。5.4计算机控制系统仿真若需要求出控制器的输入和输出信号时，则可分别求 出控制器部分及连续部分的差分方程。对于 每一个采样点，这两部分各计算一次，然后对控制器的信号进行综合，以便得到下一次计算控制器的输入。如此循环迭代下去，便可逐步求出各采样时刻的响应。以上两种算法所得到的结果是相同的，但后一种算法的运算时间要长一些。以上两种算法，最后都归结为高阶差分方程。下面介绍用计算机解这种方程的方法及其程序实现。假定高阶差分方程具有如下形式：

5.4计算机控制系统仿真

(5.4.5)式中： ——输出变量； ——输入变量。由(5.4.5)式可知，在解此差分方程时，要用到计算时刻以前若干个采样时刻的输入值和输出值，这可在内存中设置若干个单元，见数据存储起来，以便在计算时使用。下面首先讨论输出变量的处理方法。由(5.4.5)式可知，输出变量x的阶数n为，因此，需要在内存中设置n个单元以存放计算时刻以前n个采样时刻的输出量。这些单元的安排如图5.4.5(a)所示。5.4计算机控制系统仿真

图5.4.5变量存储单元的安排5.4计算机控制系统仿真在计算中，现在时刻前n个输出量 分别从，单元中取出。在每一次计算中其操作顺序都是“取出——平移——存入”。这使得“存入”，是指计算所得的现时刻的输出量由存入到第n单元中去。对于输入变量也可采用上述相似的方法进行处理。根据(5.4.5)式，E的阶数为m，在内存中设置m+1个单元存放E的数据。其安排如图5.4.5(b)所示。5.4计算机控制系统仿真在计算中，现时刻前m个输入变量 分别从单元中取出。在每一次计算中，其操作顺序都是“存入——取出——平移”。这里的“存入”，是指计算所得到的现时刻的输入存入到第m+1单元中去。差分方程仿真程序：#include<stdio.h>#include<bios.h>#include<stdlib.h>voidmain()5.4计算机控制系统仿真{ intI,j,n;intN,M,K; double*X,*E,*A,*B,e,x;clrscr();printf(“PleaseinputN,M,K,E:”);

scanf(“%d%d%d%lf”,&N,&M,*K,&e);A=malloc(N*sizeof(double));B=malloc(M*sizeof(double));X=malloc(N*sizeof(double));E=malloc(M*sizeof(double));for(i=0;i<N;i++){printf(“InputA[%d];”,i);scanf(“%lf”%A[i]);X[i]=0;}5.4计算机控制系统仿真

for(i=0;i<K;i++) {printf(“InputB[%d];”,i);scanf(“%lf”%B[i]);E[i]=0;}E[M-1]=e;for(i=0;i<K;i++){x=0for(j=0;j<K;j++)x=x+A[j]*E[M-j-1];5.4计算机控制系统仿真for(j=0;j<M;j++) x=x+B[j]*E[M-j-1];for(j=0;j<K-1;j++)X[j]=X[j+1];for(j=0;j<M-1;j++)

E[j]=E[j+1];X[K-1]=x;printf(“%d%lf\k”,i,x);}bioskey(1);}

5.4计算机控制系统仿真三、纯时延环节的数字仿真我们现在研究纯时延环节的数字仿真方法，我们将看到纯时延环节其实是一种特殊的差分方程描述的系统，因此可以利用上面讨论的方法解决纯时延环节的仿真。假设纯时延环节的输入量为，输出量为，延迟时间为，则时延环节的传递函数

(5.4.6)再假设采样周期为，且

(5.4.7)5.4计算机控制系统仿真式中，为正数部分，为小数部分。则有

(5.4.8)对(5.4.8)式作Z变换

(5.4.9)再作Z反变换，得到差分方程

(5.4.10)5.4计算机控制系统仿真若为的整数倍，则

其仿真方法同差分方程一样，如图所示

图5.4.6纯时延环节变量存储单元的安排利用BASIC语言的仿真程序为：Q(C0+1)=UY=Q(1)FORI=1TOC0Q(I)=Q(I+1)NEXTI5.4计算机控制系统仿真当时，则因，因此我们可以将的值取在和两个时刻之间，故可以利用线性插补公式求得 (5.4.11)

增设一个内存单元存放，故可取个内存单元，如如图所示

图5.4.7纯时延环节变量存储单元的安排5.4计算机控制系统仿真利用BASIC语言的仿真程序为：Q(C0+2)=UY=(1-C2)*Q(2)+C2*Q(1)FORI=1TOC0+1Q(I)=Q(I+1)NEXTI四、采样周期改变引起仿真模型的变化

为了减少计算量，加快仿真的速度，对计算机控制系统进行数字仿真时，通常会遇到这样的情况，即原来的比较小，现在希望用一个比较大的来进行仿真。这样要对离散部分的仿真模型进行修改，但当给出的离散部分仿真模型不是含时，须将原z域零极点通过映射关系对应到s域中，再采用映射到z域中找到改变为后新的零极点。5.4计算机控制系统仿真下面通过一个实际的例子来说明一下当采样周期改变时，如何改变仿真模型。假设在一个数字驾驶仪中有一个数字校正环节，它的采样周期，它的脉冲传递函数为

(5.4.12)若用来进行数字仿真，那么它的仿真模型可以用差分方程来表示，即

(5.4.13)5.4计算机控制系统仿真今要求仿真时，那么仿真模型应当怎样变化呢？假设当时，可以表示成，它可以看做是 的仿真模型，为了使时的仿真模型与 时的原模型具有相同的特性，首先要求与 两个脉冲传递函数在s平面上映射具有相同的零极点，同时稳态值相同。已知在z平面上的一个极点，同时有一个极点，它们映射到s平面上，有一个极点，一个零点。5.4计算机控制系统仿真根据零极点相配原则，当时，有

故

(5.4.14)现在还有一个增益系数要加以确定。它可以根据原型与模型具有相同的稳态值这一原则来确定。假设在的输入端加一个单位阶跃函数，则根据终值定理，输出的终值为

5.4计算机控制系统仿真而仿真模型的终值为

若要求有，则可得因此可得改变采样周期后的仿真模型为

(5.4.15)5.4计算机控制系统仿真相依的差分方程为

(5.4.16)在Matlab语言中有指令d2d()，该指令可以完成由于采样周期的改变而引起的模型变化。该指令的调用格式为：sys1=d2d(sys,Ts) sys1=d2d(sys,[],Nd)调用说明：该指令将原离散时间系统的采样时间改变为采样时间为T