是的,所以我有一个x值,y值和z值的列表(我认为我转换成数组?)我想制作一个表面图,但它不起作用。
以下是我要做的事情,您可以忽略大部分代码,因为它非常无关紧要 - 只需看看我的xdis
,ydis
和dist
我正在尝试绘制atm我正在获得ValueError:需要超过1个值来解压:(。非常感谢。
from math import *
from numpy import *
import pylab
def sweep (v,p,q,r,s):
a=.98
for i in range (1, len(v)-1):
for j in range (1, len(v)-1):
c =0.0
if i==p and j==q: c =1.0
if i==r and j==s: c= -1.0
v[i,j]=(v[i -1,j]+v[i+1,j]+v[i,j -1]+v[i,j+1]+c-a*v[i,j])/(4-a)
def main():
from mpl_toolkits.mplot3d import Axes3D
from matplotlib import cm
from matplotlib.ticker import LinearLocator, FormatStrFormatter
import matplotlib.pyplot as plt
ydis=[]
xdis=[]
resis=[]
for j in range(2,18):
for i in range(2,18):
v= zeros ((20,20),float )
p=q=9
r=i
s=j
dv =1.0e10
lastdv =0
count =0
while (fabs(dv - lastdv)>1.0e-7*fabs(dv)):
lastdv =dv
sweep(v,p,q,r,s)
dv=v[p,q]-v[r,s]
resis.append(dv)
xdis.append(r-p)
ydis.append(s-q)
X=array(xdis)
Y=array(ydis)
Z=array(resis)
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.plot_wireframe(X,Y,Z)
plt.show()
main()
答案 0 :(得分:3)
plot_wireframe期望三个2D阵列(X,Y,Z)作为输入。所以,
后:
X=np.array(xdis)
Y=np.array(ydis)
Z=np.array(resis)
添加:
X=X.reshape((-1,16))
Y=Y.reshape((-1,16))
Z=Z.reshape((-1,16))
答案 1 :(得分:0)
似乎“扫描”功能不会修改“v”,因此您将获得一个空列表。