Supongamos que $B$ es una matriz de $n\times n$ de tasas de transición de Poisson, donde $B_{ij}\geq 0$ para $i\neq j$ denota la tasa a la cual el estado $i$ transiciona al estado $j$, y $B_{ii}\leq 0$ da la tasa a la cual el estado $i$ transiciona a todos los demás estados. Cada fila de $B$ suma 0.
Entonces, si $p(t)$ denota la distribución de probabilidad en el tiempo $t$, por definición de $B$ tenemos la EDO $$\dot{p}(t) = Bp(t)$$ Sabemos cómo se ve la solución a este tipo de EDO: $p(t)=e^{Bt}p(0)$, donde $e^{Bt}$ es la exponencial de matriz de $Bt$. Entonces, si queremos que $B$ genere la matriz de transición de Markov $A$ después de $t=1$, necesitamos tener $e^B=A$.
En principio, para obtener $B$, necesitamos invertir la exponencial de matriz, tomando el logaritmo de matriz de $A$. El problema es que cada matriz tiene muchos logaritmos de matriz: el logaritmo en el espacio complejo unidimensional tiene infinitas ramas, y esto se complica cuando hablamos de matrices en un espacio $n$-dimensional. La mayoría de estos logaritmos no serán matrices de transición de Poisson satisfactorias: tal vez no serán reales, o las entradas no tendrán los signos correctos. Sin embargo, es posible que más de uno de ellos lo sea: en algunos casos hay más de un $B$ de Poisson correspondiente a un $A$ de Markov, al igual que en algunos casos no hay un $B$ de Poisson correspondiente a $A$. Es complicado.
Afortunadamente, existe una situación en la que la vida es relativamente simple, y casi con seguridad incluye tu caso: cuando todos los valores propios de $A$ son reales positivos distintos. En este caso, hay solo un logaritmo de $A$ que será real, y es fácil de calcular: simplemente diagonalizas la matriz como $A=V\Sigma V^{-1}$ y tomas el logaritmo real de los valores propios, obteniendo $B=V\Omega V^{-1}$, donde $\omega_{ii} = \log(\sigma_{ii})$. De hecho, no necesitas hacer esto tú mismo: si usas el comando $\text{logm}(A)$ en Matlab (presumiblemente también en Python), te dará precisamente este $B$.
Dado este $B$, todo lo que tienes que hacer es verificar si en realidad es una matriz de Poisson. El primer requisito, que las filas sumen cero, se cumple automáticamente debido a la construcción de $B$.** El segundo requisito, que los elementos diagonales sean negativos y los elementos fuera de la diagonal sean positivos, no siempre se cumple (creo), pero es fácil de verificar.
Para ver esto en acción, consideraré un $A$ para un proceso de Markov de 3 estados que se asemeja a un AR(1) discretizado. $$A = \begin{pmatrix}0.5 & 0.4 & 0.1 \\ 0.2 & 0.6 & 0.2 \\ 0.1 & 0.4 & 0.5\end{pmatrix}$$ Ahora, si escribo $B=\text{logm}(A)$ en Matlab, obtengo $$B = \begin{pmatrix}-0.86 & 0.80 & 0.06 \\ 0.40 & -0.80 & 0.40 \\ 0.06 & 0.80 & -0.86\end{pmatrix}$$ Esto es de hecho una matriz de transición de Poisson válida, ya que podemos verificar fácilmente que las filas suman cero y tienen los signos correctos - así que esta es nuestra respuesta.
El caso con valores propios positivos es bastante importante, ya que abarca todos los casos en los que no hay algún tipo de comportamiento oscilatorio en la cadena de Markov (que requeriría valores propios negativos o complejos), presumiblemente incluyendo tu AR(1) discretizado.
Más en general, el comando $\text{logm}$ en Matlab nos dará el logaritmo de matriz principal, un análogo del logaritmo escalar principal que hace que todos los valores propios tengan una parte imaginaria entre $-\pi$ y $\pi$. El problema es que este no es necesariamente el logaritmo que queremos, y al verlo podríamos pasar por alto un $B$ de Poisson que sí genere $A. (Por eso el caso de valores propios positivos, donde no tuvimos que preocuparnos por esto, fue tan agradable.) Sin embargo, incluso en estos otros casos no está de más intentar y ver si funciona.
Por cierto, este problema de ver si hay un $B$ que genere alguna matriz de Markov $A$ ha sido estudiado extensamente. Se le llama el problema de embebibilidad: consulta un resumen y referencias en este excelente artículo de revisión de Davies. No soy un experto en aspectos técnicos del problema, sin embargo; esta respuesta se basa más en mi propia experiencia chapucera e intuición.
Me siento obligado a cerrar secundando el comentario de ecksc y diciendo que podría haber mejores y más directas formas de convertir un AR(1) ajustado de forma discreta en un proceso continuo de tiempo finito, en lugar de simplemente tomar la matriz obtenida a través del método Tauchen y hacerla continua. ¡Pero personalmente no sé cuál es esa mejor manera!
**Explicación (aunque estoy algo oxidado): $A$ tiene un valor propio de Perron-Frobenius único de 1, y como $A$ es estocástica el autovector derecho de este valor propio es el vector unitario $e$. Este sigue siendo el autovector derecho, ahora con un valor propio de 0, cuando tomamos el logaritmo de matriz.