【精选】数值分析教案4.doc

§1 插值型数值求积公式 教学目的 1. 会求插值型数值求积公式及Gauss型数值求积公式并会讨论它们的代数精度； 2. 理解复化梯形数值求积公式及复化Simpson数值求积公式和余项的推导的基础上掌握它们； 3. 理解数值微分公式推导的基础上掌握一阶、二阶数值微分公式及余项； 4. 了解外推原理。 教学重点及难点 重点是插值型数值求积公式及Gauss型数值求积公式的求解及它们代数精度的讨论；难点是Gauss型数值求积公式节点的求解方法的推导及求解方法。 教学时数 12学时 教学过程 1．1一般求积公式及其代数精度 设是上的权函数，是上具有一定光滑度的函数。用数值方逑下积分 的最一般方法是用在节点上函数值的某种线性组合来近似 其中是独立于函数的常数，称为积分系数，而节点称为求积节点。 我们也可将（1．2）写成带余项的形式 （1．2）和（1．3）都称之为数值求积公式或机械求积公式。更一般些的求积公式还可以包含函数在某些点的低阶导数值。 在（1.3）中余项也称为求积公式的截断误差。 一个很自然的想法是数值求积公式要对低次多项式精确成立这就导出了求积公式数精度的概念。 定义1 若求积公式（1.2）对任意不高于次的代数多项式都精确成立，而对不能精确成立，则称该求积公式具有次代数精度。 一个求积公式的代数精度越高，就会对越多的代数多项式精确成立。 确定求积公式 的代数精度。 解 。 从而该求积公式的代数精度为。 对给定节点，如何选择求积系数使求积公式代数精度尽可能高，对此可用插值型求积公式来实现。 2 插值型求积公式 对给定求积节点构造求积公式的一种简单方法是利用插值多项式的准许确积分来作为数值积分值。 设是关于的Lagrange插值多项式 其中 为Lagrange基函数。取 其中 。 定义2 对给定互异求积节点，若求积系数是由（1.4）给出的，则称该求积公式是插值型的。 定理1 数值求积公式（1.2）或（1.3）是插值型的当且仅当它的代数精度。 证明 假设求积公式（1.2）是插值型的，则 上面我们假设了。从而当为次数的代数多项时必精确成立，故有 。 假设。注意到多项式的次数为，对=数值求积精确成立，从而 即其求积系数由（1．4）给出。 推论1 对给定求积节点，代精度最高的求积公式是插值型求积公式。 求插值型求积公式 并确定其代数精度。 解 。 从而求积公式为 且。 对 从而。 若我们利用Hermite插值多项式的准确积分作为数值积分值，我们可以类似地建立带有函数在某些节点导数值的插值型求积分式。 推论2 若是插值型求积公式，则有余项公式 其中 1．3 Newton-Cotes 求积公式 在[a,b]上的插值型求积公式应用最方便、最广泛，称之为Newton-Cotes求积公式。 设令则求积系数为 其中 因此，Newton-Cotes公式为 其中 由（1．6）给出。 求职系数独平于区间[a,b] 称之为Cotes系数。Cotes系数可以用（1．6）计算或查（见表4-1）给出。 n=1,2 的Newton-Cotes求积是常用公式。n=1的公式称为梯形公式，其几何意义是用直边梯形的面积来近似曲边梯形面积（图4-1）。即 表4-1 （1.8） 的Newton-Cotes公式称为Simpson公式: （1.9） Simpson公式的几何意义是用以插值抛物线为曲边的曲边梯形面积来近似为曲边的曲边梯形面积（如图4—2），因此Simpson求积公式也称为抛物线公式。 Newton—Cotes公式分别为Simpson法则（公式）和Cotes公式。 1.4 Newton—Cotes求积公式的余项 定理２　　若，则梯形公式(1.8)的余项为 　　　　　　　　　　　　　 (1.10) 证明　　由插值型求积公式的余项得 利用在上不变的号，由积分中值定理得 定理３　　若，则Simpson公式(1.9)的余项为 　　　　　　　　　　　　　　　　　　 (1.11) 证明　　由例１知Simpson公式的代数的精度为。令为的三次Hermite插值多项式，满足插值条件： 对多项式，Simpson公式精确成立，即： 从而 利用上小于等于零，由积分中值定理给出 　　　　 可以证明，对一般的，只要充分光滑，Newton—Cotes公式的余项为 　　 （1．2） 其