数值微分与数值积分初次修改稿.pptVIP

例 计算 =0.9207355 =0.9397933 =0.9445135 =0.9456909 解 先求梯形值序列 * 第六十三页，共九十一页，2022年，8月28日 2 4 1 8 0.9207355 0.9397933 0.9445135 0.9456909 0.9461459 0.9560869 0.9460833 0.9460830 0.9460831 0.9460831 ? 利用只有两三位有效数字的T1, … ,T8 经过三次外推得到7位有效数字. 可见加速的效果十分显著. 用Romberg算法计算如下 * 第六十四页，共九十一页，2022年，8月28日 ? 理论依据: 复化梯形公式的余项展开. 记 ?定理 设 则 其中系数 ?k ( k=0, 1, …)是与 h 无关的常数. ? T (h) 逼近 I 的速度是 O ( h2 )阶. * 第六十五页，共九十一页，2022年，8月28日 ? 高阶Newton-Cotes公式是不稳定的, 因此不可能通过提高阶的方法来提高精度. ? 为了提高精度通常可把积分区间分成若干个子区间, 再在子区间上用低阶求积公式. 这种方法称为复化求积法. §3 复化求积公式 * 第三十一页，共九十一页，2022年，8月28日 x0 x1 x f (x) x2 h h x3 h h x4 ? 复化梯形公式 * 第三十二页，共九十一页，2022年，8月28日 ? 将积分区间[a, b]划分为n等分, 步长 h=( b ? a )/n,分点 xk= a+kh ( k=0, 1, … , n). 在每个小区间 [xk , xk+1 ] ( k=0, 1, … , n ?1)上用梯形公式： ? 复化梯形公式 * 第三十三页，共九十一页，2022年，8月28日 = Tn( f ) ? 复化梯形公式 * 第三十四页，共九十一页，2022年，8月28日 ? 复化梯形公式的截断误差 在小区间[xk , xk+1]上, 梯形公式的误差为 当 时 * 第三十五页，共九十一页，2022年，8月28日 由连续函数介值定理知 ? ??[a, b]使得 于是复化梯形公式的余项为 ? 当 时, 可以看出误差是 h2 阶的, 且 即复化梯形公式是收敛的. * 第三十六页，共九十一页，2022年，8月28日 ? 事实上只要设 f (x) ? C[a, b], 则可得到收敛性, 当n??时, 上式右端括号内的两个和式均收敛到积分 所以复化梯形公式是收敛的. * 第三十七页，共九十一页，2022年，8月28日 设有求积公式 设计算 时产生误差 实际得到 即 记 记 若存在与? 无关的 M 0, 使得 则称该求积公式是稳定的. ? 求积公式的稳定性 * 第三十八页，共九十一页，2022年，8月28日 ? 若求积公式是稳定的, 则 f (x)的观察值的较小的误差引起的求积结果的误差也是较小的. 求积公式没有把 f (x)的误差“放大”很多. * 第三十九页，共九十一页，2022年，8月28日 证明 因此复化梯形公式是数值稳定的. 当 ? 定理 复化梯形公式是数值稳定的. * 第四十页，共九十一页，2022年，8月28日 x0 x2 x f (x) x4 h h xn?2 h xn …... h x3 x1 xn?1 ? 复化Simpson公式 分片二次多项式近似 * 第四十一页，共九十一页，2022年，8月28日 ? 将积分区间[a, b]划分为n=2m等分, 步长 h=( b ? a )/n, 分点 xk= a+kh ( k=0, 1, … , n). 在每个小区间 [x2k ?2 , x 2k ] ( k=1, … , m)上用Simpson公式： ? 复化Simpson公式 k=1, … , m * 第四十二页，共九十一页，2022年，8月28日 = Sn( f ) ? 复化Simpson公式 * 第四十三页，共九十一页，2022年，8月28日 当 f (x)在[a, b]上具有四阶连续导数时, 故得 ? 复化Simpson公式的截断误差 * 第四十四页，共九十一页，2022年，8月28日 ? 由复化Simpson公式的截断误差知, 误差阶为 h4, 收敛性是显然的, 事实上,只要 f (x)?C[a, b]则可得到收敛性, 即 ? 由于求积系数均为正, 与复化梯形公式一样的证法可得复化