中国传媒大学《MATLAB基础与应用》课件-第8章数值积分与数值微分.pptVIP

在相同的n（一定的采样条件）下取得更高的精度！！ 在相同的n时收敛的更快！！！ 利用低阶公式产生高精度的结果。 Romberg 序列 Newton-Cotes 类似的有： 龙贝格（Romberg）求积 基本思想： 梯形公式的逐次分半过程进行加速。 Richardson外推公式 ? Romberg 算法： ? ? ? ? ? ? … … … … … … ? T1 = ) 0 ( 0 T ? T8 = ) 3 ( 0 T ? T4 = ) 2 ( 0 T ? T2 = ) 1 ( 0 T ? S1 = ) 0 ( 1 T ? R1 = ) 0 ( 3 T ? S2 = ) 1 ( 1 T ? C1 = ) 0 ( 2 T ? C2 = ) 1 ( 2 T ? S4 = ) 2 ( 1 T ? ? 折半 折半 折半 直到 又叫：数值积分逐次分半加速收敛法 8.4 Gauss求积方法 将节点 x0 … xn 以及系数 A0 … An 都作为待定系数。令 f (x) = 1, x, x2, …, x2n+1 代入可求解，得到的公式具有2n+1 次代数精度。这样的节点称为Gauss 点，公式称为Gauss 型求积公式。 为使问题更具一般性,我们研究带权积分 ? (x)为权函数, Ak(k＝0,1,…,n)为不依赖于f (x)的求积系数, xk (k＝0,1,…,n)为求积节点. 要使(8.1)具有2n+1次代数精度，则需要满足 构造具有2n+1次代数精度的求积公式 (8.1) 8.5 特殊函数的积分 一、振荡函数的积分 工程问题中有时要计算如下形式的积分： 其中 。当w很大时，coswx与sinwx在区间(a,b)内与x轴有很多个交点，称其为振荡函数。相应地，当w很大时， f(x)coswx与f(x)sinwx在区间(a,b)内与x轴也有很多个交点。称上述积分为振荡函数的积分。 二、反常（广义）积分 1、无界函数的反常积分 设函数 在区间 上连续，b为奇点。若对 且 ，称极限 为无界函数 在 上的反常积分（或暇积分），记为 。 2、无穷区间上的反常积分 设对任何大于a的实数b，f (x)在[a,b)上均可积，则称极限 为f (x)在无穷区间 上的反常积分，记为 。 对于无穷区间上的反常积分通常有两种解决方法： 无穷区间逼近 变量替换 在某些情况下，可以通过变量替换将无穷区间的积分变成有限区间的积分。 做等距节点，x轴，y轴分别有： 先计算 ，将x作为常数，有 再将y作为常数，在x方向，计算上式的每一项的积分 三、重积分的计算——以二重积分的复化梯形公式为例说明。 系数：在积分区域的四个角点为1/4，4个边界为1/2，内部节点为1 8.6 数值积分的MATLAB函数求解 1、trapz()函数 MATLAB中的trapz()函数是基于复化梯形公式设计编写的，其一般调用格式为： I=trpaz(x,y,dim) 其中x,y是观测数据，x可以为行向量或列向量，y可以为向量或矩阵，y的行数应等于x向量的元素个数；dim表示按维进行求积，若dim=1（缺省值），则按行求积，若dim=2，则按列求积。 2、quad()函数 MATLAB提供的quad()函数是基于自适应辛普森法设计的，该函数的调用格式为： [q,fcnt] = quad(fun,a,b,tol,trace,p1,p2,...) 其中fun是被积函数，可以是字符表达式、内联函数、匿名函数和M函数；a,b是定积分的上限和下限；tol为指定的误差限，缺省值为；trace提供中间输出[fcnt a b-a q]，若trace=[]，则quad不提供中间输出；p1,p2,...是函数fun的附加参数。q是返回的数值积分；fcnt返回函数评估的次数。 另外，MATLAB还提供了一个新的函数quadl()。其调用格式与quad()函数完全一致，使用的算法是自适应Lobatto算法，其精度和速度均远高于quad()函数，所以在追求高精度数值解时建议使用该函数。 3、quadgk()函数 quadgk()函数是MATLAB R2007b版本提供的基于Gauss-Kronrod算法实现的数值积分函数，该函数可以用来求解振荡函数的积分、广义积分甚至是复数积分，其一般调用格式为： [q,errbnd] = quadgk(fun,a,b,param1,val1,param2,val2,...) 其中fun是被积函数，可以是字符表达式、内联函数、匿名函数和M函数；a,b是积分的上限和下限，它们可