- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
* * * * * * * * * * * * * * *我们看到,中点公式具有如下形式(8.7) 式中的系数均与步长无关。若将步长二分,则有(8.8) 则可消去误差主项,得取(8.7)与(8.8)加权平均 (8.9) 重复同样的手续,再导出下列加速公式这种加速过程还可继续下去。这种加速方法通常称作Richardson外推加速法。若令(8.10) 则进一步消去D1(h)误差主项, 有例8.2 运用加速公式和加工例8.1的结果。 表8.2 Richardson外推加速法计算结果 hD(h) 0.8 0.4 0.2 0.1 3.01765 2.79135 2.73644 2.72281 2.715917 2.718137 2.718267 2.718285 2.718276 2.71828解 计算结果见表8.2。这里,加速的效果同样是相当显著的。 注: 中点心公式及其加速方法适合用表达式表示的函数。对于列表函数,则宜使用插值方法等导出数值求导公式。 最小,步长h不宜太大,也不宜太小. 其最优步长应为它表明h 越小, 舍入误差 越大, 故它是病态的. 用中点 公式计算 的误差上界为 要使误差E(h) 则计算当 f (a+h) 及 f (a-h) 分别有舍入误差 ?1 及 ?2 时,若令 的舍入误差上界为 二、插值型的求导公式xx0x1x2…xnyy0y1y2…yn对于列表函数 y = f (x):插值多项式 y = Pn(x)作为它的近似, 我们取 统称插值型的求导公式. 的近似值, 作为 建立的数值公式依据插值余项定理,求导公式(8.11)的余项为 式中 (8.11)我们限定:求某个节点 xk 上的导数值,上面的第二项变为零,这时有余项公式 (8.12) 1.两点公式已给两节点 x0, x1 上的函数值 f (x0), f (x1),做线性插值 记 x1 – x0 = h,对上式两端求导,有 于是有下列求导公式:下面我们仅仅考察节点处的导数值.为简化讨论, 假定所给的节点是等距的. 而利用余项公式(8. 12)知,带余项的两点公式是:设已给出三节点x0, xl=x0+h, x2=x0+2h上的函数值,做二次插值 令 x = x0+th, 则 2.三点公式 (8.13) 这里撇号(’)表示对x求导.上式分别取 t = 0,1,2, 得到三种三点公式: (8.14) 而带余项的三点求导公式如下: 公式(8.6)是我们所熟悉的中点公式.在三点公式中,它由于少用了一个函数值 f(x1) 而引人注目. 例如,将式(6.5)再对t求导一次,有 于是有用插值多项式Pn(x)作为 f (x)的近似函数,还可以建立高阶数值微分公式: 而带余项的二阶三点公式如下: (8.15)例 已知函数y= ex 的下列数值:x2.52.62.72.82.9y 12.1825 13.4637 14.8797 16.4446 18.1741 试用二点、三点微分公式计算x=2.7处的一阶、二阶导数值.解: 方法1: 取h=0.1时,两点公式有两种取法当 x0=2.6, x1=2.7时,f’(2.7)[f(2.7)- f(2.6)]=[14.8797-13.4637]=14.1600 0.1 1 0.1 1当x 0= 2.7,x1= 2. 8时,f’(2.7)[f(2.8)- f(2.7)] =15.6490. 0.1 1方法2: 取h=0.2时,两点公式有两种取法当 x0=2.5, x1=2.7时,f’(2.7)[f(2.7)- f(2.5)]=13.4860 当 x0=2.7,x1=2.9时,f’(2.7)[f(2.9)- f(2.7)]=16.4720 0.2 1 0.2 1 三点公式取x0= 2.6,x1=2.7,x2=2.8,得f’(2.7)[f(2.8)- f(2.6)] =14.9045.f’’(2.7)[f(2.8)- 2f(2.7)+ f(2.6)] =14.8900. 0.12 1 2?0.1 1注: f’(2.7)和 f’’(2.7 )的真值都是 14.87973…,上面的计算表明:1)当使用两点公式时.应取步长较小的函数值;2)一般情况下,同样步长的两点公式没有三点公式准确,步长越小越精确但如果高阶导数无界或舍入误差超过截断误差时,这个结论就不一定对了. 三点公式取x0= 2.5,x1=2.7,x2=2.9,得f’(2.7)[f(2.9)- f(2.5)] =14.9790.f’’(2.7)[f(2.9)- 2f(2.7)+ f(2.5)] =14.9300. 0.22 1 2?0.2 1 附注: 与积分相比,数值微分比较困难。积分描述了一个函数的整体或宏观性质,而微分则描述
您可能关注的文档
最近下载
- 阿法拉伐分油机中文说明Instruction book1.pdf VIP
- 《Java程序设计》教案(思政版).docx VIP
- 四年级上册数学《三位数乘两位数的笔算乘法》(共18张PPT).pptx VIP
- AWS B4.0-2016 焊缝机械测试的标准方法(中文版).pdf
- 《绩效与薪酬管理》课程教学大纲(中文).docx VIP
- 国开电大学习网网络安全技术形考任务答案.pdf VIP
- 成长赛道飞机机电设备维修专业1500字.pptx VIP
- 从业人员食品安全知识培训计划及培训考核记录.docx VIP
- GB50641-2010 有色金属矿山井巷安装工程施工规范.docx VIP
- 胃癌根治术的手术护理查房2讲课文档.ppt VIP
文档评论(0)