modconcprec.m

麻省理工学院的人工智能工具箱,很珍贵,希望对大家有用!

load('..//AreaPatterns/ECMWF1/CuencaNorte/CP')
load('..//AreaPatterns/ECMWF1/CuencaNorte/SV')
CP=CP(:,1:80);dato=CP*SV(:,1:size(CP,2))';clear CP SV
VV=[];noC=15;
[VV,clases]=kmeans(dato,25);

par=5;tim=3;lvl=3; s1=19;s2=19*3;s3=19*9;
inD(1:par,1:tim,1:lvl,1:s1)=0;
for par=1:5
   for tim=1:3
      for lvl=1:3
         ini=(par-1)*s3+(tim-1)*s2+(lvl-1)*s1;
         inD(par,tim,lvl,:)=[ini+1:ini+s1];
      end
   end
end

%seleccion de campos
for par=1:5
   for lvl=1:3
      lista=squeeze(inD(par,:,lvl,:));
      V=dato(:,lista);
      [c,clases]=kmeans(V,noC);
      VV=[VV c];
   end
end

load ('..//..//RDATA/Pre0707R/RData1')
load ('..//..//RDATA/Pre0707R/RMaes1')
noFM=8;%10 clases de precipitaci髇
%se toman los datos de la cuenca de Cantabria
lista=find(MM(:,2)==num2str(1));
%se toman los datos de Parayas
%lista=120;
MM=MM(lista,:);DD=DD(:,lista);FM=0;
while max(FM)<noFM
   [FM,centers]=kmeans(pstd(DD),noFM);
end

dato=[VV(1:length(FM),:) FM];

labels={'EA','Z1000','Z850','Z500','T1000','T850','T500','U1000','U850','U500','V1000','V850','V500','H1000','H850','H500','FM'};

novar=size(dato,2);
ns = max(dato);order=1:novar;
N=abs(corrcoef(dato));
pot_par=novar-1;
max_par=3;lis=[];
[a,b]=sort(N);lis(order,:)=b(end-pot_par+1:end,:)';

%se aprende la estructura con Local-KB
[dag,score]= L_K2(dato',ns,lis,pot_par,order,max_par);
%[dag,score]= L_B(dato',ns,lis,pot_par,max_par);

draw(dag',labels);