第二章节数值微分和数值积分.ppt

? Gauss 公式的余项： /* 设P为f 的过x0 … xn的插值多项式 */ /*只要P 的阶数不大于2n+1，则下一步等式成立*/ 插值多项式的余项 Q：什么样的插值多项式在 x0 … xn 上有 2n+1 阶？ A：Hermite 多项式！ 满足 n xk Ak n xk Ak 2 0.5858864376 3.4142135623 0.8535533905 0.1464466094 5 0.2635603197 1.4134030591 3.5964257710 7.0858100058 12.6408008442 0.5217556105 0.3986668110 0.0759424497 0.0036117587 0.0000233700 3 0.4157745567 2.2942803602 6 0.7110930099 0.2785177335 0.0103892565 6 0.2228466041 1.1889321016 2.9927363260 5.7751435691 9.8374674183 15.9828739806 0.4589646793 0.4170008307 0.1133733820 0.0103991975 0.0002610172 0.0000008985 4 0.3225476896 1.7457611011 4.5366202969 9.3950709123 0.6031541043 0.3574186924 0.0388879085 0.0005392947 1、梯形公式 此处用了积分中值定理 误差 2、Simpson公式 注意到，Simpson公式有3阶代数精度，因此为了对误差有更精确地估计，我们用3次多项式估计误差 为0 一般的有 因此，N-C积分，对偶数有n＋1阶代数精度，而奇数为n阶代数精度 复化积分 数值积分公式与多项式插值有很大的关系。因此Runge现象的存在，使得我们不能用 太多的积分点计算。采用与插值时候类似，我们采用分段、低阶的方法 误差 做等距节点， 复化梯形公式 由均值定理知 可以看出，复化梯形公式是收敛的。 误差 做等距节点， 复化Simpson公式 由均值定理知 可以看出，复化Simpson公式是收敛的。 定义 　　 若一个积分公式的误差满足 且C ? 0，则称该公式是 p 阶收敛的。 ~ ~ ~ 例：计算 解： 其中 = 3.1 其中 = 3.141592502 运算量基本相同 3、积分的自适应计算 函数变化有急有缓，为了照顾变化剧烈部分的误差，我们需要加密格点。对于变化 缓慢的部分，加密格点会造成计算的浪费。以此我们介绍一种算法，可以自动在变化剧烈的地方加密格点计算，而变化缓慢的地方，则取稀疏的格点。 ①先看看事后误差估计 （不同的误差表达式，事后误差估计式是不同的） 以复化梯形公式为例 n等分区间 2n等分区间 近似有： 类似，复化Simpsom公式 ②自适应计算 记 为复化一次，2次的Simpson公式 控制 求 是 4、Romberg积分 由前面的事后误差估计式， 则， 这启发我们，可以用低阶的公式组合后称为一个高阶的公式。 类似， 记 为以步长为h的某数值积分公式，有 有如下的Euler-Maclaurin定理 若 为2m阶公式，则 Romberg积分就是不断地用如上定理组合低阶公式为高阶公式，进而计算积分 ? Romberg 算法： < ? ? < ? ? < ? ? … … … … … … ? T1 = ) 0 ( 0 T ? T8 = ) 3 ( 0 T ? T4 = ) 2 ( 0 T ? T2 = ) 1 ( 0 T ? S1 = ) 0 ( 1 T ? R1 = ) 0 ( 3 T ? S2 = ) 1 ( 1 T ? C1 = ) 0 ( 2 T ? C2 = ) 1 ( 2 T ? S4 = ) 2 ( 1 T 重积分的计算 在微积分中，二重积分的计算是用化为累次积分的方法进行的。计算二重数值积分 也同样采用累次积分的计算过程。简化起见，我们仅讨论矩形区域上的二重积分。对 非矩形区域的积分，大多可以变化为矩形区域上的累次积分。 a,b,c,d为常数，f在D上连续。将它变为化累次积分 首先来看看复化梯形公式的二重推广 做等距节点，x轴，y轴分别有： 先计算 ，将x作为常数，有 再将y作为常数，在x方向，计算上式的每一项的积分 二重积分的复化梯形公式 系数，在积分区域的四个角点为1/4，4个边界为1/2，内部节点为1 误差 类似前面有： 记 二重积分的复化Simpson公式 做等距节