Octave二次曲面绘图问题求助：代码错误排查

我来帮你排查并解决这个绘图问题～首先看你的代码，主要有两个核心问题导致绘图异常：

1. 平方根内的负数导致NaN值 #

你的表达式 zz=sqrt(-9*xx.^2+9*yy.^2) 可以简化为 zz=sqrt(9*(yy.^2 - xx.^2))，当 |x| > |y| 时，括号内的结果为负数，sqrt 会返回 NaN（非数值），这部分区域就不会被绘制出来，导致曲面残缺。

2. 采样点过少导致曲面粗糙 #

你用 x = [-3:1:3] 只生成了7个点，网格密度太低，画出来的曲面会非常生硬，完全达不到优于在线工具的效果。

修改后的完整代码 #

clear;
% 生成更密集的采样点，用linspace让分布更均匀
x = linspace(-3, 3, 50);
y = linspace(-3, 3, 50);
[xx, yy] = meshgrid(x, y);

% 计算z的平方项，避免直接在sqrt里处理负数
z_squared = 9 * (yy.^2 - xx.^2);

% 分别计算正、负z值，完整展示双曲柱面的上下两部分
zz_pos = sqrt(max(z_squared, 0));  % 负数部分取0后开方，避免NaN
zz_neg = -sqrt(max(z_squared, 0));

figure;
hold on;  % 保持画布，同时绘制两个部分
mesh(xx, yy, zz_pos);
mesh(xx, yy, zz_neg);
hold off;

% 添加标签和标题，提升图表可读性
xlabel('x轴');
ylabel('y轴');
zlabel('z轴');
title('双曲柱面 z² = 9(y² - x²)');

代码说明 #

• 用 linspace 生成50个采样点，让曲面更平滑细腻；
• 先计算 z_squared，再用 max(z_squared, 0) 处理负数情况，避免sqrt返回NaN；
• 同时绘制正、负z值的曲面部分，完整展示双曲柱面的形态；
• 添加坐标轴标签和标题，让图表更专业。

如果你的目标是绘制其他类型的二次曲面（比如单叶双曲面、椭圆抛物面等），可以调整 z_squared 的表达式，比如单叶双曲面的公式是 x²/9 + y²/4 - z²/1 = 1，只需要对应修改计算逻辑即可。

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