梯形法数值积分—— trapz

森之张卫东

函数2 trapz

功能 梯形法数值积分

格式 T = trapz(Y)       %用等距梯形法近似计算Y的积分。若Y是一向量，则trapz(Y)为Y的积分；
若Y是一矩阵，则trapz(Y)为Y的每一列的积分；若Y是一多维阵列，
*则trapz(Y)沿着Y的第一个非单元集的方向进行计算。

T = trapz(X,Y)   %用梯形法计算Y在X点上的积分。若X为一列向量，Y为矩阵，
且size(Y,1) = length(X)，则trapz(X,Y)通过Y的第一个非单元集方向进行计算。

T =trapz(…,dim)   %沿着dim指定的方向对Y进行积分。若参量中包含X，
则应有length(X)=size(Y,dim)。

例2-41

>>X = -1:.1:1;

>>Y = 1./(1+25*X.^2);

>>T = trapz(X,Y)

计算结果为：

T =

    0.5492