急、急、急毕业论文需要用到滑动T检验，跪求各位高人指点指点。

sswyj

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


我的是30年的数据，各位大神啊，第二行怎么修改呀，帮帮我吧。

dr-liu

感谢楼主，把代码都给我们免费呈现了，真是好心人呐。我最近也在研究这个，也有同样的问题。显著性水平在哪体现呢？

dr-liu

难道论坛里就没有会的吗？还是会而不愿意教呢？

dr-liu

{:5_235:}{:5_235:}

dr-liu

感谢楼主分享