MATLAB绘制Maple同款3D曲面边缘异常的解决方法咨询

你的代码核心问题出在径向网格R的生成逻辑上：当你写r=0:0.005:rho时，rho是一个与theta长度对应的数组（每个角度对应一个超椭圆半径），但MATLAB的冒号运算符只会取rho的第一个元素作为径向范围的上限，这就导致所有角度切片的径向范围都被限制在第一个角度对应的半径，自然无法正确显示完整的超椭圆边缘。

下面是修正后的完整代码，我会标注关键修改点：

clc;clear all;close all;
theta = linspace(0,2*pi,201);
rho = (cos(theta).^4 + sin(theta).^4).^(-1/4); % 把原代码里的0.1e1/0.4e1换成更易读的1/4

% 关键修正：生成与每个角度匹配的径向网格
num_r = 201; % 径向采样点数，可根据需求调整
[TH,R_normalized] = meshgrid(theta, linspace(0,1,num_r)); % 先生成0到1的归一化径向值
R = R_normalized .* rho; % 将归一化值与对应角度的超椭圆半径相乘，确保每个theta列的R都能到达正确边缘

X = R .* cos(TH);
Y = R .* sin(TH);
U = -R.^2/4 + 0.2866749079 ...
    - 2.78e-11.*R.^40.*cos(40.*TH) ...
    - 8.064151950e-11.*R.^36.*cos(36.*TH) ...
    - 4.943347436e-10.*R.^32.*cos(32.*TH) ...
    - 4.390774680e-9.*R.^28.*cos(28.*TH) ...
    - 4.536478572e-8.*R.^24.*cos(24.*TH) ...
    - 4.832799212e-7.*R.^20.*cos(20.*TH) ...
    - 5.529796642e-6.*R.^16.*cos(16.*TH) ...
    - 7.085150960e-5.*R.^12.*cos(12.*TH) ...
    - 1.134817408e-3.*R.^8.*cos(8.*TH) ...
    - 3.546317566e-2.*R.^4.*cos(4.*TH);

mesh(X,Y,U);
axis equal; % 开启等比例坐标轴，让显示效果更接近Maple的输出

关键修改说明： #

• 不再直接尝试为每个角度生成不同的径向范围（这在MATLAB的冒号语法里无法实现），而是先生成0到1的归一化径向值，再将其与对应角度的超椭圆半径相乘，这样每个角度切片的径向范围都会精确到达对应的边缘。
• 将原代码中晦涩的科学计数法（比如0.4e1）替换为更直观的数值（比如4），方便后续调试和维护。
• 添加axis equal命令，确保坐标轴比例一致，让MATLAB的绘图效果和Maple的结果视觉上更匹配。

修改完成后，MATLAB生成的3D曲面就会和Maple的结果一致，超椭圆边缘也能正确显示了。

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