Cálculos con R

Con ayuda del programa de cálculo R obtuvimos los outputs
que se solicitaban en el ejercicio. Para ello usamos el código que viene a
continuación y representamos los resultados en las gráficas correspondientes.

Funcion=function(fi,dn){
##fi latitud en grados, dn dia del año
dec=23.45*sin((2*pi*(dn+284))/365)
##declinacion en grados
ws=-acos(-tan(dec*pi/180)*tan(fi*pi/180))     ##hora del amanecer en radianes

IRad=-24/pi*1367*(1+0.033*cos(2*pi*dn/365))*(ws*sin(fi*pi/180)*sin(dec*pi/180)+cos(dec*pi/180)*cos(fi*pi/180)*sin(ws))
##Irradiancia extraatmosférica en el plano horizontal W/m2

wsg=ws*180/pi   #hora del amanecer en grados

data.frame(Declinacion=dec, Amanecer=wsg, Bo0=IRad)}

Funcion(40,240)

Funcion(40,1:365)

dn=(1:365)

dat=Funcion(40,dn)

library(lattice)

xyplot(Declinacion~dn, data=dat)

xyplot(Amanecer~dn, data=dat)

xyplot(Bo0~dn, data=dat)

Grupo 5: Elena Poza Santana; David García González; Rubén Pardo Hernangómez; Alexandru N. Cosor

Comentarios desactivados en Cálculos con R / 24 Ene 2012 por Alex Cosor
MERME /

Comentarios cerrados.