matlab梯形法求定积分完整代码

这是一个使用梯形法求定积分的 MATLAB 代码示例：

function I = trapezoidal_rule(f, a, b, n)
% 梯形法求定积分
% 输入参数：
%   f: 被积函数，必须是一个函数句柄
%   a: 积分下限
%   b: 积分上限
%   n: 将区间 [a, b] 分成 n 个小区间
% 输出参数：
%   I: 积分近似值

h = (b - a) / n;  % 计算小区间长度
x = a:h:b;  % 计算各个小区间的右端点
y = f(x);  % 计算各个小区间右端点的函数值
I = h * (sum(y) - 0.5 * (y(1) + y(n+1)));  % 计算积分值
end

你可以将被积函数 `f` 以函数句柄的形式传递给 `trapezoidal_rule` 函数。例如，如果你要求 $f(x) = x^2$ 在区间 $[0,1]$ 上的定积分，可以这样调用：

f = @(x) x.^2;  % 定义被积函数
a = 0; b = 1; n = 100;  % 积分区间和分割数
I = trapezoidal_rule(f, a, b, n);  % 使用梯形法求积分
fprintf('积分值近似为 %.4f\n', I);  % 显示积分结果

希望这个代码能对你有所帮助！