第五章 数值积分方法.pptVIP

* * 第五章 数值积分方法 5.1 插值型求积公式 5.2 复合梯形公式 5.3 其他复合求积公式 5.4 数值积分公式的代数精度与高斯求积 公式 a b f (x) 数值积分的应用背景: 1) 被积函数的原函数不能表示为初等函数 某些实际问题仅有一些离散函数值,无法给 出被积函数表达式 3) 被积函数过于复杂,难以求得其原函数 借助于被积函数在一些点的函数值,推算出满 足一定精度的定积分近似值---数值积分方法 预备知识 牛顿―莱布尼兹公式 如果函数f (x)在区间[a, b]上连续，且原函数为F(x)，则可用牛顿―莱布尼兹公式 来求定积分。 预备知识 积分中值定理 若f是[a, b]上的连续函数，则存在x∈[a, b]，使 预备知识 广义积分中值定理 若f在[a, b]上连续，g在[a, b]上可积，且g(x)在[a, b] 上不变号，存在x, x∈[a, b]，使 数值积分问题 牛顿―莱布尼兹公式 找原函数很困难，有些原函数不能用初等函数表示 原函数表达式过于复杂 f(x)是由测量或计算得到的数据表 y y=f(x) x b a o xk+1 xk xk-1 数值积分问题 5.1 插值型求积公式 f(x)在这些节点的值f(xi)，求定积分 定义 设有计算 的求积公式 如其求积系数 ，则称此求积公式为插值型求积公式. 定积分转换成被积函数的有限个函数值的线性组合，无需求被积函数的原函数. 5.1 插值型求积公式 两点公式 x0=a, x1=b, n=1 梯形公式： 5.1 插值型求积公式 一、梯形公式---两点线性插值 几何意义：用梯形面 积代替被积函数的曲 边梯形面积 梯形公式误差 5.1 插值型求积公式 广义积分中值定理 若f在[a, b]上连续，g在[a, b]上可积，且g(x)在[a, b] 上不变号，存在x, x∈[a, b]，使 利用这一定理 梯形与曲边梯形面积的对比： 正负决定 三点二次拉格朗日插值积分--辛卜生公式 x0 x2 x1 y=f(x) L2(x) 5.1 插值型求积公式 辛卜生公式: 取x0=a, x1=(a+b)/2, x2=b, n=2 辛卜生公式： 5.1 插值型求积公式 误差 精度较梯形高 y x o y= f(x) a b 5.2 复合梯形公式 分段线性插值--复合梯形法 等分求积区间，比如取步长 ，分[a, b]为n等分，分点为 ,k = 0, 1, 2,…, n 2. 在区间 [xk, xk+1]上求 3. 取和值 ，作为整个区间上的积分近似值 复合梯形公式 误差由各小区间梯形误差累加 小区间增多,误差减小→控制 … … x0 x1 x2 xk xk+1 xn-1 xn 复合梯形公式(节点加密) 复合梯形公式(节点加密) 由 递推 逐渐逼近，达到计算精度即停止。 条件成立 则终止计算并以T2n为定积分 的近似值 教材P68--例5.1 (1) 牛顿-莱布尼兹公式—0.8670 (2) 梯形公式—0.75 (3) 辛卜生公式—0.8775 (4) 复合梯形公式T4=0.8617 5.3 其它复合求积公式 借用积分中值定理 若f是[a, b]上的连续函数，则存在x∈[a, b]，使得 将其用于积分的近似计算，取ξ=b, 得 ---积分右矩形公式 复合右矩形公式 如在区间[a,b]内插入节点xj=a+jh(j=0,1,···,n), h=(b-a)/n 得到复合右矩形求积公式： 利用拉格朗日中值定理 求右矩形公式的误差估计 复合右矩形公式 复合辛卜生公式 记 每2个节点间增加一个中值节点, 节点数由n→2n. 节距变为h=(b-a)/2n. 展开, 得 利用数据表 2 2.26549 2.4600 2.8764 3.2000 3.5068 3.7647 3.93846 4 f (xk) 1 7/8 3/4 5/8 1/2 3/8 1/4 1/8 0 xk 计算积分 复合求积方法比较 取n = 8用复合梯形公式 = 取n=4, 用复合辛卜生公式 复化求积方法 定义 如果一个求