matlab 符号 积分法,MATLAB积分方法

一 相关函数

%符号积分

int(f,v)

int(f,v,a,b)

%数值积分

trapz(x,y)%梯形法沿列方向求函数Y关于自变量X的积分

cumtrapz(x,y)%梯形法沿列方向求函数Y关于自变量X的累计积分

quad(fun,a,b,tol)%采用递推自适应Simpson法计算积分

quad1(fun,a,b,tol)%采用递推自适应Lobatto法求数值积分

dbquad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol)%二重(闭型)数值积分指令

triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol)%三重(闭型)数值积分指令

二.示例:

例1:计算f(t)=exp(-t^2)在[0,1]上的定积分

本例演示:计算定积分常用方法

>>symsx

int(exp(-x^2),0,1)

ans=

1/2*erf(1)*pi^(1/2)　%erf为误差函数

>>vpa(int(exp(-x^2),0,1))

ans=

.7468241328124270

>>d=0.001;x=0:d:1;d*trapz(exp(-x.^2))

ans=

0.7468

>>quad('exp(-x.^2)',0,1,1e-8)

ans=

0.7468

例2:计算f(t)=1/log(t)在[0,x],0

注意:被积函数于x=0无义,在x-->1^-处为负无穷

本例演示:用特殊函数表示的