用Romberg方法求解积分.doc
用Romberg方法求解积分.doc
- 配套讲稿：
如PPT文件的首页显示word图标，表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制：
部分文档作品中含有的国旗、国徽等图片，仅作为作品整体效果示例展示，禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词：
- Romberg 方法 求解 积分
- 资源描述：
-
1. 用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 0 0 0.4319 0.4566 0 0.4603 0.4698 0.4707 0.4709 0.4745 0.4748 S = 0.3578 0 0 0 0.4319 0.4566 0 0 0.4603 0.4698 0.4707 0 0.4709 0.4745 0.4748 0.4748 Er = 0.0042 I = 0.4748 step = 3 I = 0.4748 step = 3 2. 设方程组试用Jacobi迭代法求解此方程，，当时终止迭代。 解：Jacobi.m文件： function Jacobi(A, b, max, eps) %max为最大迭代次数，eps为容许误差 n = length(A); x = zeros(n, 1); x1 = zeros(n, 1); k = 0; while 1 x1(1) = ( b(1) - A(1,2:n) * x(2:n,1) )/A(1,1) for i = 2:n-1 x1(i) = ( b(i) - A(i,1:i-1) * x(1:i-1,1) - A(i,i+1:n) * x(i+1:n,1))/A(i,i) end x1(n) = ( b(n) - A(n,1:n-1) * x(1:n-1,1) )/A(n,n) k = k + 1 if sum(abs(x1 - x)) < eps fprintf('number = %d\n',k) break end if k >= max fprintf('The Method is disconvergent\n') break end x = x1 end if k < max for i = 1:n fprintf( 'x[ %d ] = %f\n',i,x1(i) ) end end 命令： clear clc format short A = [5 2 1; -1 4 2; 2 -3 10]; b = [-12 20 3]'; max = 100; eps = 1e-5 Jacobi(A, b, max, eps) 计算结果： i = 1 A = 5 2 1 -1 4 2 2 -3 10 b = -12 20 3 D = 5 0 0 0 4 0 0 0 10 L = 0 0 0 1 0 0 -2 3 0 U = 0 -2 -1 0 0 -2 0 0 0 D0 = 0.2000 0 0 0 0.2500 0 0 0 0.1000 x0 = 0 0 0 B = 0 -0.4000 -0.2000 0.2500 0 -0.5000 -0.2000 0.3000 0 f = -2.4000 5.0000 0.3000 x = -2.4000 5.0000 0.3000 x0 = -2.4000 5.0000 0.3000 i = 2 x = -4.4600 4.2500 2.2800 x0 = -4.4600 4.2500 2.2800 i = 3 x = -4.5560 2.7450 2.4670 x0 = -4.5560 2.7450 2.4670 i = 4 x = -3.9914 2.6275 2.0347 x0 = -3.9914 2.6275 2.0347 i = 5 x = -3.8579 2.9848 1.8865 x0 = -3.8579 2.9848 1.8865 i = 6 x = -3.9712 3.0922 1.9670 x0 = -3.9712 3.0922 1.9670 i = 7 x = -4.0303 3.0237 2.0219 x0 = -4.0303 3.0237 2.0219 i = 8 x = -4.0139 2.9815 2.0132 x0 = -4.0139 2.9815 2.0132 i = 9 x = -3.9952 2.9900 1.9972 x0 = -3.9952 2.9900 1.9972 i = 10 x = -3.9954 3.0026 1.9960 x0 = -3.9954 3.0026 1.9960 i = 11 x = -4.0002 3.0031 1.9999 x0 = -4.0002 3.0031 1.9999 i = 12 x = -4.0012 3.0000 2.0010 x0 = -4.0012 3.0000 2.0010 i = 13 x = -4.0002 2.9992 2.0002 x0 = -4.0002 2.9992 2.0002 i = 14 x = -3.9997 2.9998 1.9998 x0 = -3.9997 2.9998 1.9998 i = 15 x = -3.9999 3.0002 1.9999 x0 = -3.9999 3.0002 1.9999 i = 16 x = -4.0000 3.0001 2.0000 x0 = -4.0000 3.0001 2.0000 i = 17 x = -4.0000 3.0000 2.0000 x0 = -4.0000 3.0000 2.0000 i = 18 x = -4.0000 3.0000 2.0000 x0 = -4.0000 3.0000 2.0000 i = 19 x = -4.0000 3.0000 2.0000 x0 = -4.0000 3.0000 2.0000 i = 20 x = -4.0000 3.0000 2.000 x = -4.0000 3.0000 2.0000 i = 20展开阅读全文
咨信网温馨提示：1、咨信平台为文档C2C交易模式，即用户上传的文档直接被用户下载，收益归上传人（含作者）所有；本站仅是提供信息存储空间和展示预览，仅对用户上传内容的表现方式做保护处理，对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿，我们不确定上传用户享有完全著作权，根据《信息网络传播权保护条例》，如果侵犯了您的版权、权益或隐私，请联系我们，核实后会尽快下架及时删除，并可随时和客服了解处理情况，尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确)，网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据，个别因单元格分列造成显示页码不一将协商解决，平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺，下载前须认真查看，确认无误后再购买，务必慎重购买；若有违法违纪将进行移交司法处理，若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传，付费前请自行鉴别，如您付费，意味着您已接受本站规则且自行承担风险，本站不进行额外附加服务，虚拟产品一经售出概不退款（未进行购买下载可退充值款），文档一经付费（服务费）、不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印，是因预览和防盗链等技术需要对页面进行转换压缩成图而已，我们并不对上传的文档进行任何编辑或修改，文档下载后都不会有水印标识（原文档上传前个别存留的除外），下载后原文更清晰；试题试卷类文档，如果标题没有明确说明有答案则都视为没有答案，请知晓；PPT和DOC文档可被视为“模板”，允许上传人保留章节、目录结构的情况下删减部份的内容；PDF文档不管是原文档转换或图片扫描而得，本站不作要求视为允许，下载前可先查看【教您几个在下载文档中可以更好的避免被坑】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用；网站提供的党政主题相关内容(国旗、国徽、党徽－－等)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
6、文档遇到问题，请及时联系平台进行协调解决，联系【微信客服】、【QQ客服】，若有其他问题请点击或扫码反馈【服务填表】；文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“【版权申诉】”，意见反馈和侵权处理邮箱：1219186828@qq.com；也可以拔打客服电话：0574-28810668；投诉电话：18658249818。
实名认证
平台协调中心【客服】
自信AI助手
微信客服
客服QQ
发送邮件
意见反馈
链接地址：https://www.zixin.com.cn/doc/10780712.html