chapt-7常微分方程数值解法幻灯片.pptVIP

§7.4 高阶常微分方程（组）的数值解法 10 subroutine rkutta(a,b,n,m,y,yy) 20 dimension hh(5),y(m),y1(m),f(m) 30 yy(mm) 40 x0=a 50 h=(b-a)/n 60 hh(1)=h/2;hh(2)=hh(1) 70 hh(3)=h;hh(4)=h;hh(5)=h/2 80 do i=1,n 90 x1=x0+(i-1)*h 100 x=x1 110 do j=1,m 120 y1(j)=y(j);yy(j)=y(j) 130 end do 140 do k=1,4 150 call fun(x,y,f) 160 x=x1+hh(k) 170 do l=1,m 180 yy(l)=yy(l)+f(l)*hh(k+1)/3 190 y(l)=y1(l)+f(l)*hh(k) 200 end do 210 end do 220 do j=1,m 230 y(j)=yy(j) 240 end do 250 end do 260 end subroutine hrkutta 练习题 7-1 取步长h=0.2再用向前Euler法和改进的Euler法以及四阶Ronge-Kutta法求解初值问题 7-2 取步长h=0.1再用四阶Ronge-Kutta法求区间[0, 1]内的解初值问题 上机实习题 7-3 取步长h=0.1，0.01和0.001，再用四阶Ronge-Kutta法求x=1处时如下解初值问题的近似解： 7-4 取步长h=0.25再用四阶Ronge-Kutta法求区间[0, 5]内的解初值问题数值解。 第7章 常微分方程及方程组的数值解法 常微分方程的定解问题一般分为初值问题和边值问题，而边值问题常常可以化为初值问题来求解，所以本章只讨论初值问题。 设在X—Y平面有： 当 求 。 在工程问题中，只有极少有解析解，大多数只能求数值解。即在区域： 求y(xi)的近似解： 。 (7-1) §7.1 常微分方程的数值法 本节研究常用的数值方法。先说明几个符号的含义： 一、Taylor级数展开法 把yi+1在yi处展开成Taylor级数，则有 因为 所以 O(h2)为h2的同阶无穷小，略去高阶项得如下的近似关系： §7.1 常微分方程的数值法 这就是常微分方程数值解法中的向前欧拉(Euler)公式。 二、数值微分法 这种方法就是利用差商近似代替导数，即 (7-2) §7.1 常微分方程的数值法 略去二阶导数项，用 代替 ，得 利用它们求问题(7-1)式数值的方法，分别称为欧拉(Euler)法，后退欧拉法，中点法。 三、数值积分法 将(7-1)式在区间 内积分得 (7-3) §7.1 常微分方程的数值法 对上式右端的积分使用左端矩形积分公式，则有： 于是(7-1)式的数值解法 但如果上式右端的积分使用的是梯形积分公式，则有 从而(7-1)式的数值解法也改为 这就是改进的欧拉公式。 (7-4) (7-5) §7.1 常微分方程的数值法 四、几个基本概念 （1）单步法：总是由前一点yi推算后一点yi+1，只需要一个点的值。除此之外，还有多步法，如Adams法。 （2）显示法：yi+1只出现等式一边。若同时出现在右边，则为隐式法。 （3）截项误差： 当h趋于零时，它是h2的同阶无穷小，（7-2-1）误差阶数为2。如果y(x)是一次式，则O(h2)?0,此时欧拉法是精确的（直线逼近直线），同时也说明欧拉法是一阶代数精度。 §7.2 一阶常微分方程的数值法 让我们研究(7-1)式形式的一阶常微分方程的数值求解方法。 7.2.1 欧拉折线法 一、算法分析 二、几何意义 参见图7-1中的折线f01f11f21f31，欧拉折线法命名由此而来。容易看出，欧拉折线法的数值解离理论曲线越来越远。折线f01f11f21f31为后面要讨论的欧拉二次法。 (7-6) §7.2 一阶常微分方程的数值法 图7-1 欧拉折线法及改进的欧拉法的几何意义 §7.2 一阶常微分方程的数值法 三、通用程序框图与程序设计 subroutine euler comp f(x0,y0) x0=a x=x0+h do i=1,n+1 y=y0+h*f(x