《开题报告-浅谈常微分方程的数值解法及其应用》由会员分享,可在线阅读,更多相关《开题报告-浅谈常微分方程的数值解法及其应用(10页珍藏版)》请在金锄头文库上搜索。
1、毕业论文开题报告信息与计算科学浅谈常微分方程的数值解法及其应用一、选题的背景、意义1、选题的背景微分方程差不多是和微积分同时先后产生的,苏格兰数学家耐普尔创立对数的时候,就讨论过微分方程的近似解. 牛顿在建立微积分的同时,对简单的微分方程用级数来求解. 后来瑞士数学家雅各布 ?贝努利、欧拉、法国数学家克雷洛、达朗贝尔、拉格朗日等人又不断地研究和丰富了微分方程的理论. 微分方程的形成及发展与力学、天文学、 物理学、生物学,以及其他科学技术的发展密切相关 .在数学学科内部的许多分支中,微分方程是常用的重要工具之一,微分方程进一步发展的需要, 有推动着其它数学分支的发展;相反, 常微分方程每一步进展
2、都离不开其他数学分支的支援 . 数学的其他分支的新发展,如复变函数、李群、组合拓扑学等,都对微分方程的发展产生了深刻的影响.当前计算机的发展更是为微分方程的应用及理论研究提供了非常有力的工具 .时至今日,可以说微分方程在所有自然科学领域和众多社会科学领域都有着广泛的应用,如自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等.只要能够列出相应的微分方程,有了解方程的方法,利用它就可以精确地表述事物变化所遵循的基本规律.从微积分理论形成以来,人们一直用微分方程来描述、解释或预见各种自然现象,不断的取得了显著的成效. 2、选题的意义微分方程的理论逐步
3、完善的时候,利用它就可以精确地表述事物变化所遵循的基本规律,只要列出相应的微分方程,有了解方程的方法. 微分方程也就成了最有生命力的数学分支. 总之, 力学、 天文学、 几何学等领域的许多问题都导致微分方程. 在当代, 甚至许多社会科学的问题亦导致微分方程,如人口发展模型、交通流模型等. 因而微分方程的研究是与人类社会密切相关的. 1“常微分方程” 是理学院数学系所有专业学生的重要专业基础课之一,也是工科、 经济等专业必学内容之一. 其重要性在于它是各种精确自然科学、社会科学中表述基本定律和各种问题的根本工具之一,换句话说, 只要根据实际背景,列出了相应的微分方程,并且能(数值地或定性地)求出
4、这种方程的解,人们就可以预见到,在已知条件下这种或那种“运动”过程将怎样进行,或者为了实现人们所希望的某种“运动” 应该怎样设计必要的装置和条件等等 . 例如,我们要设计人造卫星轨道,首先,根据力学原理,建立卫星运动的微分方程,列出初始条件, 然后求出解, 即卫星运行轨道. 随着物理科学所研究的现象在广度和深度两方面的扩展,偏微分方程的应用范围更广泛. 2从数学自身的角度看,偏微分方程的求解促使数学在函数论、变分法、级数展开、常微分方程、代数、微分几何等各方面进行发展. 从这个角度说,偏微分方程变成了数学的中心. 3总之, 微分方程从它诞生起即日益成为人类认识并进而改造自然、社会的有力工具,成
5、为数学科学联系实际的主要途径之一. 文章就常微分进行展开,对其数值解法进行简单的阐述.二、研究的基本内容与拟解决的主要问题2.1 微分方程概念介绍2.1.1 微分方程概况由一元函数得到的方程. 即:称含有自变量,未知函数及其导数的关系式22( , ,.,)0nndy d yd yF x ydxdxdx. (1)为常微分方程 . 其中出现的最高阶导数的阶数,叫做常微分方程的阶. 例如dydxx,dyydx,是一阶常微分方程. 22sin0dgdtp是二阶常微分方程. 设)(xy定义于区间J上,有直到n阶的导数,将它代入(1) ,使( 1)变成关于x的恒等式,即( )( )( , ( ),.,)0
6、,nndxdxF xxxJdxdx. 就称y( )x为( 1)的一个定义于J上的解,并称J为该解的定义区间. 4 如果一个微分方程中出现多元函数的偏导数,或者说如果未知函数和几个变量有关,而且方程中出现未知函数对几个变量的导数,那么这种微分方程就是偏微分方程.2.2 微分方程产生的历史背景微分方程差不多是和微积分同时先后产生的,苏格兰数学家耐普尔创立对数的时候,就讨论过微分方程的近似解. 牛顿在建立微积分的同时,对简单的微分方程用级数来求解。后来瑞士数学家雅各布 ?贝努利、欧拉、法国数学家克雷洛、达朗贝尔、拉格朗日等人又不断地研究和丰富了微分方程的理论。微分方程的形成与发展是和力学、天文学、
7、物理学, 以及其他科学技术的发展密切相关的. 数学的其他分支的新发展,如复变函数、李群、组合拓扑学等,都对微分方程的发展产生了深刻的影响, 当前计算机的发展更是为常微分方程的应用及理论研究提供了非常有力的工具 . 5牛顿研究天体力学和机械力学的时候,利用了微分方程这个工具,从理论上得到了行星运动规律 .后来,法国天文学家勒维烈和英国天文学家亚当斯使用微分方程各自计算出那时尚未发现的海王星的位置. 这些都使数学家更加深信微分方程在认识自然、改造自然方面的巨大力量 . 微分方程的理论逐步完善的时候,利用它就可以精确地表述事物变化所遵循的基本规律,只要列出相应的微分方程,有了解方程的方法. 微分方程
8、也就成了最有生命力的数学分支. 总之, 力学、 天文学、 几何学等领域的许多问题都导致微分方程. 在当代, 甚至许多社会科学的问题亦导致微分方程,如人口发展模型、交通流模型等. 因而微分方程的研究是与人类社会密切相关的. 6 2.3 微分方程发展现状及其基本功能在数学学科内部的许多分支中,微分方程是常用的重要工具之一,微分方程进一步发展的需要, 有推动着其它数学分支的发展;相反,微分方程每一步进展都离不开其他数学分支的支援 . 数学的其他分支的新发展,如复变函数、李群、组合拓扑学等,都对微分方程的发展产生了深刻的影响. 当前计算机的发展更是为微分方程的应用及理论研究提供了非常有力的工具 . 时
9、至今日,可以说微分方程在所有自然科学领域和众多社会科学领域都有着广泛的应用,如自动控制、 各种电子学装置的设计、弹道的计算、 飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等. 只要能够列出相应的微分方程,有了解方程的方法,利用它就可以精确地表述事物变化所遵循的基本规律. 从微积分理论形成以来,人们一直用微分方程来描述、解释或预见各种自然现象,不断的取得了显著的成效. 72.4 常微分方程的数值求解方法2.4.1 Euler 法Euler法是最简单的数值方法, , a b为求解良态初值问题( , )yf t y,0( )y ay的区间。实际上,下面的过程不是要找到满足该初值问题的可微函
10、数,而是要生成点集(,)kkty,并且将这些点作为近似解,即( )kky ty。如何构造 “近似满足微方程”的“点集”呢?首先为这些点选择横坐标,为方便起见,将区间 , a b划分为M个等距子区间,并选择网络点ktakh,k=0,1, ,M其中bahM(1)值h称为步长。然后近似解( ,)yf t y在0 ,Mtt上,00()y ty(2)设( )y t,( )y t和( )y t连续,; ;利用泰勒定理将( )y t在0tt处展开,对每个值t,存在一个0t和t之间的值1c,使得210000()()( )()()()2y ctty ty ty ttt( 3)将00( )( ,( )y tf t
11、y t和10htt代人等式( 3) ,得到1( )y t的表示:211000()( )( )( ,( )2y c hy ty thf ty t(4)如果步长h 足够小,则可以忽略2 次项(包含2h的项),得到1000( ,)yyhf ty(5)这就是欧拉近似。重复该过程,就能得到近似解曲线( )yy t的一个点序列。欧拉方法的一般步骤是1kktth,1(,)kkkkyyhf ty其中k = 0,1, ,M-18( 6)2.4.2 泰勒级数法泰勒级数法有着广泛的应用,并且是比较求解初值问题的各种不同数值方法的标准,它可设计为任意指定的精度。下面首先将泰勒定理用新的公式表示,使之适合于求解微分方程
12、。定理 9.5 (泰勒定理)设1( )Ny tC0, tb,且( )y t在不动点0, ktttb处有 N 次泰勒级数展开:1()()(, ()()NkkNkky thy thTty tO h(1)其中,()11( )(, ( )!jNjkNkkjytTty thj(2)( )1( )( , ( )jjytft y t表示函数f关 t 的(1j)次全导数。求导公式可以递归地计算:(3)22(4)23232( )( )( )2()2()33()33()(33)(2)tytyttytyyyttytyyytytttyttyyttyyyyyyytttyttyytyyyyttytyyy tfy tff
13、yff fytff yf yfyffffffff fyffyfyf yf yfy yfyfffffffffffff23()()()tyytyyytyfffffffff f(3)并且一般有()(1)( )( , ( )NNytPf t y t(4)其中P为导数算子()Pfty区间0 ,Mtt上的初值问题( )( , )y tf t y的近似数值解可由各子区间1,kktt上的公式(1)来推导。N次泰勒方法的一般步骤为323211.2!3!NNkkd hd hd hyyd hN(5)其中在各步0,1,.,1kM有()(),1,2,.,jjkdytjN。N次泰勒方法的最终全局误差是1()NO h阶的,
14、因此可选择所需大小的N,使得误差足够小。如果N是固定,则理论上可以推导出步长h,使之满足任意的最终全局误差。然而在实际运算中,通常用h和/ 2h计算两个近似结果集,然后比较其结果9。2.4.3 龙格库塔方法泰勒方法的优点是最终全局误差的阶为()NO h, 并且可以通过选择较大的 N 来得到较小的误差。然而泰勒方法的缺点是,需要先确定 N ,并且要计算高阶导数,它们可能十分复杂。每个龙格一库塔(Runge-Kutta )方法都由一个合适的泰勒方法推导而来,使得其最终全局误差为()NO h。一种折中方法是每步进行若干次函数求值,从而省去高阶导数计算。这种方法可构造任意 N 阶精度的近似公式。最常用
15、的是N= 4 的龙格一库塔方法,它适用于一般的应用,因为它非常精确、稳定,且易于编程。许多专家声称,没有必要使用更高阶的方法,因为提高的精度与增加的计算量相抵消。如果需要更高的精度, 则应该使用更小的步长或某种自适应方法。 4 阶龙格一库塔方法( RK4 )可模拟N=4 的泰勒方法的精度。它基于如下对1ky,的计算:111223344kkyyw kw kw kw k(1)其中1k,2k,3k和4k形如1211132213243415263( ,)(,)(,)(,)kkkkkkkkkhf tykhf ta h yb kkhf ta h yb kb kkhf ta h yb kb kb k(2)通
16、过与 N = 4 阶的泰勒级数方法的系数匹配,使得局部误差为5()O h,龙格和库塔得出了如下方程组:112324563123421324322221324311213babbabbbawwwww aw aw aw aw aw a(3)33321324331 3415263123431 526222313415264136141()61()81()12124w aw aw aw a bwa ba bw a a bw aa ba bw a bwa ba bw a b b该方程组有11 个方程和13 个未知量, 必须补充两个条件才可以求解。最有用的选择是112a,20b(4)其余变量的解为23134561234111,1,0,0,12221111,6336aabbbbbwwww(5)将式( 4)和( 5)中的值代入式(2)和式( 1) ,得到标准的4N阶龙格库塔方法,其描述如下。自初始点00( ,)ty开始,利用12341(22)6kkh ffffyy(6)生成近似值序列,其中10 1213243( ,)(,)22(,)22(,)kkkkkkkkff tyhhff tyfhhff tyff