Matlab 绘制椭球体的横截面等高线图（如corel图所示）

这是我使用

slice

命令获得的图像。无法获取更精细的网格网格：

此图分别显示了沿X-Y、Y-Z和X-Z平面围绕椭球体的压力等值线：

我用MATLAB绘制了这些图表。我想在一个椭球体的三维图中显示所有这些平面。可以使用MATLAB或任何其他平台完成吗

% 3D meshgrid in clast's coordinate
    xgv     = 0;           % grid vector: x'axis,a1
    ygv     = -8:0.1:8;           % grid vector: y'axis,a2
    zgv     = -6:0.1:6;    % grid vector: z'axis,a3
    [X,Y,Z] = meshgrid(xgv,ygv,zgv);



    % Exterior points
    ind   = (X./1).^2 + (Y./a(2)).^2 + (Z./a(3)).^2 > 1;
    x_ex  = X(ind);
    y_ex  = Y(ind);
    z_ex  = Z(ind);
    ep    = cat(1,x_ex',y_ex',z_ex');
    % total points number
    num   = numel(X);
    % exterior points number
    [~,n] = size(ep);


    % Exterior Fields  
    % 
      G = Ex_Gtensor(a,ep); % this is some function to calculate pressure at 
               external points
      Here I use some algorithm to calculate pressure field outside ellipsoid.
    %--------------------------------------------------
    pressure        = zeros(num,1);
    pressure(~ind)  = p1;
    pressure(ind)   = p2;
    % So this pressure variable has all the pressure values for all external points 
    pressure        = reshape(pressure,size(squeeze(Z)));


    contourf(squeeze(Y),squeeze(Z),pressure)
    hold on
    axis equal
    % Similarly I use this code for X-Y and X-Z plane.

如果需要任何其他信息，请告诉我。

---

提前谢谢

你知道吗？嗨@EBH！谢谢你的回复！我试着使用切片并得到一些图。但是，我无法得到更细网格的结果。这表明，对于更细网格，压力计算的结果是复数格式的。你有什么建议吗？（我已经编辑了问题并添加了切片结果。）你知道吗？嗨@EBH！谢谢你的回复！我试着使用切片并得到一些图。但是，我无法得到更细网格的结果。这表明，对于更细网格，压力计算的结果是复数格式的。你有什么建议吗？（我已经编辑了问题并添加了切片结果。）