我有两个值为x和y的1D阵列,还有一个值为z的2D阵列,每个点的列对应于x值,行对应于y值。有没有办法用这些数据得到一个plot_surface?当我尝试这样做时,它不会给我任何回报。这是代码:(calculate_R是我为程序制作的函数)
x=np.arange(0,10,1)
y=np.arange(0,1,0.2)
lx= len(x)
ly=len(y)
z=np.zeros((lx,ly))
for i in range(lx):
for j in range(ly):
z[i,j]=calculate_R(y[j],x[i])
fig = plt.figure()
ax = Axes3D(fig)
x, y = np.meshgrid(x, y)
ax.plot_surface(x, y, z, rstride=1, cstride=1, cmap='hot')
您忘记调用plt.show()
来显示您的绘图。
请注意,您可能能够利用numpy矢量化来加快z
:的计算
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d.axes3d import Axes3D
x = np.arange(0,10,1)
y = np.arange(0,1,0.2)
xs, ys = np.meshgrid(x, y)
# z = calculate_R(xs, ys)
zs = xs**2 + ys**2
fig = plt.figure()
ax = Axes3D(fig)
ax.plot_surface(xs, ys, zs, rstride=1, cstride=1, cmap='hot')
plt.show()
在这里,我使用了一个简单的函数,因为您没有提供一个完整的工作示例。