系统仿真技术_第2章 经典的连续系统仿真建模方法学

《系统仿真技术_第2章 经典的连续系统仿真建模方法学》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、系统仿真技术第2章经典的连续系统仿真建模方法学陈无畏合肥工业大学机械与汽车工程学院2.1离散化原理及要求问题：数字计算机在数值及时间上的离散性----被仿真系统数值及时间上的连续性?连续系统的仿真，从本质上：对原连续系统从时间、数值两个方面进行离散化并选择合适的数值计算方法来近似积分运算离散模型≈原连续模型？相似原理设系统模型为：，其中u(t)为输入变量，y(t)为系统变量；令仿真时间间隔为h，离散化后的输入变量为，系统变量为，其中表示t=nh。如果,且即,（对所有n=0,1,2,…）则可认为两

2、模型等价。u(t)hy(t)-+图2.1相似原理原连续模型仿真模型对仿真建模方法三个基本要求（1）稳定性：若原连续系统是稳定的，则离散化后得到的仿真模型也应是稳定的。（2）准确性：有不同的准确性评价准则，最基本的准则是：绝对误差准则：相对误差准则：其中规定精度的误差量。对仿真建模方法三个基本要求(续)（3）快速性：若第n步计算对应的系统时间间隔为计算机由计算需要的时间为Tn，若Tn=hn称为实时仿真，Tnhn称为超实时仿真Tnhn称为亚实时仿真，对应于离线仿真数值积分算法对，已知系统变量y

3、的初始条件要求y随时间变化的过程――初值问题计算过程：由初始点的欧拉法对任意时刻tn+1截断误差正比于数值积分算法（续）梯形法:是隐函数形式。预报－—欧拉法估计初值，校正－—用梯形法校正：校正公式预报公式反复迭代，直到满足经典的数值积分法分为两类：单步法与多步法2.2龙格库塔法2.2.1龙格-库塔法基本原理对若令：则有的数值求解：称作“右端函数”计算问题。在附近展开泰勒级数，只保留项，则有：龙格-库塔法基本原理（续）假设这个解可以写成如下形式：其中对式右端的函数展成泰勒级数，保留h项，可得：代入

4、，则有：龙格-库塔法基本原理（续）将(2)式与(1)式进行比较，可得：四个未知数但只有三个方程，因此有无穷多个解。若限定，则计算公式：其中龙格-库塔法基本原理（续）若写成一般递推形式，即为：其中（1）截断误差正比于h3，称为二阶龙格-库塔法（简称RK-2）。（2）截断误差正比于h5的四阶龙格--库塔法（简称RK-4）公式：其中：2.2.2龙格--库塔法的特点1.形式多样性例：非唯一解，可以得到许多种龙格--库塔公式：(中点公式)其中各种龙格---库塔法可以写成如下一般形式：其中龙格--库塔法的特

5、点(续）式中各系数满足以下关系s称为级数，表示每步计算右端函数f的最少次数。可以证明，1阶公式至少要计算一次，2阶公式;….；4阶公式；依此类推。有时为了某种特殊需要，可以选择的计算公式。龙格--库塔法的特点(续）2.单步法在计算时只用到，而不直接用等项。优点：存储量减小，可以自启动.3.可变步长步长h在整个计算中并不要求固定，可以根据精度要求改变，但是在一步中，为计算若干个系数，则必须用同一个步长h。龙格--库塔法的特点(续）4.速度与精度四阶方法的h可以比二阶方法的h大10倍，每步计算量仅比

6、二阶方法大一倍，高于四阶的方法由于每步计算量将增加较多，而精度提高不快。2.2.3实时龙格－库塔法实时仿真：要求仿真模型的运行速度往往与实际系统运行的速度保持一致。一般的数值积分法难以满足实时仿真的要求，这不仅仅是因为由这些方法所得到的模型的执行速度较慢，而且这些方法的机理不符合实时仿真的特点。考虑系统实时龙格－库塔法(续）RK-2公式如下：一个计算步内分两子步：tn时刻：利用当前的un，yn计算k1----计算一次右端函数f需。tn+h/2时刻：应计算k2，尽管此时yn+1/2已经得到，但un

7、+1则无法得到。(若对un+1也进行预报――加大仿真误差)。仿真执行延迟h/2――输出要迟后半个计算步距。实时龙格－库塔法(续）RK-2的计算流程实时龙格－库塔法(续）实时2阶龙格－库塔法：tn时刻，应计算k1，利用当前的un，yn，需要；tn+h/2时刻，应计算k2，此时yn+1/2已经得到，un+1/2也可得到，k2的计算就不会引入新的误差。计算一次右端函数需要，可实时输出yn+1。实时龙格－库塔法(续）实时RK-2公式计算流程2.3线性多步法2.3.1线性多步法基本原理基本原理：利用一个多

8、项式去匹配变量若干已知值和它们的导数值。设：时刻的和已知;预报：由和来计算校正：若也已知，由它们来计算线性多步法(续)采用的多项式具有以下形式(m阶)其中：是待定系数，在时刻，，可得到：（2-1）线性多步法(续)由和确定，需要m+1个独立方程。该m+1个方程可由以下等式导出：（2-2）线性多步法(续)1、预报公式令m=2k-1，从（2-2）式得到如下方程组：（2-3）将其写成矩阵形式：（2-4）其中上标p表示预报。其解为：（2-5）由于为常数阵，其逆存在，Z向量中的各元素为已知值，因而d向量的各