用Matplotlib绘制椭球体

Question

有人有绘制椭球体的样例代码吗？球体有一个
在“matplotlib”网站上，但椭球体没有。我正试图策划

x**2 + 2*y**2 + 2*z**2 = c

其中“c”是定义椭球体的常量（如10）。我试过了
meshgrid（x，y）route，修改了公式，使’z’在一边，但是
`sqrt是个问题。“matplotlib”球体示例适用于角度“u，v”，
但我不知道如何计算椭球体。

Answer 1

下面是如何通过球坐标实现的：

from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
import numpy as np

fig = plt.figure(figsize=plt.figaspect(1))  # Square figure
ax = fig.add_subplot(111, projection='3d')

coefs = (1, 2, 2)  # Coefficients in a0/c x**2 + a1/c y**2 + a2/c z**2 = 1 
# Radii corresponding to the coefficients:
rx, ry, rz = 1/np.sqrt(coefs)

# Set of all spherical angles:
u = np.linspace(0, 2 * np.pi, 100)
v = np.linspace(0, np.pi, 100)

# Cartesian coordinates that correspond to the spherical angles:
# (this is the equation of an ellipsoid):
x = rx * np.outer(np.cos(u), np.sin(v))
y = ry * np.outer(np.sin(u), np.sin(v))
z = rz * np.outer(np.ones_like(u), np.cos(v))

# Plot:
ax.plot_surface(x, y, z,  rstride=4, cstride=4, color='b')

# Adjustment of the axes, so that they all have the same span:
max_radius = max(rx, ry, rz)
for axis in 'xyz':
    getattr(ax, 'set_{}lim'.format(axis))((-max_radius, max_radius))

plt.show()

上面的程序实际上生成了一个更好看的“正方形”图形。这个解决方案的灵感来自
示例
在Matplotlib的画廊.