爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4367|回复: 0

[讨论] M-K突变检验的结果和滑动t检验的结果不一致。

[复制链接]

新浪微博达人勋

发表于 2017-4-7 17:11:03 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
这是我做出来的m-k突变检验和滑动t检验的两幅图,但是它们的突变情况不一样,请帮我解答一下,这是什么原因。谢谢各位。(下面附滑动t检验的代码)
clc,clear;
excelFile='xx';
data=xlsread(excelFile);
timeseries=data(:,1);  
dataseries=data(:,2);
datacount=length(62);
step=10;  
v=18;
ttest=2.878;
len1=step;
len2=step;
x=timeSeries(step:dataCount-step);
for i=step:dataCount-step
    n1=dataSeries(i-step+1:i);
    n2=dataSeries(i+1:i+step);
    mean1=mean(n1);
    mean2=mean(n2);
    c=(len1+len2)/(len1*len2);
    var1=1/len1*sum((n1-mean1).^2);
    var2=1/len2*sum((n2-mean2).^2);
    delta1=len1*var1+len2*var2;
    delta=delta1/(len1+len2-2);
    t(i-step+1)=(mean1-mean2)/sqrt(delta*c);
end
%% 制图
figure(1);
plot(x,t,'r-','linewidth',1.5);
xlabel('t (year)','FontName','TimesNewRoman','FontSize',12);
ylabel('统计量','FontName','TimesNewRoman','Fontsize',12);
axis([min(x),max(x),-4,4]);
hold on
plot(x,0*ones(i-step+1,1),'-.','linewidth',1);
plot(x,ttest*ones(i-step+1,1),':','linewidth',1);% 更改数字
plot(x,-ttest*ones(i-step+1,1),':','linewidth',1);% 更改数字

m-k突变检验

m-k突变检验

滑动t检验

滑动t检验
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表