《实用数值计算方法-5-数值积分.ppt》由会员分享,可在线阅读,更多相关《实用数值计算方法-5-数值积分.ppt(32页珍藏版)》请在三一办公上搜索。
1、计 算 方 法,授课老师:聂德明 仰仪北楼606,计量测试工程学院,Numerical Method,5 数值积分与微分,数值积分问题插值型求积公式复化求积公式,预备知识,牛顿莱布尼兹公式如果函数f(x)在区间a,b上连续,且原函数为F(x),则可用牛顿莱布尼兹公式,来求得定积分。,预备知识,积分中值定理若f是a,b上的连续函数,则存在a,b,使,预备知识,广义积分中值定理若f在a,b上连续,g在a,b上可积,且g(x)在a,b上不变号,存在a,b,使,预备知识,介值定理设函数y=f(x)在a,b上连续,介于f(a)与f(b)之间,则存在a,b,使 f()=,数值积分问题,牛顿莱布尼兹公式,找
2、原函数很困难,有些原函数不能用初等函数表示,原函数表达式过于复杂,f(x)是由测量或计算得到的数据表,预备知识,积分中值定理若f是a,b上的连续函数,则存在a,b,使,积分-代数,y,y=f(x),x,b,a,o,xk+1,xk,xk-1,数值积分问题,数值积分问题,xk为求积节点,Ak为求积系数,机械求积法,插值型,外推型,插值型积分公式,插值型积分公式,定义设有计算 的求积公式,如其求积系数,则称此求积公式为插值型求积公式。,插值型积分公式,两点公式 x0=a,x1=b,梯形公式:,梯形公式,x0,x1,y=f(x),L1(x),插值型积分公式,插值型积分公式,三点公式 x0=a,x1=(
3、a+b)/2,x2=b,辛普生公式:,辛普生公式,x0,x2,x1,y=f(x),L2(x),插值型积分公式,插值型积分公式,定义如果一个求积公式(a)对于次数不超过m的多项式均能准确成立,但至少对一个m+1次多项式不准确成立,则称该求积公式具有m次代数精度。,定理 对于求积公式(a)具有m次代数精度的充分必要条件为该公式对f(x)=1,x,.xm 精确成立,而对f(x)=xm+1,不精确成立。,插值型积分公式,确定以下求积公式的代数精度,插值型积分公式,插值型积分公式,n+1个求积节点的插值型求积公式的代数精度至少为n;,插值型积分公式,n+1个求积节点的插值型求积公式的代数精度至少为n;,
4、如果(a)式的代数精度为n,则(a)式必是插值型。,插值型积分公式,梯形公式的截断误差,预备知识,广义积分中值定理若f在a,b上连续,g在a,b上可积,且g(x)在a,b上不变号,存在a,b,使,插值型积分公式,梯形公式的截断误差,复化积分公式,梯形公式的截断误差,辛普生公式的截断误差,复化积分公式,复化求积方法,等分求积区间,比如取步长,分a,b为n等分,分点为,k=0,1,2,n,2.在区间 xk,xk+1上求,3.取和值,作为整个区间上的积分近似值,复化积分公式,复化梯形公式,复化积分公式的截断误差,复化梯形公式,预备知识,介值定理设函数y=f(x)在a,b上连续,介于f(a)与f(b)之间,则存在a,b,使 f()=,复化积分公式的截断误差,复化梯形公式,复化积分公式的截断误差,复化梯形公式,复化积分公式,复化辛普生公式,记,