6.数值积分-插值型.pptVIP

例4.6 求证 不是插值型的 证明: 设 x0 = -1, x1 =0, x2 =1, A0 =1/2, A1=1, A2=1/2 则以这三点为插值节点的Lagrange插值 基函数为 例4.7 给定求积公式 试确定求积系数A-1, A0 ,A1, 使其有尽可能高的代数精度，并指出其代数精度 解：令求积公式对f(x)=1, x, x2准确成立，则有 解之得 其代数精度至少为2,将f(x)=x3代入求积公式两端相等,而将f(x)=x4代入求积公式两端不相等,所以其代数精度为3次 例 4.8 确定求积公式 使其具有尽可能高的代数精度 解：不妨设a=0, b=h, b-a=h, 设所求公式的代数 精度为2,则当f(x)=1,x,x2时公式变成等式,即 其中h=b-a, 令f(x)=x3代入上式, 两端不等, 说明求积公式只有2次代数精度。 解之得： 构造插值求积公式有如下特点： 复杂函数f(x)的积分转化为计算多项式的积分 求积系数Ak只与积分区间及节点xk有关，而与被积函数f(x)无关，可以不管f(x)如何，预先算出Ak的值 n+1个节点的插值求积公式至少具有n次代数精度 求积系数之和 可用此检验计算求积系数的正确性 (1) 在积分区间[a,b]上选取节点xk (2) 求出f(xk)及利用 或解关于Ak的线性方程组求出Ak，这样 就得到了 (3) 利用f(x)=xn,…验算代数精度 构造插值求积公式的步骤 例4.9 对 构造一个至少有3次代数精度 的求积公式 解: 3次代数精度需4个节点, 在[0,3]上取0,1,2,3四个 节点构造求积公式 确定求积系数Ak(k=0,1,2,3),利用求积系数公式 因为求积公式有4个节点，所以至少具有3次代数精度，只需将f(x)=x4代入来验证其代数精度。将f(x)=x4代入两端不相等，所以只有3次代数精度 求积公式的收敛性和稳定性 第四章 数值积分（Numerical Integration） 内容提纲 数值积分的必要性 求积公式及其代数精度 插值型求积公式 Newton-Cotes公式及数值稳定性 复化求积公式及误差估计  数值积分的必要性 本章主要讨论如下形式的一元函数积分 在微积分里，按Newton-Leibniz公式求定积分 要求被积函数f(x) ? 有解析表达式； ? f(x)的原函数F(x)为初等函数． 实际问题 f(x)的原函数F(x)不能用初等函数表示 例如函数: 考虑一个实际问题: 建筑上用的一种铝制波纹瓦是用一种机器将一块平整的铝板压制而成的. 假若要求波纹瓦长4英尺,每个波纹的 高度(从中心线)为1英寸,且每个波纹以近 似2π英寸为一个周期. 求制做一块波纹 瓦所需铝板的长度L. 这个问题就是要求由函数f(x)=sin x给定 的曲线,从x=0到x=48英寸间的弧长L. 由微积分学我们知道,所求的弧长可表示为: 上述积分称为第二类椭圆积分,它不能用普通方法来计算. 2. 有些被积函数其原函数虽然可以用初等函数表示成有限形式,但表达式相当复杂,计算极不方便.例如函数 并不复杂,但它的原函数却十分复杂:  3.f(x)没有解析表达式，只有数表形式: x 1 2 3 4 5 f(x) 4 4.5 6 8 8.5 这些都说明,通过原函数来计算积分有它的 局限性,因而,研究关于积分的数值方法具有很 重要的实际意义. 求积公式及其代数精度 求积公式的概念 积分值　　　　　　　　　 在几何上可解释为由x=a, x=b, y=0和 y=f(x) 所围成的曲边梯形的面积.积分计算之所以有困难，就是因为这个曲边梯形有一条边y=f(x) 是曲的. 　　 依据积分中值定理,对于连续函数f(x) ,在 [a,b]内存在一点ξ,使得 称f(ξ)为区间[a,b]的平均高度. 问题在于 点ξ的具体位置一般是不知道的.这样,只要 对平均高度f(ξ)提供一种算法,相应地便获 得一种数值求积方法. —— 左矩形求积公式 —— 右矩形求积公式 —— 中矩形求积公式 此外,众所周知的梯形公式: I(f