上海交大数值分析课件数值分析4-1数值积分.pptVIP

第四章 数值积分与数值微分 §1 数值积分概论 一、数值求积的基本思想 二、代数精度的概念 三、插值型的求积公式 四、例题 一、数值求积的基本思想 1. 定积分及其计算 对于积分 只要找到被积函数f(x)的原函数F(x)，便有下列牛顿—莱布尼兹公式 2. 为何要进行数值积分？ 原因之一：f(x)的原函数难以求解或太复杂 例 原因之二：f(x)以离散数据点形式给出 yn xn … y1 y0 yi = f(xi) … x1 x0 xi 3. 数值积分的基本思想 角度1： 由定积分的定义 积分和式 是一些点上函数值的线性组合。 故数值积分的基本思想： 用被积函数在积分区间上某些节点xk，(k=0,…,n)处的函数值的线性组合作为定积分的近似值。即 3. 数值积分的基本思想 角度2： 由定积分的几何意义，只需计算曲边梯形的面积 x y 由此可得以下公式 梯形公式： 中矩形公式： 左矩形公式： 右矩形公式： 这些公式都可写成如下形式： 求积系数 求积节点 这类公式称为机械求积公式 注：这类公式有无穷多个 机械求积方法的特点： 将积分求值问题归结为函数值的计算，避开了牛顿—莱布尼兹公式需要寻求原函数的困难。 机械求积方法的两个要素： 求积节点xk (k=0, … , n) 求积系数Ak (k=0, … , n) 二、代数精度的概念 1. 为何引入代数精度？ 数值求积方法是近似方法，为了保证精度，我们自然希望公式能对“尽可能多”的函数准确成立，这就提出了所谓代数精度的概念。 2. 代数精度的定义 定义1 若某个求积公式对于次数≤m的多项式均能够准确成立，但对于m+1次多项式就不准确成立，则称该求积公式有m次代数精度。 指代数多项式 定义1ˊ 若某个求积公式对于1, x,…, xm 均能够准确成立，但对于xm+1就不准确成立，则称该求积公式有m次代数精度。 3. 代数精度的应用之一：判断代数精度 例 证明：梯形公式具有一次代数精度。 证明 梯形公式 当 f(x)=1 时， 左边= 右边= 故 左边=右边 公式至少具有0次代数精度 当 f(x)=x 时， 左边= 右边= 故 左边=右边 公式至少具有1次代数精度. 当 f(x)=x2 时， 左边= 右边= 故 左边≠右边 公式具有1次代数精度. 结论：中矩形公式具有一次代数精度。 3. 代数精度的应用之二：构造数值积分公式 注： (1) 从习题1(4)可以看出，数值积分公式可以设计被积函数函数值和导数值。 (2) 思考构造如下数值积分公式 从本例可知：需要求解非线性方程组，这将非常困难！ 考察书后习题1 三、插值型的求积公式 定义：插值型求积公式 设f(x)以离散数据点形式给出 yn xn … y1 y0 yi = f(xi) … x1 x0 xi 由插值条件构造插值多项式Ln(x)，并构造如下的求积公式 称作是插值型的。 插值型求积公式是机械求积公式吗？ 答： 是！因为插值型求积公式为 可看作机械求积公式的Ak取为 现在的问题是： 机械求积公式是插值型求积公式吗？ 答： 不一定！ 因为求积系数AK不一定满足 现在的问题是： 定理： 证明： 先证明充分性 由插值余项定理可知，对于插值型的求积公式，其余项为 式中ξ与变量x有关。 机械求积公式 至少有n次代数精度的充分必要条件是，它是插值型的。 若求积公式是插值型的， 则对于次数≤n的 多项式f(x)，其余项R[f]等于0，故 这时求积公式至少具有n次代数精度。 充分性得证。 再证明必要性 已知求积公式至少具有n次代数精度，则公式对于插值基函数 lk(x) 应准确成立，即有 注意到 lk(xj)=δkj，上式右端实际上等于Ak， 故求积公式是插值型的。 思考：梯形公式是插值型求积公式吗？ 试总结证明机械求积公式是插值型求积公式的方法。 四、例题 例1 判明以下两个求积公式的代数精度。 (1) (2) 解 (1)记