0

Ich baue ein Mesh-Gitter für das x-Grid und x-Vektor und auch Zeitraster ein, aber ich habe wieder ein Array für x (Position), die nur zwischen 0 und 20 liegen sollte und t (Zeit) wäre von 0 bis 1000, also um eine Wärmeleitungsgleichung zu lösen. Aber jedes Mal, wenn ich für zB möchte, mache ich die Anzahl der Schritte 10, erhalte ich eine Fehlermeldung:IndexError: Index 10 ist außerhalb der Grenzen für Achse 0 mit Größe 10

"Traceback (most recent call last): 
File "/home/universe/Desktop/Python/Heat_1.py", line 33, in <module> 
x[i] = a + i*h 
IndexError: index 10 is out of bounds for axis 0 with size 10" 

Hier mein Code:

from math import sin,pi 
import numpy 
import numpy as np 

#Constant variables 
N = int(input("Number of intervals in x (<=20):")) 
M = int(input("Number of time steps (<=1000):")) 

#Some initialised varibles 
a = 0.0 
b = 1.0 
t_min = 0.0 
t_max = 0.5 

# Array Variables 
x = np.linspace(a,b, M) 
t = np.linspace(t_min, t_max, M) 


#Some scalar variables 
n = []       # the number of x-steps 
i, s = [], []     # The position and time 

# Get the number of x-steps to use 
for n in range(0,N): 
    if n > 0 or n <= N: 
     continue 

# Get the number of time steps to use 
for m in range(0,M): 
    if m > 0 or n <= M: 
     continue 

# Set up x-grid and x-vector 
h =(b-a)/n 
for i in range(0,N+1): 
    x[i] = a + i*h 

# Set up time-grid 
k = (t_max - t_min)/m 
for s in range(0, M+1): 
    t[s] = t_min + k*s 

print(x,t) 
+0

Arrays sind Null indexiert dh 0-9 –

Antwort

0

Sie versuchen, Index außerhalb des Bereichs:

for s in range(0, M+1): 
    t[s] = t_min + k*s 

wechseln zu:

for s in range(M): 
    t[s] = t_min + k*s 

Und es funktioniert.

Sie erstellen t mit einer Länge von M:

t = np.linspace(t_min, t_max, M) 

So können Sie nur Zugriff auf M Elemente in t.

Python beginnt immer mit null zu indexieren. Deshalb:

for s in range(M): 

tun M Schleifen, während:

for s in range(0, M+1): 

M+1 Schleifen tun.

+0

Ja, es funktioniert jetzt wirklich, aber die Werte von x sollen in gleichmäßigen Abständen unterteilt werden, die von laufen würden: a = 0.0 bis b = 1.0, und es einfach nicht scheinen mir die Werte zu geben, die ich will, da es Werte über 1 gibt – Tonikami04

+0

Danke, wie ich jetzt es richtig arbeite und laufen lasse @MikeMuller – Tonikami04

+0

Großartig, dass es geholfen hat. BTW, Sie können eine Antwort (wenn Sie Ihr Problem lösen) [akzeptieren] (http://meta.stackexchange.com/questions/5234/how-does-accepting-an-answer-work). –