第2章_数值微分与数值积分 计算方法.pptVIP

2.3.3 复化积分的自动控制误差算法 在计算中常用 来估计误差。 （1） 的计算公式 对于 取分点n=1, h1=(b-a), 得： 对分(取分点n=2), h2=(b-a)/2 得： * 再对分(取分点n=4), h4=(b-a)/4 得： * ……，再对分(取分点n=2m), h2m=(b-a)/(2m), 得： 或： 其中 * （2） 计算公式的误差控制 比较 若要 只需 后者比前者便于控制 复化梯形积分的算法 1. 输入误差控制精度tol，初始分点数n，端点a,b, h=(b-a)/n 2. 建立函数文件f(x) 3. T2=(f(a)+f(b))*h/2; T1=T2+100; 6. S=Sum(S1); T2=T1/2+h*S; Clear S1 end 5. For i=1:n/2 计算 S1(i)=f(a+(2*i-1)*h) end 4. While |T2-T1|tol T1=T2 h=h/2;n=2n;%区间分半 * 7. 输出T2 计 算 方 法 课 件 §2.1 数值微分 结束 第2章 数值微分和数值积分 * 1.差商 2.1.1 差商与数值微分 导数的定义有三种： （1）向前差商 由泰勒（Taylor）展式： 向前差商的截断误差为： * （2）向后差商 向后差商的截断误差阶为： （3）中心差商 * 由泰勒（Taylor）展式： 中心差商的截断误差阶为： * 2.1.2 插值型数值微分 以 为插值点的插值多项式为 当 时 误差项： * 例2.2 给定三点 并且 求 解 过三点的插值多项式为 * 得三点公式： * 对定义在区间［a,b］上的定积分 以上公式多称为牛顿-莱布尼兹公式，F(x)为f(x)的原函数.但有时原函数不能用初等函数表示，有时原函数又十分复杂，难于求出或计算.如被积函数为: 等函数的积分都无法解决，当被积函数为一组数据时，更是无能为力. 为解决定积分的近似计算,从定积分的定义： 2.2.1 求积公式 * §2.2 数值积分 结束 * 这样就避开了求原函数的运算. 上式就叫做求积公式，Ak(k=0,1,…,n)与函数f(x)无关，叫做求积系数,显然要确定一个求积公式，要确定求积结点xk和求积系数Ak，或者说不同的求积结点和求积系数将确定不同的求积公式. 2.2.2 求积公式的余项和代数精度 一般情况下，上式两端并不相等.我们称: 为求积公式 的余项，或截断误差. 结束 * 为考查一个求积公式的误差，通常用代数精度来表示，如果一个求积公式对于不超过m次的多项式都能够精确成立(R[f]≡0)，而对m+1次以上的多项式不能精确成立，则称该求积公式的代数精度为m. 例如求积公式： 验证当f(x)=xm，m=0,1,2,3,4时，是否有R[xm]=0 结束 * 所以以上求积公式的代数精度为3. 任何一个求积公式的代数精度至少为零, 即取f(x)=1时公式应精确成立，这是求积系数应满足的起码条件，可以用它检验一个求积公式的系数的正确性. 求代数精度的另一种方法：利用余项公式 求代数精度m. 如： 若取： 则 而 有 所以求积公式的代数精度为3. 结束 * 2.2.3 插值型数值积分 由插值可知，对任一函数f(x)(包括表格形式的函数)可用一n次多项式对其插值，即 当Pn (x)为拉格朗日插值多项式时，即 结束 * 其中: 通常将公式(2.1)叫做插值型求积公式.（2.2）称为求积公式的余项或误差。 例：建立[0,2]上以x0=0,x1=0.5,x2=2为节点的 数值积分公式。 解 为便于上机计算，通常在求积