原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
NA-5-2复化求积公式[精]
第 二 节 复化求积公式 Example 1 Approximate the integral Using the Composite Trapezoidal rule and Composite Simpson’s rule x=eps:0.01:pi; y=sin(x)./x; plot(x,y); legend(f(x)=sin(x)/x); function rs= trapezoid(f,a,b,n) h = (b-a)/n; r= (feval(f,a)+feval(f,b))/2 ; for j = 1:n-1 x=a+j*h ; r= r+ feval(f,x); end rs = r* h; n=8, tp=1.84784230644461 n=16, tp=1.85091414036536 n=32, tp=1.85168137241373 n=64, tp=1.85187313510989 n=128, tp=1.85192107295303 n=256, tp=1.85193305723690 n=512, tp=1.85193605329681 n=1024, tp=1.85193680231110 n=2048, tp=1.85193698956462 n=4096, tp=1.85193703637800 n=8192, tp=1.85193704808136 function rs= simpson(s,a,b,n) h = (b-a)/n; r= feval(s,a)+feval(s,b); for j = 1:2:n-1 x=a+j*h ; r= r+ 4*feval(s,x); end for j = 2:2:n-2 x=a+j*h ; r= r+ 2*feval(s,x); end rs = r*h/3; format long f =inline( sin(x)/x ) ; a = eps; b =pi; n = 8; % eps 是Matlab最小正数 sp = simpson(f, a, b, n) n=8, sp=1.85195370083255 n=16, sp=1.85193808500560 n=32, sp=1.85193711642985 n=64, sp=1.85193705600861 n=128, sp=1.85193705223407 n=256, sp=1.85193705199819 n=512, sp=1.85193705198345 n=1024, sp=1.85193705198253 n=2048, sp=1.85193705198247 n=4096, sp=1.85193705198246 n=8192, sp=1.85193705198247 数值试验结果 用复化Simpson公式求解，当区间分成64份时，计算结果： sp=1.85193705600861 比用复化梯形公式，当区间分成8192份 时的计算结果： tp=1.85193704808136 还要精确。 启发：寻求精度更高的求解方法！ 复化梯形、复化 Simpson 公式计算结果对比 n=8, tp=1.84784230644461 n=16, tp=1.85091414036536 n=32, tp=1.851681372
文档评论（0）