2012数值分析 -lec26--常微分方程组解法.pptVIP

一阶常微分方程组的求解方法 考虑一阶方程组的初值问题 若采用向量的记号，记 则常微分方程组的初值问题可以表示为 这时四阶龙格－库塔公式具有形式 其中 这是单步法，利用节点 上的值 ， ，由上式顺序计 算 ，然后代入 即可求得 节点上的 。 朱立永 北京航空航天大学 数学与系统科学学院 Email: numerical_analysis@ Password:beihang 答疑时间：星期三下午2:00－5:00 答疑地点：主216 第二十六讲常微分方程组的解法 第七章 常微分方程初值问题的数值解法 前几节我们主要讨论了常微分方程的初值问题的数值解法。只要将 y 和f 改写为向量，那么前面提供的各种计算公式即可适用于一阶常微分方程组的初值问题。 考虑常微分方程组的初值问题： 解此方程组的Euler方法 考虑常微分方程组的初值问题： 常微分方程组解法 把F看成f，Y看成y，Y’＝F(t,Y)看成y’＝f(t,y)，前面介绍的所有的解y’＝f(t,y)的格式都可以直接推广应用于解常微分方程组Y’＝F(t,Y)； 对于隐式格式，一般需要通过非线性方程组迭代法来解得最终得数值解； 常微分方程组数值格式的相容性条件：即为格式具有至少一阶精度的条件 常微分方程组数值解的总体误差的定义形式是和对单个方程的格式是一样的，只是把对单个方程用的绝对值换成了对向量的范数； 常微分方程组解法 5. 常微分方程组数值格式的局部截断误差的定义：把常微分方程精确解带入数值格式所余下的部分称为所用格式的局部截断误差； 6. 常微分方程组数值格式的精度：如果格式的局部精度误差||Rn+1||=O(hp+1)，那么称格式是p阶精度的； 7. 常微分方程组数值格式的收敛性：其定义形式和单个方程一样； 8. 关于常微分方程组数值格式的精度和收敛性定理：和当方程的基本一样，只需把Lipschitz条件换成范数下的Lipschitz条件。 常微分方程组格式的绝对稳定性 模型方程： 绝对稳定性定义：和单个方程的定义一样。 格式的稳定性条件：为使每个方程 都稳定的条件 对常系数方程组的稳定性分析： 设A可对角化： 令 当A不可对角化时，只需对A的每个特征值作分析即可。 一般情况下： A的元素一般是t和Y的函数，因此A的特征值一般也是t和Y的函数，我们需要计算每个特征值的极值，然后讨论其在取极值下的稳定性条件。对所用特征值进行相似的讨论以得到最终的稳定性条件。 P207，例12、13 高阶常微分方程数值解法 可以把一个高阶常微分方程初值问题变成一个一阶常微分方程组。然后用解一阶常微分方程组的数值方法求解。 可以把一个高阶常微分方程组初值问题变成一个一阶常微分方程组。然后用解一阶常微分方程组的数值方法求解。 可以直接对高阶导数作差商近似来构造高阶常微分方程数值方法！ 刚性方程组的数值解法 刚性方程组 （Stiff Equations） 刚性方程组 刚性比： 定义： A-稳定区域 刚性方程组的数值解法 遗憾的是：早已证明，显式的方法不可能是A-稳定的，隐式的多步法也只有不超过2阶的方法才可能是A-稳定的，而高阶的隐式 单步法计算量又太大。 作业 课后习题23、25、27、28 Thanks for your attention! Any questions?! K (k (6 (5 (4 (3 (2 (1 (0 ( 1 2 3 4 5 6 1 1 1 1 1 1 1 -1 900 900 880 730 510 180