Question

我正在尝试使用Simpson的方法集成一组数据。我发现可以执行以下操作：result = simps（速度，时间，dx =（1/60））。我希望间隔为秒而不是分钟（因此为1/60）-但更改dx的值不会更改我的结果。如何在几秒钟内与dx集成？

# import the necessary modules
import numpy as np
import matplotlib.pyplot as plt

# reads the experimnetal data
a=np.loadtxt("VelvstUniform.txt")

# assigns the data columns
time=a[:,0]
velocity=a[:,1]

# graphs in the same figure Sunspots vs Year and LEO vs Year
plt.plot(time,velocity,markersize=10, color="green", marker="o", 
linestyle="-")
plt.xlabel("Time (min)", fontweight="bold", fontsize="large")
plt.ylabel("velocity (miles/hour)", fontweight="bold", fontsize="large", 
color="green")

from scipy.integrate import simps
result=simps(velocity, time, dx=(1/60))

print(result)

Answer 1

我认为您在那里误解了一些东西

dx更改集成块的宽度，但不更改其单位。为此，必须将值乘以所需的比例因子。而且我认为您宁愿在小时内使用time，因为您以英里/小时为单位来测量速度。

您只需更改result；它的当前单位是min * mi / h，因此您必须乘以1/60才能得到结果（以英里为单位）。

Python 3.6：使用Simpson方法

1 个答案: