- 积分
- 4568
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-3-1
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
最近在做某气象场的趋势要用到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也是好事)
|
-
-
mktrend.m
1.12 KB, 下载次数: 198, 下载积分: 金钱 -5
售价: 1 贡献 [记录]
|