Fehler mit Index. Index 0 ist außerhalb der Grenzen für Achse 0 mit Größe 0

Question

Wenn ich meinen Code ausführen ich den Fehler

IndexError: index 0 is out of bounds for axis 0 with size 0 on line 42 for v[0]=v0.

ich mit Python nicht vertraut bin und weiß nicht, wie dies zu debuggen.

""" Density at height y """ 
def rho(rh, y): 
rho = rh*math.exp(-y/7249) 
return rho 

""" Acceleration of skydiver """ 
def accl(m, g, rho, v, As, Cs, Ap, Cp): 
if y[i+1] >= 1200: 
    accl = (1/(2*m) * rho(rh, y) * v**2 * As * Cs - g) 
else: 
    accl = (1/(2*m) * rho(rh, y) * v**2 * (As * Cs + Ap * Cp) - g) 
return accl 

h0 = 4000 # initial hieght 
dt = 0.1 #timestep 
n = int(180/dt) 
#lists for displacement, time, and velocities 
v= np.array([]) 
v[0] = v0 
t= np.array([]) 
t[0] = t0 
y= np.array([]) 
y[0] = h0 

#calculation of time, velocity and displacement 
for i in range(n): 
t[i+1] = t[i] + dt 
v[i+1] = v[i] + accl(m, g, rho, v[i], As, Cs, Ap, Cp) * dt 
y[i+1] = y[i] + v[i+1] * dt 
if y[i+1] < 0: 
    #ends plot once skydiver has landed 
     break 
# Plot 
plt.plot(t, y, label = 'Position m') 
plt.plot(t, v, label = 'Speed m/s') 
plt.plot(t, accl, label = 'Acceleration m/s^2', color = 'green') 
plt.xlim([t[0], t[-1]]) 
plt.xlabel("Time (seconds)") 
plt.legend() 

Answer 1

Numpy-Arrays, die aus leeren Listen erstellt wurden, können nicht indiziert werden - sie verfügen nicht über Elemente, die indexiert werden können.

np.array([]).shape 
# (0,) 

Wenn Sie die Dimensionen v im Voraus kennen, können Sie np.empty() verwenden:

n = 10 
np.empty(n, dtype=float) 
# array([ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.]) 

Alternativ initialisieren nur v mit v0:

v = np.array([v0]) 

Answer 2

v= np.array([]) 
v[0] = v0 

Diese nicht macht Sinn: Sie versuchen, einen Wert zu setzen das erste Element eines leeren Arrays (das kein Element hat!).

Versuchen Sie stattdessen: v = np.array([v0])

, wenn Sie Ihre Array wachsen, dies zu tun:

v = [] 
# whatever logic to grow v like: 
for i in range(10): 
    v.append(2*i-1) 
# now turn it into an np.array: 
v = np.array(v) 

Answer 3

ersetzen v = np.array([]) mit v = np.zeros(n+1). Sie können nicht auf Elemente des numpy-Arrays zugreifen, die außerhalb seiner Grenzen liegen. Ihre Deklaration erstellt ein leeres numpy Array.