matlab中如何实现滑动T检验

vangkang

同求啊~~~！

怡馨

顶。。。。。。。。。。。。

高悬的瓦片

自己写个脚本吧

liuwenyan150825

正在做，正好看看

no_passion

我也需要用这个，先马了看

abc111

同问啊。。。。。。。。。。。。。。。。

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

剩下的自己改改咯

sha

好问题。 。 。同求!

藕茭白空心菜

jiajian0512 发表于 2016-4-26 09:23

给力呀...~~~~