Open ghost opened 7 years ago
Mi implementacion funciona pero es algo lenta me parece.
En el archivo proy.ipynb hay un pequeño ejemplo comparándolo con los métodos de la librería ODE
Hola
Estoy de acuerdo con lo de estandarizar las funciones, así como comentarlas. Lo de los temas no tengo muy claro cuáles serían. Por ahora lo que hice fue modificar el archivo proy.ipynb, (lo subí como proyv2.ipynb ) . Intenté hacer la gráfica en 3d, pero se necesita un único valor para la solución, mientras que el método devuelve un arreglo de dos elementos.
Te refieres a el metodo de integradorT2?
devuelve el arreglo de tiempos y uno de puntos, utilizo
x1 = map(y -> y[1], sol1)
y1 = map(y -> y[2], sol1);
Lo que hago es ejecutar la línea: (t2, sol2) = integradorT2(ec1, t0, tf, x̄, 10000, 1e-50); luego al ver cómo es sol2, se obtiene un arreglo de arreglos, por lo que sol2[1] es igual al arreglo para la x inicial. Eso es correcto, sin embargo para realizar la gráfica se necesita un valor único para la altura z (en este caso, nombré a esta altura como "solucion"), dicho valor, debe extraerse de algún modo del arreglo de arreglos, pero no sé cómo.
creo que no entiendo, que quieres graficar en 3d?
por lo que sol2[1] es igual al arreglo para la x inicial.
aqui hay un errorsillo,
(t2, sol2) = integradorT2(ec1, t0, tf, x̄, 10000, 1e-50);
entrega t2 = [tiempos] y sol2 = [x1,y1; x2, y2;....]
sol2 es una matriz de 2 por n, puedes sacar cada columna como
x = map(y -> y[1], sol2)
y = map(y -> y[2], sol2);
Ya veo. Lo que trato de hacer es obtener la gráfica (figura 15.1) de la página 377.
ok, ya te entendi, esa grafica es de la y minima en el arreglo mientras varias alpha y beta no? mas bien creo que seria crear una matriz de n*n y para cada valor de alpha y beta ir guardando la y minima en el arreglo y despues plotear este
si graficas como tal la solucion de el integrador te dara una trayectoria.
vale, trabajo en ello.
no me sale, estoy un tanto perdido en cómo hallar la y minima
Hector, en el espacio fase del archivo Proyv2.ipynb celda ln[6] no obtuve una gráfica cerrada en el espacio fase. ¿Cuál será el codigo correcto?
¿para beta=0?
Rubén comparó su método con uno de julia usando el paquete "ODE", ambos dan resultados casi iguales. Sin embargo, veré si se puede obtener y_min con ambos.
Tenemos planeada algunos temas no? como para establecerlos como puntos de trabajo y tal vez repartirlos entre nosotros, creo que así seria mas eficiente el asunto.
También creo que deberíamos estandarizar las funciones que utilizaremos.
PD: Siento reportarme tan noche, le estaré dando a esto toda la noche.