clear clc M{1}= 'Enero'; M{2}= 'Febrero'; M{3}= 'Marzo'; M{4}= 'Abril'; M{5}= 'Mayo'; M{6}= 'Junio'; M{7}= 'Julio'; M{8}= 'Agosto'; M{9}= 'Septiembre'; M{10}= 'Octubre'; M{11}= 'Noviembre'; M{12}= 'Diciembre'; mes1=3; dia1=25; dia2=21; mes2=5; dia3=23; mes3=10;dia4=28; mes4=11;dia5=12; dm = [ 31 29 31 30 31 30 31 31 30 31 30 31 ]'; raiz = '\Users\Nazaret\Desktop\PFC\'; nh = 24; HT=[]; XT=[]; GEN={'ECT3' 'ELC1','ESC6' 'ESCCC1' 'ESCCC2' 'ESCCC3' 'FOI1' 'GUA1' 'GUA2' 'LAD3' 'LAD4' 'LARES1' 'LARES2' 'LIT1' 'LIT2' 'MALA1' 'MEI1' 'NRC1' 'NRC2' 'NRC3' 'PALOS1' 'PALOS2' 'PALOS3'}; %GEN={'ABENVD1','ABO1','ABO2','ABRE011','ACAVADO','ACC1EBR','ACC2EBR','ACE3','ACE4','ADOUINT','ADOUNAC','AGENVD1','AGENVD2','AGENVD3','AGENVD4','AGUG','ALG3','ALIMA','ALL1','ALZ1','ALZ2','AMBIETA','AMEVD1','AMRE018','AMRE020','AMRE033','AMRE034','AMRE035','ARCOS1','ARCOS2','ARCOS3','ARRU1R','ARRU2R','ASC1','ASC2','ASIVD01','BES3','BES4','BES5','BRR1','CAMGI10','CAMGI20','CCO2','CCO3','CEGRB1','CEGRB2','CEGRB4','CEGVD1','CENERE1','CENERE2','CENERE4','CENERE5','CENEVD1','CENTVP3','CERE003','CERE004','CERE009','CERE022','CERE027','CERE028','CERE030','CERE031','CERE037','CERE039','CERE043','CERE064','CERE066','CERE067','CERE070','CILCOR','CNDRVP3','COF1','COG75RQ','COGEVD1','COGPLA','COL4','COM4','COM5','COPRAT2','CTGN1','CTGN2','CTGN3','CTJON1R','CTJON2R','CTJON3R','CTN3','CTN4','CTNJUC2','CTNU','DANSV03','DANSVP3','DBECV03','DBECVP3','DETCRB1','DETCRE1','DETIVD1','DEVAVD1','DUER','EBLCVD1','EBRA','EBRFEN','ECOIVD1','ECT2','ECT3','ECYRRE1','ECYRRE2','ECYRVD1','EDFTVP3','EDPSVD1','EE21VD1','EE21VD2','EFGVD1','EFRE003','EGEDRE1','EGEDRE2','EGEDRE3','EGEDRE4','EGEDRE5','EGEDRE6','EGEDVD1','EGLEV03','EGLEVD1','EGLEVP3','EGRE001','EGRE002','EGRE006','EGRE007','EGRE008','EGRE009','EGRE011','EGRE012','EGRE013','EGRE014','EGRE015','EGRE016','EGRE019','EGRE020','EGRE022','EGRE023','EGRE024','EGRE025','EGRE026','EGRE027','EGRE028','EGRE030','EGRE032','EGRE033','EGRE034','EGRE036','EGRE037','EGRE038','EGRE039','EGRE040','EGRE045','EGRE047','EGRE049','EGRE050','EGRE052','EGRE053','EGRE054','EGRE056','EGRE057','EGRE058','EGRE060','EGRE061','EGRE062','EGRE063','EGRE064','EGRE065','EGRE066','EGRE067','EGRE068','EGRE070','EGRE072','EGRE073','EGRE075','EGRE076','EGRE087','EGRE092','EGRE098','EGRE099','EGRE101','EGRE102','EGRE103','EGRE106','EGRE107','EGRE108','EGRE109','EGRE110','EGRE111','EGRE112','EGRE113','EGRE114','EGRE116','EGRE117','EGRE118','EGRE126','EGRE128','EGRE130','EGRE131','EGRE132','EGRE133','EGRE134','EGRE135','EGRE136','EGRE138','EGRE139','EGRE140','EGRE141','EGRE142','EGRE143','EGRE144','EGRE147','EGRE149','EGRE150','EGRE151','EGRE152','EGRE153','EGRE154','EGRE155','EGRE156','EGRE158','EGRE159','EGRE160','EGRE161','EGRE162','EGRE163','EGRE164','EGRE165','EGRE166','EGRE167','EGRE169','EGRE170','EGRE172','EGRE174','EGRE175','EGRE177','EGRE178','EGRE179','EGRE188','EGRE189','EGRE194','EGRE195','EGRE196','EGRE210','EGRE217','EGRE219','EGRE221','EGRE228','EGRE229','EGRE230','EGRE231','EGRE233','EGRE234','EGRE238','EGRE249','EGRE250','EGRE251','EGRE252','EGRE253','EGRE254','EGRE257','EGRE259','EGRE262','EGRE263','EGRE264','EGRE265','EGRE266','EGRE267','EGRE268','EGRE269','EGRE270','ELC1','ELYOVD1','EMORPER','ENDGVD1','ENECVD1','EONUVD1','EONVD01','EONVD02','EONVP3','ESC6','ESCCC1','ESCCC2','ESCCC3','FENVD01','FOI1','GACEVD1','GACVD1','GACVD2','GAR1','GARE009','GARE010','GARE011','GARE012','GARE013','GARE014','GARE015','GARE016','GARE017','GARE018','GARE019','GARE020','GARE021','GARE022','GARE023','GARE024','GARE025','GARE026','GARE027','GARE028','GARE029','GARE030','GARE031','GARE032','GARE033','GARE034','GARE038','GARE040','GARE041','GARE042','GARE043','GARE044','GARE045','GARE046','GARE047','GARE048','GARE049','GARE050','GARE051','GARE200','GASNRE2','GASNRE3','GASNRE4','GASNRE5','GASNRE6','GASNRE7','GASNRE8','GASNRE9','GASNVD1','GASNVD2','GASNVD3','GASNVD4','GASNVD5','GASNVD6','GASNVD7','GASNVD8','GASRE10','GASRE11','GASRE12','GASRE13','GASRE14','GDLQ','GDNA','GEDFTV1','GELTBV1','GENEVD1','GEPLAS','GEREVD1','GESTVD2','GESTVD3','GESTVD4','GHISPV1','GNRAVD1','GNSUVD1','GNSUVD2','GSEVD01','GUA1','GUA2','GUIG','HCENVD1','HCERE01','HCERE02','HCERE20','HCERE21','HCERE22','HCERE23','HCERE24','HCERE26','HCERE27','HCERE28','HCERE29','HCERE30','HCERE31','HCERE32','HCERE33','HCERE34','HCERE35','HCERE36','HCERE37','HCERE38','HCERE45','HCERE47','HCERE48','HCERE49','HCERE50','HCERE51','HCERE52','HCERE53','HCHI','HCURVD1','HELIVD1','HISVD01','HISVD02','IBEGVD1','IBEGVD2','IBEGVD4','IBEGVD7','IBEGVD8','IBEGVD9','IBMON','IBRE001','IBRE002','IBRE003','IBRE004','IBRE005','IBRE006','IBRE007','IBRE008','IBRE009','IBRE010','IBRE012','IBRE018','IBRE019','IBRE020','IBRE021','IBRE022','IBRE023','IBRE024','IBRE026','IBRE028','IBRE029','IBRE030','IBRE031','IBRE032','IBRE034','IBRE035','ICUVD02','ICUVD03','ICUVD04','ICUVD05','ICUVD06','ICUVD08','ICUVD09','ICUVD14','INTVD01','IPG','ISRE005','ISRE006','ISVD001','ISVD002','ISVD003','ISVD004','ISVD005','ISVD006','ISVD007','ISVD008','ISVD009','JUCA','LAD3','LAD4','LARES1','LARES2','LIT1','LIT2','LMERVD1','LMERVD2','LMERVD3','LMERVD4','LMRE027','MALA1','MCDRE1','MCDVD1','MCDVD2','MEI1','MEREVD1','MEREVD3','MLTG','MOREL1R','MOREL2R','MUEL','NEXVD02','NEXVD03','NEXVD05','NEXVD07','NEXVD20','NEXVD21','NEXVD22','NEXVD24','NEXVD25','NEXVD26','NEXVD28','NEXVD30','NEXVD31','NEXVD32','NEXVD33','NEXVD34','NEXVD35','NEXVD36','NRC1','NRC2','NRC3','NRENVD1','NXRE004','NXRE008','NXRE009','NXRE012','NXRE013','NXRE016','NXRE024','NXRE029','NXRE030','NXRE031','NXRE042','NXRE047','NXRE062','NXRE063','NXRE064','NXRE065','NXRE066','NXRE067','NXRE069','NXRE070','NXRE071','NXRE072','NXRE073','NXRE074','PALOS1','PALOS2','PALOS3','PAS1','PBCN1','PBCN2','PEGO3','PEGO4','PGR1','PGR2','PGR3','PGR4','PGR5','PLL1','PVENT1','PVENT2','RABID2R','REPVD01','REPVD02','RIBATE1','RIBATE2','RIBATE3','ROB1','ROB2','RPG01','RPG02','RPTRTAR','RTG01','RTG02','RTG03','SAGU1','SAGU2','SAGU3','SAMPVD1','SAMPVD2','SASRUBR','SBEU','SBO3','SEDARE1','SEDAVD1','SETUBA1','SETUBA2','SETUBA3','SETUBA4','SHEEVP3','SIL','SINES1','SINES2','SINES3','SINES4','SLTG','SMGVP3','SNICOG','SNTEVD1','SRI2','SRI3','SRI4R','SRI5R','SROQ1','SROQ2','STC4R','TAJO','TAPOWER','TARRAG','TEES','TEMON','TER1','TER2','TER3','TERE','TERVD01','TJEG','TRL1','UFBG','UFGC','UFMI','UFTA','VAN2','VIES','WMARRE1','WMARVD1','WMRE001','WMRE003','WMRE004','WMRE011','WMRE023','WMRE026','WMRE027','WMRE028','WMRE030','WMRE034','WMRE035','WMRE045','WMRE046','WMRE048','WMRE051','WMRE054','WMRE055','WMRE064','WMRE074','WMRE076','WMRE078','WMRE082','WMRE083','WMRE085','WMRE094','WMRE099','WMRE102','WMRE103','WMRE105','WMRE107','WMRE108','WMRE109','WMVD002','WMVD003','WMVD004','WMVD005','WMVD006','WMVD007','WMVD008','WMVD009','WMVD010','WMVD011','WMVD013','WMVD014','WMVD015','WMVD018','WMVD020','WMVD022','WMVD024','WMVD025','WMVD026','WMVD030','WMVD031','WMVD034','WMVD035','WMVD036','WMVD037','WMVD038','WMVD039','WMVD040','WMVD042','WMVD043','WMVD044','WMVD045','WMVD046','WMVD047','WMVD049','WMVD050','WMVD052','WMVD057','WMVD059','WMVD060','WMVD061','WMVD063','WMVD064','WMVD065','WMVD066','WMVD067','WMVD068','WMVD069','WMVD070','WMVD071','WMVD072','WMVD073','WMVD075','WMVD076','WMVD080','WMVD081','WMVD082','WMVD083','WMVD084','WMVD085','WMVD086','WMVD087','WMVD088','WMVD090','WMVD091','WMVD093'}; [x y]=size (GEN); for g=1:y ngc=GEN{g}; Precios=[]; for m = 1:12 if m < 10 nm = [ '0' num2str(m) ]; else nm = num2str(m); end for d = 1:dm(m); %Salto de los dias que no tenemos---------------------------- if (m==mes1 && d==dia1) || (m==mes1 && d==dia2) || (m==mes2 && d==dia3) || (m==mes3 && d==dia4)|| (m==mes4 && d==dia5) pause(3) continue; end %----------------------------------------------------------- carp_fichero = [ 'DATOS\' num2str(m) ' ' M{m} ' 2012\' num2str(d) '\' ]; fvc = [ raiz carp_fichero 'Ventas_Cas.txt' ]; fvcs = [ raiz carp_fichero 'Ventas_CasS.txt' ]; fr=[ raiz carp_fichero 'Resumen.txt' ]; frs=[ raiz carp_fichero 'ResumenS.txt' ]; pause(0.2) flistado=[raiz 'precios_generadores.txt']; fid=fopen(flistado,'a'); if d < 10 nd = [ '0' num2str(d) ]; else nd = num2str(d); end % Lee datos fichero casado fid1 = fopen(fvc,'r'); fid2 = fopen(fvcs,'r'); [ RE, RP ] = lee_fichero(fid1,nh); [ RES, RPS ] = lee_fichero(fid2,nh); [x0 y0]=size(RPS); % Lee fichero entrada(Ingresos mínimos) fe = [ raiz 'ENTRADAS\VENTA_2012' nm nd '.1.txt' ]; fid3 = fopen(fe,'r'); [IM nombre ng] = lee_entrada(fid3,nh); for i=1:x0 Ngen=nombre{i}; [x1 y1]=size(Ngen); [x2 y2]=size(ngc); if y1==y2 if ngc == Ngen Precios=[Precios;RPS(i,:)']; %fprintf(fid,'%s',nombre{i}); %fprintf(fid,'\t'); %fprintf(fid,'%6.2f\t',RPS(i,:)); %fprintf(fid,'\n'); end end end end fclose all end [H,X]=hist(Precios,10); figure(g) bar(X,H) title(GEN{g}) %save 'ACE3.txt' Precios -ASCII HT = [ HT H' ]; XT = [ XT X' ]; fclose all end