MATLAB实现MK突变检验

小菲菲

好吧，我在尝试一下mk检验

游泳的鱼

非常好，谢谢楼主~

18021719881

点赞，我也是搜了网上matlab程序但是和书上的不一致，这个画出来是一样的

abc111

matlab显示    未定义与 'double' 类型的输入参数相对应的函数 'MannKendall'。

是什么原因？

伽蓝鸟

abc111 发表于 2019-2-13 15:47
matlab显示    未定义与 'double' 类型的输入参数相对应的函数 'MannKendall'。

是什么原因？

这是个函数啊，你要把函数存起来，然后去调用它

zane

这个根据fortran改变的matlab代码，计算的和书本上的结果一样。但是论坛好多帖子，根据书本上的公式计算出来的，反而跟书本的结果不一致！

FeedyH09

qinkong 发表于 2015-6-2 08:38
看了你的程序后，觉得魏凤英老师的书上写的有点问题，书上写的是“按时间序列逆序再重复上述过程，同时使得 ...

符淙斌跟王强老师的《气候突变的定义和检测》里，M-K法计算反序列就是i’=N+1-i来着

zhangzncu

感谢分享，谢谢

敏子494

先收藏

drx

何必写的这么麻烦，直接编写一个函数

[m n]=size(data);
for t=2:1:m
    for j=1:1:m
      v(t-1,j)=sign(data(t-1,1)-data(j,1));
      V(t-1)=sum(v(t-1,:));
    end
end
U=cumsum(V);
loc=find(abs(U)==max(abs(U)));
K=max(abs(U));
pvalue=2*exp((-6*K^2)/(m^3+m^2));
res=[loc;pvalue];