计算方法-4.pptVIP

华长生制作     举例（一） 复合求积公式 1、复合梯形公式及余项 2.复合simpson公式及余项 举例（二） 解： 例：设 ，利用下表中的数据分别用复化梯形公式和复化simpson公式计算积分 0.936 5/8 0.909 6/8 0.990 2/8 0.977 3/8 0.841 0.877 0.954 0.997 1 f (xi ) 1.0 7/8 1/8 0 4/8 xi 与准确值0.9460837比较，复化梯形3位有效数字，复化辛普森6位 一 龙贝格公式 用梯形公式得到的积分近似值 的误差大致是 ,因此人们期望, 如果用这个误差作为对 的一种补偿,则得到 的求积公式 的代数精度会有所提高。 (1) 通过直接验证可知 (2) 也就是说 , 用梯形公式二分前后的两个积分值 与 按照公式 (1) 做线形组合,其结果正好是用抛物线公式得到的积分值 。 用抛物线公式得到的积分近似值 的误差大致是 ,因此对抛物线公式进行修正, 得到 (3) 也就是说 , 用抛物线公式二分前后的积分值 与 按照公式(3) 作线形组合, 其结果正好是用柯特斯公式得到的积分值 。 通过直接验证可知 (4) 用柯特斯公式得到的积分近似值 的误差大致是 ,因此, 对柯特斯公式进行修改,得到求积公式 (5) 为此,构造求积公式 (6) 称(6)式为龙贝格(Romberg)公式。 计算f(a) , f(b) , 算出 。 (2) 把[ a , b] 2等分, 计算 , 算出 与 。 (3) 把[ a , b] 4等分,计算 算出 与 。 龙贝格求积的计算步骤如下： 龙贝格公式是一种计算积分的方法。在变步长的求积过程中, 运用(2) , (4) , (6) 式可以将精度低的梯形值 逐步加工成精度较高的抛物线 ,柯特斯值 与龙贝格值 。 。 (5) 把[ a , b] 16等分, 计算 (4) 把[ a , b] 8 等分, 计算 算出 与 与 。 算出 与 与 , 继续重复进行, 直到 时停止计算, 就是所求的积分值. (允许误差) 用龙贝格求积法计算积分 的近似值,要求准确到小数点后第5位。 解: 例 1 (1) (2) (3) * * 用插值多项式 P(x) 代替被积函数或者被微分函数f(x) , 从而导出计算 f(x) 的定积分和微分近似值的公式。 1 Newton-Cotes公式 2 复合求积法 3 Romberg算法 4* Gauss求积法 5 数值微分 利用等距节点的Lagrange插值 多项式建立插值型求积公式 将求积区间等分为多个子区间,然后在每个子区间用低阶求积公式 利用复合梯形公式,并采用加速技术 用插值多项式代替被求微分的函数 本章要点 对于任意函数的定积分,是否都可以利用Newton-Leibniz 公式求得？ 如果已知 f(x) 的原函数 F(x), 则由 Newton-Leibniz 公式： 可求得式(1)的积分值. (1) f(x)的解析式不存在,只给出了 f(x) 的一些数值； ? f(x)的原函数 F(x) 无法求出,例如 F(x) 不是初等函数； ? f(x)的表达式结构复杂,求原函数较困难. ? 对于积分 解决方案：建立积分近似计算 —— Newton-Leibniz 公式失效的几种情况： 利用插值多项式构造数值求积公式 如何解决上述问题？ 4.1 Newton-Cotes 公式 1-1 插值型求积公式及Cotes系数 使用 Newton-Cotes 公式求解数值积分问题的步骤: 则 f(x) 可表示为 ： (2) 1.将积分区间 [a, b] n 等分, 令 , 得等距节点： 等距节点下使用 Lagrange 插值多项式建立的数值求积公式. Newton-Cotes 公式 :