Blox - Donde la lógica se vuelve natural

Nombre

graficar_lotka_volterra - Grafica la evolucion temporal de dos especies en una simulacion Lotka-Volterra.

Sinopsis

ok = graficos.graficar_lotka_volterra(trayect, titulo, archivo)

Descripcion

Recibe una trayectoria Lotka-Volterra y construye dos series contra el tiempo.

Formato general:

ejemplo.blox

TABLA trayect

trayect[1] = { t = ts, x = presas, y = depredadores }

ok = graficos.graficar_lotka_volterra(trayect, "Lotka-Volterra", "salida.png")

Se espera una tabla trayect con trayect[1].t, trayect[1].x y trayect[1].y.

Comportamiento

Interpreta t como tiempo, x como una especie y y como la otra.
Construye las dos series necesarias para la comparacion temporal.
Usa la infraestructura general de graficacion del modulo.

Funciones API relacionadas

graficos.graficar_fase_lotka_volterra()
graficos.graficar_fases_lotka_volterra()

Parámetros

Parámetro	Tipo	Descripción
trayect	TABLA	TABLA de trayectoria con `trayect[1].t`, `trayect[1].x` y `trayect[1].y` como vectores numericos. Normalmente esta estructura proviene de `numericos.lotka_volterra()`, que modela la dinamica presa-depredador mediante integracion numerica.
titulo	STRING	Titulo visible del grafico.
archivo	STRING	Nombre o ruta del archivo de entrada o salida.

Valor de retorno

VOID - No retorna valor. Genera el grafico y escribe el archivo PNG indicado.

Ejemplo

ejemplo.blox

INCLUIR "numericos.api"
INCLUIR "matrices.api"
INCLUIR "tablas.api"
INCLUIR "vectores.api"
INCLUIR "graficos.api"

FUNCION PRINCIPAL
INICIO "Lotka-Volterra (RK4)"

    STRING f, g
    NUMERICO t0, x0, y0, h, n
    VECTOR T, X, Y

    // Definimos la EDO: dy/dt = f(x(t), y(t), t)   dy/dt = g(x(t), y(t), t)
    // En este caso dy_f/dt = a1x - a2xy
    // En este caso dy_g/dt = -b1y + b2xy
    //

    f = "0.4*x-0.018*x*y"     // x(t) es la presa
    g = "-0.8*y+0.023*x*y"    // y(t) es el predador 

    n  = 220
    h = 0.1

    TABLA trayect
    t0 = 0   x0 = 30   y0 = 4
    trayect = numericos.lotka_volterra(f, g, t0, x0, y0, h, n)
    graficos.graficar_lotka_volterra(trayect, "Lotka-Volterra", "PROG41 - GRAFICO - XY - presa-predador.png")
    graficos.graficar_fase_lotka_volterra(trayect, "Lotka-Volterra - Diagrama de fase", "PROG41 - GRAFICO - FASE - presa-predador.png")


FINAL```

graficar_lotka_volterra(trayect,titulo,archivo)\texttt{graficar\_lotka\_volterra}(trayect, titulo, archivo)graficar_lotka_volterra(trayect,titulo,archivo)

Parámetros

Valor de retorno

Ejemplo

$\texttt{graficar\_lotka\_volterra}(trayect, titulo, archivo)$