如何在Matlab中绘制垂直于已知直线且过定点的平面？

已知一条由蓝色端点与黑色星型端点构成的绿色直线，需绘制一个垂直于该直线且经过红色定点的平面。目前已算出平面的系数A、B、C、D（满足平面方程 (Ax + By + Cz + D = 0)），以下是完成绘制的具体方法：

现有代码 #

plot3(x1,y1,z1,'.b','MarkerSize',10) % 蓝色端点
hold on
plot3(x2,y2,z2,'kp','MarkerSize',6,'MarkerFaceColor','k') % 黑色星型端点
plot3(point_x,point_y,point_z,'.r','MarkerSize',14) % 红色定点
plot3([x1 x2],[y1 y2],[z1 z2],'g') % 绿色直线

A = x1 - x2;
B = y1 - y2;
C = z1 - z2;

D = -A*point_x - B*point_y - C*point_z;

绘制平面的步骤与代码 #

要绘制平面，需先生成平面上的网格点，再用surf或mesh函数完成绘制：

1. 确定绘制范围：基于现有点的坐标范围设置，确保平面完整展示且不超出视图
2. 生成网格坐标：用meshgrid创建x和y的网格，通过平面方程解出对应z值；若C=0（平面平行于z轴），则换为x和z网格求解y值
3. 绘制平面：使用surf函数并设置透明度，避免遮挡原有图形

完整代码示例 #

plot3(x1,y1,z1,'.b','MarkerSize',10) % 蓝色端点
hold on
plot3(x2,y2,z2,'kp','MarkerSize',6,'MarkerFaceColor','k') % 黑色星型端点
plot3(point_x,point_y,point_z,'.r','MarkerSize',14) % 红色定点
plot3([x1 x2],[y1 y2],[z1 z2],'g') % 绿色直线

A = x1 - x2;
B = y1 - y2;
C = z1 - z2;

D = -A*point_x - B*point_y - C*point_z;

% 确定平面绘制范围，可根据实际情况调整
x_range = linspace(min(x1,x2,point_x)-2, max(x1,x2,point_x)+2, 50);
y_range = linspace(min(y1,y2,point_y)-2, max(y1,y2,point_y)+2, 50);
[X,Y] = meshgrid(x_range, y_range);

% 根据平面方程计算对应坐标，避免除以0的情况
if abs(C) > 1e-6
    Z = (-A*X - B*Y - D)/C;
else
    Z_range = linspace(min(z1,z2,point_z)-2, max(z1,z2,point_z)+2, 50);
    [X,Z] = meshgrid(x_range, Z_range);
    Y = (-A*X - C*Z - D)/B;
end

% 绘制半透明平面
surf(X,Y,Z,'FaceColor','cyan','Alpha',0.3,'EdgeColor','none');
hold off
axis equal;
xlabel('X'); ylabel('Y'); zlabel('Z');
grid on;

说明 #

• 调整linspace的范围和点数可改变平面的大小与精度
• Alpha参数设为0.3是为了让平面半透明，清晰展示后方的点与直线
• 加入了C=0的特殊情况处理，避免运行报错

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