请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 29742|回复: 33

[程序设计] matlab中如何实现滑动T检验

[复制链接]

新浪微博达人勋

发表于 2015-4-12 20:23:47 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
请问各位大师如何用MATLAB实现这个例子的滑动T检验? 无标题.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 0
发表于 2016-4-26 09:23:39 | 显示全部楼层
  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显著水平');
复制代码
密码修改失败请联系微信:mofangbao
回复 支持 6 反对 0

使用道具 举报

新浪微博达人勋

 成长值: 0
发表于 2016-4-26 09:24:28 | 显示全部楼层
剩下的自己改改咯
密码修改失败请联系微信:mofangbao
回复 支持 2 反对 0

使用道具 举报

新浪微博达人勋

 成长值: 32430
发表于 2015-4-12 20:32:53 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-4-12 20:48:21 | 显示全部楼层
二爷名声在外 发表于 2015-4-12 20:32
http://bbs.06climate.com/forum.php?mod=viewthread&tid=17633&extra=page%3D17

不好意思,我不会MATLAB,所以看不懂,请问有没有这个例子的具体操作呀?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 32430
发表于 2015-4-13 07:57:14 | 显示全部楼层
TVXQ 发表于 2015-4-12 20:48
不好意思,我不会MATLAB,所以看不懂,请问有没有这个例子的具体操作呀?

刚才你问我,我可以回答一句无可奉告,但是你又不高兴,我怎么办?
既然不会就去学啊!例子已经不错了,希望您结合原理仔细看看!再说您为什么不去问发帖人来问我,这个本身也不合适啊。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-4-13 08:21:04 | 显示全部楼层
我也要做这个
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-4-13 08:27:28 | 显示全部楼层
楼主你是想要现成答案么
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-4-14 09:08:41 | 显示全部楼层
我也要做这个
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-5-18 19:52:57 | 显示全部楼层
我也在做这个,先看看
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-7-1 17:03:50 | 显示全部楼层
这个链接的程序有问题,运行不了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-7-1 22:00:49 | 显示全部楼层
顶。。。。。。。。。
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表