MMT（滑动t)方法检验气候突变

892034114

%%%该程序用MMT（滑动t)方法检验气候突变%%%
%读数据%
M=load('global.txt');
x=M(:,1);y=M(:,2);t=zeros(134,1);y1=zeros(142,1);
%MMT方法检验%
%先设定子序列的长度（取n1=5;n2=5)%
n1=5;n2=5;
m=sqrt((1.0/n1)+(1.0/n2));
%开始计算%
for i=5:138
    averx1=mean(y(i-4:i));
    averx2=mean(y(i+1:i+5));
    s1=var(y(i-4:i));  
    s2=var(y(i+1:i+5));
    s=sqrt((n1*s1+n2*s2)/(n1+n2-2));
    t(i-4)=(averx2-averx1)/(s*m);  
end
plot(x(5:138),t,'r-');
y1(:)=3.2;                   %n=14时，显著性检验标准,99.5%
hold on;
plot(x,y1,'k:',x,-y1,'k:');
axis([1861,2002,-5,7]);
xlabel('年');ylabel('t');title('MMT(n=5)检验');

GGGGDL

话说你这里的自由度不是14吧，不应该是n1+n2-2=8？？如果是我理解的不对也请大神指点指点

yzby

GGGGDL 发表于 2017-3-30 21:35
话说你这里的自由度不是14吧，不应该是n1+n2-2=8？？如果是我理解的不对也请大神指点指点

我也觉得是，n=8

想安静的风

不管论坛里有没有！我们先感谢你的共享！

Aires

谢谢分享，可以参考

mofangbao

貌似很眼熟  谢谢啦

散落灬

感谢分享，正好学到这里，保存下

上将皮克

good!!!!!!111

kongfeng0824

谢谢分享                        

空兵营

感谢分享，学到了

永远的beyond

感谢分享，学了很多

永远的beyond

感谢分享 学到了