domingo, 8 de noviembre de 2015

Álgebra

Ecuaciones

La ecuación de Price es una ecuación de covarianza que proporciona una descripción matemática de la evolución y la selección natural. La ecuación de Price fue obtenida por George R. Price, mientras trabajaba en Londres en la reelaboración de los trabajos de W. D. Hamilton sobre selección por parentesco. La ecuación de Price tiene además aplicaciones en economía.
Supóngase que hay una población de n individuos en la cual una determinada característica presenta variaciones cuantitativas. Esos nindividuos pueden agruparse en función de la cantidad de característica que cada uno presenta. En ese caso, como máximo habra ngrupos con n valores distintos de la característica, y como mínimo un sólo grupo con un solo valor compartido de esa característica. Etiquétese cada grupo con una i de manera que el número de miembros de cada grupo es n_i y el valor de la característica compartido por todos los miembros de ese grupo es z_i. Ahora asúmase que para un valor de la característica z_i viene asociado una aptitud w_idonde el producto w_i n_i representa el número de descendientes en la siguiente generación. Desígnese este número de descendientes del grupo i como n_i' de manera que w_i=n_i'/n_i. Sea z_i' la cantidad promedio de característica presentada por la descendencia del grupo i. Se designa como \Delta z_i a la variación en la característica para el grupo i, definida por:
\Delta{z_i} \ \stackrel{\mathrm{def}}{=}\  z_i' - z_i
Considérese z como la cantidad promedio del valor de la característica en esta población, y z' el promedio del valor de la característica en la siguiente generación. Defínase el cambio en la característica promedio mediante \Delta{z}. Es decir,
\Delta{z} \ \stackrel{\mathrm{def}}{=}\  z' - z
Obsérvese que esto no es el valor promedio de \Delta{z_i}. Considérese además w como la aptitud promedio de esta población. La ecuación de Price establece que:
w\,\Delta{z}=\operatorname{cov}(w_i,z_i)+\operatorname{E}(w_i\,\Delta z_i) , \,\!
donde las funciones \operatorname{E} y \operatorname{cov} son similares a las versiones muestrales de los operadores valor esperado y de la covarianza provenientes de la probabilidad. Obsérvese que ésta es en realidad una ecuación en diferencias que relaciona el valor promedio de la característica en una generación con el valor promedio de la característica en la generación venidera. De hecho, considerando un w no nulo, a menudo es útil escribirla como
\Delta{z}=\frac{\operatorname{cov}(w_i,z_i)}{w}+\frac{\operatorname{E}(w_i\,\Delta z_i)}{w}\,
En el caso específico de que la característica z_i = w_i (es decir, que la aptitud sea la característica que nos interesa), entonces la ecuación de Price reformula el teorema fundamental de Fisher de la selección natural.
La ecuación de Price es un teorema. Es la expresión de un hecho matemático entre ciertas variables, y su valor yace en el entendimiento adquirido al asignar ciertos valores encontrados en la genética evolutiva a las variables. Por ejemplo, la afirmación "si cada par de pájaros tiene dos descendientes, entonces entre diez parejas de pájaros habrán veinte descendientes" es la expresión verbal de un teorema sencillo. En realidad no transmite ninguna información nueva sobre los pájaros sino que organiza nuestros conceptos sobre los pájaros y su descendencia. La ecuación de Price es mucho más sofisticada que la anterior expresión, pero en lo primordial, es también un teorema matemáticamente demostrable.

Price aplicó la teoría de juegos a la comprensión de los comportamientos y pudo demostrar, mediante modelos matemáticos, cómo la selección natural podía actuar sobre individuos, pero también sobre genes o sobre grupos familiares. La ecuación de Price describe cómo tendencias o rasgos pasan de una generación a la siguiente, convirtiéndose en un referente en el estudio de la evolución.






La ecuación de tiempo es la diferencia entre el tiempo solar medio (medido generalmente por un reloj) y el tiempo solar aparente (tiempo medido por un reloj de sol). Esta diferencia varía a lo largo del año y alcanza una mayor diferencia a principios de noviembre, cuando el tiempo solar medio está a más de 16 minutos por detrás del tiempo solar aparente (en concreto a 16 minutos 33 segundos cerca del 3 de noviembre), y a mediados de febrero, cuando el tiempo solar medio va más de 14 minutos por delante del aparente.
Son iguales el tiempo solar medio y el tiempo solar aparente en cuatro instantes del año: el 15 de abril14 de junio,1 de septiembre y el 25 de diciembre (coinciden aproximadamente con los equinocios y solsticios). La ecuación del tiempo se representa gráficamente con un diagrama denominado analema, que suele indicarse a veces a manera de leyenda en los globos o esferas terrestres y que tiene forma de un 8 algo asimétrico. El analema indica la misma información que la expresada a través del gráfico adjunto, por lo que también a este gráfico se le podría considerar como un analema abierto.

Etimología

En la antigua astronomía la palabra 'ecuación' ('equatio' en latín) significaba 'corrección',1 y con ello se indicaba que había que sumar algebraicamente un valor para corregirlo. De esta forma, la ecuación de tiempo es el valor que hay que añadir al valor del tiempo solar aparente para 'corregirlo' (hacerlo 'regular'). Otro ejemplo es laecuación de centro de Copérnico, muy empleada en los cálculos de mecánica celeste para 'corregir' la anomalía verdadera.

Fundamentación

Observación reiterada del Sol en periodos de 24 horas de tiempo solar medio; se puede ver dibujada el analema y el concepto de ecuación del tiempo.
El origen de este concepto se deriva de la distinta velocidad del movimiento de traslación terrestre alrededor del Sol, y también de la inclinación del eje de rotación de la Tierra con respecto al plano de su órbita. (En otras palabras, aún cuando la órbita terrestre fuera perfectamente circular y la Tierra se desplazara sin variar su distancia al Sol, habría igualmente una diferencia entre el tiempo solar y el tiempo de reloj.) La órbita terrestre se denomina eclíptica(porque es en ella donde se producen los eclipses cuando la órbita de la Luna coincide en un punto con la de la Tierra y con la del Sol) y no es circular sino elíptica, ocupando el Sol uno de los focos de la elipse. De acuerdo con las leyes de movimiento orbital formuladas por Kepler sobre los movimientos de traslación, "tiempos iguales barren espacios iguales", lo cual significa que la Tierra disminuye la velocidad de traslación cuando se encuentra más alejada del Sol (porque la atracción del mismo es menor al encontrarse más lejos) y lo acelera al acercarse.
Si no existiera esta diferencia de velocidad, la Tierra se escaparía del Sistema Solar cuando se encontrara más lejos o chocaría con el Sol al acercarse. Así pues, el movimiento de traslación terrestre es un movimiento uniformemente variado. Con la Luna sucede lo mismo: cuando el Sol se encuentra más lejos de la Tierra y la Luna está más cerca, el disco lunar puede tapar por completo el disco solar (en este caso podría producirse un eclipse total de Sol) mientras que cuando sucede lo contrario (el sol más cerca y la Luna más lejos), puede producirse un eclipse anular de Sol, en el que queda un anillo luminoso del Sol alrededor de la sombra de la Luna.

Valores de la ecuación de tiempo

Los valores de la ecuación del tiempo se suelen publicar para cada año en los almanaque náuticos, en los anuarios de los observatorios, en revistas especializadas, etc. Generalmente en el apartado de efemérides solares. El motivo de esta publicación previa es proporcionar a los astrónomos la posibilidad de planificar sus observaciones. Suele representarse en forma de tabla en la que una de las entradas es el día del año y la salida es la diferencia entre tiempo medio y el verdadero (m-v), o viceversa (v-m). En algunas fórmulas empíricas dadas por los observatorios se puede averiguar de forma analítica la ecuación del tiempo. Un ejemplo es:2

E \left(d\right) = 595^{s} \sin \left(198^o + 1^o.9713 d\right)+ 441^s \sin \left(175^o - 0^o.386564 d\right)/(15.4/1.9)=1022.5541379169550677564089194965d
Donde el valor obtenido por esta fórmula semiempírica es en segundos (v-m), siendo d el día del año (del año 1999). Esta ecuación no es precisa y puede llegar a cometer errores de medio minuto como máximo.

La ecuación del tiempo

Inclinación del eje de la Tierra

En la página titulada El giróscopo estudiamos el movimiento de precesión de una peonza
Esta descripción es aproximada, siempre que la velocidad angular ω, sea grande en comparación con la velocidad angular de precesión Ω.
Un análisis más detallado del problema indica que el sólido tiene tres movimientos:
  • De rotación alrededor de su eje principal de inercia con velocidad angular ω
  • De precesión alrededor del eje vertical Z, con velocidad angular Ω
  • De nutación o de oscilación del eje vertical entre dos círculos.
No solamente el movimiento de rotación contribuye al momento angular L, como hemos supuesto, sino también y en menor medida, los movimientos de precesión y de nutación, esto es lo que hace difícil el análisis detallado de este sistema mecánico.
Los fenómenos giroscópicos tienen muchas aplicaciones: la tendencia de un giróscopo a mantener el eje de rotación fijo en el espacio en ausencia de momento es utilizado en la estabilización de los barcos y en los pilotos automáticos de los aviones.
Otro ejemplo interesante es la precesión de los equinocios. El plano del ecuador hace un ángulo de 23.45° con el plano de la órbita terrestre o eclíptica. La intersección de los dos planos es la línea de los equinoccios. La Tierra es un giróscopo gigante cuyo eje de rotación precesa alrededor del eje perpendicular al plano de la eclíptica con un periodo de 27 725 años. La precesión de los equinocios se debe al momento de las fuerzas ejercido por el Sol y la Luna sobre la Tierra. El eje de la Tierra experimenta también una nutación con una amplitud de 9.2" y un periodo de oscilación de 19 años.
A lo largo de un intervalo de tiempo de varios años (pequeño comparado con el periodo de precesión de 27725 años) podemos considerar que el eje de la Tierra permanece fijo apuntando al mismo lugar en el cielo, la estrella Polar. Pero a medida que la Tierra gira alrededor del Sol siguiendo la órbita elíptica, el ángulo de declinación formando por el plano de la eclíptica y el plano del ecuador va cambiando. En la figura, se representa el mínimo y el máximo que corresponden a las fechas 21 de diciembre y 21 de junio denominadas respectivamente, solsticio de invierno y solsticio de verano en el hemisferio norte. Ambos planos (declinación nula) coinciden en las fechas de 20 de marzo y 22 de septiembre.
El ángulo de declinación es el ángulo formado por los rayos de sol con el plano ecuatorial de la Tierra, es también la latitud a la cual los rayos de sol inciden perpendicularmente a la superficie terrestre a mediodía. El ángulo es positivo cuando el Sol está por encima del ecuador (figura de la derecha) y negativo (figura de la izquerda) cuando el Sol está por debajo del ecuador. En la siguiente dirección se puede encontrar una tabla de ángulos de declinación para cada día de año. Table of the Declination of the Sun.
La declinación solar δ se calcula mediante la siguiente fórmula
δ=0.006918-0.399912·cos(x)+0.070257·sin(x)-0.006758·cos(2x)+0.000907·sin(2x)-0.002697·cos(3x)
+0.001480·sin(3x).
x=2π365(N1+h1224)
es el número de día y h la hora. Véase, General Solar Position calculations (NOAA)
La tabla siguiente proporciona la conversión de día del mes d a número de día N en el año
MesNúmero de día, N
Enerod
Febrerod+31
Marzod+59
Abrild+90
Mayod+120
Juniod+151
Juliod+181
Agostod+212
Septiembred+243
Octubred+273
Noviembred+304
Diciembred+334
En caso de año bisiesto se añadirá uno más a partir del mes de Marzo. Por ejemplo en la casilla Marzo se pondría d+60, en la casilla Abril d+91, etc.
Mediante el siguiente script representamos la declinacion δ en grados en función de N (1-365) el día del año. Tómese h=12.
h=12;
N=1:365;
x=2*pi*(N-1-(h-12)/24)/365;
delta=(0.006918-0.399912*cos(x)+0.070257*sin(x)-0.006758*cos(2*x)
+0.000907*sin(2*x)-0.002697*cos(3*x)+0.001480*sin(3*x))*180/pi;

plot(N,delta,'r')
line([0 365],[23.45 23.45])
line([0 365],[-23.45 -23.45])
xlim([0 365])
set(gca,'XTick',[15 46 74 105 135 166 196 227 258 288 319 349])
mes = ['Ene';'Feb';'Mar';'Abr';'May';'Jun';'Jul';'Ago';'Sep';'Oct';
'Nov';'Dec'];
set(gca,'XTickLabel',mes)
grid on
ylabel('\delta grados')
title('Declinación a lo largo del año')
Como podemos apreciar en el script se ha cambiado las divisiones del eje X por defecto, sustituyéndolas por los días 15 de cada mes.
Así pues, la inclinación del eje de la Tierra respecto de la perpendicular al plano de la eclíptica de 23.45° es la responsable de la sucesión de las estaciones. La rotación de la Tierra alrededor de su eje es la responsable de la sucesión de los días y las noches.
La Tierra se mueve alrededor del Sol en una órbita elíptica cuyo periodo es 365.25 días. La distancia media entre el Sol y la Tiera es 1.496·1011 m (una unidad astonómica UA). El plano de la órbita se denomina eclíptica.
  • La distancia más alejada del Sol se denomina afelio de 152 097 525 km y ocurre el 5 de julio
  • La distancia más cercana al Sol, se denomina perihelio y ocurre alrededor del 2 de enero cuando la Tierra está a una distancia del Sol de 147 097 984 km
  • El ángulo de declinación δ=23.45° se denomina solsticio de verano y ocurre el 21 de junio a las 07h 04m (inicio del verano)
  • El ángulo de declinación δ=-23.45° se denomina solsticio de invierno y ocurre el 21 de diciembre a las 18h 11m (inicio del invierno)
  • El ángulo de declinación δ=0° se denomina equinocios y ocurren el 20 de marzo a las 12h 02m (inicio de la primavera) y 22 de septiembre a las 22h 44m (inicio del otoño)
Véase la web del Ministerio de Fomento: "Efemérides y fenomenos astronómicos relevantes en el año 2013"

Hora local y hora solar

El paso del tiempo se mide a partir de la rotación de la Tierra alrededor de su eje, una vuelta cada día. Todos las localidades que están en un meridiano tienen la misma hora. Se toma como referencia el meridiano que pasa por Greenwich (Inglaterra). El Tiempo Universal está definido dividiendo el día en 24 horas y tomando como hora cero a medianoche.

Latitud y longitud de un lugar

La posición de un punto sobre la superficie de la Tierra se mide mediante una latitud λ y una longitudφ. La latitud va de 0 en el Ecuador a ±90°, hacia el norte las latitudes son positivas y hacia el sur negativas. La longitud de 0 a ±180° se mide tomando como referencia el meridiano que pasa por Greenwich, hacia el Este las longitudes son positivas y hacia el Oeste las longitudes son negativas.
Podemos ver en Google Maps la latitud y longitud de cualquier localidad. Por ejemplo situamos el cursor sobre Eibar (País Vasco) y pulsamos el botón derecho del ratón. En el menú que aparece seleccionamos Â¿Qué hay aquí?. En el control de edición situado en la parte superior izquierda aparece la latitud y longitud, tal como apreciamos en la figura.
Dado que Eibar está en el hemisferio norte, la latitud 43.19° es positiva y como está está al Oeste del meridiano 0 o de Greenwich la longitud es -2.49, negativa.
Introducimos las coordenadas (latitud, longitud) del lugar en dicho control de edición y Google maps nos lo muestra. Se suguiere poner Google maps en baja resolución (el cursor de la barra de resolución hacia -).

Hora local

La duración variable del día se debe a las siguientes causas, solamente las dos primeras son importantes:
  1. La Tierra describe una órbita alrededor del Sol que no es circular, de modo que la Tierra tiene una velocidad mayor cuando está cerca del Sol y una velocidad menor cuando está lejos, de acuerdo con la ley de las áreas de Kepler.
  2. El eje de la Tierra está inclinado respecto de la perpendicular al plano de la eclíptica 23.45°.
  3. La velocidad angular de rotación de la Tierra alrededor de su eje va disminuyendo con el tiempo, debido al efecto de las fuerzas de marea ejercidas por la Luna.
  4. El eje de la Tierra también describe un movimiento de nutación
La Tierra se divide en zonas de anchura 15 grados tomadas desde el meridiano cero. Cada zona tiene la misma hora y se diferencia en una hora de la siguiente. Sin embargo, esta no es una regla general y como podemos ver en la figura tomada de la web Time and Date. España peninsular tiene la misma hora que Europa central, a pesar de que el oeste de Galicia está a -9° de longitud y el este de Polonia a 23.5° de longitud por lo que la anchura de esta zona horaria es de 31.5° el doble de una standard. Inglaterra, Portugal y Canarias tienen una hora menos.
Situando el cursor sobre el punto rojo que representa Madrid obtenemos la hora, por ejemplo 23:00:28. Pinchando en el punto obtenemos información suplementaria. Que Madrid está en la zona UTC/GMT+1, además aplica un adelanto horario para ahorrar energía (Daylight saving time) por lo que la hora local de Europa central CEST (Central European Summer Time) es UTC/GMT+2.
El adelanto horario, empieza el último domingo de Marzo y acaba el último domingo de Octubre. Ello viene regulado por la Directiva 2000/84/CE del Consejo de la Comunidad Europea. El horario adelantado para el año 2013 se inicia a las 2h de la madrugada (hora peninsular) del domingo 31 de marzo, momento en que los relojes se adelantarán una hora y durará hasta las 3 h de la madrugada (hora peninsular) del domingo 27 de octubre, momento en que los relojes se retrasarán una hora.

La ecuación del tiempo

En la página titulada El problema de Kepler vimos la dificultad para calcular la posición angular θ de un satélite artificial o de un planeta en el instante t. Véase Equation of Time. Problem in Astronomy. M. Muller.
El movimiento aparente del Sol no es uniforme y la duración del día solar no es constante a lo largo del año. La diferencia entre el movimiento aparente del Sol y el movimiento medio se denomina ecuación del tiempo. La expresión tomada de General Solar Position calculations (NOAA) nos da la ecuación del tiempo en minutos
Δt=229.18·(0.00075+0.001868·cos(x)-0.032077·sin(x)-0.014615·cos(2x)-0.040849·sin(2x))
Donde x se define en función del número de día N y la hora h
x=2π365(N1+h1224)
Mediante el siguiente script representamos la ecuación del tiempo en minutos en función deN (1-365) el día del año. Tómese h=12. Cambiar las divisiones del eje X por defecto, sustituyéndolas por los días 15 de cada mes.
h=12;
N=1:365;
x=2*pi*(N-1-(h-12)/24)/365;
y=229.18*(0.00075+0.001868*cos(x)-0.032077*sin(x)-0.014615*cos(2*x)
-0.040849*sin(2*x));
plot(N,y,'r')
xlim([0 365])
set(gca,'XTick',[15 46 74 105 135 166 196 227 258 288 319 349])
mes = ['Ene';'Feb';'Mar';'Abr';'May';'Jun';'Jul';'Ago';'Sep';'Oct';
'Nov';'Dec'];
set(gca,'XTickLabel',mes)
grid on
ylabel('minutos')
title('Ecuación del tiempo')

Hora local y hora solar

La hora solar se fundamenta en que la Tierra da una vuelta completa en 24 horas y que a las 12:00 h el Sol apunta exactamente hacia el sur. La hora solar verdadera se calcula mediante las dos ecuaciones siguientes:
time_offset=eqTime-4·longitude+60·timezone
  • eqTime es Δt, la ecuación del tiempo en minutos
  • longitude, es la longitud del lugar en grados (un grado de longitud equivale a 4 minutos, ó 15 grados a una hora),
  • timezone es el tiempo en horas desde UTC.
La hora solar verdadera tst, en minutos, se obtiene mediante la siguiente ecuación
tst=hr·60+min+seg/60+time_offset
donde hr es la hora (0-23), min es el minuto (0-60) y seg es el segundo (0-60)
La utilización más importante en las aplicaciones solares es la determinación de la hora local de la salida del Sol, de la puesta del Sol o del mediodía, que por otra parte, disponemos en tablas en la web del Ministerio de Fomento.
Se denomina Orto o salida del Sol al instante que corresponde a la aparición del borde superior del Sol en un horizonte hipotético en que no se considera el relieve del horizonte real, ni obstáculos cercanos (casas, árboles), ni la presencia de nubes o niebla.
El orto real en un lugar dado se puede retrasar en varios minutos si hay elevaciones del terreno en la dirección de salida del Sol, lo cual puede dar lugar a que sobre el terreno se aprecie que el Sol no ha salido pero que ya haya tenido lugar la salida de Sol en el sentido astronómico que generalmente se aplica y que hemos considerado en nuestras tablas.
Del mismo modo, se denomina Ocaso o puesta del Sol al instante que corresponde a la desaparición del borde superior del Sol en un horizonte hipotético en que no se considera el relieve del horizonte real, ni obstáculos cercanos (casas, árboles), ni la presencia de nubes o niebla.

Angulo horario

Para describir la rotación de la Tierra alrededor de su eje se utiliza el concepto de ángulo horario ω.
Como se muestra en la figura, el ángulo horario es la distancia ente el meridiano del observador P y el meridiano cuyo plano contiene al Sol. El ángulo horario es cero al mediodía (cuando el Sol alcanza la máxima altura en el cielo). El Sol apunta hacia el sur pues el plano meridiano del observador contiene al Sol. El ángulo horario se incrementa 15 grados por cada hora.
El ángulo horario es -180° a las 0h, 0° a mediodía (12h) y 180° al final del día (24h)
Por ejemplo, tres horas después del mediodía, la hora solar es 15:00 horas y el ángulo horario es de 45°. Dos horas y 20 minutos antes del mediodía el hora solar es 9:40 y ángulo solar es -35°.

No hay comentarios:

Publicar un comentario