您将要找到将第一个轨道的第i个点连接到第二个轨道的第i个点的直线的方程式。然后，您可以使用i和z作为参数，在所有可能的值上变化以找到X和Y。
例：
import numpy as npimport matplotlibimport matplotlib.pyplot as pltfrom mpl_toolkits.mplot3d import Axes3DZ1 = 8.0Z2 = 9.0font = {'size' : 18}matplotlib.rc('font', **font)fig = plt.figure()ax = fig.add_subplot(111, projection='3d')t = np.linspace(0, 2 * np.pi, 100)x = np.cos(t)y = np.sin(2 * t)N = len(x)z = np.zeros(N)z[:,] = Z1t = np.linspace(0, 2 * np.pi, 100)x2 = 2 * np.cos(t)y2 = 2 * np.sin(2*t)N = len(x2)z2 = np.zeros(N)z2[:,] = Z2#Plot the first orbitax.plot(x, y, z, 'k-', linewidth=3.0)#Plot second orbitax.plot(x2, y2, z2, 'k-', linewidth=3.0)i, h = np.meshgrid(np.arange(len(x)), np.linspace(Z1, Z2, 10))X = (x2[i] - x[i]) / (Z2 - Z1) * (h - Z1) + x[i]Y = (y2[i] - y[i]) / (Z2 - Z1) * (h - Z1) + y[i]surf = ax.plot_surface(X, Y, h, color='m', alpha=0.3, linewidth=0)#Set axis and thingsax.set_xticks([1.0,1.5,2])ax.set_yticks([32,35,38])ax.set_ylabel('$||u||_{2}$', fontsize=26, rotation=0, labelpad = 26)ax.set_xlabel('$h$', fontsize=26)ax.set_zlabel('$mu$', fontsize=26, rotation=90)plt.tight_layout()plt.show()
