如何用Matlab的trapz函数分别计算x轴上下方的曲线面积

要分开统计x轴上方（y≥0）和下方（y≤0）的区域面积，需要拆解trapz的梯形积分逻辑，先计算每个分段梯形的面积，再按正负分类求和，具体实现步骤如下：

实现思路 #

1. 计算相邻x值的间隔（每个梯形的底边长）
2. 计算每个梯形的平均高度（相邻y值的均值）
3. 得到每个分段梯形的面积
4. 分别累加正面积（x轴上方）和负面积（x轴下方）

代码示例 #

% 示例输入
x = [0,1,1,2];
y = [7,7,-4,-4];

% 计算每个分段的x间隔
dx = diff(x);
% 计算每个分段的平均y值（梯形的高）
y_avg = (y(1:end-1) + y(2:end)) ./ 2;
% 计算每个梯形的面积
segment_areas = dx .* y_avg;

% 分别统计上下方面积
area_above = sum(segment_areas(segment_areas > 0));
area_below = sum(segment_areas(segment_areas < 0));

disp(['x轴上方面积: ', num2str(area_above)]);
disp(['x轴下方面积: ', num2str(area_below)]);

运行上述代码会输出：

x轴上方面积: 7
x轴下方面积: -4

trapz函数核心说明（中文翻译） #

trapz 基于梯形积分法计算数值积分：

• 输入 trapz(x,y) 时，函数将每对相邻的(x,y)点构成的区域视为梯形，累加所有梯形面积得到总积分。
• 若仅传入 trapz(y)，默认x为步长为1的等间隔序列。
• 处理矩阵时，函数会按列独立计算每列的积分值。

内容的提问来源于stack exchange，提问作者user2587726