Python：數值微分的索引錯誤

Question

在我的代碼中，我從時間、position 測量中提取速度和加速度，並且在執行數值微分時收到索引錯誤：

    VelocityVsTime = np.empty((2,0), float)

for i in range(1, len(PosVsTime[0])-1):
    velocity = (PosVsTime[1][i+1] - PosVsTime[1][i-1]) / (PosVsTime[0][i+1] - PosVsTime[0][i-1])
    VelocityVsTime = np.append(VelocityVsTime, [[PosVsTime[0][i]], [velocity]], axis = 1)
    
#print(VelocityVsTime)

AccelerationvsTime = np.empty((2,0), float)

for j in range(1, len(VelocityVsTime[1])-1):
    acceleration = (VelocityVsTime[1][i+1] - VelocityVsTime[1][i-1]) / (VelocityVsTime[0][i+1] - VelocityVsTime[0][i-1])
    AccelerationvsTime = np.append(AccelerationvsTime, [VelocityVsTime[0][i]], [acceleration], axis=1)
    
print(AccelerationvsTime)

錯誤是：

IndexError: index 50 is out of bounds for axis 0 with size 49

有關如何糾正此問題的任何提示？ 謝謝

Answer 1

繼承人完整代碼：錯誤發生在我聲明加速度變量的第 42 行

import numpy as np

import matplotlib.pyplot as plt


PosVsTime = np.loadtxt("balldata.txt", delimiter=",").transpose()

#print(PosVsTime[0][0])

#t_0 = PosVsTime[0][0]

#pos_0 = PosVsTime[1][0]


#print("The initial state of this system at time = 0 is ", pos_0)


VelocityVsTime = np.empty((2,0), float)

for i in range(1, len(PosVsTime[0])-1):
    velocity = (PosVsTime[1][i+1] - PosVsTime[1][i-1]) / (PosVsTime[0][i+1] - PosVsTime[0][i-1])
    VelocityVsTime = np.append(VelocityVsTime, [[PosVsTime[0][i]], [velocity]], axis = 1)
    
#print(VelocityVsTime)

#plt.errorbar(VelocityVsTime[0], VelocityVsTime[1], fmt = '--k')

AccelerationvsTime = np.empty((2,0), float)

for j in range(1, len(VelocityVsTime[0])-1):
    #acceleration = (VelocityVsTime[1][i+1] - VelocityVsTime[1][i-1]) / (VelocityVsTime[0][i+1] - VelocityVsTime[0][i-1])
    #AccelerationvsTime = np.append(AccelerationvsTime, [VelocityVsTime[0][i]], [acceleration], axis=1)
    
print(AccelerationvsTime)