第六章节常微分方程解法幻灯片.pptVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
舍去En,用yn-i代替y(xn-i),(i=-1,0,1,…,r, 得到 其中 若令 §6.8 线性多步法 则上式写为 这就是Adams-Bashforth显式公式。 §6.8 线性多步法 用泰勒展开构造 一般的线性多步公式具有形式: 其中 注意到k=-1时 中含有未知yn+1所以 时,方法是隐式, 时,是显式。 §6.8 线性多步法 记 泰勒展开 §6.8 线性多步法 代入 得 §6.8 线性多步法 欲使上式成为P阶,要求局部截断误差为O(hP+1),只要令上式与y(xn+1)泰勒展开 在hP项之前相同,为此,要求成立 §6.8 线性多步法 这是一个2r+3个未知数的p+1个线性方程组。可以证明,当p=2r+2时,存在唯一解。 §6.8 线性多步法 稳定性: 定义:绝对稳定性:对于给定的常微分方程,当用某一数值方法使用给定的步长求解时,如果由一个值yn上大小为的扰动所引起的其后值ym(mn)上的变化都不大于 ,则称方法是绝对稳定的。 对于一般方程,这种稳定性难以分析。 §6.8 线性多步法 对于欧拉法: 因此,使得yn上的扰动对后面的计算不一步一步增加的要求是 ,所以欧拉方法的绝对稳定区间是(-2,0),即 §6.8 线性多步法 §6.6 龙格-库塔方法 于是 其中P可以任意选取。这个式子叫P阶泰勒方法。 §6.6 龙格-库塔方法 这样,便得到了一种P阶显式单步法 P=1时, ,泰勒方法变为欧拉方法。 §6.6 龙格-库塔方法 使用泰勒方法时,需要计算f(2)(x,y), f(3)(x,y),…,对于复杂的f,计算较繁,所以高阶泰勒方法并不实用。 在欧拉法中, ,即 取为xn处的折线斜率。 Remarks: §6.6 龙格-库塔方法 在改进欧拉方法中, 为平均值,其中k1=f’(xn,yn), k2=f’(xn+h,y1+hk1),于是可以想象,如果取为更多个斜率k的加权平均值,有可能使阶数提高。 其中 §6.6 龙格-库塔方法 Runge-Kutta方法一般用R个斜率,称为R级。定义为: 也可以用另外一种方法构造R-K法。由微分中值定理知: §6.6 龙格-库塔方法 其中 h=xn+1-xn §6.6 龙格-库塔方法 可以看作y(x)在区间[xn,xn+1]上的平均斜率 ,所以只要对平均斜率 给出一种近似算法,即可得到一种计算y(xn+1)的近似值计算公式。 例如,取左端点斜率作为平均斜率,即 可以得到: §6.6 龙格-库塔方法 这即是一阶精度的欧拉公式。 又如:用左端点处斜率和右端点处的斜率的算术平均值作为平均斜率,可得: 这就是改进欧拉公式 §6.6 龙格-库塔方法 可以想象,如果能够在[xn,xn+1]上多估计几个点(包括端点)处的斜率值k1, k2, …, km,然后用它们的加权平均值作为平均斜率 ,那么可以构造出这样的公式: §6.6 龙格-库塔方法 其中 , 是 预估值。 这种公式叫龙格-库塔(R-K)公式。式中 都是待定参数,可通过使公式的局部截断误差关于h的阶数尽可能高来确定。 现在由 §6.6 龙格-库塔方法 来推导二、三、四阶龙格-库塔方法。 回顾: 二元函数的泰勒展开 §6.6 龙格-库塔方法 其中 二阶R-K方法 §6.6 龙格-库塔方法 为了使 是二阶的 考虑 P阶泰勒方法的 为 把k2按h作泰勒展开 §6.6 龙格-库塔方法 于是 §6.6 龙格-库塔方法 比较上式与2阶泰勒方法,使两式的前两项符合,得到c1,c2,a2,b21满足的关系式: 这是三个方程,未知数有四个,解不唯一,于是可以得到不同的二级二阶R-K方法。 §6.6 龙格-库塔方法 (1)例如,取c2=c1=1/2, a2=b21=1, 得到 这便是改进欧拉方法。 (2)又例如,取c2=1,c1=0,a2=b21=1/2, 得到 §6.6 龙格-库塔方法 即 这种方法叫中点法。

文档评论(0)

精品课件 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档