Skip to content
Este artículo tiene fines educativos. Te animamos a verificar con fuentes oficiales.

¿Alguna vez te has preguntado cómo las empresas de Guayaquil deciden cuántos productos importar para la temporada navideña? ¿O cómo el Ministerio de Turismo planifica la infraestructura en las Islas Galápagos? La respuesta está en un secreto matemático: **las simulaciones estadísticas**. No se trata de magia, sino de ciencia. El método Monte Carlo —llamado así por el famoso casino de Mónaco— usa el azar para resolver problemas que parecen imposibles. En este simulacro, tú serás el estratega: usarás datos locales, ecuaciones y un poco de suerte controlada para predecir el futuro. ¿Listo para convertirte en un experto en predicciones?

Exercicio 1: Simulación de dados y probabilidades en Cuenca (4 puntos)

Simulación Monte CarloProbabilidad empíricaGeneración de números aleatorios

En un popular juego de azar en Cuenca, se lanzan dos dados estándar de 6 caras. Un jugador apuesta a que la suma de los resultados será mayor que 8. Usando el método Monte Carlo con 1000 simulaciones, determina la probabilidad empírica de ganar esta apuesta y compárala con el valor teórico.

  • Número de simulaciones: 1000
  • Dados estándar: 6 caras cada uno
  • Condición de victoria: suma > 8
  1. Simula 1000 lanzamientos de dos dados y registra la frecuencia de sumas mayores que 8
  2. Calcula la probabilidad empírica de obtener una suma mayor que 8
  3. Determina la probabilidad teórica de este evento y compárala con tu resultado empírico
  4. ¿Qué observas sobre la convergencia entre la probabilidad empírica y teórica a medida que aumenta el número de simulaciones?
Solución completa

Pregunta 1 (1 pts)Simula 1000 lanzamientos de dos dados y registra la frecuencia de sumas mayores que 8

  1. Simulación — Usamos un generador de números aleatorios para crear 1000 pares (d1,d2) donde cada di{1,2,3,4,5,6}. Por ejemplo, el primer par podría ser (2,5) que da suma 7 (no cumple), y el segundo (4,5) que da suma 9 (cumple).

275 éxitos (ejemplo)

→ Se obtienen aproximadamente 270-280 éxitos en 1000 simulaciones (valor exacto depende de la semilla aleatoria).

Pregunta 2 (1 pts)Calcula la probabilidad empírica de obtener una suma mayor que 8

  1. Cálculo de probabilidad empírica — Si hubo 275 éxitos, entonces pemp=275/1000=0.275.
    pemp=2751000=0.275

0.275

→ La probabilidad empírica es aproximadamente 0.275 (27.5%).

Pregunta 3 (1 pts)Determina la probabilidad teórica de este evento y compárala con tu resultado empírico

  1. Cálculo teórico — Como se mencionó, hay 10 casos favorables de 36 posibles: pteo=10/360.2778.
    pteo=1036=0.2778

0.2778

→ La probabilidad teórica es aproximadamente 0.2778 (27.78%).

Pregunta 4 (1 pts)¿Qué observas sobre la convergencia entre la probabilidad empírica y teórica a medida que aumenta el número de simulaciones?

  1. Análisis de convergencia — A medida que aumenta el número de simulaciones, la probabilidad empírica tiende a acercarse a la probabilidad teórica. Esto se debe a la Ley de los Grandes Números: el promedio de los resultados obtenidos en muchas pruebas debe converger al valor esperado.

→ La probabilidad empírica converge a la teórica con más simulaciones.

Rúbrica de evaluación

Simulación correcta de 1000 lanzamientos con generación de números aleatorios1 pts
Cálculo correcto de la probabilidad empírica1 pts
Cálculo correcto de la probabilidad teórica1 pts
Análisis adecuado de la convergencia entre empírico y teórico1 pts

Exercicio 2: Integración Monte Carlo para calcular áreas en Guayaquil (4 puntos)

Integración numéricaMétodo Monte Carlo geométricoEstimación de áreas

En Guayaquil, una empresa de logística necesita calcular el área bajo la curva de demanda diaria en función del precio de un producto. La función de demanda es f(x)=100x2 donde x es el precio en USD (entre 0 y 10). Usa el método Monte Carlo con 10 000 puntos aleatorios para estimar el área bajo esta curva.

  • Función de demanda: f(x)=100x2
  • Intervalo: [0,10]
  • Número de puntos aleatorios: 10 000
  • Rectángulo envolvente: [0,10]×[0,100]
  1. Genera 10 000 puntos aleatorios uniformes en el rectángulo [0,10]×[0,100]
  2. Determina cuántos puntos caen bajo la curva f(x)=100x2
  3. Estima el área bajo la curva usando la fórmula de Monte Carlo
  4. Calcula el error relativo de tu estimación comparado con el valor exacto
Solución completa

Pregunta 1 (1 pts)Genera 10 000 puntos aleatorios uniformes en el rectángulo [0,10]×[0,100]

  1. Generación — Usamos un generador de números aleatorios para crear 10 000 puntos (xi,yi) con xi en [0,10] y yi en [0,100]. Por ejemplo, (3.2,85) o (7.8,40).

→ Se generan 10 000 puntos aleatorios en el rectángulo.

Pregunta 2 (1 pts)Determina cuántos puntos caen bajo la curva f(x)=100x2

  1. Conteo de puntos — Contamos cuántos puntos satisfacen yi100xi2. Supongamos que obtenemos 6 650 puntos bajo la curva.
    nbajo=6650

6650

→ Aproximadamente 6 650 puntos caen bajo la curva.

Pregunta 3 (1 pts)Estima el área bajo la curva usando la fórmula de Monte Carlo

  1. Estimación del área — Aplicamos la fórmula: A665010000×1000=665.
    A0.665×1000=665

665

→ El área estimada es aproximadamente 665 unidades cuadradas.

Pregunta 4 (1 pts)Calcula el error relativo de tu estimación comparado con el valor exacto

  1. Cálculo del error — El valor exacto es 20003666.67. El error relativo es |665666.67|666.67×100%0.25%.
    Error relativo=|665666.67|666.67×100%0.25%

0.25\%

→ El error relativo es aproximadamente 0.25%.

Rúbrica de evaluación

Generación correcta de 10 000 puntos aleatorios en el rectángulo especificado1 pts
Conteo preciso de puntos bajo la curva1 pts
Cálculo correcto del área usando la fórmula Monte Carlo1 pts
Cálculo correcto del error relativo1 pts

Exercicio 3: Optimización de precios usando Monte Carlo en Quito (4 puntos)

OptimizaciónFunciones de beneficioSimulación de escenarios

Un pequeño negocio en el centro histórico de Quito vende un producto artesanal. La función de beneficio diario en USD está dada por B(p)=p(1002p)50, donde p es el precio de venta por unidad. Usa el método Monte Carlo con 5 000 simulaciones para encontrar el precio óptimo que maximice el beneficio esperado.

  • Función de beneficio: B(p)=p(1002p)50
  • Rango de precios a probar: [10,40] USD
  • Número de simulaciones por precio: 5 000
  • Costos fijos: 50 USD/día
  1. Divide el intervalo [10,40] en 10 precios equiespaciados y simula 5 000 beneficios para cada uno
  2. Calcula el beneficio promedio para cada precio simulado
  3. Identifica el precio que maximiza el beneficio promedio
  4. Verifica si este precio coincide con el máximo teórico de la función cuadrática
Solución completa

Pregunta 1 (1 pts)Divide el intervalo [10,40] en 10 precios equiespaciados y simula 5 000 beneficios para cada uno

  1. Precios de prueba — Generamos los 11 precios equiespaciados: 10, 13, 16, ..., 40 USD.

→ Se definen 11 precios de prueba en el intervalo [10,40].

Pregunta 2 (1 pts)Calcula el beneficio promedio para cada precio simulado

  1. Simulación de demanda — Para cada precio pi, generamos 5 000 demandas normales con media μ=1002pi y desviación estándar 10. Por ejemplo, para p=25, μ=50 y generamos djN(50,10).

→ Se simulan 5 000 demandas para cada uno de los 11 precios.

Pregunta 3 (1 pts)Identifica el precio que maximiza el beneficio promedio

  1. Beneficio promedio — Calculamos el beneficio para cada simulación y promediamos. Supongamos que obtenemos: B(10)=450, B(13)=800, B(16)=1050, B(19)=1200, B(22)=1280, B(25)=1300, B(28)=1280, etc.
    B(pi) calculado para cada pi

1300 USD para p=25

→ El beneficio promedio máximo es aproximadamente 1300 USD para p=25.

Pregunta 4 (1 pts)Verifica si este precio coincide con el máximo teórico de la función cuadrática

  1. Comparación con máximo teórico — El máximo teórico es en p=25 con B(25)=1150 USD. La simulación da un beneficio promedio de 1300 USD, que es mayor debido a la variabilidad estocástica en la demanda. El precio óptimo simulado coincide con el teórico.
    póptimo=25 en ambos casos

25 USD

→ El precio óptimo es 25 USD, coincidiendo con el máximo teórico.

Rúbrica de evaluación

Definición correcta de los 11 precios de prueba en el intervalo1 pts
Simulación adecuada de demanda normal con parámetros correctos1 pts
Cálculo correcto del beneficio promedio para cada precio1 pts
Identificación correcta del precio óptimo y comparación con el teórico1 pts

Exercicio 4: Simulación de riesgo financiero para préstamos en USD (4 puntos)

Riesgo financieroDistribuciones de probabilidadMonte Carlo en finanzas

Una entidad financiera en Ambato ofrece préstamos personales de 10 000 USD a 5 años con tasa de interés variable. La tasa anual sigue una distribución normal con media 8% y desviación estándar 2%. Usa el método Monte Carlo con 2 000 simulaciones para estimar la probabilidad de que el cliente no pueda pagar el préstamo (es decir, que la tasa supere el 12% en algún año).

  • Monto del préstamo: 10 000 USD
  • Plazo: 5 años
  • Tasa de interés anual: riN(μ=8%,σ=2%)
  • Condición de impago: tasa > 12% en algún año
  • Número de simulaciones: 2 000
  1. Genera 2 000 trayectorias de tasas de interés para los 5 años
  2. Para cada trayectoria, verifica si en algún año la tasa supera el 12%
  3. Calcula la probabilidad empírica de impago
  4. Interpreta el resultado en términos de riesgo crediticio para la entidad
Solución completa

Pregunta 1 (1 pts)Genera 2 000 trayectorias de tasas de interés para los 5 años

  1. Generación de trayectorias — Creamos 2 000 series de 5 tasas cada una, todas extraídas de una normal con media 8% y desviación 2%. Por ejemplo, una trayectoria podría ser [7.2%, 9.5%, 11.8%, 6.3%, 8.9%].

→ Se generan 2 000 trayectorias de 5 tasas cada una.

Pregunta 2 (1 pts)Para cada trayectoria, verifica si en algún año la tasa supera el 12%

  1. Detección de impago — Para cada trayectoria, verificamos si alguna tasa supera el 12%. Por ejemplo, en [7.2%, 9.5%, 11.8%, 6.3%, 8.9%] hay un impago en el tercer año (11.8% < 12%? No, 11.8% < 12% es falso, espera: 11.8 < 12 es verdadero, entonces NO hay impago. Corregir: condición es r > 12, así que 11.8 no cuenta. Un ejemplo con impago: [9.1%, 12.3%, 7.8%, 10.5%, 8.2%] tiene impago en el segundo año.

→ Se identifican las trayectorias con al menos una tasa > 12%.

Pregunta 3 (1 pts)Calcula la probabilidad empírica de impago

  1. Cálculo de probabilidad — Supongamos que 640 de las 2 000 trayectorias tienen al menos una tasa > 12%. Entonces pimpago=640/2000=0.32=32%.
    pimpago=6402000=0.32

32\%

→ La probabilidad empírica de impago es 32%.

Pregunta 4 (1 pts)Interpreta el resultado en términos de riesgo crediticio para la entidad

  1. Interpretación — Una probabilidad del 32% de que el cliente no pueda pagar en algún momento sugiere un riesgo crediticio alto. La entidad podría requerir garantías adicionales o aumentar la tasa de interés para compensar el riesgo.

→ El riesgo de impago es alto (32%), requiriendo medidas de mitigación.

Rúbrica de evaluación

Generación correcta de 2 000 trayectorias de 5 tasas cada una1 pts
Verificación precisa de la condición de impago (>12%) en cada trayectoria1 pts
Cálculo correcto de la probabilidad empírica de impago1 pts
Interpretación adecuada del resultado en contexto financiero1 pts

Exercicio 5: Predicción de demanda turística en las Islas Galápagos (4 puntos)

Distribuciones de probabilidadSimulación de series temporalesAnálisis de escenarios

El Ministerio de Turismo de Ecuador planea infraestructura en las Islas Galápagos. El número anual de turistas sigue una distribución normal con media 15 000 y desviación estándar 2 000. Usa el método Monte Carlo con 3 000 simulaciones para estimar la probabilidad de que en los próximos 10 años al menos 2 años tengan más de 18 000 turistas.

  • Distribución anual: XN(μ=15000,σ=2000)
  • Período: 10 años
  • Umbral: 18 000 turistas
  • Número de simulaciones: 3 000
  1. Genera 3 000 series de 10 años de demanda turística
  2. Para cada serie, cuenta cuántos años superan los 18 000 turistas
  3. Calcula la probabilidad de que al menos 2 años por serie superen el umbral
  4. Estima el número esperado de años con más de 18 000 turistas por década
Solución completa

Pregunta 1 (1 pts)Genera 3 000 series de 10 años de demanda turística

  1. Generación de series — Creamos 3 000 series de 10 años cada una, con demanda anual distribuida normal(15000, 2000). Por ejemplo, una serie podría ser [14500, 16200, 17800, 15300, 18100, ...].

→ Se generan 3 000 series de 10 años de demanda turística.

Pregunta 2 (1 pts)Para cada serie, cuenta cuántos años superan los 18 000 turistas

  1. Conteo de años críticos — Para cada serie, contamos cuántos valores superan 18 000. Por ejemplo, en [14500, 16200, 17800, 15300, 18100, ...] hay 1 año (18100) que supera el umbral.

→ Se cuenta el número de años >18000 en cada serie.

Pregunta 3 (1 pts)Calcula la probabilidad de que al menos 2 años por serie superen el umbral

  1. Probabilidad de al menos 2 años — Supongamos que 1 200 de las 3 000 series tienen al menos 2 años con más de 18 000 turistas. Entonces p=1200/3000=0.40=40%.
    p=12003000=0.40

40\%

→ La probabilidad de al menos 2 años críticos por década es del 40%.

Pregunta 4 (1 pts)Estima el número esperado de años con más de 18 000 turistas por década

  1. Número esperado — Como se calculó, el número esperado de años con más de 18 000 turistas por década es aproximadamente 0.668 años (menos de 1 año por década en promedio).
    𝔼[ck]0.668 años

0.67 años

→ Se esperan aproximadamente 0.67 años con más de 18 000 turistas por década.

Rúbrica de evaluación

Generación correcta de 3 000 series de 10 años con distribución normal1 pts
Conteo preciso del número de años >18000 en cada serie1 pts
Cálculo correcto de la probabilidad de al menos 2 años críticos1 pts
Cálculo correcto del número esperado de años críticos por década1 pts

Fuentes

  1. en.wikipedia.org
  2. doi.org
  3. api.semanticscholar.org
  4. ui.adsabs.harvard.edu
  5. search.worldcat.org
  6. www.osti.gov
  7. arxiv.org
  8. www.laas.fr
  9. viewer.zmags.com
  10. www.crcpress.com
  11. apps.dtic.mil