python - Python:具有数值微分的索引错误
问题描述
在我的代码中,我从时间、位置测量中提取速度和加速度,并且在执行数值微分时收到索引错误:
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
有关如何纠正此问题的任何提示?谢谢
解决方案
继承人完整代码:错误发生在我声明加速度变量的第 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)
推荐阅读
- terraform - 是否可以在 Azure 中启用 linux 性能计数器来记录分析工作区?
- java - 在 K 次试验后,数字 N 减少到非正值的概率
- c++ - 在 QT 中使用 dll 文件
- python - 将 string.split() 作为参数传递给另一个函数时不起作用
- typescript - 在创建堆栈导航器之前如何导航到屏幕?
- sql - 使用计数选择最小值和最大值
- python - 在 Python 中实现表面模糊
- c - 有没有办法通过使用索引来迭代/递归二叉树?
- javascript - 如何在 Javascript HTML 文件中显示图像?
- java - 能够从 PDF 中删除页眉、页脚但同时删除页码 ex: 1/1 在此能够删除 1 但仍然 / 出现在文件的末尾