4数值求解常微分方程课件.pptx

1.引言 2.欧拉法 3.龙格-库塔法 4.线性多步法 5.高阶微分方程; 方程中包含未知函数的导数或者微分的方程叫做微分方程,求解微分方程必须附加某种定解条件.微分方程和定解条件一起组成定解问题,定解条件分为初始条件(初值问题)和边界条件(边值问题)两种.未知函数为一元的微分方程叫做常微分方程,未知函数为多元函数,叫做偏微分方程.微分方程中导数的最高阶叫做微分方程的阶.本章主要讨论一阶常微分方程.; 常微分方程在物理学中具有重要意义,牛顿第二定律就表示为一组二阶微分方程.前面提到的著名的三体问题就是由九个二阶常微分方程组成的方程组;在量子力学和电动力学中,由薛定谔及麦克斯韦方程分离变量得到本征值问题也是二阶常微分方程的求解问题.;现考虑上述一阶常微分方程初值问题.在常微分方程中，我们已经掌握了一些典型方程的解析解法。但解析方法只能解一些特殊类型的方程,求解从实际当中得出来的微分方程只能用数值方法求近似解.初值问题也就是由已知点(x0,y0)求解在某些点x上满足一定精度的y的近似解。;即使对有些函数可以求出解析解.但实际应用时,仍需求它在某一点的值,如果解析解的形式比较复杂,还是需要数值方法来求得.例 y=1-2xy y(0)=0 其解析解如下,要计算具体的函数值y,还需要积分的方法,不如直接用数值法求解.;2.1 欧拉法(折线法) 若将函数y(x)在点xi处的导数y’(xi)用(向前)差商来表示，即; 上式就是欧拉公式,它有很明显的几何意义:我们知道初值问题中的微分方程的解是平面上的一条曲线，这条曲线上任意点(x，y)的切线的斜率y’(x)等于函数f(x，y)在这点的值，而初值问题的解是过点(x0,y0)的一条特定的曲线。 ;过点(x0,y0)以f(x0,y0)为斜率的方程为 当x=x1时得下式,取y1作为解y(x1)的近似解 过(x1,y1)以f(x1,y1)为斜率的方程为;若已求得点(xn,yn),过这点,以f(xn,yn)为斜率的直线 当x=xn+1时得: 取y(xn+1)≈yn+1;这样每个xn求出对应数值解yn;还可使用泰勒展开证明欧拉公式： 对y(xn+1)在xn处泰勒展开有 舍去余项可得欧拉公式;例1:用欧拉法求初值问题的数值解(取h=0.1)。;12;f(x,y)=x-y*y; float y=0; float h=0.1; float a=0.0; float b=1.0;;隐式欧拉公式 若将函数y(x)在点xi处的导数y’(xi)用(向后)差商来表示，即 欧拉公式变为:;两步欧拉公式 为了改善精度,将函数y(x)在点xi处的导数y’(xi)用中心差商来表示，即 欧拉公式变为:误差正比于h3,是二阶精度;2.2 改进的欧拉法 欧拉法虽然形式简单，计算方便，但比较粗糙，精度也低。特别当y=y(x)的曲线曲率较大时，欧拉法的效果更差。为了达到较高精度的计算公式，对欧拉法进行改进，将在一点(xi，yi)的切线斜率f(xi，yi)用两点的平均斜率来代替，即 ; 代入欧拉公式得 ; 不难发现，欧拉公式是关于yi+1的显式，只要已知yi，经一次计算可立即得到yi+1的值;而改进的欧拉公式中的yi+1以隐式给出，且yi+1含在函数f(xi+1 ,yi+1)中，因此，通常用迭代法求解。具体做法是：先用欧拉公式求出一个y(0)i+1作为初始近似，然后再用改进的欧拉公式进行迭代，即; 直到满足 则 ;2.3 预估校正法 改进的欧拉公式在实际计算时要进行多次迭代，因而计算量较大。所谓预估校正法，就是先用欧拉公式算出yi+1的预估值y(p)i+1，然后再用改进的欧拉公式式进行一次迭代便得到校正值y(c)i+1，即;预估： 校正： 并取 ;在实际计算时，还常常将上式写成下列形式： ;23; 例2:求解初值问题 ; 其解析解为 ;26;27;28; 具有二阶以上精度的方法称为龙格-库塔法， 做为一种训练，我们给出一种直接推导法，再给出一种比较容易推广到高阶的推导法，希望同学们自己体会这里介绍的方法和思路！ ;30;31;32;33;34;35;由四阶龙格库塔公式计算得到的例2的结果;; 我们已经知道，初值问题等价于积分方程，即; 为此，我们自然可以想到，若用更高次的插值多项式来代替f(x，y)，则所得公式的精度会更高。这就是线性多步法的起源思想。 本章前面介绍的方法称为单步法，因为在计算yi+1时，只用到前面yi的值。而对于线性多步法是要利用前面已经算出的若干个值yi-k，…，yi-1，yi来求yi+1。;