高等数值算法与应用十-wenjianyu.pdfVIP

高等数值算法与应高等数值算法与应用 （十十） AdvancedAdvanced NumericalNumerical Alggorithms Apppplications 计算机科学与技术系计算机科学与技术系 喻文健喻文健 SummarySummary -- OrdinaryOrdinary differentialdifferential equationsequations 求解未知函数求解未知函数 （（单自变量单自变量tt）） 阶数－－方程中未知函数求导阶数 高阶高阶ODEODE →→一阶阶ODEODE形成的方程组形成的方程组 （（一阶向量阶向量ODEODE ）） 初值问题 （IVP） StabilityStability ofof problemproblem （（初值的扰动对准确解的影响初值的扰动对准确解的影响）） Stable, Asymptotically stable, Unstable 简单的例子简单的例子： Eig(A)≤0 Im(hλ) 一般的非线性问题 局部稳定局部稳定：：Eig(g(J )) ≤0 ff Euler’’s method 整体截断误差 Difference Eqn Re(hλ) 局部截断误差局部截断误差 -22StabilityStability regionregion 00 Order p accuracy: StabilityStability, stabilitystability regionregion 1阶准确度（一般性） = y −u (t ) = (1 + hλ)y −y e λh k k k -1 k k -1 k -1 Wenjian Yu 2 OrdinaryOrdinary DifferentialDifferential EquationsEquations (II)(II) Wenjian Yu 3 内容概要内容概要 常微分方程数值解法－初值问题 向后欧拉法向后欧拉法, 梯形法梯形法 稳定性、准确度(考察模型问题, 一般问题) Runge-Kutta法 2级方法的推导级方法的推导,, 经典经典4阶方法阶方法，，自动自动R-K算法算法 刚性ODE-IVP 问题 两个例子两个例子、一般概念般概念 多步法 待定系数法推导, Adams公式, BDF公式思想 MatlabMatlab与应用与应用 Matlab命令, Lorenz吸引子 Wenjian Yu 4 优点：显格式，因此计算