数值分析8-01.ppt

湘潭大学数学与计算科学学院1第八章常微分方程数值解§1引论湘潭大学数学与计算科学学院2解析法求解常微分方程的初值问题很多时候解析解求不出来,如(1.1)湘潭大学数学与计算科学学院3常微分方程的初值问题(1.1)湘潭大学数学与计算科学学院4若存在常数L0使得不等式则称f(x,y)在G上关于y满足Lipschitz条件，而式中的常数L称为Lipschitz常数.一切在G上关于y满足Lipschitz条件的连续映射f所构成的集合记为，而相应的初值问题（1.1）构成的问题类记为.湘潭大学数学与计算科学学院5定理1中的任何初值问题在[a,b]上有连续可微的解存在并且惟一.定义1初值问题(1.1)称为在[a,b]上是适定的，当时初值问题湘潭大学数学与计算科学学院6定理2中的任何初值问题在[a,b]上是适定的.以上各定理的证明在常微分方程的教材上都已经给出.定理3(Bellman不等式)设上的非负连续函数，则当时,必有(1.2)湘潭大学数学与计算科学学院7证明先设并记所以得得湘潭大学数学与计算科学学院8从而得到我们证明我们有引用上面已证明的结果，得到故不等式(1.2)仍成立.因此，命题得证.湘潭大学数学与计算科学学院9(1.3)定理4(离散的Bellman不等式)设则必有证明令湘潭大学数学与计算科学学院10由假设知，因而有或即由此即得□湘潭大学数学与计算科学学院11§2Euler方法一、Euler方法二、误差分析三、Euler方法的收敛性和稳定性湘潭大学数学与计算科学学院121、Euler方法记：因为：(等距剖分)(积分方程)令：有：湘潭大学数学与计算科学学院13——Euler方法又称Euler折线法.可得递推关系：湘潭大学数学与计算科学学院14欧拉方法的几何意义：h步长Euler方法的几何意义湘潭大学数学与计算科学学院152、误差分析有：精确值时，湘潭大学数学与计算科学学院16其中：湘潭大学数学与计算科学学院17几何分析：Euler公式的误差湘潭大学数学与计算科学学院18整体截断误差：由：湘潭大学数学与计算科学学院19从而有：有：湘潭大学数学与计算科学学院20于是便得Euler方法的整体截断误差界(*)湘潭大学数学与计算科学学院21定理5设f(x,y)属于F且关于x满足Lipschitz条件，其Lipschitz常数为K，则成立.并且有估计式(*)湘潭大学数学与计算科学学院22说明Euler方法的整体截断误差与h同阶。湘潭大学数学与计算科学学院23注意:对于Euler方法隐式Euler方法等价于积分方程：微分方程：有：称为隐式Euler方法.则：湘潭大学数学与计算科学学院27改进的Euler方法等价于积分方程：微分方程：湘潭大学数学与计算科学学院28有：湘潭大学数学与计算科学学院29称为改进的Euler方法.误差分析：仍记注意：则：湘潭大学数学与计算科学学院30于是：整体截断误差的阶由局部截断误差的阶来决定.可见改进的Euler方法误差比Euler方法要高一阶.则有湘潭大学数学与计算科学学院31三、Euler方法的收敛性和稳定性结论:注意:收敛性湘潭大学数学与计算科学学院32稳定性定义湘潭大学数学与计算科学学院33结论:湘潭大学数学与计算科学学院34计算问题：隐式计算格式由迭代法去完成.将上式变形为记湘潭大学数学与计算科学学院35总结通过对Euler方法的讨论可以看到,微分方程数值方法的研究应包括以下方面1.数值计算公式的构造;2.方法稳定性,收敛性的研究;3.方法的误差估计;4.方法的实现等.湘潭大学数学与计算科学学院36§3龙格-库塔法斜率一定取K1K2的平均值吗？步长一定是一个h吗？湘潭大学数学与计算科学学院37Step1:将K2在(xi,yi)点作Taylor展开湘潭大学数学与计算科学学院38Step3:将yi+1与y(xi+1)在xi点的泰勒展开作比较这里有个未知数，个方程。32存在无穷多个解。所有满足上式的格式统称为2阶龙格-库塔格式。Q:为获得更高的精度，应该如何进一步推广？湘潭大学数学与计算科学学院39其中i(i=1,…,m)，i(i=2,…,m)和ij(i=2,…,m;j=1,…,i1)均为待定系数，确定这些系数的步骤与前面相似。最常用为四级4阶经典龙格-库塔法：湘潭大学数学与计算科学学院40由于龙格-库塔法的导出基于泰勒展开，故精度主要受解函数的光滑性影响。对于光滑性不太好的解，最好采用低阶算法而将步长h取小。湘潭大学数学与计算科学学院41§4收敛性与稳定性收敛性解：该问题的精确解为欧拉公式为湘潭大学数学与计算科学学院42稳定性1.00002.00004.00008.00001.60001013.20001011.00002.50001016.25001021.562510