Python编写M-K突变检验

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
plt.rcParams[‘font.sans-serif’] = [‘SimHei’] #用来正常显示中文标签
plt.rcParams[‘axes.unicode_minus’] = False #用来正常显示负号

input_data = pd.read_csv(r"D:\PycharmProjects\shuiziyuan\M-K检验\runoff.csv") # 输入文件的绝对路径，数据样式如下图

print(input_data.head())
print(len(input_data))

正序列计算

Sk_list = [0] # 定义累计量序列Sk，初始值=0
UFk_list = [0] # 定义统计量UFk，初始值 =0
Sk = 0 # 定义Sk序列元素s，初始值 =0
Exp_Sk_value_list = [0] # 均值
Var_Sk_value_list = [0] # 方差

i从1开始，因为根据统计量UFk公式，i=0时，Sk(0)、E(0)、Var(0)均为0

此时UFk无意义，因此公式中，令UFk(0)=0

for i in range(1, len(i