matlab数值计算的优点,Matlab中几个数值积分函数的比较和优缺点

具体参数的意义我这里不详细说明，大家可以看帮助系统，我们这里只是讨论各大函数的区别比较和注意

一、Z = trapz(X,Y,dim)

梯形数值积分，通过已知参数x,y按dim维使用梯形公式进行积分

例1 计算int(sin(x),0,pi)

>>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'

意味着在积分计算时࿰