第八章基于MATLAB的科学计算—常微分方程数值解法.doc
《第八章基于MATLAB的科学计算—常微分方程数值解法.doc》由会员分享,可在线阅读,更多相关《第八章基于MATLAB的科学计算—常微分方程数值解法.doc(12页珍藏版)》请在咨信网上搜索。
1、科学计算理论、方法及其基于MATLAB的程序实现与分析微分方程(组)数值解法1 常微分方程初值问题的数值解法微分方程(组)是科学研究和工程应用中最常用的数学模型之一。如揭示质点运动规律的Newton第二定律:()和刻画回路电流或电压变化规律的基尔霍夫回路定律等,但是,只有一些简单的和特殊的常微分方程及常微分方程组,可以求得用公式给出的所谓“解析解”或“公式解,如一阶线性微分方程的初值问题:()的解为:()但是,绝大多数在实际中遇到的常微分方程和常微分方程组得不到“解析解”,因此,基于如下的事实:、绝大多数的常微分方程和常微分方程组得不到(有限形式的)解析解;、实际应用中往往只需要知道常微分方程
2、(组)的解在(人们所关心的)某些点处的函数值(可以是满足一定精度要求的近似值);如果只需要常微分方程(组)的解在某些点处的函数值,则没有必要非得通过求得公式解,然后再计算出函数值不可,事实上,我们可以采用下面将介绍的常微分方程(组)的初值问题的数值解法,就可以达到这一目的.一般的一阶常微分方程(组)的初值问题是指如下的一阶常微分方程(组)的定解问题:()其中() ()常微分方程(组)的初值问题通常是对一动态过程(动态系统、动力系统)演化规律的描述,求解常微分方程(组)的初值问题就是要了解和掌握动态过程演化规律。1 常微分方程(组)的Cauch问题数值解法概论假设要求在点(时刻),处初值问题()
3、的解的(近似)值,如果已求得时刻的值或它的近似值(如时刻的值),那么将式()的两端在区间上积分(10)可得 (11)或(12)显然,为了利用式(11)或(12)求得的精确值(近似值),必须计算右端的积分,这是问题的关键也是难点所在,如前所述,一般得不到精确的公式解,因此需要采用数值积分的方法求其近似解,可以说,不同的式值积分方法将给出不同的Cauch问题的数值解法。1。2最简单的数值解法Euler 方法假设要求在点(时刻),处初值问题()的解的近似值。首先对式()的两端积分,得(13)对于式(13)的右边,如果用积分下限处的函数值代替被积函数作积分(从几何上的角度看,是用矩形面积代替曲边梯形面
4、积),则有(14)进而得到下式给出的递推算法Euler 方法(15)例用Euler 方法解如下初值问题,取,解:由(15)得结果如下:open Euler_Method.m 如果取,其结果如下图所示:Euler_Method 1。3改进的Euler 方法 对于(15)的右边,如果被积函数用积分限和处的函数值的算术平均值代替(几何上,是用梯形面积代替曲边梯形面积),则有(16)进而得到下式给出的递推算法:(17)通常算法(17)比Euler 方法(15)的精度高,但是,按算法(17)求时要解(非线性)方程(组),这是算法(17)不如Euler 方法的方面,为了) 尽可能地保持算法(17)精度高的
5、优点;) 尽可能地利用Euler 方法计算简单的长处;人们采取了如下的称之为改进的Euler 方法的折衷方案:预测 (18)修正 (19) 例Euler 方法与改进的Euler 方法的比较下图是当时比较的结果:open Improved_Euler_Method。m 1.4Euler 方法和改进的Euler 方法的误差分析由Taylor 公式 (19)说明Euler 方法的截断误差是,类似地,由 (20) (21)以及 (22)让式(20)的两端减式(21)的两端,可得 (23)从上述推导Euler 方法、改进的Euler 方法的过程以及例、例容易看出,改进的Euler 方法Euler 方法的
6、精度高,其原因在于:1 在推导Euler 方法时,我们是用待求解函数在一点处的变化率代替在区间上的平均变化率: (24)2 而在推导改进的Euler 方法时,我们是用待求解函数在两点处变化率的平均值代替在区间上的平均变化率;显然,通常比更接近于在区间上的平均变化率。由此启发人们:适当地选取区间上函数若干点处的变化率,用它们加权平均值代替在区间上的平均变化率,近似解的精度应更高。下面将要介绍的Runge-Kutta法就是基于上述想法得到的。2RungeKutta法RungeKutta法是按选取区间上函数变化率的个数的多少和截断误差的阶数来区分的一系列方法,如 二阶的RungeKutta法(改进的
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第八 基于 MATLAB 科学 计算 微分方程 数值 解法
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【天****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【天****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。