2 de 4
Agregar imagen

La respuesta corta es: "Se debe al juego de datos que posees".

La respuesta larga, es mejor con un ejemplo... imagínate por un momento que tienes otro juego de datos. Si estas de acuerdo podrían ser:

datos = [
    30.05251300, 19.14849600, 25.31769200, 27.59143700,
    32.07645600, 23.48796100, 28.47594000, 35.12375300,
    36.83848500, 25.00701700, 30.72223000, 28.69375900,
    36.64098600, 23.82460900, 29.31168300, 31.77030900,
    35.17787700, 19.77524400, 29.60175000, 34.53884200,
    41.27359900, 26.65586200, 28.27985900, 35.19115300,
    42.20566386, 24.64917133, 32.66733514, 37.25735401,
    45.24246027, 29.35048127, 36.34420728, 41.78208136,
    49.27659843, 31.27540139, 37.85062549, 38.83704413,
    51.23690034, 31.83855162, 41.32342126, 42.79900337,
    55.70835836, 33.40714492, 42.31663797, 45.15712257,
    59.57607996, 34.83733016, 44.84168072, 46.97124960,
    60.01903094, 38.37117851, 46.97586413, 50.73379646,
    61.64687319, 39.29956937, 52.67120908, 54.33231689,
    66.83435838, 40.87118847, 51.82853579, 57.49190993,
    65.25146985, 43.06120822, 54.76075713, 59.83447494,
    73.25702747, 47.69662373, 61.09776802, 66.05576122]

indice = pd.date_range("2020-10-31 23:57:00", periods=len(datos), freq="T")

datos = pd.Series(data=datos, index=indice)
        
datos = datos.asfreq(freq='T')

De esta manera tendrías un codigo similar a esto:

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
from statsmodels.tsa.exponential_smoothing.ets import ETSModel
    
datos = [
        30.05251300, 19.14849600, 25.31769200, 27.59143700,
        32.07645600, 23.48796100, 28.47594000, 35.12375300,
        36.83848500, 25.00701700, 30.72223000, 28.69375900,
        36.64098600, 23.82460900, 29.31168300, 31.77030900,
        35.17787700, 19.77524400, 29.60175000, 34.53884200,
        41.27359900, 26.65586200, 28.27985900, 35.19115300,
        42.20566386, 24.64917133, 32.66733514, 37.25735401,
        45.24246027, 29.35048127, 36.34420728, 41.78208136,
        49.27659843, 31.27540139, 37.85062549, 38.83704413,
        51.23690034, 31.83855162, 41.32342126, 42.79900337,
        55.70835836, 33.40714492, 42.31663797, 45.15712257,
        59.57607996, 34.83733016, 44.84168072, 46.97124960,
        60.01903094, 38.37117851, 46.97586413, 50.73379646,
        61.64687319, 39.29956937, 52.67120908, 54.33231689,
        66.83435838, 40.87118847, 51.82853579, 57.49190993,
        65.25146985, 43.06120822, 54.76075713, 59.83447494,
        73.25702747, 47.69662373, 61.09776802, 66.05576122]
    
indice = pd.date_range("2020-10-31 23:57:00", periods=len(datos), freq="T")
    
datos = pd.Series(data=datos, index=indice)
    
datos = datos.asfreq(freq='T')
          
          
pasado = datos[:48]
futuro = datos[47:]

              
modelo = ETSModel(datos, error="add", trend="add", seasonal="add",
                    damped_trend=True, seasonal_periods=4)
#modelo_fit = modelo.fit(maxiter=10000)
fit = modelo.fit()
    
print(fit.summary())
    
pred = fit.get_prediction(start='2020-11-01 00:44:00', end='2020-11-01 01:04:00')
    
df = pred.summary_frame(alpha=0.05)
    
    
simulated = fit.simulate(anchor="end", nsimulations=10, repetitions=100)

for i in range(simulated.shape[1]):
  simulated.iloc[:,i].plot(label='_', color='gray', alpha=0.1)
      
df["mean"].plot(label='mean prediction')
df["pi_lower"].plot(linestyle='--', color='tab:cyan', label='95% interval')
df["pi_upper"].plot(linestyle='--', color='tab:cyan', label='_')

pred.endog.plot(label='data')
plt.legend()
plt.show()

Obtendrías un resultado de este tipo:

Ejemplo: simulación_1

Tus datos estan representados en color anaranjado. El modelo ETS, estima un promedio de los datos en azul, y un rango en el que los datos pueden variar basado en la media (que son las lineas intermitentes de color cian). Después (en la predicción) el modelo realiza una simulación tratando de pronosticar, 10 pasos adelante, y realiza 100 intentos (son las lineas grises).

En este caso particular, el modelo se ajusta muy bien a los datos... ¡pero claro!, es un ejemplo de libro de texto, así que funcionara a la perfección.

En tu caso los datos, son a falta de otra palabra estocásticos a los ojos del modelo y este modelo en particular, le cuesta generar o decidirse hacia donde pueden ir los datos en el futuro, así que estima una media, unos rangos superiores e inferiores en los cuales pueden estar los datos en el futuro.

Tomemos el mismo código, y solo variemos los datos, tendrías algo asi:

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
from statsmodels.tsa.exponential_smoothing.ets import ETSModel

    
pasado = datos[:100]
futuro = datos[99:]
print(futuro)
        
modelo = ETSModel(datos, error="add", trend="add", seasonal="add",
              damped_trend=True, seasonal_periods=4)
#modelo_fit = modelo.fit(maxiter=10000)
fit = modelo.fit()

print(fit.summary())

#prediccion = modelo_fit.get_prediction(start='2020-11-01 01:37:00', end='2020-11-01 01:57:00')
pred = fit.get_prediction(start='2020-11-01 01:36:00', end='2020-11-01 01:56:00')

df = pred.summary_frame(alpha=0.05)




simulated = fit.simulate(anchor="end", nsimulations=20, repetitions=100)
for i in range(simulated.shape[1]):
  simulated.iloc[:,i].plot(label='_', color='gray', alpha=0.1)


df["mean"].plot(label='mean prediction')
df["pi_lower"].plot(linestyle='--', color='tab:cyan', label='95% interval')
df["pi_upper"].plot(linestyle='--', color='tab:cyan', label='_')
pred.endog.plot(label='data')

pasado.plot(label='Pasado')
futuro.plot(label='Futuro')



plt.legend()
plt.show()

Simulación

Después de los datos de entrenamiento (en color verde), sigue una especie de burbuja (lo contenido entre las lineas punteadas de color cian), que es una estimación (segun el modelo), de el lugar donde podrían estar lo datos en el futuro, así que la linea que te aparece generalmente con el mismo valor, es el promedio estimado de los valores futuros que augura el modelo. O sea, según los datos, el modelo no se logra ajustar con precisión a los datos que tienes en la variable futuro.

introducir la descripción de la imagen aquí Futuros

Un modelo que se puede ajustar con mayor precisión podría ser SARIMA o SARIMAX, lo mejor es buscar (para los casos anteriores) algún mecanismo/libreria que ajuste los valores order=(p,d,q) y seasonal_order=(P,D,Q,s) de manera automática (aunque el coste computacional puede empezar a elevarse).