Teorema de Gauss-Márkov, formulado por Carl Friedrich Gauss y Andréi Márkov, establece que en un modelo linealgeneral (MLG) en el que se establezcan los siguientes supuestos:
- Correcta especificación: el MLG ha de ser una combinación lineal de los parámetros () y no necesariamente de las variables:
- Muestreo aleatorio simple: la muestra de observaciones del vector es una muestra aleatoria simple y, por lo tanto, el vector es independiente del vector
- Esperanza condicionada de las perturbaciones nula:
- Correcta identificación: la matriz de regresoras (X) ha de tener rango completo: rg(X)=K<=N
- Homocedasticidad:
el estimador mínimo cuadrático ordinario (MCO) de B es el estimador lineal e insesgado óptimo (ELIO o BLUE: best linear unbiased estimator), es decir, el estimador MCO es el estimador eficiente dentro de la clase de estimadores lineales e insesgados.
Dicho teorema se basa en 10 supuestos, denominados, Supuestos de Gauss Márkov; que sirven como hipótesis a la demostración del mismo:
- El modelo esta correctamente especificado.
- Debe ser lineal en los parámetros.
- El valor de la media condicional es cero.
- Hay homocedasticidad.
- No existe correlación entre las perturbaciones.
- La covarianza entre y es cero.
- El número de observaciones es mayor que el de parámetros.
- Existe variabilidad entre los .
- No hay multicolinealidad perfecta.
- Las son no estocásticas, es decir, son fijas en muestras repetidas.
teorema de Gelfond-Schneider es un resultado que establece la trascendencia de una gran clase de números. Fue probado originalmente por Alexander Gelfond en 1934 y de nuevo de forma independiente por Theodor Schneider, en 1935. El teorema Gelfond–Schneider es una respuesta parcial al séptimo problema de Hilbert.Si y son números algebraicos en el cuerpo de los números complejos (siendo ), y si no es un número racional, entonces cualquier valor de αβ es un número trascendente.
- En general, es multivaluada, donde "log" es el logaritmo complejo. Ésta es la razón de la expresión "cualquier valor de" en el enunciado.
- La siguiente es una formulación equivalente del teorema: si y son números algebraicos diferentes de cero, y , entonces es (real) racional o trascendente.
- Si se elimina la restricción de que sea algebraica, el enunciado no será cierto en el caso general (escójanse y , que es trascendente, y , que es algebraico). No se conoce una caracterización de los valores de α y β que produzca un αβ trascendente.
Se deriva inmediatamente del teorema la trascendencia de los siguientes números:
- Los números (la constante de Gelfond-Schneider) y (véase demostración no constructiva).
- El número (constante de Gelfond), dado que es uno de los valores de .
teorema de Gerschgorin es utilizado en álgebra lineal para encontrar una cota de los autovalores de una matriz compleja (esto incluye también a las reales) de orden nxn. Fue publicado por el matemático soviético S.Gerschgorin en 1931. 1
Dada una matriz se definen los círculos con centro en y radio . El teorema afirma que los autovalores de la matriz A se encuentran en la unión de los n círculos. Además, cada componente conexa de esa unión contiene tantos autovalores como círculos haya en ella, donde círculos y autovalores son contados con multiplicidad.
Teorema de Gerschgorin, condición mecánica y forma final de la relación de dispersión.
Ahora analizamos una forma, para encontrar el rango de eigenvalores de una matriz y una condición mecánica para el modelo empleado. Supongamos que es una matriz arbitraria y construyamos el problema de eigenvalores siguiente en forma matricial, esto es:
nos interesa una componente, que explícitamente será:
Pasando al lado derecho el término resulta:
Como siempre existe por lo menos un eigenvector distinto de cero supongamos que es el , por lo tanto la ecuación anterior resulta:
Supongamos además, que la componente -ésima de este eigenvector tiene el módulo máximo, llamémosla por lo tanto resulta, tomando además el valor absoluto en ambos miembros que:
Como es posible suponer que por lo tanto resulta finalmente que:
Este resultado es conocido como ``teorema de Gerschgorin'' (Ref. [12]) que nos dice: Que en el plano complejo cada eigenvalor de la matriz de elementos complejos , se encuentra localizado en al menos uno de los círculos con centro y radio (que es la suma de los elementos de un renglón).
Como no sabemos cual de las componentes es la máxima, entonces se construyen los discos para cada una de las componentes y si no es el disco que interesa, no lo tomamos en cuenta, pero si lo es, entonces va a incluir el eigenvalor en cuestión; aunque es posible que al cambiar el eigenvalor cambie el elemento donde se alcanza el máximo, pero no importa, ya que se tienen todos los discos y en algunos de ellos se encuentra el eigenvalor.
Por lo tanto lo anterior nos da una información de los eigenvalores en particular, en el caso que estamos tratando, nos dice que se tienen un círculo con centro y radio . Naturalmente que para nuestro problema únicamente nos interesa la línea real, ya que se trata de una matriz hermitiana con eigenvalores reales, o sea que el intervalo de eigenvalores para el modelo que estamos analizando es:
Ahora en el caso de vibraciones mecánicas, se tiene la siguiente restricción . Esto quiere decir que la suma de los elementos a través de un renglón tiene que ser cero. Esta condición es debida a que un desplazamiento, en una dirección produce una fuerza la cual produce un desplazamiento en otra dirección y debido a que el efecto de estas fuerzas es producido por resortes que se encuentran a ambos lados de las partículas, entonces son los resortes los que dan lugar a las constantes que intervienen en la condición mecánica. Esta condición mecánica se puede visualizar mejor, por medio de un diagrama de fuerzas, esto es
Como vemos en este diagrama de fuerzas, si movemos una partícula desde su posición de equilibrio, el efecto de la interacción en ambos lados debe ser el mismo debido a la simetría de la cadena, por lo tanto va a interaccionar con las partículas en ambos lados y es precisamente la constante del resorte la que influye en el equilibrio de fuerzas.
La restricción es gráficamente razonable de postular, únicamente para nuestro caso que es un sistema vibracional, ya que en otros sistemas no es válida, como por ejemplo en problemas de mecánica cuántica, esto es sistema de electrones .
Además si no fuera físicamente razonable postular que y la empleáramos, el resultado sería únicamente una traslación de todos los eigenvalores, como lo podemos ver del siguiente hecho.
Tenemos la restricción en la forma:
y del teorema de Gershgorin dado por la relación
(4.4) |
lo cual nos permite escribir esto como
Sustituyendo el valor de resulta finalmente:
Analizando este resultado, vemos una justificación más de que las frecuencias resultan imaginarias, con la posibilidad de que sean cero, pero nunca tendremos el caso . Otro resultado que nos da es la justificación del denominador de la relación (4.2).
Por lo tanto de (4.5) y (4.2) resulta:
Al menos en el rango en que y ya que debido a la estabilidad de las redes, puede suceder que o que sea negativa también y puede darse el caso en que estas dos constantes sean iguales (como se verá al analizar modelos especiales), por lo tanto es posible obtener interacciones negativas.
Ahora lo que nos interesa es calcular los eigenvalores de la matriz en términos de los parámetros (2.13), esto es y por lo tanto seguimos tratando una cadena homogénea, en este caso la ecuación característica de , como vimos anteriormente (Ec. 2.21) es
también vimos que las raíces ocurren en parejas recíprocas, supongamos que son: por lo tanto podemos definir el siguiente arreglo:
y debido precisamente a la forma de la ecuación característica, o lo que es lo mismo de la matriz de transferencia, es posible escribir el polinomio característico como
factorizando y reagrupando términos resulta finalmente que
(4.8) | |||
(4.9) |
Resolviendo este sistema de ecuaciones resulta:
(4.10) |
Además estamos interesados en regiones en que ``'' es complejo o real puro (si es real puro, entonces debe ser mayor o menor que uno en valor absoluto), así como analizar el efecto de si las mismas raíces son puramente imaginarias o puramente reales.
La condición para que sea real puro es que el descriminante de la ecuación sea cero, esto es:
y utilizando la condición mecánica, esto es
resulta finalmente:
Que es exactamente la misma relación que se obtiene al considerar únicamente la condición mecánica, después de ciertos artificios matemáticos.
Otro resultado importante que se puede deducir de las ecuaciones anteriores, es la relación entre las constantes de interacción a primeros y segundos vecinos, esto lo podemos ver de la siguiente manera.
Usando el resultado de que
factorizando y agrupando términos resulta
Analizando esta última desigualdad vemos que
lo cual quiere decir que en cuyo caso, hay raíces complejas y por lo tanto habrá factores de fase compleja; para existirán amortiguaciones puras, o soluciones puramente oscilatorias, en el valor extremo de ; el caso en que es una razón crítica.
Anteriormente mencionamos como una relación de dispersión puede ser deducida, a partir de la ecuación característica, la obtención es directa como podemos apreciarlo de la siguiente forma: Usando la ecuación (4.6), junto con los parámetros definidos en (2.13) resulta:
despejando resulta:
(4.12) |
Que es la relación de dispersión, ya que nos da una relación de frecuencia espacial, con la temporal.
La gráfica de esta relación de dispersión es la siguiente:Otra relación de dispersión que nos puede dar bastante información se obtiene de la siguiente forma.
Pongamos los elementos de la matriz de transferencia y en función de y de , esto es, de (4.1) y (4.2) resulta:
Para solucionar este sistema de ecuaciones, despejamos de (4.3) y sustituyámoslo en (4.14) resultando
Por lo tanto las soluciones para ``'' están dadas por
Esta ecuación es una relación de dispersión porque relaciona a ``'' como función de y y de las constantes del modelo, además se tienen dos alternativas para tratar esta ecuación; una es solucionar ó dependiendo de la conveniencia.
Es conveniente poner la ecuación (4.15) en una forma más explícita anteriormente definimos
por lo tanto
Por lo tanto la ecuación (4.15) queda finalmente como:
Que es la relación de dispersión para el interior de cualquier cadena por lo tanto es conveniente analizarla en forma detallada.
Lo que interesa en principio es analizar cuándo la expresión es real o compleja, ya que si es complejo necesariamente que lo es también.
Por lo tanto la frontera entre real y complejo es cuando el discriminante de (4.16) es igual a cero, esto es
despejando resulta
(4.17) |
Esta ecuación la podemos graficar, ya que por lo tanto es necesario analizarla en la forma siguiente: tiene un polo en .
Cuando entonces y si entonces .
Para obtener un extremo calculemos la derivada, esto es
Como nos interesa un extremo, hagamos , entonces
Por lo tanto en . Evaluando en los puntos resulta
Y es fácil comprobar que ocurre un mínimo en . También es posible ver cuando la expresión es positiva o negativa si y entonces lo cual nos indica que la expresión es negativa.
Además, cuando entonces
Con todos estos datos podemos construir la gráfica de la frontera entre la parte real y compleja , gráfica IV-II.
Como se puede ver, por combinaciones de y en la región compleja es complejo, lo cual quiere decir que es complejo, es decir, la frecuencia en esta región vamos a ver una onda oscilatoria con factor exponencial.
Analicemos las diferentes regiones correspondientes a los distintos valores que puede tomar el parámetro definido como
Caso Real. real
Caso Imaginario
Caso en que existe un cambio de fase. es real, entonces resulta
Ahora estamos interesados en contornos en los que pasa, por estos puntos, para lo cual es necesario hacer algunas manipulaciones algebraícas del , en la forma siguiente: De (4.16) tenemos que:
Despejando resulta
Que es finalmente la relación de dispersión para el interior de cualquier cadena, que nos interesa graficar e interpretar físicamente, además graficaremos ya que se desean encontrar contornos con fijo.
Analizando la ecuación (4.19), vemos que es necesario tomar valores seleccionados de para poder interpretar la gráfica, por lo tanto, si hacemos , ó obtendremos las fronteras de las curvas que separan las partes real e imaginaria pura, esto es:
Si entonces:
Si resulta:
(4.20) |
Debido a que , entonces puede ser que la dada por (4.19) se escape de los límites, por lo tanto es necesario analizar su comportamiento gráficamente. De (4.19) resulta:
Si entonces
Si entonces
En resulta que que es un polo.
Por lo tanto con estos elementos podemos graficar la ecuación (4.19) esto es:
Prosiguiendo el análisis de la ecuación (4.18). Si entonces
Cuando , entonces como necesariamente , pero cuando , las condiciones varían ligeramente.
En existe una línea crítica ya que
En la relación entre y es:
lo cual quiere decir que hay siempre dos valores de para cualquier que esté en el intervalo , tal que está en el intervalo . En otra región, por supuesto que también dos valores reales para , pero necesariamente ambas estarán situadas en el intervalo (esto es ).
Por lo tanto, como las gráficas (IV.II) y (IV.III) pertenecen a la misma ecuación, deducidas únicamente con análisis diferentes, entonces podemos superponerlas y formar una sola, juntamente con los resultados del análisis de dicha ecuación (4.18), esto esAnalizando esta gráfica, junto con las ecuaciones (4.18) y (2.20) resultan las siguientes regiones características:
- Región I.
- Se encuentra en , esta es una región en donde el número de onda () es un número complejo con parte real cero, en este caso el movimiento de la cadena es ondulatorio, ya que resulta una función senosoidal o cosenoidal.
- Región II.
- Está limitada por dos curvas asintóticas, descritas anteriormente (gráfica IV.II) correspondiente a la ecuación . Es una región de complejo, por lo tanto tenemos una onda modulada por un factor exponencial que depende de la parte real de , este factor no es otra cosa que la amplitud, que nos va a modular las funciones seno o coseno.
- Región III.
- Esta región corresponde a valores de , en esta región el número de onda es real positivo o negativo, por lo tanto el movimiento de las partículas en la cadena estará descrito por una exponencial creciente o decreciente dependiendo de que ó .
- Región IV.
- Se encuentra localizada en donde toma valores reales negativos (en la gráfica corresponde a ) este caso se cumple cuando dando como resultado que por lo tanto en este caso se observarán oscilaciones con inversiones de fase, producidas por un desfasamiento de , este caso tiene que ver con el hecho de que exista cooperación o antagonismo entre primeros y segundos vecinos.
No hay comentarios:
Publicar un comentario