- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数值分析方法第7章
? 考虑一阶常微分方程的初值问题 /* Initial-Value Problem */:;例如:;这是微积分的发明者之一Leibniz在1686年曾经让当时
数学界人士求解的一阶微分方程式,吸引了许多数学家
的注意,大约经过150年的探索
到1838年, 刘维尔(Liouville)在理论上证明了这个微分方程不能用初等积分法求解,得借助于数值方法
;数值方法的基本思想是:在解的存在区间上取n + 1个节点
;步进式:
根据已知的或已求出的节点上的函数值计算当前节点上的函数值,一步一步向前推进。
因此只需建立由已知的或已求出的节点上的函数值求当前节点函数值的递推公式即可。;§1 欧拉方法 /* Euler’s Method */;也称欧拉折线法.;例 用欧拉法求初值问题
当h = 0.02时在区间[0, 0.1]上的数值解。;具体计算结果;表中y(xn),是初值问题的真解
在xn上的值。
为近似值yn的误差。
从表中可以看出,随着n的增大,误差也在增大,所以说,欧拉法计算简便,对一些问题有较大的使用价值,但是,它的误差较大,所得的数值解精确度不高。;定义;Euler’s Method;§ Euler’s Method;§ Euler’s Method;例 ;?
xn
;? 中点欧拉公式 /* midpoint formula */;方 法;? 改进欧拉法 /* modified Euler’s method */;§2 龙格 - 库塔法 /* Runge-Kutta Method */;§2 Runge-Kutta Method;§2 Runge-Kutta Method;其中?i ( i = 1, …, m ),?i ( i = 2, …, m ) 和 ?ij ( i = 2, …, m; j = 1, …, i?1 ) 均为待定系数,确定这些系数的步骤与前面相似。
文档评论(0)