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

《连续系统离散化处理基本方法》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

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

2、方面对原系统进行离散化，并选择合适的数值计算方法来近似积分运算，由此得到的离散模型来近似原连续模型。如何保证离散模型的计算结果从原理上确能代表原系统的行为，这是连续系统数字仿真首先必须解决的问题。设系统模型为：.V=f(y,u,t)9其中为输入变量，y(0为系统变量；令仿真吋间间隔为//，离散化后的输入变量为列~)，系统变量为y(G)，其中G表示狄。如,y(tk)«y(D，即)-u(tk)^09ey(tk)=y(tk)-y(tk)^o(对所有h0，/，2,...),则可认为两模型等价，这称为相似原理(参见图2.1)。原连续模堺L/.仿真模彻

3、图2.1实际上，要完全保证4(~)=0，~(~)=0是很困难的。进一步分析离散化引的误差，随着计算机技术的发展，rt!计算机字长引入的舍入误差可以忽略，关键是数值积分算法，也称为仿真建模方法。相似原理用于仿真时，对仿真建模方法有三个基本要求：(1)稳定性：若原连续系统是稳定的，则离散化后得到的仿真模型也应是稳定的。关于稳定性的详细讨论将在2.4节中进行。(2)准确性：有不同的准确性评价准则，最基本的准则是：绝对误差准则：，么)

4、=

5、y么)—＞，么)

6、s5相对误差准则:

7、，即由某一初始值出发，逐步计算，得到yOjv(G),…,7(G),每一步计算所需时间决定了仿真速度。若第々步计算对应的系统时间间隔为么=Q+1—计算机由＞，匕)计算jv(g+i)需要的时间为7＞则，若称为实时仿真，7；＜~称为超实时仿真，而大多数情况是7＞么，对应于离线仿真。图2.2数值积连续系统数字仿真屮离散化最基本的算法是数值积分算法。对于形如y=/Cv，w，z)的系统，已知系统变量:v的初始条件=yu，现在要求:v随时间变化的过程计算过程可以这样考虑(参见图2.2):首先求出初始点＞，(/(＞)=凡的/(/9,y0),微分方程可以写作

8、：(2.1)图2.2所示曲线下的而积就是由于难以得到f(y，u，t)积分的数值表达式，人们对数值积分方法进行了长期探索，其中欧拉法是最经典的近似方法。欧拉法用矩形面积近似表示积分结果，也就是当时，＞，(^)的近似值为％:wh)三)’()+△，•/(，()，y())(2.2)重S上述作法，当z=z2时y2=y(t2)=yi+(t2-tx)-f(ti9yx)所以，对任意时刻&+/，有：人+1=y(~+1)三人+GK)(2.3)令~+1一称为第K步的计算步距。若积分过程中步距不变可以证明，欧拉法的截断误差正比于/z2。为进一步提高计算精度，人

9、们提山了“梯形法”。梯形法近似积分形式如式(2.4)所示，令：tK+l-tK=hK=h己知：，=~吋)’(~)的近似值yK，那么：1yK+i=X^+i)=八+y[，(，K，))+，(，K+1’yK+1)](2-4)可见，梯形法是隐函数形式。采用这种积分方法最简单的预报一校正方法是用欧拉法估计初值，用梯形法校正，即：(2.5)y^i}=y+去M/(，K，yK)+/(u；^)]式(2.6)称作预报公式，采用欧拉法，式(2.5)为校正公式，采用梯形法。用欧拉法估计一次的值，代入校正公式得到yK+1的校正值^了。设e是规定的足够小正整数，称作允

10、许误差，若/=0,/+1=1称作第一次校正；z=l,/+1=2称作第二次校正；通过反复迭代，直到满足这时是满足误差要求的校正值。上述方法是针对(2.3)式所示的微分方程在已知初值情况下进行求解，因此也称为微分方程初侪问题数值计算法，为统一起见，本书屮称为数值积分法。连续系统数字仿真的离散化方法有两类，它们是数值积分方法和离散相似方法，本文讨论数值积分法。数值积分方法采用递推方式进行运算，而采用不同的积分方法会引进不同的计算误差，为了提髙计算精度，往往会增加运算S。就同一种积分算法而言，为提高计算精度，减小积分步距计算量增大，影响系统运算速

11、度。因此，计算精度和速度是连续系统仿真屮常迂到的一对矛盾，也是数字仿真屮要求解决的问题之一。也就是说，选择合适的算法、合适的软、硬件环境，在保证计算精度的前提下，考虑怎样提岛仿真