Зачем вызывать вызов дважды, показывая ошибку

Question

У меня возникла проблема с вызовом функции дважды. Если я прокомментирую мои последние 3 строки и сохраняю команду show(), то у меня не будет никакой ошибки, и что-то произойдет, поскольку они предполагают .. но если я dnot комментарий не как я пытаюсь вызвать последнюю функцию снова показывает эту ошибку

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() 

Answer 1

Как объяснены в других ответах, вы не можете вычесть регулярные списки Python. Так почему он работает в первый раз и не работает второй? Давайте посмотрим на код.

t, hp, hn = np.loadtxt(...) 
... 
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) 
... 
domega = c_d(len(t), t, z) 
... 

При первом вызова c_d параметры x и y кажется numpy массив (по крайней мере, t и phase являются результатом numpy вызовов функциев), а также для тех, - является юридическим действием. Но внутри c_d, вы создаете и, наконец, вернуть регулярный список Python, dy, поэтому, когда вы затем вызвать c_dснова с результатом первого вызова, как y, эта часть y[2:]-y[0:-2] потерпит неудачу, поскольку y сейчас обычный список ,

Убедитесь, что dy является numpy массив, тоже, т.е. dy = np.array([0.0] *order) или просто dy = np.zeros(order), то он должен работать.

Answer 2

проблема очень четко поясняется в сообщении об ошибке: Оператор '-' не применяется для списка типов.

(y[2:]-y[0:-2])/(x[2:]-x[0:-2]) 

y [2:] разрезает список и возвращает список. y [0: -2] разрезает список и возвращает список. Итак, у вас есть 2 списка.

y[2:] (a list) -(your operator) y[0:-2] (a list) 

И список - список не определен (нет синтаксиса: 'ListObject' - 'ListObject').

НО: оператор + определен для списков (например):

l = ["ja"] 
m = ["nein"] 

n = l + m 

print n 
# output: ['ja', 'nein'] 

Посмотрите здесь для этих различных видов возможных операторов: https://docs.python.org/2/library/stdtypes.html

Answer 3

Как заявил 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])