数值积分与微分方法.docx

数值积分与微分 摘要本文首先列举了一些常用的数值求积方法，一是插值型求积公式，以公式为代表，并分析了复合型的公式；另一个是求积公式，并给出几个常用的求积公式。其次，本文对数值微分方法进行分析，主要是差分型数值微分和插值型数值微分，都给出了几种常用的微分方法。然后，本文比较了数值积分与微分的关系，发现数值积分与微分都与插值或拟合密不可分。本文在每个方法时都分析了误差余项，并且在最后都给出了MATLAB的调用程序。关键词：插值型积分 差分数值微分 插值型数值微分 MATLAB常用的积分方法计算积分时，根据公式，但如果碰到以下几种情况：1）被积函数以一组数据形式表示；2）被积函数过于特殊或者原函数无法用初等函数表示3）原函数十分复杂难以计算这些现象中，公式很难发挥作用，只能建立积分的近似计算方法，数值积分是常用的近似计算的方法。插值型积分公式积分中的一个常用方法是利用插值多项式来构造数值求积公式，具体的步骤如下：在积分区间上上取一组节点：。已知的函数值，作的次插值多项式，则其中，为次插值基函数，则得公式写成一般形式：其中，显然，当被积函数为次数小于等于的多项式时，其相应的插值型求积公式为准确公式，即：求积公式的代数精度定义：求积公式对于任何次数不大于的代数多项式均精确成立，而对于不精确成立，则称求积公式具有次代数精度。定理：含有个节点的插值型求积公式的代数精度至少为。公式公式将积分区间等分，并取分点为求积公式，这样构造出来的插值型求积公式就是公式。其中，且系数满足重要的关系式：时，求积公式为梯形公式（两点公式）：梯形公式具有1阶代数精度，余项为：=2时，求积公式为公式(三点公式)：公式具有3阶代数精度，余项为：=4时，求积公式为公式（五点公式）：其中，公式具有5次代数精度，余项为：复合公式当积分区间过大时，直接使用公式所得的积分的近似值很难得到保证，因此在实际应用中为了既能够提高结果的精度，又使得算法简便且容易在计算机上实现，往往采用复合求积的方法。所谓复合求积，就是先将积分区间分成几个小区间，并从每个小区间上用低阶公式计算积分的近似值，然后对这些近似值求和，从而得到所求积分的近似值，由此得到一些具有更大实用价值的数值求积公式，统称为复合求积公式。将区间等分，记分点为，其中，称为步长，然后在每个小区间内利用梯形公式，即可导出复合梯形公式：若将所得积分近似值记为，并注意到，则复合梯形公式为：其余项为：类似可得复合公式：其中，.其余项为：公式在MATLAB中的实现1）复合梯形数值积分：调用形式：Z=trapz(X,Y)其中，X,Y分别代表数目相同的向量或者数值，Y与X的关系可以是函数形态或者离散形态；Z代表返回的积分值。2）自适应公式基本调用格式：q=quad（fun，a，b，tol，trace，p1，p2）其中：fun代表被积函数；a，b为积分的上下限；q为积分结果；tol为默认误差限，默认了1.e-6；trace表示取0表示不用图形显示积分过程，非0表示用图形显示积分过程；p1，p2为直接传递给函数fun的参数3）自适应Lobatto法数值积分：quadl（） Quadl是高阶的自适应数值积分法函数，比quad函数更有效，精度更高，使用方法与quad完全相同。求积公式1、精度较高公式（1）多项式。以点为零点的n次多项式： 上式称为多项式（2）求积公式。以多项式的n个实根为节点的插值求积公式为求积公式。考虑在上求积公式的构造1）一个节点2）两个节点二次正交多项式所以两点的求积公式为：对于一般区间的积分，可以用将区间转化为，即然后用相应的求积公式计算。（3）一般形式的求积公式为：其中是一个权重函数，为系数，为横坐标上的节点。因为，所以，一个n点的求积公式具有如下形式：其中，是函数在节点处的值，节点是正交多项式的根。给出x和A的表格：n正交多项式102213,42、在MATLAB中的实现MATLAB没有提供的有关计算函数，此处给出一部分的编程代码：function q=gaussL（f，a，b，x，A）N=length（x）；T=zeros（1，N）；T=（a+b）/2+((b-a)/2)*x;q=((b-a)/2)*sum(A.*feval(‘f’.T));其中，f为被积函数；x和A的值可有上表查到。数值微分数值微分的建立常用的三种思路：直接从微分的定义出发，通过近似的处理（泰勒展开），得到数值微分的近似公式；利用插值的基本思想，采用插值近似公式，对插值公式的近似求导得到原数值微分的近似公式根据已知数据，利用最小二乘拟合的方法，得到近似的函数，然后对此近似函数求微分就可以得到数值微分的近似公式。差分法近似微分1、计算公式在微积分中，一阶微分的计算可以在相邻点和间函数取得极限求得。所以给出下列差分近似式子：一阶向前差分：一阶向后差分：精