复合梯形公式(compound trapezoidal formula)是1993年公布的数学名词。
将积分区间分为若干份, 在每一个“小区间”上用低阶梯形求积公式可得
复合梯形公式的收敛阶为2阶。
matlab程序
function I = ftrapz(fun,a,b,n)%fun,a,b,n分别为被积分函数、积分下限、积分上限、积分区间数目h = (b-a)/n;%求积分区间长度x = linspace(a,b,n+1);%bailinspace(x1,x2,N)用于产生x1,x2之间的N点行线性的矢量;%其中x1、x2、N分别为起始值、终止值、元素个数。%其中x1、x2、N分别为起始值、终止值、元素个数。y = feval(fun,x);I = h * (0.5*y(1) + sum(y(2:n)) + 0.5*y(n+1) );
被积函数f(x) Matlab表示
f=@(x)(sin(1./x)); %@表示参数;%(x)表示一个参数;%(sin(1./x))表示被积分函数;
求解
f=@(x)(sin(1./x)); m=ftrapz(f,1,1.5,5);format