浅谈常微分方程的数值解法及其应用开题报告文档格式.docx
- 文档编号:22748193
- 上传时间:2023-02-05
- 格式:DOCX
- 页数:10
- 大小:90.07KB
浅谈常微分方程的数值解法及其应用开题报告文档格式.docx
《浅谈常微分方程的数值解法及其应用开题报告文档格式.docx》由会员分享,可在线阅读,更多相关《浅谈常微分方程的数值解法及其应用开题报告文档格式.docx(10页珍藏版)》请在冰豆网上搜索。
(1)
为常微分方程.其中出现的最高阶导数的阶数,叫做常微分方程的阶.例如?
?
?
是一阶常微分方程.是二阶常微分方程.设定义于
区间上,有直到阶的导数,将它代入
(1),使
(1)变成关于的恒等式,即就称=为
(1)的一个定义于上的解,并称为该解的定义区间.[4]
如果一个微分方程中出现多元函数的偏导数,或者说如果未知函数和几个变量有关,而且方程中出现未知函数对几个变量的导数,那么这种微分方程就是偏微分方程.
2.2微分方程产生的历史背景
微分方程差不多是和微积分同时先后产生的,苏格兰数学家耐普尔创立对数的时候,就讨论过微分方程的近似解.牛顿在建立微积分的同时,对简单的微分方程用级数来求解。
后来瑞士数学家雅各布?
贝努利、欧拉、法国数学家克雷洛、达朗贝尔、拉格朗日等人又不断地研究和丰富了微分方程的理论。
微分方程的形成与发展是和力学、天文学、物理学,以及其他科学技术的发展密切相关的.数学的其他分支的新发展,如复变函数、李群、组合拓扑学等,都对微分方程的发展产生了深刻的影响,当前计算机的发展更是为常微分方程的应用及理论研究提供了非常有力的工具.[5]
牛顿研究天体力学和机械力学的时候,利用了微分方程这个工具,从理论上得到了行星运动规律.后来,法国天文学家勒维烈和英国天文学家亚当斯使用微分方程各自计算出那时尚未发现的海王星的位置.这些都使数学家更加深信微分方程在认识自然、改造自然方面的巨大力量.
微分方程的理论逐步完善的时候,利用它就可以精确地表述事物变化所遵循的基本规律,只要列出相应的微分方程,有了解方程的方法.微分方程也就成了最有生命力的数学分支.总之,力学、天文学、几何学等领域的许多问题都导致微分方程.在当代,甚至许多社会科学的问题亦导致微分方程,如人口发展模型、交通流模型等.因而微分方程的研究是与人类社会密切相关的.[6]
2.3微分方程发展现状及其基本功能
在数学学科内部的许多分支中,微分方程是常用的重要工具之一,微分方程进一步发展的需要,有推动着其它数学分支的发展;
相反,微分方程每一步进展都离不开其他数学分支的支援.数学的其他分支的新发展,如复变函数、李群、组合拓扑学等,都对微分方程的发展产生了深刻的影响.当前计算机的发展更是为微分方程的应用及理论研究提供了非常有力的工具.时至今日,可以说微分方程在所有自然科学领域和众多社会科学领域都有着广泛的应用,如自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等.只要能够列出相应的微分方程,有了解方程的方法,利用它就可以精确地表述事物变化所遵循的基本规律.从微积分理论形成以来,人们一直用微分方程来描述、解释或预见各种自然现象,不断的取得了显著的成效.[7]
2.4常微分方程的数值求解方法
2.4.1Euler法
Euler法是最简单的数值方法,为求解良态初值问题,的区间。
实际上,下面的过程不是要找到满足该初值问题的可微函数,而是要生成点集,并且将这些点作为近似解,即。
如何构造“近似满足微方程”的“点集”呢?
首先为这些点选择横坐标,为方便起见,将区间划分为个等距子区间,并选择网络点
k0,1,……,其中
(1)
值称为步长。
然后近似解
在上,
(2)
设,和连续,;
;
利用泰勒定理将在处展开,对每个值,存在一个和之间的值,使得
(3)
将和代人等式(3),得到的表示:
(4)
如果步长h足够小,则可以忽略2次项(包含的项),得到
(5)
这就是欧拉近似。
重复该过程,就能得到近似解曲线的一个点序列。
欧拉方法的一般步骤是
其中k0,1,……,M-1[8](6)
2.4.2泰勒级数法
泰勒级数法有着广泛的应用,并且是比较求解初值问题的各种不同数值方法的标准,它可设计为任意指定的精度。
下面首先将泰勒定理用新的公式表示,使之适合于求解微分方程。
定理9.5(泰勒定理)设,且在不动点处有N次泰勒级数展开:
其中,
(2)
表示函数关t的()次全导数。
求导公式可以递归地计算:
(3)并且一般有
其中为导数算子
区间上的初值问题的近似数值解可由各子区间上的公式
(1)来推导。
次泰勒方法的一般步骤为(5)
其中在各步有。
次泰勒方法的最终全局误差是阶的,因此可选择所需大小的,使得误差足够小。
如果是固定,则理论上可以推导出步长,使之满足任意的最终全局误差。
然而在实际运算中,通常用和计算两个近似结果集,然后比较其结果[9]。
2.4.3龙格?
库塔方法
泰勒方法的优点是最终全局误差的阶为,并且可以通过选择较大的N来得到较小的误差。
然而泰勒方法的缺点是,需要先确定N,并且要计算高阶导数,它们可能十分复杂。
每个龙格一库塔(Runge-Kutta)方法都由一个合适的泰勒方法推导而来,使得其最终全局误差为。
一种折中方法是每步进行若干次函数求值,从而省去高阶导数计算。
这种方法可构造任意N阶精度的近似公式。
最常用的是N4的龙格一库塔方法,它适用于一般的应用,因为它非常精确、稳定,且易于编程。
许多专家声称,没有必要使用更高阶的方法,因为提高的精度与增加的计算量相抵消。
如果需要更高的精度,则应该使用更小的步长或某种自适应方法。
4阶龙格一库塔方法(RK4)可模拟N4的泰勒方法的精度。
它基于如下对,的计算:
其中,,和形如
(2)
通过与N4阶的泰勒级数方法的系数匹配,使得局部误差为,龙格和库塔得出了如下方程组:
该方程组有11个方程和13个未知量,必须补充两个条件才可以求解。
最有用的选择是,(4)
其余变量的解为
将式(4)和(5)中的值代入式
(2)和式
(1),得到标准的阶龙格?
库塔方法,其描述如下。
自初始点开始,利用
(6)
生成近似值序列,其中[10]
(7)
2.4.4预报?
校正方法
欧拉方法、休恩方法、泰勒方法以及龙格一库塔方法都称为单步长方法,因为它们只利用前一个点的信息来计算下一个点,即计算时只使用了初始点。
一般地,只有用来。
当计算出若干个点之后,就可以利用几个已计算出的点来计算下一个点。
以亚当斯一巴什福斯4步法的推导为例,计算需要,,和。
该方法不是自启动的,要生成点,必须先给出其4个初始点,,,(可用前面各节中的方法完成)。
多步法的一个优点是,可以确定它的局部截断误差(localtruncationerror,简称L.T.E.),并可以包含一个校正项,用于在每一步计算中改善解的精确度。
该方法还可以确定步长是否小到能得到的精确值,同时又大到能够免除不必要的和费时的计算。
使用预报子和校正子的组合在每一步只需要进行两次函数求值[11]。
2.4.4.1亚当斯一巴什福斯一莫尔顿方法
亚当斯一巴什福斯一莫尔顿方法(Adams?
Bashforth?
Moulton)是由基本微积分定理推导出的多步法:
1
预报子使用基于点和的的拉格朗日多项式逼近值,并在区间上对式
(1)积分,这个过程产生亚当斯一巴什福斯预报子:
校正子的推导类似。
这时可以实用刚刚计算出的值。
基于点,和新的点构造的一个新的拉格朗日多项式逼近,然后在区间上对该多项式积分,即可得到亚当斯一莫尔顿校正子[12]:
2.4.4.2米尔恩?
辛普森方法
米尔恩?
辛普森方法是预报子基于区间上的对的积分:
预报子使用基于和的拉格朗日多项式逼近,在区间上对它积分,得到米尔恩预报子:
5
校正子的推导类似。
此时值已知,基于点,和新点构造的新的拉格朗日多项式,然后在区间上对该多项式积分,结果为大家所熟悉的辛普森公式[13]:
三、研究的方法与技术路线、研究难点,预期达到的目标
1.研究内容
(1)阐述常微分方程研究的现状,了解微分方程的形成,发展,以及微分方程在描述客观世界中的作用;
(2)掌握一些常见微分方程的数值求解方法;
(3)理解如何用常微分方程解决实际问题,能够做到将理论知识与实际问题相结合,利用常微分方程知识解决实际生活中遇到的几个问题.
2.研究方法及技术路线
通过阅读有关常微分方程方面的论著及文献,了解常微分方程研究的现状.采取了从大量阅读已有的数据资料?
然后对这些内容进行总结?
最后运用相关知识进行分析.
3.研究难点
(1)常微分方程数值求解方法深入研究具有一定的难度;
(2)由于论题比较广泛,很难有独创或新颖之处;
(3)常微分方程应用领域太广,很难研究到多方面.
4.预期达到的目标
了解常微分方程的形成,发展,以及常微分方程在描述客观世界中的作用,掌握一些常微分方程的数值求解方法,理解如何用常微分方程解决实际问题,能够做到将理论知识与实际问题相结合.
四、论文详细工作进度和安排
第一阶段(2010年11月24日?
2010年12月11日):
熟悉毕业论文题目,查阅文献,收集信息、材料并进行加工整理,完成毕业论文文献检索、开题报告、文献综述及外文文献翻译初稿.
第二阶段(2010年12月12日?
2010年1月10日):
完成毕业论文开题报告、文献检索及外文文献翻译,交指导老师.
第三阶段(2011年1月11日?
2011年2月28日):
完成毕业论文的数据收集、论文初稿.
第四阶段(2011年2月29日?
2011年5月3日):
1.进入实习单位进行实习,对论文进行修改;
2.完成毕业实习返校,并递交毕业实习报告,进一步完善毕业论文.
第五阶段(2011年5月16日?
2011年5月23日):
准备毕业论文答辩.
五、主要参考文献:
[1]张良勇,董晓芳.常微分方程的起源与发展[J].高等函授学报自然科学版,20063:
34-38.
[2]黄赛.常微分方程发展的主要历史沿革[J].教育与职业.20062:
36-41.
[3]黄焕福.常微分方程课程建设初探[J].教育与职业,200714:
146-148.
[4]林建平.常微分方程早期发展概观[J].南京工程学院学报,20012:
1-4.
[5]周仲旺.几类特殊的常微分方程[J].潍坊学院学报.20036:
7.
[6]邵晓锋,徐卫卫,李龙星.试析常微分方程模型的归结方法[J].黄冈职业技术学院报,20071:
37-39.
[7]张同斌,管军军,焦万堂.基本微分方程的应用研究[J].河南科学,20076:
888-890.
[8]胡建伟,汤怀民.微分方程数值解方法[M].科学出版社,20034.
[9]JohnH.Mathews,Kurtis//.ericalMethodsUsingMATLABFourthEdition[M].电子工业大学出版社.200512.
[10]王高雄,周之铭等.常微分方程[M].高等教育出版社,20081.
[11]JohnH.Mathews,Kurtis//.ericalMethodsUsingMATLABFourthEdition[M].电子工业大学出版社.200512.
[12]//.mas.偏微分方程的数值解法[M].世界图书出版公司,19958.
[13]郭晓梅.常微分方程的数值解法[J].网络财富,20098.
[14]姜启源,谢金星,叶俊.数学模型[M].高等教育出版社.20038.
[15]JamesStewart.Calculus[M].FifthEdition.Beijing:
高等教育出版社.20047:
1141-1161.
[16]欧阳瑞,孙要伟.常微分方程在数学建模中的应用[J].宿州教育学院学报,20082:
[17]JamesStewart.Calculus[M].FifthEdition.Beijing:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 浅谈 微分方程 数值 解法 及其 应用 开题 报告