一阶常微分程的初值问题.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-