滑动窗口的t检验分析Python代码
发布时间
阅读量：
阅读量
近期在进行时间序列突变点检测相关工作时，经过长时间查找，发现目前缺乏对应的Python实现代码，因此根据相关公式自行编写了一段程序。若其中存在不妥之处，恳请各位专家和同行提出宝贵意见与建议。
代码如下：
import numpy as np
import matplotlib.pyplot as plt
def huaT(inputdata):
inputdata = np.array(inputdata)
n = inputdata.shape[0]
t = np.zeros(90)
n1 = 5.0 #把n1, n2设置成float型
n2 = 5.0
n11 = 1 / n1
n22 = 1 / n2
m = np.sqrt(n11 + n22)
k=[]
for i in range (5, 85):
x1_mean = np.mean(inputdata[i-4 : i]) #均值
x2_mean = np.mean(inputdata[i+1 : i+5])全部评论 (0)
还没有任何评论哟~