Octave绘制r=sin²θ旋转曲面异常：期望甜甜圈却出现尖刺

核心问题 #

你的脚本存在两个关键错误，导致生成的曲面出现尖刺：

1. 未定义变量R #

脚本中计算3D坐标X和Y时直接使用了R，但这个变量从未定义过。虽然你之前定义了二维曲线的半径f = sin(theta).^2，但没有将其扩展为网格形式以匹配meshgrid生成的U和T矩阵。

2. Z坐标参数化错误 #

你注释的Z坐标定义Z = repmat(theta', 1, length(u))是问题的核心：

• 原二维曲线r=sin²θ在theta=±π时，r=0，对应二维点为原点；
• 当你把Z设为theta时，这些r=0的点会映射到Z=±π处的(0,0,±π)，形成尖刺；
• 同时，这种参数化方式本身就不是甜甜圈（环面）的正确参数化逻辑。

修正后的脚本 #

如果想要生成无尖刺的曲面，可选择以下两种方案：

方案1：修正现有脚本（基于原曲线消除尖刺） #

theta=linspace(-pi,pi,200);
f=sin(theta).^2;
x=f.*cos(theta);
y=f.*sin(theta);

figure;
plot(x,y);
axis equal;

% 生成旋转曲面的点
u = linspace(-pi, pi, 200);     % 旋转角度
[U, T] = meshgrid(u, theta);    % 创建网格坐标

% 正确定义半径网格
R = sin(T).^2;  

% 绕x轴旋转的参数化（避免原曲线零点集中在Z轴）
X = repmat(x, 1, length(u));    % X保持原曲线的x值
Y = R .* sin(U);                % Y为旋转后的y分量
Z = R .* cos(U);                % Z为旋转后的z分量

% 绘制曲面
figure;
surf(X, Y, Z);
xlabel('X'); ylabel('Y'); zlabel('Z');
title('Surface de révolution: r = sin^2(theta)');
shading interp;
axis equal;

方案2：生成标准甜甜圈（环面） #

如果你确实想要标准的无尖刺甜甜圈，需要使用无零点的母线，比如r = 1 + 0.5*cos(theta)（外半径1.5，内半径0.5）：

theta=linspace(-pi,pi,200);
u = linspace(-pi, pi, 200);
[U, T] = meshgrid(u, theta);

% 标准环面参数化
R = 1 + 0.5*cos(T);  % 母线半径随theta变化
X = R .* cos(U);
Y = R .* sin(U);
Z = 0.5*sin(T);

figure;
surf(X, Y, Z);
xlabel('X'); ylabel('Y'); zlabel('Z');
title('标准甜甜圈曲面');
shading interp;
axis equal;

关键说明 #

• 原曲线r=sin²θ本身在theta=±π处有零点，若直接绕z轴旋转并将Z设为theta，必然会产生尖刺；
• 标准甜甜圈的母线是一个不经过原点的圆，因此旋转后不会出现奇点。

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