correlacion_ea.m

来自「麻省理工学院的人工智能工具箱,很珍贵,希望对大家有用!」· M 代码 · 共 46 行

close all
clear all
noclas=9;
load('..//AreaPatterns/ECMWF1/Nacional10\CP');
CP=CP(:,1:25);
c=kmeans(CP,noclas);
load ('..//AreaStations/completas/RdataPp')
load ('..//AreaStations/completas/RMaesPp')
data=RdataPp;
data(find(data==-999))=NaN;lista=sum(isnan(data));
data=data(1:size(CP,1),:);
data=data(:,1:2:end);Maes=Maes(1:2:end,:);cucu=[]
%data=data(:,1:14);Maes=Maes(1:14,:);cucu=[]
for i=0:9
   cucu=[cucu;max(find(Maes(:,1)==num2str(i)))+1];
end

%global
A=[];
[RK,MN,DV]=nanpstd(data);
[Cx,Cor] = rnancov(RK);
A=[A Cor(:,3)];
round(Cor*100);
pcolor(Cor)
colormap(flipud(gray))
figure
%para cada una de las clases
for i=1:noclas
   [RK,MN,DV]=nanpstd(data(find(c==i),:));
   [Cx,Cor] = rnancov(RK);
	A=[A Cor(:,3)];
   round(Cor*100);
   subplot(3,3,i)
   pcolor(Cor)
   colormap(flipud(gray))
end

round(A*100)


%break
%RK=fgauss(RK,Cx);
%recuperacion de los datos
%RK=RK.*repmat(DV,size(RK,1),1)+repmat(MN,size(RK,1),1);
%RK(find(RK<0))=0;data=round(RK);