用Romberg方法求解积分.doc

下载文档 降价啦

29
0
约5.44千字
约 14页
2016-08-18 发布于湖北
举报
版权申诉
保障服务

用Romberg方法求解积分.doc

1、本文档共14页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

用Romberg方法求解积分

用Romberg方法求解积分，要求误差不超过解：Romberg.m文件： function [I, step] = Romberg(f, a, b,EPS) % Romberg.m 是用龙贝格公式求积分 % f 为被积函数 % EPS 为积分结果精度 % a,b 为积分区间的上下限 % I 为积分结果；step 为积分的子区间数 m = 1 k = 0 Er = 0.1 H =b-a S = zeros(1, 1) S(1, 1) = (H/2) * (subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b)) while Er EPS k = k + 1 f1 = 0 H = H/2 for i = 1:m x = a +H*(2*i-1) f1 = f1 + subs(sym(f),findsym(sym(f)),x) end S(k+1, 1) = S(k, 1)/2 + H*f1 m = 2 * m for n = 1:k S(k+1, n+1) = S(k+1, n) + (S(k+1, n)-S(k, n))/(4^n-1) end Er = abs(S(k+1, n+1)-S(k, n)) end I = S(k+1, k+1) step = k 命令： clear clc format short a = 0; b = 0.8; EPS = 1e-2; [I, step] = Romberg(x^(1/2), a, b, EPS) 计算结果： m = 1 k = 0 Er = 0.1000 H = 0.8000 S = 0 S = 0.3578 k = 1 f1 = 0 H = 0.4000 x = 0.4000 f1 = 0.6325 S = 0.3578 0.4319 m = 2 S = 0.3578 0 0.4319 0.4566 Er = 0.0988 k = 2 f1 = 0 H = 0.2000 x = 0.2000 f1 = 0.4472 x = 0.6000 f1 = 1.2218 S = 0.3578 0 0.4319 0.4566 0.4603 0 m = 4 S = 0.3578 0 0.4319 0.4566 0.4603 0.4698 S = 0.3578 0 0 0.4319 0.4566 0 0.4603 0.4698 0.4707 Er = 0.0141 k = 3 f1 = 0 H = 0.1000 x = 0.1000 f1 = 0.3162 x = 0.3000 f1 = 0.8640 x = 0.5000 f1 = 1.5711 x = 0.7000 f1 = 2.4077 S = 0.3578 0 0 0.4319 0.4566 0 0.4603 0.4698 0.4707 0.4709 0 0 m = 8 S = 0.3578 0 0 0.4319 0.4566 0 0.4603 0.4698 0.4707 0.4709 0.4745 0 S = 0.3578