爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4100|回复: 1

[程序设计] MATLAB处理nc数据T检验问题

[复制链接]

新浪微博达人勋

发表于 2014-5-3 20:34:15 | 显示全部楼层 |阅读模式

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

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

x

我在做一个臭氧对应紫外辐射高低值年(已知2001为高值年,2008为低值年)的求平均,然后进行T检验的工作。现有20012008年臭氧nc数据,其中o3
信息为 o3.png
编了求平均和T检验的程序如下:
%读取高值年2001年的数据
o3_01=ncread('o3_2001_150-1hPa.nc','o3');
level_01=ncread('o3_2001_150-1hPa.nc','level');
long_01=ncread('o3_2001_150-1hPa.nc','longitude');
lat_01=ncread('o3_2001_150-1hPa.nc','latitude');
time_01=ncread('o3_2001_150-1hPa.nc','time');
%读取低值年2008年的数据
o3_08=ncread('o3_2008_150-1hPa.nc','o3');
level_08=ncread('o3_2008_150-1hPa.nc','level');
long_08=ncread('o3_2008_150-1hPa.nc','longitude');
lat_08=ncread('o3_2008_150-1hPa.nc','latitude');
time_08=ncread('o3_2008_150-1hPa.nc','time');
x=nanmean(o3_01,4);
y=nanmean(o3_08,4);
%作图
for k=1:13
   figure;
   contour(lat_01,long_01,o3_01(:,:,k));
   contour(lat_08,long_08,o3_08(:,:,k));
end
for i=1:241
for j=1:480
x1=squeeze(x(j,i,:));
y1=squeeze(y(j,i,:));
s=ttest2(x1,y1);
ss(j,i)=s;
end
end
%T检验图
figure;
contour(ss)
然后只出了13张臭氧图,最后的T检验图上什么都没有 问题2.jpg

并且显示 问题1.jpg
而且打开变量ss后发现全是0,T检验该怎么弄?




密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-12-25 11:34:17 | 显示全部楼层
全部是零,说明都通过了检验
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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