滑动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)检验');

我这是1971-2000年30年的气象资料，请各位大神帮我看看第二行括号里的冒号那怎么改，还有检验水平怎么看，非常感谢！！！

dr-liu

良心帖，把代码都贴出来了，不像有的人还要积分。。我也是遇到这样的问题，第二行x,y括号里的东西表示什么？？？

dr-liu

{:5_217:}{:5_217:}

dr-liu

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

小英

这个代码就是个数为30的序列，可以把第第2列 t=zeros(30,1)改成t=zeros(21,1),

小英

dr-liu 发表于 2016-7-24 15:14
良心帖，把代码都贴出来了，不像有的人还要积分。。我也是遇到这样的问题，第二行x,y括号里的东西表示什么 ...

x表示文件M第一列的全部数据；y表示第二列的全部数据

kangkangstudy

感谢，试试您的代码

dr-liu

小英 发表于 2016-12-9 16:04
这个代码就是个数为30的序列，可以把第第2列 t=zeros(30,1)改成t=zeros(21,1),

非常感谢！

xdeamond

感谢楼主分享

gxx123

EXCEL一啦就出来了