matlab中如何实现滑动T检验

TVXQ

请问各位大师如何用MATLAB实现这个例子的滑动T检验?

jiajian0512

1. clc;clear;
   
2. %% 读取数据
   
3. myData=xlsread('H:\极端气温变化\年最高气温.xls'); % 读取数据
   
4. timeSeries=myData(:,1);  % 时间序列数据
   
5. dataSeries=myData(:,2);  
   
6. dataCount=length(dataSeries);
   
7. %% 设置步长与检验值
   
8. step=10;  % 步长
   
9. v=step+step-2; % 计算自由度
   
10. ttest=2.878; % 查表得t检验值，修改
    
11. len1=step;
    
12. len2=step;
    
13. x=timeSeries(step:dataCount-step);
    
14. for i=step:dataCount-step
    
15.     n1=dataSeries(i-step+1:i);
    
16.     n2=dataSeries(i+1:i+step);
    
17.     mean1=mean(n1);
    
18.     mean2=mean(n2);
    
19.     c=(len1+len2)/(len1*len2);
    
20.     var1=1/len1*sum((n1-mean1).^2);
    
21.     var2=1/len2*sum((n2-mean2).^2);
    
22.     delta1=len1*var1+len2*var2;
    
23.     delta=delta1/(len1+len2-2);
    
24.     t(i-step+1)=(mean1-mean2)/sqrt(delta*c);
    
25. end
    
26. %% 制图
    
27. figure(1);
    
28. plot(x,t,'r-','linewidth',1.5);
    
29. xlabel('t (year)','FontName','TimesNewRoman','FontSize',12);
    
30. ylabel('统计量','FontName','TimesNewRoman','Fontsize',12);
    
31. axis([min(x),max(x),-4,4]);
    
32. hold on
    
33. plot(x,0*ones(i-step+1,1),'-.','linewidth',1);
    
34. plot(x,ttest*ones(i-step+1,1),':','linewidth',1);% 更改数字
    
35. plot(x,-ttest*ones(i-step+1,1),':','linewidth',1);% 更改数字
    
36. legend('t统计量','0.01显著水平');

复制代码

jiajian0512

剩下的自己改改咯

二爷名声在外

http://bbs.06climate.com/forum.p ... amp;extra=page%3D17

TVXQ

二爷名声在外 发表于 2015-4-12 20:32
http://bbs.06climate.com/forum.php?mod=viewthread&tid=17633&extra=page%3D17

不好意思，我不会MATLAB，所以看不懂，请问有没有这个例子的具体操作呀？

二爷名声在外

TVXQ 发表于 2015-4-12 20:48
不好意思，我不会MATLAB，所以看不懂，请问有没有这个例子的具体操作呀？

刚才你问我，我可以回答一句无可奉告，但是你又不高兴，我怎么办？
既然不会就去学啊！例子已经不错了，希望您结合原理仔细看看！再说您为什么不去问发帖人来问我，这个本身也不合适啊。

气象男孩

我也要做这个

微末1314

楼主你是想要现成答案么

气象男孩

我也要做这个

气象小黑

我也在做这个，先看看

michellegoo

这个链接的程序有问题，运行不了

人在江湖

顶。。。。。。。。。