python如何实现MK突变检验方法,代码复制修改可用（python如何导入模块soundfile）干货满满

🌐 经济型：买域名、轻量云服务器、用途：游戏网站等《腾讯云》特点：特价机便宜适合初学者用点我优惠购买

🚀 拓展型：买域名、轻量云服务器、用途：游戏网站等《阿里云》特点：中档服务器便宜域名备案事多点我优惠购买

🛡️ 稳定型：买域名、轻量云服务器、用途：游戏网站等《西部数码》特点：比上两家略贵但是稳定性超好事也少点我优惠购买

import matplotlib.dates as mdates #處理日期
import matplotlib.pyplot as plt
import numpy as np
from pylab import mpl
from matplotlib.pyplot import MultipleLocator
mpl.rcParams[‘font.sans-serif’]=[‘SimHei’] #防止标题出现乱码。
plt.rcParams[‘axes.unicode_minus’]=False #防止出现图上的负数为方框。
# y值和x值分别输入六个站点的最大冻土深度值，将值以列表的方式导入
a=[150,150,114,109,96,95,83,76,109,80,115,80,94,86,133,91,110,116,114,128,172,172,
162,121,175,151,110,92,116,156,134,110,89,97,109,157,153,105,76,87,122,78,97,93,141,162,
123,133,161,128,138,104,133,102,140,109,118,86,126,92,121,149,116] #这个部分值可以替换成为要检验的气温、水文等值
x_values=list(range(1961,2022))
uf,ub=mktest(a)
plt.figure(figsize=(8,4)) #图片的大小
plt.plot(uf,’r’,label=’UFk’)
plt.plot(ub,’b’,label=’UBk’)
plt.xticks([0,5,10,15,20,25,30,35,40,45,50,55,60],[‘1960′,’1965′,’1970′,’1975′,’1980′,’1985′,’1990′,’1995′,’2000′,’2005′,’2010′,’2015′,’2020’,])
#将默认的x轴数值替换为年份的X轴，默认是0-61，一共62个值，代表X轴内容。
# 0.01显著性检验
plt.legend()
plt.axhline(1.96)
plt.axhline(-1.96)
#设置图片的标签（标题）
plt.title(“富蕴点最大冻土深度突变检验结果”)#x轴上的名字
plt.xlabel(“年份（1960年-2022年）”)#x轴上的名字
plt.ylabel(“突变值波动参数”)#y轴上的名字
plt.grid() #形成网格线输出
x_major_locator=MultipleLocator(5)
plt.show()