第二章常微分方程的初值问题教材.ppt

练习：二阶Runge-Kutta法处理二阶常微分方程 2.4 算法的稳定性 以欧拉法的一个简单扩展为例 对微分方程进行积分时，一个首要的考虑是所用算法的数值稳定性，也就是说，舍入误差或数值计算中的其它误差能被放大的程度。 为了启动上面的递推关系，还需要 y1 的值，这可以由泰勒展开获得 当把这一方法应用到下述问题 这个问题的解析解是 y=e-x 设解为指数形式 y=Arn，代入递推关系得 正根略小于1，它对应于我们要找的按指数减小的解，但是负根略小于-1，因而它对应于一个虚假的解 上述方程的解为 其大小随n增大并且在格点上逐点发生振荡。 线性差分方程的通解正是这两个指数解的一个线性组合。虽然可以精心安排初值 y0 和 y1， 使得当 x 值小时只呈现指数衰减的解，但是在递推过程中的舍入误差将会掺入一个小的“坏”解，它最后将增长到在解中占压倒地位。 一个好的经验规则是，每当积分一个随着迭代过程急剧衰减的解时，就应当小心不稳定性和舍入误差。 例 1. 强迫钟摆 一根长度为 l 钟摆被限制在一个垂直的平面内，在强迫外力 fd 和阻力 fr 的作用下振荡运动。钟摆的运动可以通过 Newton 方程来描述 其中 fg = -mg sinθ 是重力在运动方向的分力，a=l d2 θ/dt2 是沿切线方向的加速度， θ 是杆和垂线的夹角。 2.5 动力学中的有序和混沌 那么 Newton 运动方程就可以写成如下形式 假设强迫外 力为 , 阻力为 并且取 (l/g)1/2 为时间单位。 其中 该运动方程就可以化为一阶方程组 设 q=0.5,b=0.9, ω0 =2/3, t=100。 在这 种情况下,钟摆运动是一个有序的周期运动 杆和垂线的夹角、角速度随时间演化过程 角速度与夹角的轨迹 q=0.5, b=1.15, ω0 =2/3, t=1000，在这 种情况下钟摆运动是一个无序运动,呈现出分形的特征 杆和垂线的夹角、角速度随时间演化过程 角速度与夹角的轨迹 例 2. 自激振动——范德波尔方程 范德波尔描述非线性有阻尼的自激动振动系统 其中 μ 是一个正的小量 自激系统能将非振动的能源通过系统本身的反馈调节吸收进来，以补充被损耗的能量。自激系统的例子包括心脏等。 VDP方程不能准确的表示心脏振荡图形的细节，为了更准确的表示心脏的波动，FitzHugh对VDP方程进行了修改，提出了所谓的BVP方程，形式为 例 3. 化学振荡——BZ反应 CA, CB， CC， CD 分别为四种化合物浓度。 例 4. 二维粒子的运动 考虑一个单位质量的粒子,它在一个位势 V 中作二维运动。 设在 t 时刻,粒子的平面坐标为(x, y), 它的共轭动量为(px , py), 则 Hamilton 量的形式为 粒子的轨迹就由坐标和动量随时间的演化由如下 Hamilton 方程规定 该约束条件把粒子运动的轨道限制在四维相空间的三维流形上. 它是四个耦合的一阶微分方程组.这些方程使能量 E 守恒,即满足约束条件 在可分离变量的系统中,位势 V 具有 (x,y) 的可分离形式,即 V(x,y)=Vx(x)+Vy(y)，此时 Hamilton 方程可写为 可分离变量的系统 和 由此可见, x 方向和 y 方向的运动是相互不耦合的，每一个 Hamilton 量单独都是一个运动常数，其中 作为一个例子,考虑两个坐标方向的运动都是简谐运动的情况,即位势 对应的 Hamilton 方程为 和 该问题是精确可解的。 我们可以通过直接积分,求得满足初始条件 (x0, y0 , px0 , py0 ) 的解为 (x, px)平面和 (y, py) 平面的轨迹图 (x, y) 平面的运动轨迹 在 Henon-Heiles 位势系统中，位势的形式为 Henon-Heiles 系统 这个位势原来是由 Henon 和 Heiles 在研究恒星穿过星系的轨道时引入的。它具有三重对称性，并且在原点处位势为零,当坐标值大时位势无限增大。对于小于 1/6 的能量，轨道被限制在一个等边三角形之 内 Henon-Heiles 位势所确定的 Hamilton 方程为 该问题不能精确求解，所以必须用数值方法研究. 二维粒子的运动轨迹 (x, px ) 平面和 (y, py ) 平面的轨迹图 * 单击此处编辑母版标题样式 单击此处编辑母版标题样式 第二章 常微分方程的初值问题 本章要研究的物理问题：经典动力学中的有序和混沌 本章内容 4 简单方法 1 2 3 多步法和隐式法 龙格库塔法 稳定性问题 4 动力学中的有序和混沌 5 2.0 引子 常微分方程是物理学中经常碰到，以一维运动粒子为例 一般形式为M个藕合的一阶方程 实际问题往往涉及不止一类问题，例如对偏微分方程分离变量时