复合梯形公式matlab

复合梯形公式是一种数值积分方法，用于计算函数在一定区间上的定积分。其 matlab 代码实现如下：

function [I] = composite_trapezoidal(f, a, b, n)
% 输入：f-被积函数，a,b-积分区间，n-子区间数目
% 输出：I-积分近似值

h = (b - a) / n; % 计算子区间宽度
x = linspace(a, b, n+1); % 划分子区间
y = f(x); % 计算函数值
I = h * (sum(y) - (y(1) + y(n+1)) / 2); % 计算积分近似值
end

其中，`f`是被积函数，`a`和`b`是积分区间的上下界，`n`是子区间数目，`h`是子区间宽度，`x`是子区间的端点，`y`是对应的函数值，`I`是积分的近似值。