如何用Octave绘制两个梯形模糊隶属度函数的复合曲面？

你的问题核心是数据格式不匹配：cartprod生成的是一维点序列（N×2矩阵），但绘制曲面需要二维网格形式的X、Y、Z矩阵——这就是mesh无法正常工作、plot3只画出线条的原因。

修正步骤与代码 #

用meshgrid生成二维网格，配合向量化计算就能得到正确的曲面效果：

% MAX复合梯形模糊隶属度曲面绘制
pkg load linear-algebra

% 生成二维网格（替代cartprod，满足曲面绘图的数据格式要求）
tx = linspace(-10, 10, 100);
ty = linspace(-10, 10, 100);
[X, Y] = meshgrid(tx, ty);

% 梯形隶属度函数参数
a = -6; b = -2; c = 2; d = 6;

% 向量化计算X方向的梯形隶属度（保持矩阵维度）
trape_x = max(min(min(ones(size(X)), (X - a)/(b - a)), -(X - d)/(d - c)), 0);
% 向量化计算Y方向的梯形隶属度
trape_y = max(min(min(ones(size(Y)), (Y - a)/(b - a)), -(Y - d)/(d - c)), 0);

% 计算MAX复合（逐元素取最大值，对应模糊集合的MAX运算）
Z = max(trape_x, trape_y);

% 绘制曲面（mesh展示网格，surf带颜色填充，按需选择）
mesh(X, Y, Z);
xlabel('X'); ylabel('Y'); zlabel('Membership Degree');
title('MAX Composite Trapezoidal MF');
grid on;

关键说明 #

• meshgrid(tx, ty)生成的X和Y都是100×100的矩阵，每个位置对应平面上的一个网格点，后续计算的隶属度和复合结果也保持同维度，完全适配曲面绘图函数的要求。
• 向量化计算既避免了循环，又保证了矩阵维度一致，max/min默认逐元素操作，正好对应模糊集合的MAX复合逻辑。
• 用mesh或surf替代plot3：plot3是绘制三维线条，而mesh/surf是绘制曲面网格/填充曲面，能得到你期望的类似uLT9T.png的效果。

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