微分方程常用的两种数值解法:欧拉方法与龙格—库塔法.doc
《微分方程常用的两种数值解法:欧拉方法与龙格—库塔法.doc》由会员分享,可在线阅读,更多相关《微分方程常用的两种数值解法:欧拉方法与龙格—库塔法.doc(33页珍藏版)》请在咨信网上搜索。
1、四川师范大学本科毕业论文微分方程常用的两种数值解法:欧拉方法与龙格库塔法学生姓名XXX院系名称数学与软件科学学院专业名称信息与计算科学班 级2006级 4 班学 号20060640XX指导教师Xxx四川师范大学教务处二一年五月微分方程常用的两种数值解法:欧拉方法与龙格库塔法学生姓名:xxx 指导教师:xx【内容摘要】微分方程是最有生命力的数学分支,在自然科学的许多领域中,都会遇到常微分方程的求解问题。当前计算机的发展为常微分方程的应用及理论研究提供了非常有力的工具,利用计算机解微分方程主要使用数值方法,欧拉方法和龙格库塔方法是求解微分方程最典型常用的数值方法。本文详细研究了这两类数值计算方法的
2、构造过程,分析了它们的优缺点,以及它们的收敛性,相容性,及稳定性。讨论了步长的变化对数值方法的影响和系数不同的同阶龙格库塔方法的差别。通过编制C程序在计算机上实现这两类方法及对一些典型算例的结果分析比较,能更深切体会它们的功能,优缺点及适用场合,从而在实际应用中能对不同类型和不同要求的常微分方程会选取适当的求解方法。关键词:显式单步法 欧拉(Euler)方法 龙格库塔(RungeKutta)方法 截断误差 收敛性Two commonly used numerical solution of differential equations:Euler method and Runge - Kutt
3、a methodStudent Name: Xiong Shiying Tutor:Zhang Li【Abstract】The differential equation is the most vitality branch in mathematics. In many domains of natural science, we can meet the ordinary differential equation solution question. Currently, the development of computer has provided the extremely po
4、werful tool for the ordinary differential equation application and the fundamental research, the computer solving differential equation mainly uses value method. The Euler method and the RungeKutta method are the most typical commonly value method to solve the differential equation. This article dis
5、sects the structure process of these two kinds of values commonly value method to solve the analyses their good and bad points, to their astringency, the compatibility, and the stability has made the proof. At the same time, the article discuss the length of stride to the numerical method changing i
6、nfluence and the difference of the coefficient different same step Rungekutta method. Through establishing C program on the computer can realize these two kind of methods, Anglicizing some models of calculate example result can sincerely realize their function, the advantage and disadvantage points
7、and the suitable situation, thus the suitable solution method can be selected to solve the different type and the different request ordinary differential equation in the practical application .Keywords: Explicit single-step process Euler method RungeKutta method truncation error convergence目 录微分方程常用
8、的两种数值解法:欧拉方法与龙格库塔法前言 常微分方程的形成与发展是和力学、天文学、物理学以及其他科学技术的发展密切相关的。数学其他分支的新发展,如复变函数、群、组合拓扑学等,都对常微分方程的发展产生了深刻的影响,当前计算机的发展更是为常微分方程的应用及理论研究提供了非常有力的工具。牛顿在研究天体力学和机械力学的时候,利用了微分方程这个工具,从理论上得到了行星运动规律。后来,法国天文学家勒维烈和英国天文学家亚当斯使用微分方程各自计算出那时尚未发现的海王星的位置。这些都使数学家更加深信微分方程在认识自然、改造自然方面的巨大力量,微分方程也就成了最有生命力的数学分支。然而,我们知道,只有少数十分简单
9、的微分方程能够用初等方法求得它们的解,多数情形只能利用近似方法求解。在常微分方程课中的级数解法,逐步逼近法等就是近似解法。这些方法可以给出解的近似表达式,通常称为近似解析方法。还有一类近似方法称为数值方法,它可以给出解在一些离散点上的近似值,利用计算机解微分方程主要使用数值方法。本文主要讨论一阶常微分方程初值问题 (1.1)在区间上的数值解法,其中为关于,的已知函数,为给定的初始值,将上述问题的精确解记为。该问题常用的数值解法有:欧拉(Euler)方法、龙格库塔(RungeKutta)方法及一些常用的线性多步法。本文重点介绍欧拉(Euler)方法和龙格库塔(RungeKutta)方法。并对这两
10、种方法编制程序,体会它们的功能、优缺点及适用场合,对不同类型常微分方程会选取适当的求解方法。 基本概念和准备知识 一阶常微分方程初值问题是:其中是平面上某个区域上的连续函数,式(1.1.1)的微分方程一般有无穷多个解,式(1.1.2)是确定解的初始条件,如果一元函数对一切满足(1);(2);(3)存在而且;则称是初值问题(1.1)在区间上的解。误差:假定在计算时,用到的前一步的值是准确的即,把用计算得到的近似值记为,估计误差:= y(xn+1),这种误差称为局部截断误差。如果不作这一假定,在每一步计算时除局部截断误差以外,还有由于前一步不准确而引起的误差,称为总体截断误差。收敛性:对于解初值问
11、题的数值方法,我们希望它产生的数值解收敛于初值问题的准确解,“收敛性”的一般定义为:对于所有满足引理1.1条件的初值问题(1.1),如果有一种显式单步法:产生的近似解,对于任意固定的,均有,则称该显式单步法是收敛的。相容性:显式单步法(1.2.1)称为与原微分方程相容,如果 (1.2.3)成立。并称式(1.2.3)为相容性条件。稳定性:在实际计算中,一方面初始值不一定精确,往往带有一定误差,同时由于计算机字长有限,在计算过程中有舍入误差,而且这种误差在式(1.2.1)的递推过程中会传递下去,对以后的结果产生影响。因此要考虑舍入误差的积累是否会得到控制,也即要考虑数值方法的稳定性。当时,若舍入误
12、差引起的后果是有限的,则可以认为该方法是数值稳定的。2 欧拉方法2.1欧拉方法简介对常微分方程初值问题(1.1)用数值方法求解时,我们总是认为(1.1)的解存在且唯一。欧拉方法是解初值问题的最简单、最原始的数值方法,它是显式单步法。下面介绍几种导出欧拉法的途径,每个途径皆可以推导出更为有效的数值方法。(1)Taylor展开在点将作Taylor展开得: 当充分小时略去误差项,并注意到,得,以近似代替,以近似代替,且用“=”代替“”得差分方程初值问题: , (2.1.1)用式(2.1.1)求解初值问题(1.1)的方法称为欧拉方法。(2)数值微分由导数的定义知,对于充分小的整理得,对此作相似的处理也
13、可以得到欧拉方法(2.1.1)。(3)数值积分在区间对积分得 (2.1.2)用数值积分的左矩形公式计算式(2.1.2)右端的积分,得,于是同样可以得到欧拉方法(2.1.1)。(4)多项式插值利用解和其导数在点的值,作一次埃尔米特插值,得到关于的插值多项式:,用近似代替就得到欧拉方法。(5)待定系数法在第步,已知和,利用这两个值估计出下一步的,将已知的值与估计值作线性组合:,其中, 为待定系数。为确定这两个参数,要求这个估计值对和(为常数)精确成立。如果,则,得到方程:,得。如果,则,这样有:,说明,这样估计值为:,即为欧拉方法。欧拉方法的几何意义:由点斜式得切线方程等步长为,则,可由切线方程算
14、出:,逐步计算出在点的值:,用分段的折线逼近函数为“折线法”而非“切线法”,除第一个点是曲线上的切线,其它都不是。0图1欧拉方法的几何意义2.2欧拉方法的截断误差,收敛性,相容性,稳定性设,把在处展开成泰勒级数,即 再由欧拉方法: 两式相减得欧拉方法的局部截断误差为:,若在上充分光滑,且令,则,故欧拉方法是一阶方法或具有一阶精度。欧拉方法的增量函数就是,由引理1.3、引理1.4知当满足Lipschitz条件时欧拉方法是收敛的而且是相容的。用欧拉方法求解典型方程(1.2.4)的计算公式为:,有 。要让,必须有,因此欧拉方法的绝对稳定域为 :,当为实数时,绝对稳定区间为。在复平面上,是以1为半径、
15、以为圆心的内部。3 龙格库塔法3.1 龙格库塔法的基本思想为了导出龙格库塔法的一般公式,我们取如下的线性组合形式: (3.2.3)其中 (3.2.4)即 ,;a21,a31除c1=0外均为待定系数。显然用公式(3.2.3)每计算一个新值要计算函数的值s次,又因每个都能以一种明显的方式由,计算出来,故将公式(3.2.3)称为s级显式龙格库塔法。s级显式龙格库塔法又可以写成下面既简洁又直观的阵列形式: 0 3.2 二、三、四级龙格库塔法 常见的二级二阶龙格库塔法有:中点方法(取)二阶Heun方法(取)取s = 3,完全仿照上述方法推导出三阶龙格库塔公式。这时参数满足下列条件这个方程组解也不唯一,从
16、而可以得到不同的三级三阶龙格库塔法。两个较为常见的三级三阶龙格库塔法是:Kutta方法(取=, = 1,a21 = ,a31 = ,a32 = 2, = , = , = )将它代入(3.2.3)得(3.3.2)三阶Heun方法(取=, = ,a21 = ,a31 = 0,a32 = , = , = , =)将它代入(3.2.3)得通常人们所说的龙格库塔法是指四阶而言的。取s=4,我们同样可以仿照二阶的情形推导出此公式。这样就得到如下含12个未知数但仅有10个方程的方程组 此方程组的解同样不是唯一的,从而有许多不同的四级四阶龙格库塔法,最常用的两个四阶公式是:经典龙格库塔法 (3.3.3)RKG
17、(Runge-Kutta-Gill)方法 (3.3.4)经典龙格库塔方法的几何解释:将(3.3.3)式中的改记为,从而得到经典龙格库塔方法的计算公式为只考虑区间上的解曲线。由(3.3.3)式可知,曲线在处的斜率,在处的斜率,而和则是在中点处的斜率的两个近似值,如图2:图2另一方面有:,用Simpson积分公式逼近此式右端的积分,得上式右端要计算的三个值。由前面的结果,可取,见图3:图3代入后得到 (3.3.5)另外,从(3.3.3)式容易看出,当与无关时,公式(3.3.5)实际上就是用标准的Simpson积分公式计算积分得到的,因此,可以将经典龙格库塔方法视为是用Simpson积分公式的推广形
18、式计算积分而得到的数值方法。经典龙格库塔法的稳定性:增量函数:,其中代入后得于是有=从而得出龙格库塔方法的绝对稳定区间为 1经典龙格库塔方法的算法框图:开 始读入输出=结 束图4 经典龙格库塔法算法框图4 应用举例建模及其分析4.1 欧拉方法解题及其数学模型4.1.1 问题提出 例4.1.1 假定某公司的净资产因资产本身产生了利息而以4%的年利率增长,同时,该公司以每年100万的数额支付职工工资。净资产的微分方程为 (t以年为单位)分别以初始值问题:用Euler公式预测公司24后的净资产趋势。4.1.2 模型建立 分析:这是求微分方程的数值积分,为的是预测公司24年后的净资产趋势。确定变量:设
19、净资产是时间的微分函数,不妨设变量t为时间(以年为单位)。设为第n年的净资产,为第n+1年的净资产,利息以每年4%的速度增长,且公司每年支付职工工资为100万,则第n+1年的净资产增长数额为,由于已得第n年的净资产为,则第n年的净资产加上第n+1年的净资产增长数额就得到第n+1年的净资产。归纳公式: (4-1)确定其微分方程的标准形式,这就是例4.1.1的微分方程模型。4.1.3 解决问题分别以代入,x表示利息赢利低于工资支出的数额,y表示利息赢利与工资支出平衡的数额,z表示利息赢利高于工资支出的数额,计算结果见表1.表1nn11460.25003540.13834.92625004165.0
20、721418.425003581.614768.32425004231.6831375.1425003624.8615699.05625004300.9441330.1425003669.8616627.01925004372.9851283.3525003716.6517552.125004447.961234.6825003765.3218474.18325004525.8271184.0725003815.9319393.15125004606.8581131.4325003868.5720308.87725004691.1291076.6925003923.3121221.232250
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 微分方程 常用 种数 解法 方法 库塔法
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【人****来】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【人****来】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。