如何在MATLAB中绘制分段函数f(x,y)与垂直平面的3D组合图

没问题，这事儿在MATLAB里很好实现！我给你一步步拆解，保证你能画出想要的图：

1. 生成网格数据 #

首先得创建x和y的网格矩阵，这是MATLAB 3D绘图的基础。你可以根据自己的需求调整x、y的范围和精度：

% 定义x和y的取值范围，这里以[-2, 2]为例，100个采样点保证平滑
x = linspace(-2, 2, 100);
y = linspace(-2, 2, 100);
[X, Y] = meshgrid(x, y); % 创建网格矩阵

2. 计算分段函数的Z值 #

用逻辑索引来给不同区间的Z赋值，这样比循环高效得多。这里假设你已经定义了h(x,y)、g(x,y)、k(x,y)三个函数（注意要支持向量化运算，用点运算符.*、./、.^等）：

% 初始化Z矩阵，和X/Y同尺寸
Z = zeros(size(X));

% 第一个区间：x - y < -0.5
mask1 = X - Y < -0.5;
Z(mask1) = h(X(mask1), Y(mask1)); % 替换成你的h函数，比如h=@(x,y) x.^2 + y.^2;

% 第二个区间：-0.5 ≤ x - y < 0
mask2 = (X - Y >= -0.5) & (X - Y < 0);
Z(mask2) = g(X(mask2), Y(mask2)); % 替换成你的g函数，比如g=@(x,y) x.*y;

% 第三个区间：x - y ≥ 0
mask3 = X - Y >= 0;
Z(mask3) = k(X(mask3), Y(mask3)); % 替换成你的k函数，比如k=@(x,y) sin(x) + cos(y);

3. 绘制分段函数的3D曲面 #

用surf或mesh绘制函数曲面，记得打开hold on来叠加后续的垂直平面：

figure;
surf(X, Y, Z, 'EdgeColor', 'none'); % 关闭边缘线让曲面更干净
hold on; % 保持当前图窗，方便叠加其他图形

4. 绘制两个垂直平面 #

垂直平面x-y=-0.5和x-y=0可以通过构造平面的网格来绘制，设置半透明效果能让平面不遮挡后面的函数图像：

% 先获取Z的取值范围，让平面覆盖函数的高度
z_min = min(Z(:));
z_max = max(Z(:));
x_plane = linspace(-2, 2, 100);
z_plane = linspace(z_min, z_max, 100);

% 绘制平面x - y = -0.5（即y = x + 0.5）
[X_plane1, Z_plane1] = meshgrid(x_plane, z_plane);
Y_plane1 = X_plane1 + 0.5;
mesh(X_plane1, Y_plane1, Z_plane1, 'EdgeColor', 'r', 'FaceAlpha', 0.3); % 红色半透明

% 绘制平面x - y = 0（即y = x）
[X_plane2, Z_plane2] = meshgrid(x_plane, z_plane);
Y_plane2 = X_plane2;
mesh(X_plane2, Y_plane2, Z_plane2, 'EdgeColor', 'b', 'FaceAlpha', 0.3); % 蓝色半透明

5. 美化图形（可选但推荐） #

添加标签、图例、网格，让图的可读性更强：

xlabel('X');
ylabel('Y');
zlabel('Z');
title('Piecewise Function f(x,y) with Vertical Planes');
legend('f(x,y)', 'x-y=-0.5', 'x-y=0', 'Location', 'best');
grid on;
hold off;

小提示 #

• 如果你的h/g/k函数不支持向量化，可以用arrayfun来计算，但效率会低一些，比如Z(mask1) = arrayfun(@(x,y) h(x,y), X(mask1), Y(mask1));
• 可以调整linspace的采样点数，点数越多图形越平滑，但计算时间也会增加
• 若想更换绘图风格，比如用mesh代替surf，或者调整曲面颜色，直接修改对应的绘图函数参数即可

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