Se describe un nuevo modelo estadístico de distribución espacial de galaxias observadas. Las correlaciones estadísticas se incluyen mediante ensambles de cadenas de Markov, cuyos parámetros se extraen del espectro de potencia observable mediante la adopción del ansatz de Uchaikin-Zolotarev. Las trayectorias de las cadenas de Markov con longitudes de paso distribuidas de manera Lévy-Feldheim forman el conjunto de nodos que imitan las posiciones de las galaxias. El modelo reproduce plausiblemente las funciones de correlación a dos puntos, los datos de recuento de celdas y algunas otras propiedades importantes. Puede utilizarse de manera efectiva en el postprocesamiento de datos astronómicos para estudios cosmológicos.