- 积分
- 7779
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-9-12
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
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;
|
-
第一第二方差贡献率
|