系统仿真技术第2章 经典的连续系统仿真建模方法学  陈无畏  合肥工业大学机械与汽车工程学院  	 2.1 离散化原理及要求  问题：数字计算机在数值及时间上的离散性----被仿真系统数值及时间上的连续性? 连续系统的仿真，从本质上：对原连续系统从时间、数值两个方面进行离散化并选择合适的数值计算方法来近似积分运算   离散模型≈原连续模型？   相似原理      设系统模型为：                 ，其中u(t)为输 入变量，y(t)为系统变量；令仿真时间间隔为 h，离散化后的输入变量为       ，系统变量 为      ，其中    表示t=nh。       如果                 , 且 即                             ,  （对所有n=0,1,2,…）  则可认为两模型等价。  对仿真建模方法三个基本要求     （1）稳定性：若原连续系统是稳定的，则离散化后 得到的仿真模型也应是稳定的。     （2）准确性：有不同的准确性评价准则，最基本的 准则是：  绝对误差准则：  相对误差准则：   其中? 规定精度的误差量。  对仿真建模方法三个基本要求(续)       （3）快速性：若第n步计算对应的系统时间间隔 为               计算机由计算需要的时间为Tn，若        Tn=hn   称为实时仿真，Tn?hn称为超实时仿真         Tn?hn  称为亚实时仿真，对应于离线仿真  数值积分算法       对                    ，已知系统变量y的初始条件 要求y随时间变化的过程――初值问题       计算过程：由初始点                 的   欧拉法 对任意时刻tn+1  截断误差正比于    数值积分算法（续）        梯形法:  是隐函数形式。预报－—欧拉法估计初值，校 正－—用梯形法校正：   校正公式    预报公式       反复迭代，直到满足 经典的数值积分法分为两类：单步法与多步法 2.2 龙格库塔法  max.book118.com-库塔法基本原理        对       若令：       则有           的数值求解：称作“右端函数”计算问题。 在   附近展开泰勒级数，只保留   项，则有： 龙格-库塔法基本原理（续）       假设这个解可以写成如下形式：   其中  对    式右端的函数展成泰勒级数，保留h项， 可得：  代入，则有：   龙格-库塔法基本原理（续）       将(2)式与(1)式进行比较，可得：  四个未知数                            但只有三个方程， 因此有无穷多个解。       若限定             ，则 计算公式：  其中     龙格-库塔法基本原理（续）       若写成一般递推形式，即为：  其中     （1）截断误差正比于h3，称为二阶龙格-库塔法 （简称RK-2）。       （2）截断误差正比于h5的四阶龙格--库塔法（简称 RK-4）公式： 其中：   max.book118.com--库塔法的特点  1.形式多样性       例：                   非唯一解，可以得到许多 种龙格--库塔公式：                         (中点公式)  其中        各种龙格---库塔法可以写成如下一般形式：  其中  龙格--库塔法的特点(续）        式中各系数满足以下关系            s称为级数，表示每步计算右端函数f的最少次数。 可以证明，1阶公式至少要计算一次，2阶公式          ; ….；4阶公式           ；依此类推。有时为了某种特殊 需要，可以选择            的计算公式。  龙格--库塔法的特点(续） 2.单步法      在计算       时只用到    ，而不直接用             等 项。优点：存储量减小，可以自启动. 3.可变步长        步长h在整个计算中并不要求固定，可以根据精度 要求改变 ，但是在一步中，为计算若干个系数    ，则 必须用同一个步长h。   龙格--库塔法的特点(续） 4.速度与精度         四阶方法的h可以比二阶方法的h大10倍，每步计 算量仅比二阶方法大一倍，高于四阶的方法由于每步 计算量将增加较多，而精度提高不快。  max.book118.com  实时龙格－库塔法       实时仿真：要求仿真模型的运行速度往往与实际 系统运行的速度保持一致。        一般的数值积分法难以满足实时仿真的要求，这 不仅仅是因为由这些方法所得到的模型的执行速度较 慢，而且这些方法的机理不符合实时仿真的特点。        考虑系统  实时龙格－库塔法(续）      RK-2公式如下：       一个计算步内分两子步：  tn时刻：利用当前的un，yn计算k1----计算一次右端函数f需         。 tn+h/2时刻：应计算k2，尽管此时yn +1/2已经得到，但un +1则无法得到。(若对un +1也进行预报――加大仿真误差)。仿真执行延迟h/2――输出要迟后半个计算步距。 实时龙格－库塔法(续） 实时龙格－库塔法(续）     实时2阶龙格－库塔法：         tn时刻，应计算k1，利用当前的un，yn，需要     ； tn+h/2时刻，应计算k2，此时yn +1/2已经得到，un +1/2也 可得到，k2的计算就不会引入新的误差。计算一次右 端函数   需要      ，可实时输出yn +1。  实时龙格－库塔法(续） 2.3 线性多步法  max.book118.com  线性多步法基本原理        基本原理：利用一个多项式去匹配变量若干已知 值和它们的导数值。        设：                 时刻的                   和                已知;       预报：由                  和                   来计算       校正：若     也已知，由它们来计算   线性多步法(续)       采用的多项式具有以下形式(m阶)     其中：   是待定系数， 在     时刻，    ，可得到：                                                                       （2-1）   线性多步法(续)       由                 和                确定                 ， 需要m+1个独立方程。该m+1个方程可由以下 等式导出：                                                                            （2-2）                                                             线性多步法(续) 1、预报公式       令m=2k-1，从（2-
系统仿真技术_第2章 经典的连续系统仿真建模方法学.ppt