微分方程数值解法.ppt
- 文档编号:18745021
- 上传时间:2023-10-26
- 格式:PPT
- 页数:48
- 大小:2.40MB
微分方程数值解法.ppt
《微分方程数值解法.ppt》由会员分享,可在线阅读,更多相关《微分方程数值解法.ppt(48页珍藏版)》请在冰点文库上搜索。
,第6章常微分方程数值解法,6.1引言,6.2欧拉方法,6.3龙格库塔方法,6.1引言,微分方程数值解一般可分为:
常微分方程数值解和偏微分方程数值解。
自然界与工程技术中的许多现象,其数学表达式可归结为常微分方程(组)的定解问题。
一些偏微分方程问题也可以转化为常微分方程问题来(近似)求解。
Newton最早采用数学方法研究二体问题,其中需要求解的运动方程就是常微分方程。
许多著名的数学家,如Bernoulli(家族),Euler、Gauss、Lagrange和Laplace等,都遵循历史传统,研究重要的力学问题的数学模型,在这些问题中,许多是常微分方程的求解。
作为科学史上的一段佳话,海王星的发现就是通过对常微分方程的近似计算得到的。
本章主要介绍常微分方程数值解的若干方法。
1、常微分方程与解,为n阶常微分方程。
如果函数在区间a,b内n阶可导,称方程,为方程满足定解条件的解。
一、初值问题的数值解法,解的图示,本教材重点讨论定解问题(初值问题),定解条件(初始条件),是否能够找到定解问题的解取决于,仅有极少数的方程可以通过“常数变易法”、“可分离变量法”等特殊方法求得初等函数形式的解,绝大部分方程至今无法理论求解。
2、数值解的思想,
(1)将连续变量离散为,
(2)用代数的方法求出解函数在点的近似值,如果找不到解函数数学界还关注:
解的存在性解的唯一性解的光滑性解的振动性解的周期性解的稳定性解的混沌性,求函数y(x)在一系列节点a=x0x1xn=b处的近似值的方法称为微分方程的数值解法。
称节点间距为步长,通常采用等距节点,即取hi=h(常数)。
称为微分方程的数值解。
所谓数值解法:
称在区域D上对满足Lipschitz条件是指:
记,3、相关定义,
(2)一般构造方法:
离散点函数值集合+线性组合结构近似公式,4、迭代格式的构造,
(1)构造思想:
将连续的微分方程及初值条件离散为线性方程组加以求解。
由于离散化的出发点不同,产生出各种不同的数值方法。
基本方法有:
有限差分法(数值微分)、有限体积法(数值积分)、有限元法(函数插值)等等。
(3)如何保证迭代公式的稳定性与收敛性?
5、微分方程的数值解法需要解决的主要问题,
(1)如何将微分方程离散化,并建立求其数值解的迭代公式?
(2)如何估计迭代公式的局部截断误差与整体误差?
二、初值问题解的存在唯一性,考虑一阶常微分方程的初值问题/*Initial-ValueProblem*/:
则上述IVP存在唯一解。
只要在上连续,且关于y满足Lipschitz条件,,即存在与无关的常数L使,对任意定义在上的都成立,,三、初值问题的离散化方法,离散化方法的基本特点是依照某一递推公式,,如果计算需用到前r步的值,则称这类方法为r步方法。
6.2Euler方法,第一步:
连续变量离散化,第二步:
用直线步进,1、Euler格式,18世纪最杰出的数学家之一,13岁时入读巴塞尔大学,15岁大学毕业,16岁获得硕士学位。
1727年-1741年(20岁-34岁)在彼得堡科学院从事研究工作,在分析学、数论、力学方面均有出色成就,并应俄国政府要求,解决了不少地图学、造船业等实际问题。
24岁晋升物理学教授。
1735年(28岁)右眼失明。
1741年-1766(34岁-59岁)任德国科学院物理数学所所长,任职25年。
在行星运动、刚体运动、热力学、弹道学、人口学、微分方程、曲面微分几何等研究领域均有开创性的工作。
1766年应沙皇礼聘重回彼得堡,在1771年(64岁)左眼失明。
Euler是数学史上最多产的数学家,平均以每年800页的速度写出创造性论文。
他去世后,人们用35年整理出他的研究成果74卷。
在假设yi=y(xi),即第i步计算是精确的前提下,考虑的截断误差Ri=y(xi+1)yi+1称为局部截断误差/*localtruncationerror*/。
定义2.2,若某算法的局部截断误差为O(hp+1),则称该算法有p阶精度。
定义2.1,2、欧拉法的局部截断误差,欧拉法的局部截断误差:
Ri的主项/*leadingterm*/,欧拉法具有1阶精度。
例1:
用欧拉公式求解初值问题,取步长。
解:
应用Euler公式于题给初值问题的具体形式为:
其中。
计算结果列于下表:
可用来检验近似解的准确程度。
进行计算,数值解已达到了一定的精度。
这个初值问题的准确解为,,从上表最后一列,我们看到取步长,3、欧拉公式的改进:
隐式欧拉法/*implicitEulermethod*/,向后差商近似导数,由于未知数yi+1同时出现在等式的两边,不能直接得到,故称为隐式/*implicit*/欧拉公式,而前者称为显式/*explicit*/欧拉公式。
一般先用显式计算一个初值,再迭代求解。
隐式欧拉法的局部截断误差:
即隐式欧拉公式具有1阶精度。
梯形公式/*trapezoidformula*/,显、隐式两种算法的平均,注:
梯形公式的局部截断误差,,即梯形公式具有2阶精度,比欧拉方法有了进步。
但注意到该公式是隐式公式,计算时不得不用到,迭代法,其迭代收敛性与欧拉公式相似。
中点欧拉公式/*midpointformula*/,中心差商近似导数,假设,则可以导出即中点公式具有2阶精度。
简单,精度低,稳定性最好,精度低,计算量大,精度提高,计算量大,精度提高,显式,多一个初值,可能影响精度,改进欧拉法/*modifiedEulersmethod*/,Step1:
先用显式欧拉公式作预测,算出,注:
此法亦称为预测-校正法/*predictor-correctormethod*/,可以证明该算法具有2阶精度,同时可以看到它,是个单步递推格式,比隐式公式的迭代求解过程,简单。
后面将看到,它的稳定性高于显式欧拉法。
改进的欧拉法,在实际计算时,可将欧拉法与梯形法则相结合,计算公式为,应用改进欧拉法,如果序列收敛,它的极限便满足方程,改进欧拉法的截断误差,因此,改进欧拉法公式具有2阶精度,解:
对此初值问题采用改进Euler公式,其具体形式为,计算结果列于下表:
例1:
用欧拉公式求解初值问题,改进的Euler法,Euler法,通过计算结果的比较可以看出,改进的Euler方法,的计算精度比Euler方法要高。
欧拉法误差概述,6.3龙格库塔方法,对许多实际问题来说,欧拉公式与改进欧拉公式精度还不能满足要求,为此从另一个角度来分析这两个公式的特点,从而探索一条构造高精度方法的途径.,改进欧拉法,三阶龙格-库塔方法,三阶龙格-库塔方法是用三个值k1,k2,k3的线性组合,要使三阶龙格-库塔方法具有三阶精度,必须使其局部截断误差为O(h4)将k1,k2,k3代入yn+1的表达式中,在(xn,yn)处用二元泰勒公式展开,与y(xn+1)在xn处的泰勒展开式比较,类似二阶龙格-库塔方法的推导过程,8个待定系数c1,c2,c3,a2,a3,b21,b31,b32应满足:
8个未知参数,6个方程,有无穷多组解,三阶龙格库塔公式,四阶Runge-Kutta方法,附注:
二阶Runge-Kutta方法的局部截断误差只能达到,五阶Runge-Kutta方法的局部截断误差只能达到,四阶Runge-Kutta方法的局部截断误差只能达到,三阶Runge-Kutta方法的局部截断误差只能达到,附注:
龙格-库塔法的主要运算在于计算的值,即计算的值。
Butcher于1965年给出了计算量与可达到的最高精度阶数的关系:
由于龙格-库塔法的导出基于泰勒展开,故精,太好的解,最好采用低阶算法而将步长h取小。
度主要受解函数的光滑性影响。
对于光滑性不,感谢您的听讲!
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 微分方程 数值 解法