Відео

¡Muchas gracias por tu retroalimentación! =D

En breve subo uno con Octave, un software muy ameno y libre

¡Hola! Una forma en la que puedes graficar el plano fase es extrayendo los valores para diferentes tiempos mediante la función Table, luego utilizando la función ListPlot. Enseguida te dejo el código que sirve para el ejemplo del video. Muchas gracias por tu pregunta. CÓDIGO TPP = Table[{x[t] /. sol[[1]], y[t] /. sol[[1]]}, {t, ti, tf, 0.01}]; (*Aquí se tabula y se guarda en TTP*) ListPlot[TPP, AxesLabel -> {"x", "y"}] (*Aquí se grafican los pares ordenados guardados en TTP*)

Hola Liz 😊 En este otro video te mostramos como implementar Runge Kutta de 4to orden. ua-cam.com/video/qw7onY6YRG0/v-deo.html ¿Te refieres al ode45? En el blog que aquí te dejo puedes encontrar una buena discusión sobre el método utilizado. En particular... ode45 is a six-stage, fifth-order, Runge-Kutta method. blogs.mathworks.com/cleve/2014/05/26/ordinary-differential-equation-solvers-ode23-and-ode45/#:~:text=ode45%20is%20a%20six%2Dstage,order%2C%20Runge%2DKutta%20method.

Que bueno iniciar un curso de álgebra, una de las ramas más importante de la matemática. Tienes el potencial de enseñar 👍

Claro, los estaré haciendo con regularidad antes del examen... Gracias por comentar, like y comparte! 😃🙌

Ignora mi comentario de dónde está la función. Ya vi que ODE45 es de MATLAB. Saludos.

¡Muchas gracias!

No comprendo bien tu pregunta. Si tienes datos en Excel y los quieres abrir en Matlab, lo que puedes hacer es dar clic en... HOME -> Import Data -> Eliges el archivo -> Seleccionas lo que quieres importar -> Import Selection. Ahora bien, si lo que tienes son datos en formato hh:mm:ss, recomendaría que los pasaras a una sola unidad y así tratarlos más fácilmente. Si posteriormente vas a graficar, puedes elegir una escala para que la presentación sea elocuente. Si es necesario ese formato (hh:mm:ss), te dejo un link con la función que te permitirá regresar a dicho formato. www.mathworks.com/matlabcentral/answers/375392-convert-seconds-to-time Saludos.

Ya que se resuelven numéricamente, lo que se obtiene es la matriz de dos columnas que contiene los valores de y1 y y2 en el intervalo de integración. Con éste método no se obtienen las funciones explícitas y1(t), y2(t).

Hola Me parece que no esta bien el problema como lo has escrito, por favor, corrobora si es: dx/dt = 0.08x ; x(0)=0.5.

Son dos ejercicios por separado, ¿cierto? Ambos ejercicios se pueden resolver en forma analítica y son relativamente sencillos, si los quires resolver numéricamente te dejo aquí el código para el primero: clc; clear; %Tiempo integración xspan = (0:0.1:1); %Condiciones Iniciales y0 = [0]; %Solver [x,y] = ode45(@miODE,xspan,y0); %Graficar plot(x,y(:,1),'r'); hold on; legend('y'); xlabel('tiempo'); ylabel('y'); %EDO function dy = miODE(x,y) dy(1) = 1 + y(1)^2; dy = [dy(1)]; end

@@matesenminutos8351 Sí, son 2 ejercicios pero no los enumere. Disculpa. Ya corrí el programa y me da la gráfica con las solución. Muchísimas gracias, llevaba atorado en esa parte varios días y no entendía la sintaxis cómo lo planteé porque normalmente era dy/dt.

Tienes que guardar la función en otro archivo Guarda lo siguiente como miODE.m %EDO function dy = miODE(t,y) dy(1) = y(2)+1; dy(2) = y(1)*y(2)-6; dy = [dy(1);dy(2)]; end Luego, En otro script que guardes en la misma carpeta pero con otro nombre, podría ser Ejemplo.m escribe lo siguiente clc; clear; %Tiempo integración tspan = (0:0.1:5); %Condiciones Iniciales y0 = [1,0]; %Solver [t,y] = ode45(@miODE,tspan,y0); %Graficar plot(t,y(:,1),'r'); hold on; plot(t,y(:,2),'b'); legend('y1','y2'); xlabel('tiempo'); ylabel('y');

@@matesenminutos8351 excelente, más tarde lo intento. Muchas gracias de antemano.

¿Te funcionó?

Hola hola, da doble clic en la variable que quieras ver en el workspace y vas a poder observar la tabla con los valores

Hola hola, dale clic sobre el punto cuyo valor quieras observar (sobre la gráfica) y aparecen los valores de x y y

Gracias por avisar sobre el link, ya deberías poder ver el PDF con el código

Sería: yExacta = (1/4.4817)*exp((3/2)*x.^2); err = (abs(yExacta-yRK4)./yExacta).*100; En nuestra página web matesenminutos.wixsite.com/memin , subí un PDF llamado "RK4 error relativo porcentual" donde puedes ver el código para calcularlo e incluye la gráfica del error relativo porcentual para este ejercicio.

@@matesenminutos8351 muchas gracias

Hola, el valor vi en el código es el valor de x al tiempo cero, de forma que tendrías que poner vi=0.

@@matesenminutos8351 muchas gracias

¡Excelente! Próximamente subiremos más material

De nada, espero te haya ayudado :)