IA Math
IA Math
interiores de la tierra, mejor conocidas como placas tectónicas “cuando las placas se
caracterizado por ser una de las zonas geográficas con mayor movimiento sísmico en
convergencia entre Nazca y Sudamérica” (CSN, 2010) lo cual hace que la sismicidad en
esta zona sea mayormente intensa y como efecto se produce que haya muchos
terremotos en Chile. Elegí este tema porque en noveno en clase de física realizamos
una investigación sobre diferentes tipos de ondas que podemos ver en la naturaleza,
entre esas las ondas sísmicas. Por supuesto, yo seleccioné ese tema y me terminó
encantando el poder comprender más allá del funcionamiento de las placas tectónicas
y la función de las ondas dentro del movimiento sísmico en la tierra. Asimismo, previo a Commented [LHEG2]: Criterio C: Muestra su interés
en el tema, y la motivación de poder profundizar en el
esa investigación en clase yo ya había indagado un poco sobre los sismos, porque me conocimiento del mismo.
llamó mucho la atención la particularidad con Chile. El cual, al ser un país tan cercano
la escala de Richter, y están preparados para enfrentar cualquier efecto que esto les
pueda causar.
2
y tsunamis que han llegado a causar a la muerte y desaparición de personas y animales,
falta de necesidades básicas, entre otras afecciones provocadas por este fenómeno
natural. Es por eso, que por medio de esta investigación, desarrollaré un modelo
de sismos en Chile para 2020, 2021 y 2022. Esto consiste en la recopilación de datos Commented [LHEG3]: Criterio A: Establece un
objetivo claro, en el que es evidente lo que va a hacer,
desde 2011 a 2019 del número de sismos que han ocurrido en cada uno de estos 10 y el para qué de su modelo.
años; después de esto, se procederá con la elaboración de tablas y gráficas con los
datos obtenidos que permitan realizan un ajuste de mínimos cuadrados. Tras obtener
resultados, se podría establecer una función matemática con la que se puede realizar
una aproximación a futuro de la cantidad de sismos que pueden llegar a haber entre
MARCO TEORICO
una regresión polinomial de tercer grado, también conocida como ajuste cúbico de
mínimos cuadrados. Esta técnica matemática, también llamada regresión lineal consiste
general de los datos, sin coincidir necesariamente en todos los puntos.” (Chapra &
Canale, 2006), es por eso que se termina considerando como una aplicación donde se
optimiza el uso de las matemáticas que ayuden a determinar la función más apropiada
que se aproxime a los datos reales. El ajuste de mínimos cuadrados nos proporciona
una representación visual entre las coordenadas de los datos y la línea de tendencia es
proporcionada por el gráfico de dispersión, el cual nos muestra con mayor facilidad la
menor distancia entre los puntos del modelo y de la línea de tendencia del modelo que
tiene una mayor aproximación. Su creación se debe al matemático alemán Carl Friedrich
3
matemático francés Andrien-Marie Legendre en 1805, quien no copió la idea de Gauss,
Pearson, que oscila entre un valor de 0 a 1, siendo 1 el modelo perfecto donde todos los Commented [LHEG4]: Criterio B: Define términos que
posteriormente van a ser fundamentales en el
datos del modelo se ajustan a los datos reales. Esto nos permite evaluar el porcentaje desarrollo de su modelo.
entre el modelo y las observaciones” (Chapra & Canale, 2006), donde esta diferencia es
de valor real y del modelo que se predijo con la suma de 𝑎0 y 𝑎1 𝑥. Ahora el siguiente
paso para tener el mejor ajuste con base a lo previamente mencionado, se puede por
medio de la línea minimizar la suma de los errores residuales de todos los datos, o
continuación:
𝑛 𝑛 𝑛
Commented [LHEG5]: Criterio B: Explicó el
𝑆𝑟 = ∑ 𝑒𝑖2 = ∑(𝑦𝑖, 𝑚𝑒𝑑𝑖𝑑𝑎 − 𝑦𝑖,𝑚𝑜𝑑𝑒𝑙𝑜 )² = ∑(𝑦𝑖 − 𝑎0 − 𝑎1 𝑥𝑖 ) 1.1significado de algunas de las variables, pero todas no
𝑖=1 𝑖=1 𝑖=1 están claramente definidas.
En este caso, el ajuste que más conviene para una mejor línea es el criterio
“minimax” (Chapra & Canale, 2006), la cual toma como referencia datos de una función
simple a una compleja, es decir con números complejos. La cual “consiste en minimizar
modelo lineal.” (Chapra & Canale, 2006) como vemos a continuación, donde se derivan
4
la cual al ser una función con diferentes variables, nos indica la función (𝑑𝑓) cuando se
En los 3 casos se consideran las derivadas de cada una de las variables (𝑥, 𝑦, 𝑧) y se
surge la relación entre la derivada parcial y la ecuación 1.1, puesto que al tener 3
variables 𝑆𝑟 :
𝜕𝑆𝑟
= −2 ∑(𝑦𝑖 − 𝑎0 − 𝑎1 𝑥𝑖 )
𝜕𝑎0
1.2
𝜕𝑆𝑟
= −2 ∑[(𝑦𝑖 − 𝑎0 − 𝑎1 𝑥𝑖 )𝑥𝑖 ]
𝜕𝑎1
𝑛𝑎0 , que se reflejan como un conjunto de dos ecuaciones lineales simultaneas con dos
continuación:
𝑛𝑎0 + (∑ 𝑥𝑖 ) 𝑎1 = ∑ 𝑦𝑖
1.3
(∑ 𝑥𝑖 ) 𝑎0 + (∑ 𝑥𝑖2 ) 𝑎1 = ∑ 𝑥1 𝑦𝑖
5
Como se ve acá, este es el origen del resto de ajustes o regresiones lineales o
polinomiales. De hecho, para esta investigación se utilizará la regresión polinomial de Commented [LHEG6]: Criterio E: Considera
elementos teóricos valiosos para llegar a una
tercer grado (o cubica), que al igual de la de segundo grado, se basan de la expresión comprensión sólida en la creación de su modelo
utilizando la técnica de mínimos cuadrados, sin
embargo no profundiza o explica más detalladamente
general lineal, pero se ve de esta manera: 𝑦 = 𝑎0 + 𝑎1 𝑥 + 𝑎2 𝑥 2 + 𝑎3 𝑥 3 , esta la algunos elementos, por ejemplo los mencionados
respecto a funciones complejas o derivadas parciales.
consideráremos como la ecuación general u origen del resto de regresiones, mientras
Por lo que, en este caso “tenemos una sola variable explicativa continua, x, pero
𝑥, para describir diversos tipos de curvatura en la relación 𝑦 𝑥.” (Vinuesa, 2016). En este
𝑛𝑎0 + (∑ 𝑥 ) 𝑎1 + (∑ 𝑥 2 ) 𝑎2 + (∑ 𝑥 3 ) 𝑎3 = ∑ 𝑦
(∑ 𝑥 ) 𝑎0 + (∑ 𝑥 2 ) 𝑎1 + (∑ 𝑥 3 ) 𝑎2 + (∑ 𝑥 4 ) 𝑎3 = ∑ 𝑥𝑦 2.1
(∑ 𝑥 2 ) 𝑎0 + (∑ 𝑥 3 ) 𝑎1 + (∑ 𝑥 4 ) 𝑎2 + (∑ 𝑥 5 ) 𝑎3 = ∑ 𝑥 2 𝑦
(∑ 𝑥 3 ) 𝑎0 + (∑ 𝑥 4 ) 𝑎1 + (∑ 𝑥 5 ) 𝑎2 + (∑ 𝑥 6 ) 𝑎3 = ∑ 𝑥 3 𝑦
de sismos que pueden llegar a ocurrir entre los años de 2020 a 2022 en Chile, sin
importar su magnitud, intensidad o nivel de afectación. Es por eso, que para poder
realizar este análisis se utilizaron datos reales del Centro Nacional de Sismología de
Chile partiendo del número de sismos que ocurrieron en 2010, año 1, mismo año en el
cual ocurrió el terremoto de 2010 de magnitud de 8,8; hasta 2019, año 10. Aquí podemos
ver como esta cantidad varia en estos 10 años, en los cuales los años 7, 8 y 9 (2016,
6
2017 y 2018 respectivamente) son los que presentan un menor movimiento sísmico, a
𝟏 7733
2 7079
3 8094
4 6437
5 6499
6 6637
7 5119
8 4816
9 5447
10 7906
Viendo esta tabla, se podría desarrollar con los datos tomados un modelo
aparece al comparar los datos con el ajuste, mientras más cercano el 𝑅² esté de 1, más
7
apropiado va a ser la ecuación que resulte del ajuste para la investigación. En este caso,
𝑦 = 𝑎3 𝑥 3 + 𝑎2 𝑥 2 + 𝑎1 𝑥 + 𝑎0 3.1
esté más cerca de 1, mayor el ajuste. Aquí ninguno llega 1 pero todos se
Cuadrático 0,4437
Cubico 0,7676
Logarítmico 0,2872
Exponencial 0,2403
Potencial 0,2944
anualmente en Chile desde 2010 hasta 2019. En la gráfica 1, podemos ver la línea de
tendencia del modelo cúbico. Claramente, se puede ver que hay datos que coinciden
8
Grafica 1: Grafica de dispersión de la cantidad de sismos anuales en Chile. Es
es la gráfica del modelo cúbico y los puntos son los datos reales en Chile.
el ajuste de mínimos cuadrados del modelo cubico. En este caso, primero se organizan
la variable dependiente (𝑦) y la variable independiente (𝑥) en una tabla con ayuda de
sumatoria de cada una de las columnas. Como se puede ver en la tabla 3 a continuación.
9
Tabla 3: Ajuste cúbico de mínimos cuadrados con los datos reales. Se presentan
V. V. 𝒙𝟐 𝒙𝟑 𝒙𝟒 𝒙𝟓 𝒙𝟔 𝒙𝒚 𝒙𝟐 𝒚 𝒙𝟑 𝒚
independ. depend.
(𝒙) (𝒚)
1 7733 1 1 1 1 1 7733 7733 7733
2 7079 4 8 16 32 64 14158 28316 56632
3 8094 9 27 81 243 729 24282 72846 218538
4 6437 16 64 256 1024 4096 25748 102992 411968
5 6499 25 125 625 3125 15625 32495 162475 812375
6 6637 36 216 1296 7776 46656 39822 238932 1433592
7 5119 49 343 2401 16807 117649 35833 250831 1755817
8 4816 64 512 4096 32768 262144 38528 308224 2465792
9 5447 81 729 6561 59049 531441 49023 441207 3970863
10 7906 100 1000 10000 100000 1000000 79060 790600 7906000
= 19039310
= 1978405
= 2404156
= 220825
= 346682
= 25333
= 3025
= 65767
= 𝟓𝟓
𝒙
𝑥3
𝑥2 𝑦
𝑥4
𝑥𝑦
= 385
𝑥3𝑦
𝑥5
∑
𝑥2
𝑥6
𝑦
∑
∑
𝑎1 𝑥 + 𝑎0 , se deben hallar los valores de cada uno de los coeficientes (𝑎3 , 𝑎2 , 𝑎1 , 𝑎0) que
dependiente (y que equivale a la cantidad de sismos) es usada para calcular el ajuste Commented [LHEG7]: Criterio B: Define las variables
que usará en su modelo.
cubico de mínimos cuadrados por medio de un sistema de ecuaciones 4 × 4. Para
𝑛𝑎0 + (∑ 𝑥 ) 𝑎1 + (∑ 𝑥 2 ) 𝑎2 + (∑ 𝑥 3 ) 𝑎3 = ∑ 𝑦
3.2/
2 3 4 2.1
(∑ 𝑥 ) 𝑎0 + (∑ 𝑥 ) 𝑎1 + (∑ 𝑥 ) 𝑎2 + (∑ 𝑥 ) 𝑎3 = ∑ 𝑥𝑦
(∑ 𝑥 2 ) 𝑎0 + (∑ 𝑥 3 ) 𝑎1 + (∑ 𝑥 4 ) 𝑎2 + (∑ 𝑥 5 ) 𝑎3 = ∑ 𝑥 2 𝑦
(∑ 𝑥 3 ) 𝑎0 + (∑ 𝑥 4 ) 𝑎1 + (∑ 𝑥 5 ) 𝑎2 + (∑ 𝑥 6 ) 𝑎3 = ∑ 𝑥 3 𝑦
10
Para poder definir la ecuación, la n se debe reemplazar por el número de datos
en la variable independiente, en este caso 10, que corresponden a los años entre 2010
a 2019. Después de esto, se reemplazaron todas las ∑ 𝑥 por sus respectivas sumatorias
donde se digitó el sistema de ecuaciones, solo que en lugar de poner “𝑎”, se sustituyó
por una “𝑥”. Al realizar esto en la ventana de “Ecuación simultánea con 4 incógnitas”,
nos dieron los siguientes datos para la variable dependiente como se ve en las figuras
Figura 2: Tabla con resultados del sistema de ecuaciones del ajuste cúbico
11
Figura 3: Resultado y variables de la ecuación del modelo cúbico
matemáticamente confiable. Con esta ecuación ya es posible poder trazar una gráfica y
determinar los datos que se acomodan al ajuste. En este caso, para poder obtener los
datos que se ven en la tabla 4, se tuvo que digitar la ecuación en la calculadora Casio
fx-9750GII, en la sección “tabla”, la cual nos brinda los datos del modelo que se utilizan
12
Tabla 4: La cantidad de sismos anuales en Chile usando el modelo cúbico. En la Commented [LHEG11]: Criterio B – (3): La
comunicación matemática es pertinente, adecuada y,
columna de la izquierda están los años que se tomaron para realizar la en su mayor parte, coherente, sería importante haber
tenido en cuenta métodos deductivos que den cuenta
de algunos de sus resultados, la notación utilizada ha
investigación, enumerados del 1 al 10, siendo el año 1, 2010 y el año 10, 2019. sido apropiada, definiendo la mayoría de términos
clave y las variables, y la exploración apoya su
Cantidad de sismos anuales en Chile usando el desarrollo de forma apropiada en múltiples formas de
representación matemática.
modelo cúbico
Año (2010 - 2019) Datos con el modelo
1 7350
2 7799
3 7641
4 7089
5 6357
6 5660
7 5212
8 5228
9 5923
10 7510
con el modelo cúbico en el plano cartesiano, en una gráfica que refleja los datos por
medio de una curva. En este caso el modelo, al ser comparado con los datos originales,
como podemos ver en la gráfica 1, donde el ajuste con el modelo cubico es una línea
curva punteada, deja ver que no se ajusta completamente a los datos. Por eso la
particularidad de la forma de la gráfica, porque, aunque hay puntos que llegan a ser
similares, hay datos muy diferentes como en el año 6, donde en la gráfica de los datos
originales hay 6637, con el modelo cúbico son 5660. Esto lo podemos ver a continuación.
13
Grafica 2: Cantidad de sismos anualmente en Chile con el modelo cúbico. Se
presenta una curva que inicia en las coordenadas (𝟏, 𝟕𝟑𝟓𝟎) y finaliza en
(𝟏𝟎, 𝟕𝟓𝟏𝟎).
minimicen las discrepancias entre los puntos y la curva, aunque estos no tengan una
completa relación con la realidad. La ecuación desarrollada con el modelo puede ser
comparada con los datos reales para comprobar su aproximación con la realidad.
Inicialmente, tras las diferencias visibles entre la línea de tendencia del modelo y los
investigación, dado a que no se aproxima completamente a la realidad, cuando la idea Commented [LHEG13]: Criterio D: Hace referencia a
los objetivos de la exploración en el análisis plateado
es que sea útil para poder determinar la cantidad de sismos anuales en Chile a futuro. frente a la validez del modelo.
Por eso en la gráfica 4, podemos ver de cerca la comparación entre el modelo cubico
que fue desarrollado previamente y los datos reales de la primera tabla. La línea en color
azul es la que contiene los datos reales, por eso podemos observar que hay algunos
picos como en el año 3, 6 y 10, que corresponden a los años con mayor movimiento
14
sísmico en Chile. Mientras que la línea de color rojo, es la del modelo cúbico, en esta
no se observan picos o desniveles entre cada punto, por el contrario los puntos forman
periodicidad del movimiento sísmico en Chile. En lugar de picos, la curva también refleja
puntos altos y bajos que no se alejan de los datos uniformes del modelo. Asimismo, la
comparación realizada nos deja ver que si hay puntos en común entre los datos reales
y el modelo cubico, así como lo podemos ver en los años 5 y 7, permitiéndonos ver la
corresponde a los datos reales y la curva en rojo es del modelo cúbico. El eje Y
Ahora, si nos planteamos la proyección a futuro que se quiere realizar con esta
investigación por medio del modelo, se deben tener en cuenta algunas limitaciones que
Commented [LHEG14]: Criterio D: Considera
hacen que este modelo cúbico no sea 100% cercano a la realidad. En primer lugar, la limitaciones de su modelo en torno al comportamiento
de los datos y el contraste con otros modelos.
15
forma curvilínea no se acomoda por completo a lo que realmente ocurre anualmente,
este modelo supone que la cantidad anual de sismos en Chile es en cierto modo algo
continuo y que no deben haber picos al comparar datos, algo que en la realidad no es
posible, ya que por el choque de las placas tectónicas, hay años en los cuales se
verificar que hayan otros ajustes, no solo el cubico, para garantizar la veracidad
2020 y 2021, plantearé los puntos a favor de este modelo. Por un lado, los datos son los Commented [LHEG15]: Criterio D: Discute los puntos
fuertes y débiles de su modelo, desde un enfoque
más cercanos a la realidad aunque la gráfica de una curva, y al ser una aproximación gráfico y numérico de aproximación.
inicial, que luego puede tener sus modificaciones, es lo suficientemente acertado como
para establecer una cantidad anual de sismos análoga a la realidad, lo que permite que
se tenga conciencia de los posibles sismos a lo largo de un año, sin importar que esto
correspondiente a ese año. Por el lado de 2021, también se debe reemplazar, pero por
16
𝑦 = 61777,728 − 74710,08 + 21079,2 + 6075,4
El único dato que se tiene hasta el día de hoy es la cantidad de sismos que hubo en
2020, puesto que es posible calcular el porcentaje de error de una de las 3 predicciones
debe al 𝑅² de 0,7676 porque no estuvo cercano a 1, para garantizar una mayor precisión
en las predicciones.
proceso del desarrollo del modelo cúbico para calcular la cantidad de sismos anuales
en Chile, me dejó ver la relación matemática que existe en todas las cosas que nos
17
rodean, me permitió comprender de una mejor manera los conceptos que he visto en
clase como la distribución normal, todo alrededor del desarrollo de ecuaciones que sean
y no subjetiva. Ya por el lado del desarrollo del modelo, encontré limitaciones no frente Commented [LHEG19]: Criterio D: Menciona
aspectos matemáticos que ha aprendido o en los
a los datos pero al ajuste del modelo, porque al tener un 𝑟 2 por debajo de 0,90, era cuales ha mejorado su comprensión a partir de la
exploración.
complicado comparar la forma de la función con los valores reales. Considero importante
poder encontrar otro modelo que no sea lineal, cuadrático, potencial, exponencial, etc.,
porque tienen menores 𝑟 2 del cúbico, pero si comprender las particularidades de las
ondas sísmicas dentro del movimiento sísmico para poder encontrar un modelo
matemático mucho más aproximado a la realidad. Aun así, el ajuste cúbico, permitió que
se identificaran cantidades aproximadas del número de sismos anuales; así como una
modelo y los datos de la vida real, que al ser análogos se pudo comprender el
todos. Al final se tuvieron resultados mayores a los de la realidad, lo que nos haría creer
que el desarrollo de todo el ajuste cúbico es erróneo, sin embargo, este aumento
18
de la densidad de la población en zonas de riesgo como lo es Chile. Esto nos ayuda a Commented [LHEG21]: Criterio D: Considera
diferentes perspectivas que pueden llegar a explicar el
comprender, que a pesar de que el 𝑟 2 no alcanzará a ser 1, si se aproximó a la realidad, comportamiento de modelo en el contexto de los
sismos en Chile.
pero a raíz del cálculo del porcentaje de error, pudimos descubrir que esta diferencia
entre los datos reales y las predicciones hacen que la ecuación desarrollada no sea
realidad, si bien se muestra un aumento de los sismos por año, no aumentan Commented [LHEG22]: Criterio C – (2): El tópico
elegido es del interés personal de la estudiante,
exponencialmente pero si explica matemáticamente el porqué del aumento de sismos realizando y comprobando predicciones de su modelo,
y considerando el tema desde otras perspectivas
(como la densidad de la población en zonas de riesgo),
por año en Chile. Esta predicción y aclaración resultan útiles para comprender el sin embargo, no fue evidente el pensamiento creativo.
como un caso de estudio, pero también en el resto del mundo, aun así, es pertinente
que se realice una investigación más extensiva donde el modelo se acerque cada vez
REFERENCIAS BIBLIOGÁFICAS
BBC. (2010). BBC News Mundo. Obtenido de
https://www.bbc.com/mundo/ciencia_tecnologia/2010/03/100303_terremotos_fr
ecuencia_men
Chapra, S., & Canale, R. (2006). Capitulo 17: Regresión de Mínimos Cuadrados. En S.
http://www.sismologia.cl/pdf/difusion/001_terremotos_y_sismicidad_chile.pdf
mas-8-mil-sismos-2017/
https://www.universoformulas.com/matematicas/analisis/derivadas-parciales/
19
Gonzalez, E. (n.d). Departamento de Geografía y Ordenación del Territorio UCML.
Obtenido de https://previa.uclm.es/profesorado/egcardenas/subduccion.htm
https://www.meganoticias.cl/nacional/324653-temblor-chile-aumento-temblores-
ano-2020-antartica-jrx07.html
https://www.significados.com/sismo/
https://www.oei.es/historico/divulgacioncientifica/noticias_566.htm
Uchile, D. (s.f.). 27/F: Los cambios que generó el terremoto de 2010. Diario Uchile,
Santiago.
https://www.ccg.unam.mx/~vinuesa/R4biosciences/docs/Tema9_regresion.html
20