用Rom berg方法求解积分.docVIP

下载本文档

37
0
约5.19千字
约 10页
2020-11-15 发布于江西
举报
版权申诉

用Rom berg方法求解积分.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

用Romberg方法求解积分,要求误差不超过解:Rombｅrg、m文件: function [I, steｐ] =　Rombeｒg(f, a,　ｂ,EＰS) ％ Rｏmberｇ。ｍ就是用龙贝格公式求积分 %　f 为被积函数％ EPS 为积分结果精度％ a,ｂ为积分区间得上下限％　I 为积分结果;ｓtep 为积分得子区间数 m = 1 k ＝ 0 Er = ０、1 Ｈ　＝b-a Ｓ＝ｚeroｓ(1,　1) S(1,　1) ＝ (H/２)　＊　(suｂs(ｓym(f),findsyｍ(syｍ(f)),ａ)＋subs(ｓym(f),findｓym(sym(f)),b)) whilｅＥr　＞ EＰS k = k + 1 　　　 f1 ＝０ H = Ｈ／2 foｒｉ = 1:ｍ　 x = a +Ｈ＊(2*i-１) 　　　　ｆ1 = f1　+ subs(ｓym(f),finｄsym(sym(f)),x) 　　eｎd 　 S(ｋ＋１, 1)　= S(k,　1)/２　+ H*f1 　 m　= 2 *　ｍ　 for n = 1:k 　S(k+１, n+1) = Ｓ(ｋ＋1,　n) + (S(k+1, n)-S(k,　n))/(4＾n—1) enｄ Eｒ = abｓ(S(k+1,　ｎ+1)—Ｓ(ｋ,　n)) end I = Ｓ(ｋ+１, k+1) step = k 命令: ｃleaｒｃlｃ fｏrmat ｓhort a ＝　0; b　= 0、8; EＰS ＝ 1e—２; [I,　ｓtep] =　Ｒｏmbeｒg(＇ｘ^(１／２),　a,　b, ＥＰS) 计算结果: m = 　　1 k = 　　0 Er　= 　0、１000 H　= 　　０。8000 S　= 　 0 Ｓ　= 　0、35７８ｋ＝ 1 f1 = 　 0 Ｈ　＝　　0.4000 ｘ = ０、４0０0 ｆ1 ＝　０.6３25 S = ０。３５78 0、43１9 m = 　 2 S　= 0、3578 　　　 0 　　0.4319 0。4５6６ Er = 　０、0９8８ k　＝　　　　2 f１ = 0 H = 　 0。２０00 x = 　 0.2000 f１ = 0、4472 ｘ＝　　0、60００ｆ1　＝　１.2218 S = ０.3578 　　０　　0、4319 0、4566 　０.4603 0 ｍ = 　　　4 Ｓ = 　　0、3578 　 0 　　 0.４３19 　　 0。4５66 　　　 0。4603 　　0。469８ S = 0、３578　　　　０　　　 0 　　0、43１9　　 0。4５66 　　 0 　　　0。４603 　０.4６98 　　 0。4707 Er = 0.0141 k　＝　 3 ｆ１＝ 0 H　= 　　 0、100０ｘ = 　　 0。10０0 ｆ1 = 　0。3162 x = 0.３0０0 f1　= 　　 0。8640 x = 0、50０0 f1　＝　 1、5711 ｘ = 　　0。７0００ f1 = 　 2、407７Ｓ＝ 0.３578　　 0 　　0 　０。4319　　 0。4566 　０ 0。460３　　０.4698 　０。4７0７　　 0。４70９　　 0 　　　　０ｍ = 8 Ｓ　＝０、３５78 　　　０　　　　０　　　0、4319 　　　0、45６6 　　0 　 0、46０3 0.469８　0。4７07 　０。4709 　０、４74５　　 0 S ＝　 0.３578 　　０　　０　　 0、43１9 0、45６6 　　　 0 　 0、460３　 0。4698 　0.470７ 0、4709 0。4７45 　　０.4748 S = 　０、3578 　 0 　　 0 　　　　 0 ０。４31９　0、45６6 　　　０　　 0 　0、4603