trapz数组

使用方法如下：

I=trapz(x,y)

其中x和y分别是自变量和对应的值，例如有函数y=x^3-2x-3，为了计算在[0,1]上的积分，可以这么做：

>> format compact

>> x=0:0.05:1;

>> y=x.^3-2.*x-3;

>> I=trapz(x,y)

I =

-3.7494

这个函数是可以直接使用经典积分理论计算的，精确值为 -15/4=-3.75，误差为0.016%。

扩展资料：

注意事项

MATLAB中的trapz()函数是基于复化梯形公式设计编写的，其一般调用格式为：

I=trpaz(x,y,dim)

其中x,y是观测数据，x可以为行向量或列向量，y可以为向量或矩阵，y的行数应等于x向量的元素个数；dim表示按维进行求积，若dim=1（缺省值），则按行求积，若dim=2，则按列求积。

如：计算函数y=x^3-2x-3，为了计算在[0,1]上的积分

x=0:0.05:1;

y=x.^3-2.*x-3;

trapz(x,y)

ans =

-3.7494