******所以程序都是matlab版本**********
1，mk是趋势性检验程序，程序中没有考虑结（就是序列中相等的值）
2，mkabrpt是突变点分析的程序，mkabrpt会调用smk；smk是计算秩的程序。
3，上海气温.mat是测试数据。你可以试试。
4，程序代码里有简单的参数说明
特别说明：
*上述程序均反复测试绝对没问题。
* mk是我从网上收集的。mkabrpt和smk是我自己写的。
* 我的mkabrpt程序和魏凤英编著的《现代气候统计诊断与预测技术》书上的结果可能不一样。我仔细分析的结论是我的是对的。
smk程序的13行-15行：
if x(i)>x(j)
m(i)=m(i)+1;
end
上面是我的写法。
如果改成:
if x(i)>=x(j)
m(i)=m(i)+1;
end
结果会和那本书上结果一致。
按我的理解不应该是 x(i)>=x(j)，因为如果那样的话，一个等值序列也会检测到显著增长，那是不合理的。
请自己斟酌。
*欢迎大家和我交流 QQ：82033959
评论30