Matlab

Tarea n°1 de Series Meteorológicas

Nombre: Cinthya Bravo
Profesor: Pablo Bilbao

1.
% ejercicio numero 1
ano=1957:2002; mes=1:12; tt=[];
for i=1:length(ano)
dano=diasenano(ano(i));
dia=1:dano;
for j=1:dano;
index=find(temp_850(:,1)==ano(i) & temp_850(:,3)==dia(j));
for jj=5:9
dd(jj-4)=mean(temp_850(index,jj));
endaux=[ano(i) dia(j) dd];
tt=[tt;aux];
end
end
subplot(2,3,1)
plot(tt(:,3)-273.15,'r')
title('promedios diarios desde 1/9/1957 - 31/8/2002 para longitud 300°','fontsize',9,'color','blue')
xlabel('numero de dias desde 1/9/1957-31/8/2002','fontsize',8,'color','blue')
ylabel('temperatura promedio del dia en °C','fontsize',8,'color','blue')
subplot(2,3,2)
plot(tt(:,4)-273.15,'r')title('promedios diarios desde 1/9/1957 - 31/8/2002 para longitud 302.5°','fontsize',9,'color','blue')
xlabel('numero de dias desde 1/9/1957-31/8/2002','fontsize',8,'color','blue')
ylabel('temperatura promedio del dia en °C','fontsize',8,'color','blue')
subplot(2,3,3)
plot(tt(:,5)-273.15,'r')
title('promedios diarios desde 1/9/1957 - 31/8/2002 para longitud 305°','fontsize',9,'color','blue')xlabel('numero de dias desde 1/9/1957-31/8/2002','fontsize',8,'color','blue')
ylabel('temperatura promedio del dia en °C','fontsize',8,'color','blue')
subplot(2,3,4)
plot(tt(:,6)-273.15,'r')
title('promedios diarios desde 1/9/1957 - 31/8/2002 para longitud 307°','fontsize',9,'color','blue')
xlabel('numero de dias desde 1/9/1957-31/8/2002','fontsize',8,'color','blue')
ylabel('temperatura promedio deldia en °C','fontsize',8,'color','blue')
subplot(2,3,5)
plot(tt(:,7)-273.15,'r')
title('promedios mensuales desde 1/9/1957 - 31/8/2002 para longitud 310°','fontsize',9,'color','blue')
xlabel('numero de meses desde 1/9/1957-31/8/2002','fontsize',8,'color','blue')
ylabel('temperatura promedio del dia en °C','fontsize',8,'color','blue')

2.
%ejercicio numero 2
ano=1957:2000;
mes=1:12;aux=[];
tt=[];
for i=1:length(ano)

for ii=1:length(mes)
idex=find(temp_850(:,1)==ano(i) & temp_850(:,2)==mes(ii));
for j=5:9
mm(j-4)=mean(temp_850(idex,j));
end
aux=[ano(i) mes(ii) mm];
tt=[tt;aux];
end
end
subplot(2,3,1)
plot(tt(:,3)-273.15,'g')
title('promedios mensuales desde 1/9/1957 - 31/8/2002 para longitud300°','fontsize',9,'color','blue')
xlabel('numero de meses desde 1/9/1957-31/8/2002','fontsize',8,'color','blue')
ylabel('temperatura promedio del mes en °C','fontsize',8,'color','blue')
subplot(2,3,2)
plot(tt(:,4)-273.15,'g')
title('promedios mensuales desde 1/9/1957 - 31/8/2002 para longitud 302.5°','fontsize',9,'color','blue')
xlabel('numero de meses desde 1/9/1957-31/8/2002','fontsize',8,'color','blue')ylabel('temperatura promedio del mes en °C','fontsize',8,'color','blue')
subplot(2,3,3)
plot(tt(:,5)-273.15,'g')
title('promedios mensuales desde 1/9/1957 - 31/8/2002 para longitud 305°','fontsize',9,'color','blue')
xlabel('numero de meses desde 1/9/1957-31/8/2002','fontsize',8,'color','blue')
ylabel('temperatura promedio del mes en °C','fontsize',8,'color','blue')
subplot(2,3,4)plot(tt(:,6)-273.15,'g')
title('promedios mensuales desde 1/9/1957 - 31/8/2002 para longitud 307°','fontsize',9,'color','blue')
xlabel('numero de meses desde 1/9/1957-31/8/2002','fontsize',8,'color','blue')
ylabel('temperatura promedio del mes en °C','fontsize',8,'color','blue')
subplot(2,3,5)
plot(tt(:,7)-273.15,'g')
title('promedios mensuales desde 1/9/1957 - 31/8/2002 para longitud310°','fontsize',9,'color','blue')
xlabel('numero de meses desde 1/9/1957-31/8/2002','fontsize',8,'color','blue')
ylabel('temperatura promedio del mes en °C','fontsize',8,'color','blue')

3.
>> %media de la estacion ubicada en 300°
>> nanmean(tt(:,3))

ans =

274.8469

>> %media de la estacion ubicada en 302.5°
>> nanmean(tt(:,4))

ans =

275.1768

>> %media de la estacion ubicada en...