电子电子仿真第三章.pptx

第3章 系统仿真算法 ;大 纲;大 纲;第3章 系统仿真算法;3.1 数值积分法; 若被仿真系统可表示成一阶微分方程组或状态方程的形式。一阶向量微分方程及初值为 （3-1） 则上式在 处的形式上的连续解为 通过上式的演化，连续系统的数值解就转化为相邻两个时间点上的数值积分问题。 ;所谓数值解法，就是寻求初值问题式（3-1）的真解在一系列离散点 上的近似数值解 （3-1） 相邻两个时间离散点的间隔hn = tn+1 - tn，称为计算步距或步长，通常取hn = h为定值。 数值积分法的主要问题归结为对函数F(t，Y)的数值积分问题，即如何求该函数定积分的近似解。 ;初值问题 根据已知的条件y0，逐步地递推计算后续时刻的数值解yi（i = 1, 2, …） 共同特点：步进式 采用不同的递推算法，就出现不同的数值积分方法 采用不同的数值积分公式，对系统的求解精度、速度和稳定性均有不同的影响 ;;Euler法 已知微分方程 假定y(t)为其解析解，将展成泰勒级数： 则有 写成差分方程的形式 只要给定初始条件y0及步长h，就可根据 f(t0, y0)算出y1，再由y1算出y2，如此递推算出y3，y4，y5…。;Euler法的特点 方法简单，导出清晰，几何意义明显 计算量小 计算精度比较低 在实际中很少采用 【例一】 已知系统方程 ， 。设 ， 试用Euler法求其数值解。 ;2. 改进Euler法（预估 — 校正法） 在Euler法中，是用矩形面积来代替曲边梯形面积。 在改进Euler法中，改用直边梯形面积代替曲边梯形面积，则可提高精度。 曲边梯形的面积为 直边梯形的面积为 ;当h比较小时，以直边梯形面积取代曲边梯形的面积，可得 即 上式称为梯形公式或二阶隐式Adams公式。 公式右端隐含有待求量，故梯形法不能自行起步运算，而需有其他算法的帮助。 为了提高计算精度常常需要迭代运算。 ;为减小计算量，常常迭代一次就求得近似解。这样就可以得到改进的Euler法公式： 通常称这类方法为预估— 校正法。 实质上是采用了两点斜率平均值的结果。 将这种思想引申出加权方法，可得到龙格 — 库塔 （Runge-Kutta）法 。 ;3. 几个基本概念 （1）算法自启动 如Euler法 （2）单步法与多步法 单步法，求解yn+1只需要知道yn，故可自启动 多步法，求解yn+1不仅需要知道yn，还需要知道yn-1、yn-2 等，不可自启动 （3）显式与隐式 显式算法，求解yn+1不需要tn+1及其后时刻的数据 隐式算法，计算yn+1 的公式含有未知量yn+1;3. 几个基本概念 （4）截断误差 基于泰勒级数的精度分析 Euler法—泰勒展开式前2项，具有一阶精度 改进Euler法—泰勒展开式前3项，具有二阶精度 （5）舍入误差 由计算机字长引起的误差 与步长h成反比，步长越小，舍入误差越大 （6）初始误差 给定的初始值与真实值之间的差异;Runge-Kutta（简称RK）法的基本思想：用几个点上的y(t)的一阶导函数值的线性组合来近似代替在某一点的各阶导数，然后用Taylor级数展开式确定线性组合中各加权系数。 优点：既可避免计算高阶导数，又可提高数值积分的精度; 三阶RK法公式： 四阶RK法公式： ;四阶RK法程序框图 四阶RK法的特点 计算量大； 精度较高。 ;【例3-2】已知系统方程 取步长 h=0.1，试用RK4公式计算t=0.1，t=0.2时y的值。 解：将原系统方程转化成下列方程组： （1）计算所有变量的第一个RK系数k1(k11, k21) ;（2）计算所有变量的第二个RK系数k2(k12, k22) ;（3）计算所有变量的第三个RK系数k3(k13, k23) ;（4）计算所有变量的第四个RK系数k4(k14, k24) ;（5）计算时的近似值 用同样的方法，在y11，y21的基础上可以求得y12和y22的值。即为原系统中在t=0.2时的近似值。 ;单