Python如何实现MK检验
发表于:2025-01-31 作者:千家信息网编辑
千家信息网最后更新 2025年01月31日,本篇内容介绍了"Python如何实现MK检验"的有关知识,在实际案例的操作过程中,不少人都会遇到这样的困境,接下来就让小编带领大家学习一下如何处理这些情况吧!希望大家仔细阅读,能够学有所成!MK检验:
千家信息网最后更新 2025年01月31日Python如何实现MK检验
本篇内容介绍了"Python如何实现MK检验"的有关知识,在实际案例的操作过程中,不少人都会遇到这样的困境,接下来就让小编带领大家学习一下如何处理这些情况吧!希望大家仔细阅读,能够学有所成!
MK检验:时间序列进行检测,并找出突变点,本文参考网上的matlab程序改写为python代码如下:
import numpy as npimport pandas as pdimport matplotlib.pyplot as plt#读取时间序列数据data = pd.read_csv('')#定义时间和径流数据列x = list(range(len(data)))y = data.to_list#获取样本数据n = len(y)#正序计算#定义累计量序列Sk,长度n,初始值为0Sk = np.zeros(n)UFk = np.zeros(n)#定义Sk序列元素ss = 0#i从2开始,根据统计量UFk公式,i=1时,Sk(1)、E(1)、Var(1)均为0,此时UFk无意义,因此公式中,令UFk(1) = 0for i in range(2,n): for j in range(1,i): if y[i]>y[j]: s += 1 Sk[i] = s E = i * (i - 1)/4 Var = i * (i - 1) * (2 * i + 5)/72 UFk[i] = (Sk[i] - E)/np.sqrt(Var)#逆序计算,构造逆序列y2,长度为n,初值为0#定义逆累计量序列Sk2,长度n,初始值为0#定义逆统计量序列Sk2,长度n,初始值为0y2 = np.zeros(n)Sk2 = np.zeros(n)UBk = np.zeros(n)#s归零s = 0#按时间序列逆转样本y2 = y[::-1]# i从2开始,根据统计量UBk公式,i=1时,Sk(1)、E(1)、Var(1)均为0,此时UFk无意义,因此公式中,令UBk(1) = 0for i in range(2, n): for j in range(1, i): if y2[i] > y2[j]: s += 1 Sk2[i] = s E = i * (i - 1) / 4 Var = i * (i - 1) * (2 * i + 5) / 72 UBk[i] = -(Sk2[i] - E) / np.sqrt(Var)#逆转逆序列UBk2 = UBk[::-1]#画图plt.figure(figsize=(10,5))plt.plot(range(1 ,n+1),UFk,label = 'UFk',color = 'orange')plt.plot(range(1 ,n+1),UBk2,label = 'UBk',color = 'cornflowerblue')plt.ylabel('UFk-UBk')x_lim = plt.xlim()plt.plot(x_lim,[-1.96,-1.96],'m--',color = 'r')plt.plot(x_lim, [0,0],'m--')plt.plot(x_lim,[1.96,1.96],'m--',color = 'r')plt.show()
"Python如何实现MK检验"的内容就介绍到这里了,感谢大家的阅读。如果想了解更多行业相关的知识可以关注网站,小编将为大家输出更多高质量的实用文章!
序列
公式
时间
长度
检验
数据
时间序列
统计
内容
更多
样本
知识
实用
学有所成
接下来
代码
元素
初值
困境
实际
数据库的安全要保护哪些东西
数据库安全各自的含义是什么
生产安全数据库录入
数据库的安全性及管理
数据库安全策略包含哪些
海淀数据库安全审计系统
建立农村房屋安全信息数据库
易用的数据库客户端支持安全管理
连接数据库失败ssl安全错误
数据库的锁怎样保障安全
软件开发属于成果转化
什么服务器玩倩女幽魂的返利最大
九曳供应链软件开发咋样
平顶山bim软件开发工程
关于网络安全的设计论文
计算机网络技术职业发展规划
怎样制作表格数据库
施云波无线传感网络技术概论
网络安全法合集
实用办公软件开发工具
gps修改定位软件开发
软件数据库连接
网页无法连接到SQL数据库
江西本地软件开发预算
全国网络安全等级测评
南宁网络安全
济南中智互联网科技
r420服务器u盘
金融公司数据库恢复
sql数据库语句怎么备注
ros 拨号服务器
电脑的网络安全密码如何获得
mfc怎么连接数据库
软件开发设计阶段任务
张家港智能服务器货源充足
服务器链接繁忙怎么办
web属于网络安全技术吗
网络安全设计实验题目
航空公司软件开发岗笔试
软件开发中系统描述指啥