Matlab中几个数值积分函数的比较和优缺点.pdfVIP

Matlab 中几个数值积分函数的比较和优缺点 一、Z = trapz(X,Y,dim) 梯形数值积分，通过已知参数x,y 按dim 维使用梯形公式进行积分 。若不知道函数表达式，而是一组离散的数据 ，matlab 的函数库中只能使用trapz 来求积。梯形方法精度不高 ，可考虑自己动手编写高阶精度的方法。 例1 计算int(sin(x),0,pi) %by dynamic %all rights reserved by x=0:pi/100:2*pi; y=sin(x); z=trapz(x,y)%或者说使用z = pi/100*trapz(y) z = 1.0300e-017 z = pi/100*trapz(y) 二、[q,fcnt]= quad(fun,a,b,tol,trace,p1,p2...) 自适应simpson 公式数值积分，适用于精度要求低，被积函数平滑性 较差的数值积分 注意事项： 1.被积函数fun 必须是函数句柄 2.积分限[a,b]必须是有限的，因此不能为inf 3.p1 为其他需要传递的参数，一般是数值 可能警告： 1.Minimum step size reached 意味着子区间的长度与计算机舍入误差相当，无法继续计算了。 原因可能是有不可积的奇点 2.Maximum function count exceeded 意味着积分递归计算超过了10000 次。原因可能是有不可积的奇 点 3.Infinite or Not-a-Number function value encountered 意味着在积分计算时，区间内出现了浮点数溢出或者被零除。 例2 计算积分1/(x^3-2*x-p) ，其中参数p=5 ，积分区间为[0,2] %by dynamic %all rights reserved by F = @(x,n)1./(x.^3-2*x-n); Q = quad(@(x)F(x,5),0,2)%或者使用 quad(F,0,2,[],[],5) 效果是一样的，只是前者使用的函数嵌套 Q = -0.4605 quad(F,0,2,[],[],5) ans = -0.4605 三、[q,fcnt] = quadl(fun,a,b,tol,trace,p1,p2...) 自适应Lobatto 数值积分，适用于精度要求高，被积函数曲线比较平 滑的数值积分 注意事项： 同quad 可能警告： 同quad 例3 计算积分1/(x^3-2*x-p) ，其中参数p=5 ，积分区间为[0,2] %by dynamic %all rights reserved by F=@(x,p)1./(x.^3-2*x-p); Q = quadl(F,0,2,[],[],5)%或者Q = quadl(@(x)F(x,5),0,2) Q = -0.4605 四、[q,errbnd] = quadgk (fun,a,b,param1,val1,param2,val2,...) 自适应Gauss-Kronrod 数值积分，适用于高精度和震荡数值积分，支 持无穷区间，并且能够处理端点包含奇点的情况，同时还支持沿着不 连续函数积分，复数域线性路径的围道积分法 注意事项： 1.积分限[a,b]可以是[-inf,inf]，但必须快速衰减 2.被积函数在端点可以有奇点，如果区间内部有奇点，将以奇点 区间划分成多个，也就是说奇点只能出现在端点上 3.被积函数可以剧烈震荡 4.可以计算不连续积分，此时需要用到Waypoints参数，Waypoints 中的点必须严格单调 5.可以计算围道积分，此时需要用到Waypoints参数，并且为复数， 各点之间使用直线连接 6.param,val 为函数的其它控制参数，比如上面的waypoints就是，具 体看帮助 出现错误： 1.Reached the limit on the maximum number of intervals in use 2.Infinite or Not-a-Number function value encountered 例4 计算有奇点积分int(exp(x)*log(x),0,1) %by dynamic %all rights reserved by F=@(x)exp(x).*log(x);%奇点必须在端点上，否则请先进行区间划 分 Q = quadgk(F,0,1) Q = -1.3179 例5 计算半无限震荡积分int(x^5*exp(-x)*sin(