版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
欧拉方法在常微分方程第1页,共15页,2023年,2月20日,星期五第8章常微分方法的数值解法
科学技术与工程问题常常需要建立微分方程形式的数学模型,下面是这类问题的例子。设N(t)为某物种的数量,为该物种的的出生率与死亡率之差,为生物的食物供给及它们所占空间的限制,描述该物种增长率的数学模型是设Q是电容器上的带电量,C为电容,R为电阻,E为电源的电动势,描述该电容器充电过程的数学模型是第2页,共15页,2023年,2月20日,星期五以上两个例子是常微分方程初值问题,下面是一个两点边值问题的例子。设一跟长为L的矩形截面的梁,两端固定。E是弹性模量,S是端点作用力,I(x)是惯性矩,q是均匀荷载强度,梁的桡度y(x)满足如下方程针对实际问题建立的数学模型,要找出模型解的解析表达式往往是困难的,甚至是不可能的。因此,需要研究和掌握微分方程的数值解法,即计算解域内离散点上的近似值的方法。本章讨论常微分方程数值解的基本方法和理论。第3页,共15页,2023年,2月20日,星期五8.1Euler方法8.1.1Euler方法及其有关的方法考虑一阶常微分方程初值的问题:设f(x,y)是连续函数,对y满足Lipschitz条件,这样初值问题的解是存在唯一的,而且连续依赖于初始条件。为了求得离散点上的函数值,将微分方程的连续问题(8.1.1)进行离散化。一般是引入点列{},这里为步长,经常考虑定长的情形,即。记为初始问题(8.1.1)的问题准确解在处的值,用均差近似代替(8.1.1)的导数得第4页,共15页,2023年,2月20日,星期五令为的近似值,将上面两个近似写成等式,整理后得(8.1.2)(8.1.3)从处的初值开始,按(8.1.2)可逐步计算以后各点上的值。称(8.1.2)式为显式Euler。由于(8.1.3)式的右端隐含有待求函数值,不能逐步显式计算,称(8.1.3)式为隐式Euler公式或后退Euler公式。如果将(8.1.2)和(8.1.3)两式作算术平均,就得梯形公式。第5页,共15页,2023年,2月20日,星期五梯形公式也是隐式公式。以上公式都是由去计算,故称它们为单步法。例8.1取h=0.1,用Euler方法、隐式Euler方法和梯形方法解解本题有如果用Euler方法,由(8.1.2)并代入h=0.1得同理,用隐式Euler方法有(8.1.4)第6页,共15页,2023年,2月20日,星期五用梯形公式有三种方法及准确解的数值结果如表8-1所示。从表中看到,在处,Euler方法和隐式Euler方法的误差分别是和,而梯形方法的误差却是。
在例8.1中,由于f(x,y)对y是线性的,所以对隐式公式也可以方便地计算。但是,当f(x,y)是y的非线性函数时,如,其隐式Euler公式为。显然,它是的非线性方程,可以选择非线性方程求根的迭代求解。以梯形公式为例,可用显式Euler公式提供迭代初值,用公式第7页,共15页,2023年,2月20日,星期五表8-1Euler方法隐式Euler方法梯形法准确解011110.11.0000001.0090911.0047621.0048370.21.0100001.0264461.0185491.0187310.31.0290001.0513151.0406331.0408180.41.0561001.0830131.0700961.0703200.51.0904901.1209211.1062781.106531第8页,共15页,2023年,2月20日,星期五反复迭式,直到其中,步长h成为迭代参数,它需要满足一定的条件,才能收敛。若将(8.1.4)式减去该迭代公式,得假设f(x,y)关于y满足Lipschiz条件,则有第9页,共15页,2023年,2月20日,星期五这里,L是Lipschiz常数。当hL/2<1即h<2/L时,迭代序列收敛。对于隐式公式,通常采用估计-校正技术,即先用显式公式计算,得到预估值,然后以预估值作为隐式公式的迭代初值,用隐式公式迭代一次得到校正值,称为预估-校正技术。例如,用显式Euler公式作预估,用梯形公式作校正,即称该公式为改进的Euler公式。它显然等价于显式公式为,(8.1.6)第10页,共15页,2023年,2月20日,星期五也可以表示为下列平均化的形式例8.2取h=0.1,用改进的Euler方法解解按(8.1.5),改进的Euler方法解第11页,共15页,2023年,2月20日,星期五由得计算结果如表8-2。该初值问题的准确解为。表8-20.10.20.30.40.50.60.70.81.09591.18411.26621.34341.41641.48601.55251.61531.09541.18321.26491.34161.41421.48321.54921.6165第12页,共15页,2023年,2月20日,星期五8.1.2局部误差和方法的阶
初值问题(8.1.1)的单步法可以写成如下统一形式(8.1.7)其中与有关。若中不含,则方法是显式的,否则是隐式的,所以一般显式单步法表示为(8.1.8)例如,Euler方法中,有对于不同的方法,计算值与准确解的误差各不相同。所以有必要讨论方法的截断误差。我们称为某一方法在点的整体截断误差。显然,不单与这步的计算有关,它与以前各步的计算也有关,所以误差被称为整体的。分析和估计整体截断误差是复杂的。为此,我们假设处的没有误差,即,考虑从到这一步的误差,这就是如下的局部误差的概念。第13页,共15页,2023年,2月20日,星期五定义8.1设是初值问题(8.1.1)的准确解,则称为单步法(8.1.7)的局部截断误差。
定义8.2如果给定方法的局部截断误差,其中为整数,则称该方法是p阶的,或具有p阶精度。若一个p阶单步法的局部截断误差为则称其第一个非零项为该方法的局部截断误差的主项。对于Euler方法,有Taylor展开有第14页,共15页,2023年,2月20日,
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 《信息产业》课件
- 证券结构化产品协议三篇
- 《球墨铸铁直埋热水管道技术规程》公示稿
- 校园美术作品长廊建设规划计划
- 典当服务相关行业投资规划报告范本
- 工具台车相关项目投资计划书
- 情感教育与道德认知的结合计划
- 增强幼儿园团队建设的策略计划
- 青少年犯罪预防的保安策略计划
- 理财规划师课件(综合案例分析)
- 期末综合复习(试题)-2023-2024学年五年级上册数学人教版
- 充电桩竞争格局分析
- 线性规划课后题答案
- 压裂作业质量、安全、环保管理体系与措施
- JGJ120-2012建筑基坑支护技术规程
- 国防教育基地建设项目运营方案
- 维生素B市场研究报告
- 四川省广安市2023年九年级上学期期末化学试题附答案
- 2022年云南省昆明市中考语文真题答案详解
- 五年级信息技术上册期末试卷及答案
- 工厂空气重度污染应急预案
评论
0/150
提交评论