[工学]数值分析第四章.pptVIP

上一页 下一页 第四章 数值积分与数值微分  在微积分里，按Newton-Leibniz公式求定积分 要求被积函数f(x) ? 有解析表达式； ? f(x)的原函数F(x)为初等函数．  ? f(x)没有解析表达式，只有数表形式 ? f(x)有表达式，但原函数不是初等函数 ，　　　　　　　　 8.5 5 8 4 6 3 4.5 2 4 1 f(x) x 定积分定义： 几何意义：曲边梯形的面积. f(x) a b f(a) f(b) 梯形公式 中矩形公式 辛甫生（Simpson）公式 （1） 希望恰当的选取 ， ，使得求积公式的 精确度很高。 可以验证：梯形公式、中矩形公式的代数精度均 为1，Simpson公式代数精度为3。 如果求积公式（1）对所有次数 的多项式 是精确的，但对 次多项式不精确，则称（1） 具有 次代数精度。 定义1 近似计算 思路 利用插值多项式 则积分易算。 ? 在[a, b]上取 a ? x0 x1 … xn ? b，做 f 的 n 次插值多项式 ，即得到 Ak 由 决定， 与 无关。 节点 f (x) 误差 插值型积分公式 例：对于[a, b]上1次插值，有 考察其代数精度。 f(x) a b f(a) f(b) 梯形公式 解：逐次检查公式是否精确成立 代入 P0 = 1： = 代入 P1 = x ： = 代入 P2 = x2 ： ? 代数精度 = 1 注：形如 的求积公式至少有 n 次代数精度 ? 该公式为插值型（即： ） ? 当节点等距分布时： 令 Cotes系数 注：Cotes 系数仅取决于 n 和 i，可查表得到。与 f (x) 及区间[a, b]均无关。 ? §2 Newton-Cotes 公式 6 5 4 3 2 1 一阶和二阶Newton-Cotes公式分别是梯形公式和辛甫生公式. n = 1: 梯形公式 /* 令 x = a+th, h = b?a, 用中值定理 */ 代数精度 = 1 n = 2: Simpson公式 代数精度 = 3 n = 3: Simpson’s 3/8公式, 代数精度 = 3, n = 4: Cotes 公式, 代数精度 = 5, n 为偶数阶的Newton-Cotes 公式至少有 n+1 次代数精度。 Newton-Cotes公式的代数精度为： ? 收敛性问题 是否对所有可积函数 都有 ？ 注意到多项式插值的Runge现象： 它的（取等距剖分点为节点）Lagrange插值多项式不收敛. 例如,数值求积过程 是发散的（见表）： 用Newton-Cotes公式计算的例子 Newton-Cotes公式在实际应用中，常用的是 的情形，特别是： 的情形.即梯形公式，Simpson公式，Cotes公式. 5.4902 3.5956 10 1.9411 8 3.3288 6 2.2776 4 2 §3 复化求积公式 高次插值有Runge 现象，故采用分段低次插值 ? 分段低次合成的 Newton-Cotes 复合求积公式。 ? 复化梯形公式： 在每个 上用梯形公式： = Tn /*中值定理*/ ? 复化 Simpson 公式： 4 4 4 4 4 = Sn 注：为方便编程，可采用另一记法：令 n’ = 2n 为偶数， 这时 ，有 ? 收敛速度与误差估计： 定义 　　 若一个积分公式的误差满足 且C ? 0，则称该公式是 p 阶收敛的。 ~ ~ ~ 例：计算 解： 其中 = 3.138988494 其中 = 3.141592502 运算量基本相同 Q: 给定精度 ?，如何取 n ? 例如：要求 ，如何判断 n = ？ ? 上例中若要求 ，则 即：取 n = 409 通常采取将区间不断对分的方法，即取 n = 2k 上例中2k ? 409