Tema: brújula solar

Como ahora casi todo los móviles vienen con internet, se me ocurrió hacer una aplicación para saber donde esta el sur utilizando el sol; así me ahorro tener que llevar una brújula y saco el sur verdadero, no el magnético. El enlace esta en Brujula solar

Muy sencillo pero simpático, me gusta.

Muy bueno rodri, muy bueno (tengo que cubrir 20 caracteres)

Iniciado por wenner

Muy bueno rodri, muy bueno (tengo que cubrir 20 caracteres)

gracias, pero tiene muchos fallos que ire limando.
por cierto alguien sabe como se puede obtener automáticamente la diferencia horaria de un punto respecto a la hora GMT?

Iniciado por rodri

por cierto alguien sabe como se puede obtener automáticamente la diferencia horaria de un punto respecto a la hora GMT?

La oficial o la solar?

la oficial, debe haber alguna base de datos pero no la he encontrado
la solar se saca facil 15º=1hora

Iniciado por rodri

la oficial, debe haber alguna base de datos pero no la he encontrado

Leyendo de base de datos es posible, aunque desconozco si existe y dónde. Por coordenadas del lugar es más complicado puesto que se suelen seguir fronteras políticas que nada tienen que ver con esas coordenadas.

la solar se saca facil 15º=1hora

Más la ecuación del tiempo.

Iniciado por Breitling

Más la ecuación del tiempo.

edito:es correcto

estoy modificando el programa, para la declinación solar utilizo la formula de Spencer

$declinacion = 0.006918 - (0.3999 * cos($ang)) + (0.070257 * sin($ang)) - (0.006758 * cos(2 * $ang)) + (0.000907 * sin(2 * $ang)) - (0.002697 * cos(3 * $ang)) + (0.00148 * sin(3 * $ang))

donde $ang = (2 * $pi * ($dia-1)/365)


sin embargo el valor de esta declinación es diferente del que se puede obtener en almanaques, creo que es por ser la formula "antigua" de 1971 y al pasar los años el error se ha hecho mas grande, conocéis alguna formula para la declinación solar que este mas actualizada?

Prueba esta:

d = arcsin(0.3978 * sin(j’ - 1.4 + 0.0355 * sin(j’ - 0.0489)))

Donde:

j' = 2 pi j / 365,25
j = día del año

A día de hoy me da:

-0.221956597 radianes
-12.7171762 grados

Iniciado por rodri

estoy modificando el programa, para la declinación solar utilizo la formula de Spencer

$declinacion = 0.006918 - (0.399912 * cos($ang)) + (0.070257 * sin($ang)) - (0.006758 * cos(2 * $ang)) + (0.000907 * sin(2 * $ang)) - (0.002697 * cos(3 * $ang)) + (0.00148 * sin(3 * $ang))

donde $ang = (2 * $pi * ($dia-1)/365)

Tengo la misma, el error máximo se estima en 0,0006 radianes

gracias por la formula,

las ecuaciones que utilizo son:

$HS=$hora-$USO-$AD-($et/60)+($LON/15);

$angulohorario=($HS-12)*15*$pi/180;

$altura =asin( (sin($declinacion) * sin(deg2rad($LAT))) + (cos(deg2rad($LAT)) * cos($declinacion) * cos($angulohorario)));

$azimut=asin(sin($angulohorario)* cos($declinacion)/cos($altura));

$LON es la longitud
$LAT es la latitud
$USO es el uso horario
$AD es el adelanto si estamos en horario de invierno o de verano.
$et es la ecuación de tiempo

Con estas ecuaciones estoy viendo que uno de los problemas a la hora de resolver las la posición solar respeto a la hora legal es el tema de que los valores de azimut, vienen dados por una formula que incluye un arcoseno, y el ordenador siempre toma las soluciones con valores entre 0 y 90º.
Esto hace que para días de verano a primeras y ultimas horas (cuando el sol esta detras del este y del oeste) el ordenador interprete que el sol esta en valores entre -90 y 90.
Gráficamente se puede ver en este gráfico que representa la posición del sol durante las horas legales en punto, y puede verse que este sistema no distingue valores mayores por detrás del este y el oeste :

Como resolvéis esto vosotros?
Saludos

Realmente el azimut del sol es su ángulo horario, pero como está referido al cénit tiene valores negativos por la mañana y positivos por la tarde. Yo lo que hago es sumar 180º y referirlo así al antimeridiano, con lo que la lectura es directa y siempre positiva. Me va bien:

Iniciado por Breitling

Realmente el azimut del sol es su ángulo horario, pero como está referido al cénit tiene valores negativos por la mañana y positivos por la tarde. Yo lo que hago es sumar 180º y referirlo así al antimeridiano, con lo que la lectura es directa y siempre positiva. Me va bien:

hola breitling lo he hecho como dices, suponiendo el azimut como el angulo horario. lo de los 180 o cero en realidad de igual. Sin embargo en verano la gráfica es muy triangular. a ti te queda asi también?
saludos
[/QUOTE]

Te refieres a una curva más picuda, no? Sí, sale así, es por efecto de proyectar sobre plano una curva que se dibuja en una esfera. De hecho tus curvas de radiación (creo que las mides en tu instalación, verdad?) te saldrán más picudas en verano y más aplanadas en invierno. Además ahora ya ves que tus analemas salen verticales, como debe de ser, adelantando o atrasando la hora solar en cada meridiano con la ecuación del tiempo.

Yo pongo 180 porque así la gráfica me da rumbos reales, no tengo que pensar si son 30 más o 27 menos que el sur.

precisamente por eso lo digo, no son tan picudas, ademas en los pct del idae, en el censolar etc siempre he visto las curvas estas como mas aplanadas por eso como que extraña..

mira con las tools del pvsyst las curvas son mas aplanadas, y las analemas no son planas..

Que os parecen estos recorridos?

Iniciado por rodri

precisamente por eso lo digo, no son tan picudas, ademas en los pct del idae, en el censolar etc siempre he visto las curvas estas como mas aplanadas por eso como que extraña..

Lo miro mañana despacio y te cuento.

Iniciado por Titugorlov

Que os parecen estos recorridos?

muy buenos, que fórmulas has utilizado?

Muy buenas,

Para la declinación solar he utilizado la siguiente expresión:

decl = 23,45º·[sen((360/365)·(dn+284))]

donde dn es el dia del año.

Cuando tenga un momentito te pongo el resto de fórmulas.

Salut.

Sigo con el tema, las formulas "clasicas" que se pueden encontrar en internet fallan cuando el sol esta detras del este u oeste, o cuando en latitudes tropicales el ecuador "supera" los 90º de altura solar.
Aplicando alguna correción se eliminan parte de los fallos de las formulas clasicas. Sin embargo sigo teniendo algún problema en los puntos de las 5 horas solares centrales que tienen un azimut mas alla del este u oeste.

Rodri, para la elevación ($altura) utilizo la misma equación, pero para la orientación ($azimut) he utilizado esta:

orientación= acos((sin ELE·sin LAT-sin DECL)/(cos ELE·cos LAT)

Os dejo aqui la hoja de cálculo del seguidor solar, aunque es un poco cutre y hay cosas no muy bien echas, disfruté haciéndolo y sacié un poco la curiosidad.

Salut!

Iniciado por Titugorlov

Rodri, para la elevación ($altura) utilizo la misma equación, pero para la orientación ($azimut) he utilizado esta:

orientación= acos((sin ELE·sin LAT-sin DECL)/(cos ELE·cos LAT)

Os dejo aqui la hoja de cálculo del seguidor solar, aunque es un poco cutre y hay cosas no muy bien echas, disfruté haciéndolo y sacié un poco la curiosidad.

Salut!

con esa formula va bastante mejor!!

gracias!!