Series de tiempo
Publicado por Miguel Angel (1 intervención) el 25/05/2016 18:03:05
Que tal, busco ayuda para realizar un programa para una serie de tiempo de un año, donde tengo doce archivos referentes a cada mes, cargo cada mes mediante un ciclo, y del mismo modo, me salto los archivos fallidos.
Entonces, necesito obtener vectores de 12 elementos para ese año, donde cada elemento es un promedio mensual, pero como elimino dos meses, solo tengo 10/12 datos. la otra pregunta es que realizo promedios diarios, pero no todos los meses tienen 30,31 o 28 dias, de manera errónea realizo lo siguiente:
Ojalá puedan orientarme
saludos
Entonces, necesito obtener vectores de 12 elementos para ese año, donde cada elemento es un promedio mensual, pero como elimino dos meses, solo tengo 10/12 datos. la otra pregunta es que realizo promedios diarios, pero no todos los meses tienen 30,31 o 28 dias, de manera errónea realizo lo siguiente:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
MDirvhp=[];
MRapvhp=[];
MTemphp=[];
MPreshp=[];
MDirvdp=[];
MRapvdp=[];
MPresdp=[];
MTempdp=[];
MTemppm=[];
MDirvpm=[];
MRapvpm=[];
MPrespm=[];
MHrhp=[];
MHrdp=[];
MHrpm=[];
MRhp=[];
MRdp=[];
MRpm=[];
PMT=[];
PMP=[];
PMD=[];
PMWS=[];
PMR=[];
MMaxT=[];MMinT=[];
MMaxP=[];MMinP=[];
MMaxD=[];MMinD=[];
MMaxWS=[];MMinWS=[];
MMaxR=[];MMinR=[];
MTmin=[]; MTmax=[];
MDirmin=[]; MDirmax=[];
MW_Smin=[]; MW_Smax=[];
MPmin=[]; MPmax=[];
MHrmin=[]; MHrmax=[];
MRmin=[]; MTmax=[];
for i = 201301:201312 % archivos de cada mes
if i==201310|i==201311%archivos fallidos
continue
end
archivo = [ num2str( i ) ,'pa' '.xls' ]; %lectura de datos
[num,text]=xlsread(archivo);
x=1:1:length(num(2:end,1));
Dirv=num(2:end,2);cd=2;
Rapv=num(2:end,4);cr=4;
Temp=num(2:end,6);ct=6;
Pres=num(2:end,8);cp=8;
Hr=num(2:end,7);
R=num(2:end,9);
%Promedios cada hora
i_nv=1:6:length(Dirv);
a=1:6:length(Dirv);
b=a(2)-1:6:length(Dirv);
for k=1:length(i_nv);
Dirvh=Dirv(a(k):b(k),1);
Dirvhp(k)=nanmean(Dirvh);
Rapvh=Rapv(a(k):b(k),1);
Rapvhp(k)=nanmean(Rapvh);
Temph=Temp(a(k):b(k),1);
Temphp(k)=nanmean(Temph);
Presh=Pres(a(k):b(k),1);
Preshp(k)=nanmean(Presh);
Dirvhp=[Dirvhp];
Hrh=Hr(a(k):b(k),1);
Hrhp(k)=nanmean(Hrh);
Rh=R(a(k):b(k),1);
Rhp(k)=nanmean(Rh);
Rh=R(a(k):b(k),1);
Rhp(k)=nanmean(Rh);
end
MDirvhp=[MDirvhp Dirvhp];
MRapvhp=[MRapvhp Rapvhp];
MTemphp=[MTemphp Temphp];
MPreshp=[MPreshp Preshp];
MHrhp=[MHrhp Hrhp];
MRhp=[MRhp Rhp];
%=============================================Promedios diarios y mensuales
i_nvdd=1:24:length(Dirvhp);
Dirvd=zeros(length(i_nv),1);
ad=1:24:length(Dirvhp);
bd=ad(2)-1:24:length(Dirvhp);
meses=[];
for ii=1:length(i_nvdd);
Dirvd=Dirvhp(1,ad(ii):bd(ii));
Dirvdp(ii)=nanmean(Dirvd);
Rapvd=Rapvhp(1,ad(ii):bd(ii));
Rapvdp(ii)=nanmean(Rapvd);
Tempd=Temphp(1,ad(ii):bd(ii));
Tempdp(ii)=nanmean(Tempd);
Hrd=Hrhp(1,ad(ii):bd(ii));
Hrdp(ii)=nanmean(Hrd);
Rd=Rhp(1,ad(ii):bd(ii));
Rdp(ii)=nanmean(Rd);
Presd=Preshp(1,ad(ii):bd(ii));
Presdp(ii)=nanmean(Presd);
Temppm=nanmean(Tempdp);
Dirvpm=nanmean(Dirvdp);
Rapvpm=nanmean(Rapvdp);
Prespm=nanmean(Presdp);
Hrpm=nanmean(Hrdp);
Rpm=nanmean(Rdp);
Tmin=min(Tempdp); Tmax=max(Tempdp);
Dirmin=min(Dirvdp); Dirmax=max(Dirvdp);
W_Smin=min(Rapvdp); W_Smax=max(Rapvdp);
Pmin=min(Presdp); Pmax=max(Presdp);
Hrmin=min(Hrdp); Hrmax=max(Hrdp);
Rmin=min(Rdp); Rmax=max(Rdp);
end
Ojalá puedan orientarme
saludos
Valora esta pregunta
0