第4节 连续系统的数字仿真通用算法——3.pdf

计算机仿真技术基础 李丽 lili.shu@shu.edu.cn 电机楼218室 连续系统的数字仿真通用算法 常微分方程求解介绍 常用数值积分方法☆ • 欧拉方法（Euler Method）☆ • 龙格-库塔方法（Runge-Kutta）☆ • 阿达姆斯方法（Adams) ☆ 步长的自动选择 数值积分法算法稳定性☆及收敛性 方程组及高阶方程的数值解法☆ 边值问题的数值解法 2 常微分方程求解 常微分方程(ordinary differential equation，ODE) ——是未知函数只含有一个自变量的微分方程 解析解：只有少数较简单和典型的常微分方程 （例如线性常系数常 微分方程等）可求出其解析解。 • 变系数常微分方程的解析求解就比较困难； • 一般的非线性常微分方程就更难求解析解。 近似解：在大多数情况下，常微分方程只能用近似法求解。 近似解法可分为两大类： 一类是近似解析法，如级数解法、逐次逼近法等； 另一类则是数值解法，它给出方程在一些离散点上的近似解。 3 常微分方程求解 在具体求解微分方程时，需要具备某种定解条件，微 分方程和定解条件合在一起组成定解问题。 定解条件有两种：一种是给出积分曲线在初始点的状 态，称为初始条件，相应的定解问题称为初值问题 ； 另一种是给出积分曲线首尾两端的状态，称为边界条 件 ，相应的定解问题则称为边值问题。 4 常微分方程求解 例1： 弹簧-质量系统的振动问题，如图1所示 x m x c o 图1 5 常微分方程求解 例1： 作一定的简化后，可表示为一个二阶常微分方程 d 2x k  x 0 2 dt m 式中，x是质量m离平衡位置（0点）的距离； t 是 时间 ；k是弹簧常数。 6 常微分方程求解 当弹簧在振动开始时刻 = 时的初始位置 = 和初速度 dx   x (t ) x dt 0 0 t t0 确定时，弹簧的