用Matlab进行MK检验
用Matlab进行Mann-Kendall分析
一、M-K趋势分析法与M-K突变检验的原理介绍
1、Mann-Kendall趋势分析法
非参数检验,又称任意分布检验,它不对变量的分布做严格假定,检验不针对特定的参数,而是模糊地对变量分布的中心位置或分布状态做检验,由于其不对总体分布做严格假定,因而适用性强[12]。因此,本文采用非参数的Mann-Kendall检验法对昌马河流域近50年的气候水文要素时间序列显著性检验,定量反映变化趋势的显著性。
计算公式如下:
(1.1)
在公式中,
(1.2)
(1.3)
(1.4)
式中:xk,xi为连续的气候、水文等数据序列,n为数据集合的总长度,t为每个单位的宽度,Σ表示所有单位的总和。
衡量趋势大小的指标为:
(1.5)
式中:1
Mann-Kendall检验如下所示:
零假设H0:β=0
当 ,拒绝H0假设。
式中: 为标准正态方差,α为显著性检验水平。
2、Mann-kendall突变检验
气候系统变化是一个不稳定且不连续的变化过程,而检验其变化的常用方法之一就是Mann-kendall突变检验方法[13],该方法对于变化要素从一个相对稳定状态变化到另一个状态的变化检验非常有效。且广泛应用于水文,气候,化学,矿物成分检验等各个方面。Mann-kendall突变检验方法如下: