mk趋势检验（非突变检验） matlab实现代码

lycwin

最近在做某气象场的趋势要用到Mann-kendall 趋势检验方法（m-k test）,在家园搜了一圈发现要么从网上荡的国外的，估计自己也没搞清楚怎么回事，要么是Fortran的程序，所以综合了下自己写了个函数挂上来，关于原理这位前辈已经讲的很透彻了http://bbs.06climate.com/forum.p ... p;extra=&page=1，函数在他的基础上改进了下，在此谢过

预览：
function [slope,zc,sign]=mktrend(x)
%此为MannKendall趋势检验的程序
% x: 输入变量，需要检验变化趋势的统计量,注意是一维数组，如果是多维写个循环就好了（比如144×98×9（经纬度×时间序列））
% zc: 检验统计量，zc>0表示上升，zc<0表示下降，绝对值达到1.65/1.96/2.58，分别表示通过90/95/99%显著性检验(具体希望多少设置在33行)
% sign为显著性结果，1显著，0不显著（对应上条显著性水平）
% slope斜率（采用Theil-Sen方法计算(ncl官网mk函数有详细的介绍)）



传送门：http://www.ncl.ucar.edu/Document ... /trend_manken.shtml（ncl还是很好用的，直接内置函数就解决了，原理讲的也很透彻，现在并给pathon也是好事）

暖气团

len=size(x,3);
如果x是一列数据的话，显示结果是NAN呢

lycwin

暖气团 发表于 2019-3-31 18:22
len=size(x,3);
如果x是一列数据的话，显示结果是NAN呢

如果一列数据直接len=size(x);吧，3代表第三维

暖气团

OK，就是求证一下，已经这么操作了。谢谢

kongmch15

请问ncl里面那个trend_manken函数得到的显著性怎么看 啊  怎么和90,95,99%的对应上

kongmch15

楼主  看到回复一下  非常需要你的帮助

lycwin

kongmch15 发表于 2019-4-11 09:38
楼主  看到回复一下  非常需要你的帮助

首先放下官网的例子：x(NT)            : p=trend_manken(x,False,0) => p(2) where p(0)=probability,  p(1)=trend
我理解是p(0)=probability就是显著性水平，若你要95%以上的设个判断p（0）>0.95的筛选出来即可，p（1）就是算的斜率k,本人现在转战matlab,ncl没有特别研究深入，讲的不对的地方请大家指出来哈~

kongmch15

lycwin 发表于 2019-4-11 13:42
首先放下官网的例子：x(NT)            : p=trend_manken(x,False,0) => p(2) where p(0)=probability,   ...

谢谢！！

不懂爱的秋刀鱼

楼主，如果是一列数据不能直接改为len=size(x);而是应该修改为len=length(x);才可以运行成功。

lycwin

不懂爱的秋刀鱼 发表于 2019-5-30 23:20
楼主，如果是一列数据不能直接改为len=size(x);而是应该修改为len=length(x);才可以运行成功。

我感觉用size还是length可能和行数据或者是列数据有关吧