Trapz和数据集集成

Question

我正在尝试根据数据集的平均值积分定理找到平均值。我的尝试一直无济于事，最终以平均值打印nan。是否有任何代码建议可以对其进行修改并获得平均值？进口商品为scipy，numpy和pandas。我想在某个部分而不是整个数据集上获得MV。

def meanvalue():
         b = float(input("What is the upper limit? "))
         a = float(input("What is the lower limit? "))
         integral = np.trapz(n_array,thz_array, axis = 0)
         mean_value = float(integral/(b-a))
return b, a, mean_value
data_1 = np.genfromtxt('C:/Users/Light_Wisdom/KeanDataCSV/index_test1.txt', delimiter=",")
thz_array = data_1[:,0] 
phase_array = data_1[:,1]-data_1[:,2]
thick=35 #micron
thick=thick*10**-6 #m
c_m=sp.constants.c
w_array= 2 * np.pi*(thz_array*10**12)
n_array = 1 + c_m * phase_array / (w_array *thick)
b, a, mean_value = meanvalue()
ax.hlines(mean_value, xmin = a, xmax = b)