我在调用一个函数两次时遇到了问题。如果我注释最后3行并保留show()
,我不会收到任何错误,事情就会按预期进行。但是,如果我不注释掉它们,再次调用最后一个函数会产生以下错误:
Traceback (most recent call last):
File "second_1.py", line 29, in <module>
domega=c_d(len(t),t,z)
File "second_1.py", line 25, in c_d
dy[1:-1]=(y[2:]-y[0:-2])/(x[2:]-x[0:-2])
TypeError: unsupported operand type(s) for -: 'list' and 'list'
下面是函数:
import numpy as np
from pylab import *
import time
t_initial=time.time()
clf()
t,hp,hn= np.loadtxt("Richardson.dat", usecols=(0,1,2),comments='#', unpack=True) # to select just a few columns
print(time.time()-t_initial)
def phi(y,x):
return(np.arctan(y/x))
phase=[0.0]*len(t)
phase=phi(hp[0:],hn[0:])
#plot(t,phase)
#show()
def c_d(order,x,y):
dy=[0.0]*order
dy[0]=(y[1]-y[0])/(x[1]-x[0])
dy[-1]=(y[-1]-y[-2])/(x[-1]-x[-2])
dy[1:-1]=(y[2:]-y[0:-2])/(x[2:]-x[0:-2])
return(dy);
z=c_d(len(t),t,phase);
plot(t,z)
print(len(z)-len(t))
domega=c_d(len(t),t,z)
plot(t,domega)
show()
发布于 2015-01-20 20:59:32
这个问题在错误消息中解释得非常清楚:'-‘操作数不适用于类型列表。
(y[2:]-y[0:-2])/(x[2:]-x[0:-2])
y2:对列表进行切片并返回列表。y0:-2也对列表进行切片并返回列表。所以你有两个列表。
y[2:] (a list) -(your operator) y[0:-2] (a list)
并且没有定义list - list (没有语法:'listObject‘- 'listObject')。
但是:+运算符是为列表定义的(示例):
l = ["ja"]
m = ["nein"]
n = l + m
print n
# output: ['ja', 'nein']
在这里看看这些不同类型的可能运算符:https://docs.python.org/2/library/stdtypes.html
发布于 2015-01-20 21:07:40
正如Cyber和ProgrammingIsAwsome所说的,错误是在线的
(y[2:]-y[0:-2])/(x[2:]-x[0:-2])
在那里你实际尝试减去列表。
你可以明确地写:
for i in range(1, order - 1):
dy[i]=(y[i+1]-y[i-1])/(x[i+1]-x[1-1])
https://stackoverflow.com/questions/28045804
复制相似问题