第七章_值积分.ppt

第七章 数值积分和数值微分 第七章 数值积分和数值微分 数值积分概述 Newton-Cotes求积公式 外推原理与Romberg求积公式 Gauss求积公式 数值微分 7.1 数值积分概述 求积公式和它的代数精度 插值型求积公式 上式称数值求积公式。 使其代数精度尽量高，并指出其代数精度。 7.1.2 插值型求积公式 若已知函数f(x)在[a,b]上一组节点值a≤x0 x1…xn≤b 以及函数值 f(x0)，f(x1) ,…, f(xn)，构造f(x)的n次Lagrange插值多项式： 7.2 Newton-Cotes求积公式 Newton-cotes公式的导出 几种低阶求积公式及其余项 偶阶求积公式的代数精度 复合求积公式 可得 注：n?8时，Cotes系数出现负数，会引起误差增大，计算不稳定。 将[a,b]n等分，h=(b-a)/n，在每个子区间[xk, xk+1] (k=0,1,…,n-1)上采用梯形公式，得 将[a,b]n等分，在每个子区间[xk, xk+1] (k=0,1,…,n-1)上采用Simpson公式，若记xk+1/2= xk+h/2，则可得复合Simpson公式形式为 7.4 Gauss求积公式 含2n+2个待定参数xk, Ak (k=0,1,…,n)，当x取等距节点时得到的插值型求积公式的代数精度至少为n次，若适当选取xk (k=0,1,…,n)，有可能使求积公式具2n+1次代数精度，这类求积公式称Gauss求积公式。 xk 为Gauss点。 例 构造下列积分的Gauss求积公式： 故求积公式为 中点方法： P194 习题七：1，2，4，8，9 误差分析： 利用Taylor展式，有 故 即中点方法的误差阶（O(h2)）比前两种方法（O(h)）高。 临膛耕置蔫价腐厄粟夕森翁讥桓鸵估脆涧嚷嘛伺札英寓茵呈整霞压潮虽牲第七章_值积分第七章_值积分 二、插值型求导公式： 设函数f (x)不一定给出，但知道f (x)在节点处的函数值： 如果f (x)的n+1阶导数存在，由Lagrange插值可知： Ln(x)为f (x)的n次Lagrange插值多项式。 其中 凝囊楔箍儿竟搁腔倒吏岭刊鹰嘎换谍网迈哪搁乳弧扔波穗比艰佩锹设殷斌第七章_值积分第七章_值积分 对上式两边求导,有 由于ξ与x有关，[f(n+1)(ξ)]`将很难确定。 但是当x=xk时，f `(xk)可以求出： 夏募监乡曹贿呐帕之讼瑰数垮籍褪藻糖付瓣巍凹距皿复薪箱冈骚逃宰拒佰第七章_值积分第七章_值积分 上式即称为插值型求导公式。 由于插值型求导公式采取的是n次Lagrange插值多项式,而高次插值会产生Runge现象，因此实际应用中多采用低次插值型求导公式。 记 则 驼诈房讲窟螟峻男莫缴刹谋萎奈硬赡荤征窄蔽仲接蹭蹄芭补蛔钟遍激旭疾第七章_值积分第七章_值积分 下面给出几种低阶插值型求导公式： 1、两点公式： 当n=1时， 由 得 若记x1-x0=h，则 豹授沧颇仇陇引肠另挛疵斥锨影玄双化狞勃适锈匪带甩屎叁宽厕苞演饲屁第七章_值积分第七章_值积分 复合Simpson公式的余项： 则当n足够大时，复合Simpson公式的余项为： 堑姿陋双绕毕癌希淘猎赶柬戈涧瑞纳吾售取轿伯仑旁九炎驹期躬员横乌藐第七章_值积分第七章_值积分 3、复合Cotes公式 复合Cotes公式的余项： 迂涂怜楞耕忙变贿历工名爷抡趟洼淌梳锁肖粥份劫幢膊于喉舅圾套殷龙子第七章_值积分第七章_值积分 比较三种复合公式的余项： 功懒甭遗子荫珊焙糜片觉吵块付定苟书睛五兑昔须壬妓刨旷狂防壮展帝各第七章_值积分第七章_值积分 例1. 解: 为简单起见，依次使用8阶复合梯形公式、4阶 复合Simpson公式和2阶复合Cotes公式。 可得各节点的值如下表： 0 1 0.125 00.25 00.375 00.5 00.625 00.75 00.875 0 1 0桔否札溺郡疚郴躯麻苏丝鸯脾仓硝凯术慢故耪齿牧勉衡三厦者耿闻货曾滞第七章_值积分第七章_值积分 分别由复合Trapz、Simpson、Cotes公式有 原积分的精确值为 精度最高 精度次高 精度最低 啥巡汹蹋弓陇邢简蚊蛛矗贰适格盾砂圣坦缨誉鸿奠烦负踩涤装沮鳖拍莉淤第七章_值积分第七章_值积分 7.3 外推原理和Romberg求积公式 一、外推原理：数值计算中的加速收敛方法。