数值积分方法.pptVIP

利用数据表xk01/81/43/81/25/83/47/81f(xk)43.938463.76473.50683.20002.87642.46002.265492计算积分复合求积方法比较取n=8用复合梯形公式=第23页，课件共36页，创作于2023年2月取n=4,用复合辛卜生公式复化求积方法第24页，课件共36页，创作于2023年2月定义如果一个求积公式(a)对于次数不超过m的多项式均能准确成立，但至少对一个m+1次多项式不准确成立，则称该求积公式具有m次代数精度。定理对于求积公式(a)具有m次代数精度的充分必要条件为该公式对f(x)=1,x,…..xm精确成立，而对f(x)=xm+1，不精确成立。5.4数值积分公式的代数精度第25页，课件共36页，创作于2023年2月求代数精度的阶数--确定以下求积公式的代数精度5.4数值积分公式的代数精度?阶代数精度1阶代数精度?阶代数精度第26页，课件共36页，创作于2023年2月关于数值积分方法第1页，课件共36页，创作于2023年2月abf(x)数值积分的应用背景:1)被积函数的原函数不能表示为初等函数某些实际问题仅有一些离散函数值,无法给出被积函数表达式3)被积函数过于复杂,难以求得其原函数借助于被积函数在一些点的函数值,推算出满足一定精度的定积分近似值---数值积分方法第2页，课件共36页，创作于2023年2月预备知识牛顿―莱布尼兹公式如果函数f(x)在区间[a,b]上连续，且原函数为F(x)，则可用牛顿―莱布尼兹公式来求定积分。第3页，课件共36页，创作于2023年2月预备知识积分中值定理若f是[a,b]上的连续函数，则存在x∈[a,b]，使第4页，课件共36页，创作于2023年2月预备知识广义积分中值定理若f在[a,b]上连续，g在[a,b]上可积，且g(x)在[a,b]上不变号，存在x,x∈[a,b]，使第5页，课件共36页，创作于2023年2月数值积分问题牛顿―莱布尼兹公式找原函数很困难，有些原函数不能用初等函数表示原函数表达式过于复杂f(x)是由测量或计算得到的数据表第6页，课件共36页，创作于2023年2月yy=f(x)xbaoxk+1xkxk-1数值积分问题第7页，课件共36页，创作于2023年2月5.1插值型求积公式f(x)在这些节点的值f(xi)，求定积分第8页，课件共36页，创作于2023年2月定义设有计算的求积公式如其求积系数，则称此求积公式为插值型求积公式.定积分转换成被积函数的有限个函数值的线性组合，无需求被积函数的原函数.5.1插值型求积公式第9页，课件共36页，创作于2023年2月两点公式x0=a,x1=b,n=1梯形公式：5.1插值型求积公式一、梯形公式---两点线性插值几何意义：用梯形面积代替被积函数的曲边梯形面积第10页，课件共36页，创作于2023年2月梯形公式误差5.1插值型求积公式广义积分中值定理若f在[a,b]上连续，g在[a,b]上可积，且g(x)在[a,b]上不变号，存在x,x∈[a,b]，使利用这一定理梯形与曲边梯形面积的对比：正负决定第11页，课件共36页，创作于2023年2月三点二次拉格朗日插值积分--辛卜生公式x0x2x1y=f(x)L2(x)5.1插值型求积公式第12页，课件共36页，创作于2023年2月辛卜生公式:取x0=a,x1=(a+b)/2,x2=b,n=2辛卜生公式：5.1插值型求积公式误差精度较梯形高第13页，课件共36页，创作于2023年2月yxoy=f(x) ab5.2复合梯形公式第14页，课件共36页，创作于2023年2月分段线性插值--复合梯形法等分求积区间，比如取步长，分[a,b]为n等分，分点为 ,k=0,