Introducción al Modelado Continuo

Materia de la Licenciatura en Ciencia de Datos.
Departamento de Matemática, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires.

Julián Fernández Bonder · Notas de clase, notebooks que reproducen las figuras y laboratorios en Colab

Las notas

Un curso de un cuatrimestre en tres partes: dinámica no lineal (ecuaciones diferenciales ordinarias como modelos: equilibrios, estabilidad, ciclos límite, bifurcaciones), análisis de Fourier (series, transformada, DFT y sus aplicaciones al análisis de señales) y ecuaciones en derivadas parciales (difusión, Laplace y ondas, con el énfasis puesto en la convergencia al equilibrio). El objetivo es que quien lo termine pueda tomar un problema concreto, decidir si conviene modelarlo con ecuaciones diferenciales, construir el modelo, analizarlo y volver al problema con respuestas. Cada parte se organiza alrededor de uno o dos problemas conductores, planteados al comienzo y retomados cada vez que la teoría permite decir algo nuevo; los laboratorios los llevan hasta los datos.

Parte I · una epidemia y un circuito que oscila solo

El modelo SIR, con la epidemia de gripe de un internado inglés (1978) como dato, y el oscilador de van der Pol: umbrales, equilibrios, ciclos límite y bifurcaciones.

Parte II · una señal de temperatura

Tres años de temperatura horaria en Buenos Aires: qué frecuencias contiene, cómo separar los ciclos del ruido, qué se pierde al submuestrear y cuánto se puede comprimir.

Parte III · el suelo y la cuerda

La temperatura del suelo a varias profundidades (difusión con forzado periódico) y la cuerda de una guitarra (modos normales, espectro real, inarmonicidad).

Capítulos

Parte I · Dinámica no lineal

Capítulo 1

Introducción

Del problema al modelo: ¿cuándo alcanza una EDO? Los dos problemas conductores: la epidemia (SIR) y el circuito de van der Pol.

notas

Capítulo 2

Modelos poblacionales

Malthus y logística, la recta de fase, Lotka–Volterra con y sin capacidad de carga, ciclos límite, el modelo SIR y el número reproductivo básico.

notas · notebook

Capítulo 3

Sistemas mecánicos

Leyes de Newton, oscilador armónico, fricción, péndulo, circuitos eléctricos, reducción a sistemas de primer orden y la energía como herramienta.

notas · notebook

Capítulo 4

Flujo y diagramas de fase

Existencia y unicidad, el flujo de una EDO, órbitas y retratos de fase.

notas · notebook

Capítulo 5

El Teorema de Hartman–Grobman

Conjugación, clasificación de los sistemas lineales en el plano (traza y determinante), linealización en equilibrios hiperbólicos; vuelta a los modelos poblacionales.

notas · notebook

Capítulo 6

Sistemas conservativos y funciones de Lyapunov

Cantidades conservadas y potenciales, funciones de Lyapunov, el principio de LaSalle; el circuito de van der Pol.

notas · notebook

Capítulo 7

Dinámica global en dimensión dos

Conjuntos límite, el teorema de Poincaré–Bendixson, el estudio gráfico con nulclinas: competencia entre especies.

notas · notebook

Capítulo 8

Bifurcaciones

Estabilidad estructural, bifurcaciones silla-nodo, transcrítica y de horquilla, bifurcaciones en el plano, la bifurcación de Hopf (y el caso degenerado de van der Pol); la epidemia con nacimientos.

notas · notebook

Capítulo 9

Ejercicios de la Parte I

Ejercicios teóricos (con las series guiadas sobre los problemas conductores) y ejercicios de laboratorio.

notas · laboratorios

Parte II · Análisis de Fourier

Capítulo 10

Introducción

Motivación y el problema conductor: una serie temporal de temperatura.

notas

Capítulo 11

Series de Fourier

Ondas fundamentales, superposición, ortogonalidad, la serie de Fourier compleja y real.

notas · notebook

Capítulo 12

Convergencia de la serie de Fourier

Propiedades de los coeficientes, convergencia puntual (Dini), uniforme y en media cuadrática, el fenómeno de Gibbs.

notas · notebook

Capítulo 13

Transformada de Fourier

De las series a la transformada, propiedades, inversión, convolución y el teorema de Plancherel.

notas · notebook

Capítulo 14

La transformada de Fourier discreta

La DFT y su representación matricial, muestreo, la grilla de frecuencias y el aliasing, la transformada rápida (FFT).

notas · notebook

Capítulo 15

Aplicaciones del análisis de Fourier

Filtrado de ruido en señales, la señal de temperatura, compresión de imágenes.

notas · notebook

Capítulo 16

Ejercicios de la Parte II

Ejercicios teóricos y de laboratorio.

notas · laboratorios

Parte III · Ecuaciones en derivadas parciales

Capítulo 17

Introducción

Los dos problemas conductores: la temperatura del suelo y la cuerda vibrante.

notas

Capítulo 18

La ecuación de difusión

Conservación, movimiento browniano y paseo al azar, la ecuación del calor en un intervalo, en la recta y en dimensión mayor; convergencia al equilibrio y la desigualdad de Poincaré; la temperatura del suelo.

notas · notebook

Capítulo 19

La ecuación de Laplace/Poisson

Superficies mínimas, juegos aleatorios y valor esperado, principio del máximo y valor medio, solución fundamental y series de Fourier en un rectángulo.

notas · notebook

Capítulo 20

La ecuación de ondas

Tipos de ondas, la cuerda vibrante, energía, series de Fourier y modos normales, la fórmula de d'Alembert y las dos descripciones de la misma solución.

notas · notebook

Capítulo 21

Ejercicios de la Parte III

Ejercicios teóricos y de laboratorio.

notas · laboratorios

Notebooks del texto

Once notebooks de Python, uno por capítulo o par de capítulos, que reproducen las figuras de las notas (salvo unos pocos esquemas dibujados a mano y algunas figuras que salen de los scripts tools/fig_*.py) y permiten cambiar los parámetros de cada modelo. No requieren instalar nada: el botón abre una copia en Google Colab (hace falta una cuenta de Google para ejecutarla) y la primera celda instala el módulo imc del repositorio. Cada uno termina con una sección para experimentar.

Modelos poblacionales

Capítulo 2

Logística y recta de fase, Lotka–Volterra con nulclinas, el modelo de bacterias y su ciclo, el SIR en el plano (s, i).

Abrir en Colab ver en GitHub

Sistemas mecánicos y eléctricos

Capítulo 3

Oscilador armónico, fricción, péndulo, la corriente del circuito de van der Pol y su energía no monótona.

Abrir en Colab ver en GitHub

Sistemas lineales y Hartman–Grobman

Capítulos 4–5

Los retratos lineales (silla, nodo, foco, centro), el plano traza–determinante, linealización de los modelos poblacionales.

Abrir en Colab ver en GitHub

Lyapunov y dinámica global

Capítulos 6–7

Energía del péndulo, funciones de Lyapunov, van der Pol, Poincaré–Bendixson, competencia entre especies, Lorenz.

Abrir en Colab ver en GitHub

Bifurcaciones

Capítulo 8

Silla-nodo, transcrítica, horquilla, un ejemplo en el plano, Hopf genérica contra van der Pol, el SIR con nacimientos.

Abrir en Colab ver en GitHub

Series de Fourier y convergencia

Capítulos 11–12

Sumas parciales, núcleo de Dirichlet, decaimiento de los coeficientes, Gibbs, la serie de Fourier de la temperatura.

Abrir en Colab ver en GitHub

Transformada, DFT y FFT

Capítulos 13–14

Transformadas clásicas, la matriz de Fourier, aliasing, la grilla de frecuencias, Cooley–Tukey y su costo, la temperatura submuestreada.

Abrir en Colab ver en GitHub

Filtrado y compresión

Capítulo 15

Filtro por umbral, la temperatura regular y su residuo, compresión de la serie y de una imagen.

Abrir en Colab ver en GitHub

La ecuación de difusión

Capítulo 18

Paseo al azar, calor en un intervalo y en el cuadrado, tasa de decaimiento, la solución periódica del suelo y el ajuste de la difusividad.

Abrir en Colab ver en GitHub

Laplace y Poisson

Capítulo 19

Superficies mínimas, el juego aleatorio y Monte Carlo, valor medio, principio del máximo, solución fundamental, series en el rectángulo.

Abrir en Colab ver en GitHub

La ecuación de ondas

Capítulo 20

Catálogo de ondas, modos de la cuerda, energía, d'Alembert, dominios de dependencia, reflexión de un pulso.

Abrir en Colab ver en GitHub

Laboratorios

Ocho laboratorios, uno por problema conductor y uno por cada bloque de métodos numéricos (que se explican acá y no en las notas). Cada uno viene en dos versiones generadas por el mismo script: la de estudiantes, con la explicación de los métodos, las consignas, esqueletos con # TODO y celdas de verificación, y la guía del docente, con las soluciones completas ejecutadas y respuestas modelo. Todos terminan con preguntas de interpretación escrita, que es lo que se evalúa. Cada uno está pensado para una sesión de laboratorio de cuatro horas (los más largos, para dos).

Resolución numérica de EDOs

Parte I · métodos

Euler desde cero, orden de convergencia, estabilidad y rigidez, Runge–Kutta, sistemas y energía, eventos con solve_ivp, SIR y SEIR; opcional: un termostato con histéresis.

EstudiantesGuía del docente

Ajuste del modelo SIR a datos

Parte I · epidemia

La gripe del internado de 1978: mínimos cuadrados, identificabilidad, validación con media serie, SEIR, vacunación y aislamiento; opcional: la logística y la población de varios países.

EstudiantesGuía del docente

La forma de onda del circuito

Parte I · van der Pol

Régimen permanente, período y amplitud con eventos, espectro, y el problema inverso: estimar λ a partir de un registro con ruido.

EstudiantesGuía del docente

La señal de temperatura

Parte II · señal

Espectro de datos reales, ventanas, amplitud y fase de los ciclos, aliasing, filtrado y promedios móviles, compresión, espectrograma.

EstudiantesGuía del docente

Diferencias finitas para EDPs

Parte III · métodos

Matrices de diferenciación, Poisson 1D, capa límite, Neumann, calor explícito e implícito (r ≤ 1/2), leapfrog y CFL.

EstudiantesGuía del docente

La temperatura del suelo

Parte III · difusión

Amplitud y fase por capa, cuatro estimaciones de la difusividad, simulación con la superficie como borde, residuos, la cañería.

EstudiantesGuía del docente

El espectro de una cuerda real

Parte III · ondas

Fundamental y armónicos de una grabación, inarmonicidad por rigidez, decaimiento por armónico, timbre, simulación con leapfrog y fricción.

EstudiantesGuía del docente

Laplace en 2D e inpainting

Parte III · Laplace

Cinco puntos contra Monte Carlo, costo y error, métodos iterativos, el juego con fuente, rellenar regiones de una imagen.

EstudiantesGuía del docente

Datos

Los notebooks descargan los datos solos desde el repositorio (datos.obtener("archivo")). Fuentes y licencias:

ArchivoContenidoFuente y licencia
gripe_internado_1978.csvAlumnos en cama por día durante una epidemia de gripe en un internado inglés (763 alumnos, 14 días).Anónimo, Influenza in a boarding school, BMJ, 1978; reproducido en Murray, Mathematical Biology I.
temperatura_horaria.csvTemperatura del aire cada hora, 2021–2023, en el punto de Aeroparque (Buenos Aires).Reanálisis ERA5 vía Open-Meteo, CC BY 4.0 (Copernicus). Reanálisis, no medición directa.
suelo_horaria.csvTemperatura del aire y del suelo en cuatro capas (0–7, 7–28, 28–100, 100–255 cm), 2023, mismo punto.Reanálisis ERA5-Land vía Open-Meteo, CC BY 4.0. Promedios por capa.
cuerda_guitarra.wavCuerda de guitarra (Sol) pulsada, 6 segundos, mono 44,1 kHz.Freesound, sonido 8403 de speedy, CC0.
vanderpol_registro.csvCorriente simulada con ruido de un circuito de van der Pol con λ desconocido.Propia (generada por el laboratorio), CC BY 4.0.
poblacion_banco_mundial.csvPoblación total de Argentina, Nigeria, Japón, Corea del Sur, China y el mundo, 1960–2021.Banco Mundial, World Development Indicators (SP.POP.TOTL), CC BY 4.0.
astronauta.pngImagen en escala de grises para el ejemplo de compresión y el inpainting.NASA, dominio público (distribuida con scikit-image).

Cómo usar el material

Para estudiar por cuenta propia, el orden natural es el de los capítulos: leer el capítulo, correr el notebook del texto cambiando parámetros, hacer los ejercicios teóricos (están al final de cada parte, con la misma numeración que en las notas) y, cuando toque, el laboratorio correspondiente. Los notebooks del texto son cortos y sirven para ver: reproducen cada figura de las notas desde una celda marcada con su nombre. Los laboratorios son para hacer: traen la explicación de lo que hace falta, consignas con lo que se espera y celdas para completar; la guía del docente tiene las soluciones. Todo el código común vive en el módulo imc del repositorio (retratos de fase, espectros, esquemas numéricos, acceso a los datos), que se instala solo en Colab.

Correcciones y sugerencias: abrir un issue o un pull request en el repositorio, o escribir a jfbonder@dm.uba.ar.