如何在Matlab中绘制三元函数？求助椭球面的正确绘图函数

首先你代码里有两个明显问题：

1. meshgrid(1:0.5:10,1:20,1:20)生成了三维网格数组，X、Y、Z都是3维的，计算出来的F自然也是3维数组，但surf()要求输入的高度值必须是二维数组，这就是报错的直接原因。
2. F的计算公式写错了，应该是sqrt(X.^2 + Y.^2 + Z.^2)，你写成了两个Y.^2，这会导致结果错误。

你要画的是隐式曲面（$F=x^2+y2+z^2$，本质是球面族，固定F值就是球面），surf()只适合显式曲面（比如$z=f(x,y)$），这里推荐两种可行的绘图方法：

方法一：用isosurface绘制隐式曲面 #

isosurface专门用来提取三维数据中等值面，适合画这类隐式定义的曲面，示例代码：

% 生成三维网格，调整范围更合理（比如-10到10）
[X,Y,Z] = meshgrid(-10:0.5:10, -10:0.5:10, -10:0.5:10);
% 定义隐式函数，比如画半径为10的球面：x²+y²+z²=100
V = X.^2 + Y.^2 + Z.^2 - 100;
% 提取等值面（V=0的面）
p = patch(isosurface(X,Y,Z,V,0));
% 设置曲面颜色和光照，让效果更清晰
set(p, 'FaceColor', 'blue', 'EdgeColor', 'none');
axis equal;
camlight; 
lighting gouraud;

方法二：参数化绘制球面 #

用球面的参数方程（θ∈[0,π]，φ∈[0,2π]）生成二维参数网格对应的三维坐标，这样就能用surf()正常绘制：

% 定义参数范围
theta = linspace(0, pi, 50);
phi = linspace(0, 2*pi, 50);
[Theta, Phi] = meshgrid(theta, phi);
% 球面参数方程，半径r=10
r = 10;
X = r.*sin(Theta).*cos(Phi);
Y = r.*sin(Theta).*sin(Phi);
Z = r.*cos(Theta);
% 用surf绘制
surf(X,Y,Z);
axis equal;
shading interp;

如果只是想画显式曲面（比如固定z，画x-y平面的等高线，或者把z表示成x,y的函数），那需要把meshgrid改成二维的，比如[X,Y] = meshgrid(1:0.5:10,1:20)，然后计算Z=f(X,Y)，再用surf(X,Y,Z)，但这不适用于你要的三维球面/椭球面。

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