连续系统离散化处理基本方法

1、在数字计算机上对连续系统进行仿真时，首先遇到的问题是如何解决数字计算机 在数值及时间上的离散性与被仿真系统数值及时间上的连续性这一基本问题。从根本意义上讲，数字计算机所进行的数值计算仅仅是“数字”计算，它表示数 值的精度受限于字长，这将引入舍入误差；另一方面，这种计算是按指令一步一步进 行的，因而，还必须将时间离散化，这样就只能得到离散时间点上系统性能。用数字 仿真的方法对微分方程的数值积分是通过某种数值计算方法来实现的。任何一种计算方法都只能是原积分的一种近似。因此，连续系统仿真，从本质上是对原连续系统从 时间、数值两个方面对原系统进行离散化，并选择合适的数值计算方法来近似积分运 算，由此得

2、到的离散模型来近似原连续模型。如何保证离散模型的计算结果从原理上确能代表原系统的行为，这是连续系统数字仿真首先必须解决的问题。设系统模型为：y f(y,u,t),其中u(t)为输入变量，y(t)为系统变量；令仿真 时间间隔为h,离散化后的输入变量为u?h),系统变量为y?(tk),其中tk表示t=khc 如 果 u?(tk) u(tk) ,W(tk) y (tk), 即 eu(tk) U(tk) u(tk) 0 ,ey(tk) y?(tk) y(tk) 0 (对所有k=0,1,2, 九 则可认为两模型等价，这称为相似 原理(参见图)。ey (tk)0图相实际上，要完全保证eu(tk) 0,ey

3、(tk) 0是很困难的。进一步分析离散化引的误差，随着计算机技术的发展，由计算机字长引入的舍入误差可以忽略，关键是数值积分算法，也称为仿真建模方法。相似原理用于仿真时，对仿真建模方法有三个基本要求：(1)稳定性：若原连续系统是稳定的，则离散化后得到的仿真模型也应是稳定的。关于稳定性的详细讨论将在节中进行。(2)准确性：有不同的准确性评价准则，最基本的准则是：绝对误差准则：ey(tk)y?(tk) y(tk)相对误差准则：ey (tk)y?(tk) y(tk)y?(tk)其中规定精度的误差量。(3)快速性：如前所述，数字仿真是一步一步推进的，即由某一初始值y(t0)出发，逐步计算，得到y (tj

4、y (t2)， ，y (tk),每一步计算所需时间决定了仿真速度。若第 k步计算对应的系 统时间间隔为hk tk i tk，计算机由y (tk)计算y (tk i)需要的时间为Tk，则，若Tk=hk称为实 时仿真，Tk hk称为超实时仿真，而大多数情况是 Tk hk ,对应于离线仿真。连续系统数字仿真中离散化最基本的算法是数值积分算法。对于形如y f (y ,u,t)的系统，已知系统变量y的初始条件y(to) y0,现在要求y随时间变化的过程 y(t)。计算过程可以这样 考虑(参见图)：首先求出初始点y(t0) y0的f(t0, y0),微分方程可以写作：y(t) y0；0 f(t,y)出 (

5、)ff(Ly)图所示曲线下的面积就是y(t),由于难以得到f(y,u,t) 积分的数值表达式，人们对数值积分方法进行了长期探索，其中欧拉法是最经典的近似方法。欧拉法用矩形面积近似表示积分结果，也就是当t=t i时，y(t1)的近似值为y1 ：重复上述作法，当t t2时y2y(t2)yi& ti) f(ti, yi)所以，对任意时刻tk+i,有：令t iy i y(t i) y (t i t ) f(t , y )图数值积分法yiy(ti) yo t f(t, y)t h称为第 步的计算步距。若积分过程中步距不变h h,可以证明，欧拉法的截断误差正比于 h2。梯形法近似积分形式如式所示，令:为进

6、一步提高计算精度，人们提出了 “梯形法”t 1t h =h已知：t t时y(t )的近似值y，那么： i . y i y(t i) y 2hf(t ,y ) f(t i,y i)可见，梯形法是隐函数形式。采用这种积分方法最简单的预报校正方法是用欧拉法估计初值，用梯形法校正，即：i，:、y( i)y 2hf(t ,y ) f(t i,y()i)yiy h f(t ,y )式称作预报公式，采用欧拉法，式为校正公式，采用梯形法。用欧拉法估计一次 丫的值， 代入校正公式得到 y1的校正值丫(1)。设 是规定的足够小正整数，称作允许误差，若i=0, i +i=i称作第一次校正；i 1, i 1 2称作第二次校正；通过反复迭代，直到满足yi 11 yi 1,这时y。；是满足误差要求的校正值。上述方法是针对式所示的微分方程在已知初值情况下进行求解，因此也称为微分方程初值问题数值计算法，为统一起见，本书中称为数值积分法。连续系统数字仿真的离散化方法有两类，它们是数值积分方法和离散相似方法，本文讨论数值积分法。数值积分方法采用递推方式进行运算，而采用不同的积分方法会引进不同的计算误差，为了提高计算精度，往往会增加运算量。就同一种积分算法而言， 为提高计算精度，减小积分步距h