滑动t检测

小猫伸懒腰

%%%该程序用MMT（滑动t)方法检验气候突变%%%
%读数据%
M=load('C:\Users\gt5fvbhy6jmnihy6jmn\Desktop\sy3.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)检验');




为什么总说这行plot(x,y1,'k:',x,-y1,'k:');
错误使用plot 矢量长度必须相同
为什么呢

lovjiangyp

{:eb502:}{:eb502:}

ttttt

我也在写论文，看不懂啊亲

desperado

顶顶顶顶顶

诗遇上歌

急需matlab程序，要用滑动T

leizigoal

有没有excel版本的