- 1、本文档共53页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
(2) 将区间 [ 0, 1] 4等分, 步长 采用复化 Simpson 公式计算, 仍然利用原来 9 个分点处的函数值, 求得 这两种方法计算量基本相同, 但所得到的结果与真值π= 3.1415926 ... 比较可以看出复化 Simpson 公式求得的结果要精确得多. 6.4 变步长的梯形公式与Romberg算法 1.变步长的梯形公式 复化求积公式称为定步长的求积公式,它对提高精度是行之有效的。但对于给定的精度,要确定一个合适的步长往往难以办到。因此实际上一般常采用变步长的求积公式。即让步长逐次折半的过程中,反复使用复化求积公式进行计算,直到相邻两次计算结果之差的绝对值小于允许精度ε的要求时终止计算,这种方法称为变步长的求积方法。 例如: 对于积分 采用变步长的梯形公式进行计算. 将区间 [ a , b] n 等分 , 步长 , 按复化梯形公式 计算时, 需调用 n +1 个函数值。 现在将 h 折半, 再将上述每个区间 [ xk , xk+1 ] 对分一次, 分点增至 2n + 1 个, 设上述小子区间的中点为 在[ xk , xk+1 ] 上用复化梯形公式并求和得 上式称为变步长的梯形公式. 即在求 T2n 时, 可以利用前面已求出的结果Tn , 剩下的仅仅需要求出 n 个新分点处的函数值. 注意:h = xk+1-xk 变步长的梯形公式的算法 Step 1. 给定精度 ? 0,m为正整数,步长h =(b - a)/2m。 即将积分区间分割成2m等份。 Step 2. 计算 这里 Step 3. 计算 这里 将每一个小子区间二等分,即步长折半。 Step 4. 如果 ,则停止,输出值 , 否则,置 m = m+1,h : = h/2 ,转到Step 3。 例6.4.1 用变步长的梯形公式计算积分 解: 对于 , 定义 f (0) = 1, 首先在区间 [0 , 1] 上用梯形公式(即步长 h = 1),求得 将 [0 , 1] 对分, 它的中点函数值 , 则有 如果 不成立,则 h := h/2 = 1/2 ,计算 (精确到10-6) 如此继续下去, 计算结果如下表 如果 不成立,则 h := h/2 = 1/4 , 继续计算 。 k k T 2k T 2k 0 1 2 3 4 5 6 7 8 9 10 0.920 735 5 0.939 793 3 0.944 513 5 0.945 690 9 0.945 985 0 0.946 059 6 0.946 076 9 0.946 081 5 0.946 082 7 0.946 083 0 0.946 083 1 从上表可看出, 将积分区间对分了10次, 求得 I 的近似值为0.9460831 (积分精确值为0.9460831. . . ), 可见收敛速度比较缓慢。 为了加速变步长梯形公式的收敛速度,我们采用外推策略。 2. Richardson外推算法 若用一个步长为 h 的函数 I1( h ) 去逼近问题 I , 设其 截断误差可表示为 为了提高逼近的精度,选取 q 为满足 的 正数, 将上式(1)中的 h 换为 qh , 则有 其中 是与h无关的常数,并且 , (1) 由 (1) 可知I1(h)逼近I的误差为 。 (2) (2) 式减上式 , 得 式(1)两端同乘以 得 (1) (2) 则 I2(h) 逼近I 误差降为 令 其中 是与h无关的常数,则有 , 如此继续。 一般地, 选取 q 为满足 的正数, 由此得到序列 则 Im+1 ( h ) 逼近 I 的误差由下面的定理给出。 定理 6.4.1 设 I1 ( h ) 逼近 I 的截断误差由下式给出 则 Im+1 ( h )逼近 I 的截断误差为 其中 是与 h 无关的常数。 这种利用序列{Im+1(h)} 逐步加速去逼近 I 的方法 称为Richardson外推算法 Richardson外推公式 3. Romberg 算法 Romberg 算法是利用变步长的梯形求积序列 外推加速来逼近积分真值的算法. 考虑积分 由复化梯
您可能关注的文档
最近下载
- 2024年昆明滇池国家旅游度假区国有资产投资经营管理(集团)有限责任公司人员招聘笔试备考题库及答案解析.docx
- 【教案】胚胎工程的理论基础教学设计高二下学期生物人教版(2019)选择性必修3.docx VIP
- 2024年新教材七年级上册道德与法治1.2《规划初中生活》教案.docx
- 小学五年级体育《蹲踞式跳远:助跑与起跳》教学设计.doc
- 残疾儿童送教上门一生一案资料.doc
- 铅火法冶炼工(稿).doc
- 与时代同行--中国电视文艺晚会的发展趋势初探.pdf
- 中国-象棋基础知识(适合儿童教学用)课件108页.ppt
- 安川伺服驱动器使用说明书..docx
- 第九章 马利冷却塔NC8400用户手册.pdf
文档评论(0)