滑动T检验相关，请大佬指点一下

安静多好

附上程序：              因为程序也是我在论坛上找到的，由于是小白，所以就根据自己需要处理的数据，稍微改了一下相关数值。   我做的是1961年-2008年的径流突变检测。序列长度选了5，改了一些数值，还有一些不确定怎么改，希望大佬指点一下。

M=load('数据.txt');
x=M(:,1);y=M(:,2);t=zeros(48,1);y1=zeros(48,1);
%MMT方法检验%
%先设定子序列的长度（取n1=5;n2=5)%
n1=5;n2=5;
m=sqrt((1.0/n1)+(1.0/n2));
%开始计算%
for i=5:43
    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:48),t,'r-');
y1(:)=2;                   %n=53-2时，显著性检验标准,99.5%
hold on;
plot(x,y1,'k:',x,-y1,'k:');
axis([1961,2008,-3.5,3]);
xlabel('年');ylabel('t');title('MMT(n=48)检验');

xiaoweizou

想问一下，你的数据第一列是时间吗？下面是我根据自己的数据改的，我也是刚学的，不知道是不是正确，如果有错误请大佬指点一下：
M=xlsread('jiangxue.xlsx');
x=M(:,1);y=M(:,2);t=zeros(40,1);y1=zeros(40,1);
%MMT方法检验%
%先设定子序列的长度（取n1=5;n2=5)%
n1=5;n2=5;%这是选择的步长
m=sqrt((1.0/n1)+(1.0/n2));
%开始计算%
for i=5:35
    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,t,'r-');
y1(:)=2.31;              %自由度为两个序列步长相加减去2，然后查t检验的表，我用的是显著性水平为0.05
hold on;
plot(x,y1,'k:',x,-y1,'k:');
axis([1979,2018,-3.5,3]);
xlabel('年');ylabel('t');title('MMT(n=40)检验');

xiaoweizou

想问一下，你的数据第一列是时间吗？下面是我根据自己的数据改的，我也是刚学的，不知道是不是正确，如果有错误请大佬指点一下：
M=xlsread('jiangxue.xlsx');
x=M(:,1);y=M(:,2);t=zeros(40,1);y1=zeros(40,1);
%MMT方法检验%
%先设定子序列的长度（取n1=5;n2=5)%
n1=5;n2=5;%这是选择的步长
m=sqrt((1.0/n1)+(1.0/n2));
%开始计算%
for i=5:35
    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,t,'r-');
y1(:)=2.31;              %自由度为两个序列步长相加减去2，然后查t检验的表，我用的是显著性水平为0.05
hold on;
plot(x,y1,'k:',x,-y1,'k:');
axis([1979,2018,-3.5,3]);
xlabel('年');ylabel('t');title('MMT(n=40)检验');