复合梯形公式代码matlab

复合梯形公式是一种数值积分方法，用于计算定积分。下面是使用 MATLAB 实现复合梯形公式的代码：

function I = composite_trapezoidal(f, a, b, n)
% f: 被积函数
% a, b: 积分区间 [a, b]
% n: 分割数

h = (b - a) / n; % 计算步长
x = a:h:b; % 等距节点
y = f(x); % 计算节点处函数值
I = h * (sum(y) - (y(1) + y(end)) / 2); % 计算积分值
end

使用时，可以先定义被积函数 `f`，然后调用 `composite_trapezoidal(f, a, b, n)` 函数计算积分值，其中 `a` 和 `b` 分别为积分区间的左右端点，`n` 为分割数。例如，要计算 $f(x) = x^2$ 在区间 $[0,1]$ 上的积分，可以这样调用函数：

f = @(x) x.^2; % 定义被积函数
I = composite_trapezoidal(f, 0, 1, 100); % 计算积分值
disp(I); % 输出积分值

这里使用了匿名函数 `@(x) x.^2` 来定义被积函数，`.^` 表示对向量的每个元素进行平方操作。`composite_trapezoidal` 函数的输出结果 `I` 即为所求积分值。