数值微分和数值积分.ppt

由均值定理知 可以看出，复化Simpson公式是收敛的。 定义 　　 若一个积分公式的误差满足 且C ? 0，则称该公式是 p 阶收敛的。 ~ ~ ~ 例：计算 解： 其中 = 3.138988494 其中 = 3.141592502 运算量基本相同 Lab03 复化积分 1.分别编写用复化Simpson积分公式和复化梯形积  分公式计算积分的通用程序 2.用如上程序计算积分 取节点{xi , i=0,…N}，N 为 {2k,k=0,1,…,12} ，并计算误差，同时给出误差阶 3.简单分析你得到的数据 误差阶： 记步长为h时的误差为e,步长为h/k时的误差为ek 则，相应的误差阶为： Sample Output (? represents a space) 复化梯形积分，误差和误差阶为 k=0?,?0.244934066848e00 k=1?,?0.534607244904 , 1.90 ... 复化Simpson积分，误差和误差阶为 k=1?,?0.244934066848e00 k=2?,?0.534607244904e-01 , 4.01 ... 函数变化有急有缓，为了照顾变化剧烈部分的误差，我们需要加密格点。对于变化缓慢的部分，加密格点会造成计算的浪费。以此我们介绍一种算法，可以自动在变化剧烈的地方加密格点计算，而变化缓慢的地方，则取稀疏的格点。 积分的自适应计算 ①先看看事后误差估计 以复化梯形公式为例 n等分区间 2n等分区间 近似有： 类似，复化Simpson公式 ②自适应计算 记 为复化一次，2次的Simpson公式 控制 求 是 由前面的事后误差估计式， 则， 这启发我们，可以用低阶的公式组合后成为一个高阶的公式。 类似， 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公式 做等距节点，x轴，y轴分别有： m，n为偶数 数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 第二章 数值微分和数值积分 数值微分 函数f(x)以离散点列给出时，而要求我们给出导数值， 函数f(x)过于复杂 这两种情况都要求我们用数值的方法求函数的导数值 微积分中，关于导数的定义如下： 自然，而又简单的方法就是，取极限的近似值，即差商 向前差商 x0 x0+h 由Taylor展开 因此，有误差 向后差商 x0-h x0 由Taylor展开 因此，有误差 中心差商 x0-h x0 x0+h 由Taylor展开 因此，有误差 f(x)=exp(x) h f’(1.15) R(x) h f’(1.15) R(x) 0.10 3.1630 -0.0048 0.05 3.1590 -0.0008 0.09 3.1622 -0.0040 0.04 3.1588 -0.0006 0.08 3.1613 -0.0031 0.03 3.1583 -0.0001 0.07 3.1607 -0.0025 0.02 3.1575 -0.0007 0.06 3.1600 -0.0018 0.01 3.1550 -0.0032 例： 由误差表达式，