- 积分
 - 1070
 
	- 贡献
 -  
 
	- 精华
 
	- 在线时间
 -  小时
 
	- 注册时间
 - 2017-3-4
 
	- 最后登录
 - 1970-1-1
 
 
 
 
 
 
 | 
	
 
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册 
 
 
 
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突变检验 
 
 
 
- 
滑动t检验 
 
 
 
 
 
 
 
 |