第5讲常微分方程数值解ppt课件.ppt

《第5讲常微分方程数值解ppt课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、常微分方程数值解《数值分析》第五讲第五章：常微分方程数值解§5.1引言1、常微分方程与解为n阶常微分方程。如果函数在区间[a,b]内n阶可导，称方程满足方程的函数称为微分方程的解。则如为任意常数）为方程的解,一般称为方程的通解。如果则有为方程满足定解条件的解。第五章：常微分方程数值解方程的通解满足定解条件的解微分关系（方程）解的图示第五章：常微分方程数值解本教材重点讨论定解问题(初值问题）定解条件（初始条件）是否能够找到定解问题的解取决于仅有极少数的方程可以通过“常数变易法”、“可分离变量法”等特殊方法求得初等函数形式的解，绝大部分方程至今无法理论求解。如等等2、数值解的思想第五

2、章：常微分方程数值解（1）将连续变量离散为（2）用代数的方法求出解函数在点的近似值*数学界关注工程师关注如果找不到解函数数学界还关注：解的存在性解的唯一性解的光滑性解的振动性解的周期性解的稳定性解的混沌性解的分岔性……§5.2Euler方法第五章：常微分方程数值解第一步：连续变量离散化第二步：用直线步进·····1、Euler格式Euler格式例P106第五章：常微分方程数值解初值问题Bernoulli型方程常数变易法第五章：常微分方程数值解第五章：常微分方程数值解令将代入Euler格式计算得x精确值Euler方法Euler方法误差01.00000001.00000000.000

3、00000.11.09544511.10000000.00455490.21.18321601.19181820.00860220.31.26491111.27743780.01252680.41.34164081.35821260.01657180.51.41421361.43513290.02091940.61.48323971.50896630.02572660.71.54919331.58033820.03114490.81.61245151.64978340.03733190.91.67332011.71777930.044459311.73205081.7847708

4、0.0527200第五章：常微分方程数值解Euler值2、Euler格式的误差分析第五章：常微分方程数值解事实上Euler格式的每一步都存在误差，为了方便讨论算法的好坏，假定第n步准确的前提下分析第n+1步的误差，称为局部截断误差。定义1定义2第五章：常微分方程数值解Euler格式的误差即Euler格式具有一阶精度将在点Taylor展开的计算格式第五章：常微分方程数值解3、Euler格式与后退Euler格式Euler格式的值后退Euler格式的值隐式格式第五章：常微分方程数值解4、后退格式的精度局部误差分析的要求第五章：常微分方程数值解对照即后退Euler格式具有1阶精度第五章：

5、常微分方程数值解3、梯形格式得梯形格式Euler格式的值后退Euler格式的值第五章：常微分方程数值解梯形格式几何解释第五章：常微分方程数值解4、梯形格式的精度局部误差分析的要求第五章：常微分方程数值解梯形格式5、改进的Euler格式预测校正第五章：常微分方程数值解隐式格式为方便计算，一般用以下改进格式计算用改进格式计算例5.1的结果为x精确值Euler方法改进Euler方法Euler方法误差改进Euler误差01.00000001.00000001.00000000.00000000.00000000.11.09544511.10000001.09590910.00455490

6、.00046400.21.18321601.19181821.18409660.00860220.00088060.31.26491111.27743781.26620140.01252680.00129030.41.34164081.35821261.34336020.01657180.00171940.51.41421361.43513291.41640190.02091940.00218840.61.48323971.50896631.48595560.02572660.00271590.71.54919331.58033821.55251410.03114490.0033

7、2080.81.61245151.64978341.61647480.03733190.00402320.91.67332011.71777931.67816640.04445930.004846311.73205081.78477081.73786740.05272000.0058166第五章：常微分方程数值解第五章：常微分方程数值解6、预测精度的改进—两步Euler格式第五章：常微分方程数值解如果令第五章：常微分方程数值解整理得两步Euler格式具有2阶精度。因此可得具有2阶