数值分析数值积分.pptVIP

求积公式的设计试设计求积公式第29页，共46页，星期日，2025年，2月5日第30页，共46页，星期日，2025年，2月5日第31页，共46页，星期日，2025年，2月5日第32页，共46页，星期日，2025年，2月5日例题2第33页，共46页，星期日，2025年，2月5日第1页，共46页，星期日，2025年，2月5日数值积分引言计算定积分微积分基本公式：(2)f(x)表达式未知，只有通过测量或实验得来的数据表。但是在许多实际计算问题中(1)f(x)表达式较复杂，原函数难求！甚至有时不能用初等函数表示。如此时需要利用数值方法来近似计算定积分。第2页，共46页，星期日，2025年，2月5日数值积分的几何意义第3页，共46页，星期日，2025年，2月5日数值求积的基本思想依据积分中值定理，对于连续函数f（x），在［a，b］内存在一点ξ，成立就是说，底为b-a而高为f（ξ）的矩形面积恰等于所求曲边梯形的面积I.问题在于点ξ的具体位置一般是不知道的，因而难以准确地算出f（ξ）的值.我们称f（ξ）为区间［a，b］上的平均高度.这样，只要对平均高度f（ξ）提供一种算法，相应地便获得一种数值求积方法.第4页，共46页，星期日，2025年，2月5日数值求积的基本思想分别用f(a)，f(b)和近似f(?)可得左矩形公式右矩形公式中矩形公式第5页，共46页，星期日，2025年，2月5日求积公式的基本思想若用f(a)和f(b)的算术平均值近似f(?)，则可得梯形公式若用f(a),f([a+b]/2)和f(b)的加权平均值近似f(?)，则可得辛甫生公式第6页，共46页，星期日，2025年，2月5日一般求积公式更一般地，可以用f(x)在[a,b]上的一些离散点上的值加权平均作为f(?)的近似值，从而构造出求积节点求积系数机械求积法：求积系数仅仅与结点xk的选取有关，而不依赖于被积函数f(x)的具体形式第7页，共46页，星期日，2025年，2月5日机械求积的问题描述已知n+1个x以及在这些x上的函数值求解此函数在某个区间的积分值如何衡量这个公式的好坏？第8页，共46页，星期日，2025年，2月5日代数精度定义如果对于所有次数不超过m的多项式f(x)，公式精确成立，但对于某一次数为m+1的多项式不精确成立，则称该求积公式的代数精度为m次。要验证一个求积公式具有m次代数精度，只需验证对f(x)＝1,x,x2,…,xm精确成立，但对f(x)＝xm+1不精确成立即可，即：(k=0,1,…,m)第9页，共46页，星期日，2025年，2月5日已知：求积公式对于xk（k=0，1，…，m）均能准确成立求证：求积公式对于对于次数不超过m的多项式均能准确成立证明：由已知条件知（k=0，1，…，m）证明两种说法的等价性第10页，共46页，星期日，2025年，2月5日则即：求积公式对于对于次数不超过m的多项式均能准确成立第11页，共46页，星期日，2025年，2月5日举例（一）例：试确定系数?i，使得下面的求积公式具有尽可能高的代数精度，并求出此求积公式的代数精度。解：将f(x)＝1,x,x2代入求积公式，使其精确成立得解得?0=1/3,?1=4/3,?2=1/3，所以求积公式为易验证该公式对f(x)＝x3也精确成立，但对f(x)＝x4不精确成立，所以此求积公式具有3次代数精度。第12页，共46页，星期日，2025年，2月5日矩形和梯形公式的代数精度容易验证：左矩形公式和右矩形公式具有零次代数精度中矩形公式和梯形公式具有一次代数精度特别地，具有m(?0)次代数精度的求积公式满足：辛甫生公式具有三次代数精度第13页，共46页，星期日，2025年，2月5日如何求解求积公式我们可以用代数精度作为标准来构造求积公式.譬如两点公式式中含有两个待定参数A0，A1，令它对于f（x）=1，f（x）=x准确成立，有第14页，共46页，星期日，2025年，2月5日解之得A0=A1=（b-a）/2.这说明，形如（5）且具有一次代数精度的求积公式必为梯形公式（1）.这一论断从几何角度来看是十分明显的.如何求解求积公式第15页，共46页，星期日，2025年，2月5日如何