一阶常微方程的初值问题课件.pptVIP

* Ordinary Differential Equations ODE 一阶常微分方程的初值问题： 节点：x1x2 … xn 步长 为常数 一 欧拉方法（折线法） yi+1=yi+h f(xi,yi) (i =0,1, …, n-1) 优点：计算简单。 缺点：一阶精度。 二 改进的欧拉方法 改进的欧拉公式可改写为 它每一步计算f(x,y)两次，截断误差为O(h3) 精确解： function [t,y] = Heun(ode,tspan,h,y0) t = (tspan(1):h:tspan(end)); n = length(t); y = y0*ones(n,1); for i=2:n k1 = feval(ode,t(i-1),y(i-1)); k2 = feval(ode,t(i),y(i-1)+h*k1); y(i) = y(i-1)+h*(k1+k2)/2; end 三 龙格—库塔法(Runge-Kutta) 欧拉公式可改写为 它每一步计算 f (xi,yi) 一次，截断误差为O(h2) 标准四阶龙格—库塔公式 每一步计算 f (x, y) 四次，截断误差为O(h5) 1/6 2/6 2/6 1/6