我在数值上为x-grid和x-vector以及时间网格设置了网格,但是我还是为x
(位置)设置了一个数组,该数组只能在0到20之间,而t
(时间)则是从因此,从0到1000直到求解热方程。但是每次我想要例如,使步骤数达到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"
这是我的代码:
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)
最佳答案
您尝试索引超出范围:
for s in range(0, M+1):
t[s] = t_min + k*s
改成:
for s in range(M):
t[s] = t_min + k*s
而且有效。
您创建的
t
的长度为M
:t = np.linspace(t_min, t_max, M)
因此,您只能访问
M
中的t
元素。Python总是从零开始索引。所以:
for s in range(M):
将执行
M
循环,而:for s in range(0, M+1):
将执行
M+1
循环。关于arrays - IndexError : index 10 is out of bounds for axis 0 with size 10,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/36489042/