数值分析 第7讲.pptVIP

故由介值定理,一定在(a,b)有一点 使 Remark：若 ，则有误差估计式 证毕 2.复化辛浦生公式 将[a,b]等份成N个子区间[xk，xk+1](k＝0，1，…,N-1)，子区间长度 由 3.复化柯特斯公式 将[a,b]等份成N个子区间[x4k，x4k+4](k＝0，1，…,N-1)，子区间长度 由 例子1 若取9个节点，用复化梯形公式、复化辛浦生公式和复化柯特斯公式计算积分，其步长以及与9个节点所对应的求积系数分别是多少？ 解：复化梯形公式：N=8，h=(b-a)/8，对应的求积系数为1、2、2、2、2、2、2、2、1。 复化辛浦生公式：N=4，h=(b-a)/4，对应的求积系数为1、4、2、4、2、4、2、4、1。 复化柯特斯公式：N=2，h=(b-a)/2，对应的求积系数为7、32、12、32、14、32、12、32、7。 # 用积分 计算ln2,要使所得积分近似值具有5位有效数字。问用复化梯形公式，复化Simpson公式时，至少要取多少个节点? 例子2 解:由 且 故,计算ln2时，要使误差不超过 也即计算2ln2，其误差不超过 。 即 其中 故区间应取671个，节点至少应取672。 其中 由 故区间N应取22，即45个节点 。 # 二、区间逐次分半求积法（变步长） ? 1.误差的事后估计法 复化求积公式是提高精度的一种有效方法，但在使用复化求积公式之前，必须根据复化求积公式的余项进行先验估计，以确定节点数目，从而确定合适的等分步长。因为余项表达式中包含了被积函数的导数，而估计各阶导数的最大值往往是很困难的，且估计的误差上界往往偏大。所以实际中，常常使用“事后估计误差”的方法，通过区间的逐次分半，在步长逐次分半的过程中，反复利用复化求积公式进行计算，查看相继两次计算结果的差值是否达到要求，直到所求得的积分值满足精度要求。 该法也称为步长自动选择的变步长求积法。 对梯形公式，假定区间分为N等份时，由公式 算出的积分近似值为TN ,因而有 再把各个小区间分别对分，得积分的近似值为T2N，则积分值为 其中 假定f??(x)有[a,b]上变化不大，即有 ， 则 上式可改写为 计算时只需检验 是否满足？若不满足，则再把区间分半进行计算，直到满足要求为止 。 类似的，还可以得到下面的结论： 对于辛浦生公式，假定 在[a,b]上变化不大，则有 对于Cotes公式，假定 在[a,b]上变化不大，则有 ? 2.区间逐次分半的梯形公式（P212） 据此我们得到复化梯形公式区间逐次分半时的递推计算公式： 计算时只需检验 是否满足？若不满足，则再把区间分半进行计算，直到满足要求为止 。 §7.3 Romberg求积法 一、对低精度公式经过组合构造高精度公式 从SN及T2N,TN的计算公式可验证得到： 事实上， （1） （2） （3） 类似地，可以证明： 这个公式（3）称为Romberg公式。 由（1）（2）（3）组成的方法称为Romberg算法。 序列{TN },{SN },{CN }和{RN }分别称为梯形序列，Simpson序列，Cotes序列和Romberg序列。 上述用若干个积分近似值推算出更为精确的积分近似值的方法，称为外推算法。得到Romberg序列后还可以继续外推，得到新的求积序列，称为Richardson外推算法。但由于在新的求积序列中，其线性组合的系数分别为： 因此，新的求积序列与前一个序列结果相差不大，故通常外推到Romberg序列为止。 可以证明，梯形序列，Simpson序列，Cotes序列和Romberg序列均收敛到积分值，且每次外推可使误差阶提高二阶。 二、Romberg算法的实现 T数表： R4 C8 S16 T32 32 R2 C4 S8 T16 16 R1 C2 S4 T8 8 C1 S2 T4 4 S1 T2 2 T1 1 R2k-3 C2k -2 S2k-1 T2k 区间等分数 n=2k 对上面的T数表作计算，一直到Romberg序列中前后两项之差的绝对值不超过给定的误差限为止。 Remark:Romberg算法具有占用内存少，精确度高的优点，是实际中最常用的算法之一。 §7.4 Gauss型求积公式 我们能否通过节点的选择将求积公式的代数精度从n 或者n+1提高到2n+1？ 问题： 若求