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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 10512|回复: 18

[程序设计] EOF第一模态的方差贡献率好大啊 我用的是距平场啊

[复制链接]

新浪微博达人勋

发表于 2018-2-8 09:49:00 | 显示全部楼层 |阅读模式

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

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

x
loe=112;
low=90;
lan=34;
las=20;

nclon=ncread('D:\sst.mon.mean.nc','lon');
nclat=ncread('D:\sst.mon.mean.nc','lat');
nclon=double(nclon);
nclat=double(nclat);
real_lon=find(nclon>low & nclon<loe);
real_lat=find(nclat>las & nclat<lan);
lon_length=length(real_lon);
lat_length=length(real_lat);

sst=ncread('D:\sst.mon.mean.nc','sst',[real_lon(1) real_lat(1)...
1402],[lon_length lat_length 360],[1 1 1]);
sst=double(sst);
reshape_sst=reshape(sst,[22*14,360]);
reshape_sst(reshape_sst>1000)=nan;
reshape_sst(isnan(reshape_sst(:,1)),:)=[];
juping=reshape_sst-mean(reshape_sst,2)*ones(1,360);
x=juping;

%EOF 分析
[eof,e]=eig(x'*x);
eof=fliplr(eof);%特征向量
e=fliplr(flipud(e));%特征值
pc=x*eof;%主成分
diagonal=diag(e);
sum_e=sum(diagonal); %总特征值
eof1=pc(:,1)/sqrt(e(1,1));
e1=(diagonal(1)/sum_e)*100;
e2=(diagonal(2)/sum_e)*100;


第一第二方差贡献率

第一第二方差贡献率
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-2-8 09:59:19 | 显示全部楼层
还有一个问题,我reshape以后已经把NAN值剔除了出去,但等到我画图的时候应该怎么把NAN值再添加回去呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 0
发表于 2018-2-9 10:40:41 | 显示全部楼层
逐月资料,还是什么?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-2-9 10:48:48 | 显示全部楼层
不错 谢谢了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-2-9 10:57:47 | 显示全部楼层
言深深 发表于 2018-2-9 10:40
逐月资料,还是什么?

是逐月的海温资料
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 0
发表于 2018-2-9 11:22:07 | 显示全部楼层
liuyg123 发表于 2018-2-9 10:57
是逐月的海温资料

所以,问题来了,你是怎么计算距平的呢?整个序列一起计算,还是逐月距平的呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-2-9 12:36:42 | 显示全部楼层
有启发,谢谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-2-11 09:56:43 | 显示全部楼层
言深深 发表于 2018-2-9 11:22
所以,问题来了,你是怎么计算距平的呢?整个序列一起计算,还是逐月距平的呢?

明白了!我应该求逐月的距平 谢谢深深大大!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-2-11 09:57:46 | 显示全部楼层
言深深 发表于 2018-2-9 11:22
所以,问题来了,你是怎么计算距平的呢?整个序列一起计算,还是逐月距平的呢?

另外深深大大能否再解答一下我第二个问题,这个问题也困惑了我好久。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-2-11 09:58:18 | 显示全部楼层

我也是小白 大家一起学习
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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