6 votos

Logaritmo de una matriz de Markov

Comienza con una matriz de Markov $\mathbf{M}$, cuyos elementos están todos entre $0 \le \mathbf{M}_{ij} \le 1$ y cada fila suma uno. Existe una conexión natural entre esta matriz y la matriz de tasas $\mathbf{W}$ en la Ecuación Maestra

$$ \mathbf{M} = \exp( t \mathbf{W} ) $$

Aquí, dado $\mathbf{W}$, el cálculo de $\mathbf{M}$ es único con $t$ ya que la exponencial de matriz es única y converge por la expansión de Taylor. ¿Qué tal en la otra dirección?

$$ t \mathbf{W} = \log( \mathbf{M} ) $$

¿Las propiedades de la matriz de Markov garantizan que esto sea único y que la suma alterna en la serie de Taylor de log converja?

¡Por favor, proporcione una referencia donde se discuta esto si es posible!

Motivación (por solicitud)

He estado estudiando la dinámica asociada con un sistema tipo modelo de Ising bajo dinámicas de Glauber de flip de un solo espín. Las dinámicas de Glauber dan, esencialmente, la matriz $\mathbf{W}$. Si uno observara el sistema durante un tiempo finito, se podría hacer una aproximación a $\mathbf{M}$. Me interesaba saber cuándo era permisible convertir entre las dos. En la referencia proporcionada por una de las respuestas, la pregunta se reduce a:

En términos probabilísticos, una matriz de Markov A es incrustable si se obtiene tomando una instantánea en un momento particular de un proceso de Markov finito autónomo que se desarrolla continuamente en el tiempo. Por otro lado, una matriz de Markov podría no ser incrustable si describe los cambios anuales en una población que tiene un patrón de reproducción fuertemente estacional; en tales casos, se podría construir un modelo más elaborado que incorpora las variaciones estacionales. La incrustabilidad también puede fallar porque las entradas de la matriz no son precisas; en tales casos, una técnica de regularización podría producir una matriz de Markov muy similar que sea incrustable;

6voto

Did Puntos 1

Este es un tema antiguo que se remonta al menos a Elfving en los años 30. Kingman se interesó en él en los años 60 y más recientemente Israel, Rosenthal y Wu lo estudiaron con un ojo en el aspecto práctico del problema.

Las matrices de Markov $M$ que pueden escribirse como $M=\exp(W)$ para una matriz generadora $W$ se llaman incrustables porque entonces, existe un semigrupo completo $(M_t)$ de matrices de Markov tales que $M=M_1$ (simplemente define $M_t=\exp(tW)$ para cada número real no negativo $t$). Recordemos que $W$ es un generador si la suma de sus filas es cero y si cada elemento no diagonal es no negativo.

Existen algunas obstrucciones obvias para que una matriz de Markov dada $M$ sea incrustable. Por ejemplo, el conjunto de pares de estados $(x,y)$ tal que $(M^n)_{xy}=0$ no puede depender de $n\ge1$. También puede suceder que la matriz $V=\log(M)$ definida por la expansión en serie usual de $x\mapsto\log(1+x)$ en $x=0$ aplicada a la matriz $I+(M-I)$, no sea un generador a pesar de que $M=\exp(V)$ (parece ser que esto incluso puede suceder cuando $M$ es incrustable).

A excepción de la dimensión $2$ (y quizás la dimensión $3` pero no recuerdo), no está disponible una caracterización completa de la incrustabilidad, pero se conoce mucho al respecto. Sugiero comenzar con el artículo reciente Matrices de Markov Incrustables de E. B. Davies. Está disponible en el arXiv y proporciona referencias útiles a algunos trabajos anteriores, incluidos los que mencioné en esta publicación.

Un problema relacionado es saber cuándo la ecuación $M=\exp(W)$ tiene varias soluciones admisibles $W`, sobre esto consultar el artículo de Speakman en la lista de referencias del artículo de Davies.

1voto

user8268 Puntos 13913

En general, no hay $W$ con coeficientes reales. Aquí hay un ejemplo al azar de aquí: $M=\pmatrix{1/2&1/4&1/4\\ 1/6&1/6&2/3 \\ 1/3&1/3&1/3}$ tiene autovalores $1$ y $\pm1/\sqrt{24}$. Si $M=\exp W$ entonces los autovalores de $W$ deben ser módulo $2\pi i\mathbb{Z}$ iguales a $0$, $-\log(24)/2$, $-\log(24)/2+\pi i$. Si $W$ es real entonces sus autovalores son o reales o vienen en pares conjugados complejos, y aquí no es el caso.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X