第二章数字仿真的计算方法.pptVIP

电力系统计算机仿真 福州大学 电气工程与自动化学院电力系 邵振国 第二章 数字仿真的计算方法 2.1 数值积分方法 2.11 问题描述 2.12 单步法 2.13 多步法 2.14 误差分析 2.2 相关软件的应用 2.3 离散相似法 2.4 面向结构图的数字仿真 2.1.1 问题描述 一阶常微分方程的初值问题 要计算出解函数 y(x) 在一系列节点 a = t0 t1… tn= b 上的近似值 yi≈ y(ti) 基本概念 节点 ti 步长 h = ti+1 – ti 误差 e = yi – y(ti) 2.1.2 单步法 从初值y0开始，依次求出y1 、 y2 ，后一步的值 yi+1 只依靠前一步的 yi； yi+1= yi+h·Q(ti，yi，h)； 增量函数：Q(ti，yi，h) Euler（欧拉）方法 基本思路：用切线分段逼近曲线 y(ti+1) = y( ti) +h·f (ti，yi) y(ti+1)：真实值； yi+1：数值计算值； 基本算法：yi+1= yi +h·f (ti，yi) 2.1.2 单步法－欧拉法 e1 e2 yi+1= yi +h·f (ti，yi) t y h 从数值微分构造欧拉法 欧拉法例题 例：求解常微分方程初值问题，取步长h = 0.1，计算到 t = 0.5。 解：yi+1= yi +h·f (ti，yi)= yi +0.1×