- 积分
- 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检验
|