Ecuaciones Diferenciales
Ecuaciones Diferenciales
Ecuaciones Diferenciales
NUMÉRICOS
Ing. Christian Flores
ECUACIONES DIFERENCIALES
Una ecuación diferencial es una expresión algebraica que contiene las derivadas de una o mas
variables dependientes con respecto a una o mas variables independientes. El objetivo de
resolver las ecuaciones diferenciales es hallar la función primitiva involucrada que depende de
la variable independiente. A continuación se muestran diferentes formas de representación de
las ecuaciones diferenciales.
Ejemplos:
𝑑𝑦
=3 𝑥 +5
𝑑𝑥
𝑥 𝑦 ′ =8 𝑦 +12
𝑦 ′ ′ ′ +2 𝑦 ′ ′ + 𝑦 ′ −5=0
Ecuaciones diferenciales ordinarias de primer grado
Las ecuaciones ordinarias se caracterizan por tener una sola variable independiente.
Método de Euler
𝑦 𝑖 +1= 𝑦 𝑖 + 𝑓 ( 𝑥 𝑖 , 𝑦 𝑖 ) h
𝑏− 𝑎
h=
𝑛
function [x,y]=Euler_orden1(f,a,b,y0,N)
%Método de Euler de PVI orden 1 con M sistemas de ecuaciones
h=(b-a)/N;
x=a:h:b; % crea el vector como fila
x=x(:); % : cambia de fila a columna para trabajar
y=zeros(N+1,length(y0));
y(1,:)=y0; % primera posición del vector
for k=1:N
K1=feval(f,x(k),y(k,:))';
y(k+1,:)=y(k,:)+h*K1;
end
Ejemplos:
Los métodos de Runge Kuta, tienen la exactitud de un esquema de una serie de Taylor, sin
necesidad del calculo de derivas superiores. Se ajusta a la siguiente forma general de
ecuación:
1
𝑦 𝑖 +1= 𝑦 𝑖 + ( 𝑘 +2 𝑘2+2 𝑘3 +𝑘4 ) h
6 1
𝑘 1= 𝑓 ( 𝑥 𝑖 , 𝑦 𝑖 )
1 1
𝑘2= 𝑓 ( 𝑥 𝑖 + h , 𝑦 𝑖 + h 𝑘1 )
2 2
1 1
𝑘3= 𝑓 (𝑥 𝑖 + h , 𝑦 𝑖 + h 𝑘 2)
2 2
𝑘4 = 𝑓 ( 𝑥𝑖 + h , 𝑦 𝑖 +h 𝑘 3 )
function [t,y]=RK_orden4_NM(f,a,b,y0,N)
−𝑦 𝑒𝑥
𝑥 ∈ [ 0 ;5 ] 𝑦 ( 0 )= −1𝑛=10 𝑖𝑛𝑑𝑖𝑐𝑎𝑟 𝑙𝑜𝑠 𝑣𝑎𝑙𝑜𝑟𝑒𝑠 𝑝𝑎𝑟𝑎 {0,1,4.5,5 }
′
𝑦 = + 3
3 3𝑦
Sistema de Ecuaciones Diferenciales de primer orden
𝑑 𝑦1
=𝐹 1 (𝑡 , 𝑦 1 , 𝑦 2 , … … , 𝑦 𝑛)
𝑑𝑡
𝑑 𝑦2
=𝐹 2 (𝑡 , 𝑦 1 , 𝑦 2 ,… … , 𝑦 𝑛 )
𝑑𝑡
⋮
𝑑 𝑦𝑛
=𝐹 𝑛 (𝑡 , 𝑦 1 , 𝑦 2 , … … , 𝑦 𝑛 )
𝑑𝑡
Ejemplos:
Resolver el sistema para con un h=0.1.Indicar los valores de x(t) y(t) para t=[0, 0.2, 0.4, 0.6,
0.8, 1], Las condiciones son: x(0)=1;y(0)=1
Resolver el sistema para con un h=0.1.Indicar los valores de x(t) y(t) z(t) para t=[0, 0.6, 1, 1.4,
1.8, 2], Las condiciones son: x(0)=1;y(0)=0;z(0)=1
LECCIÓN
𝑦′ ( 𝑥 ) =𝑦2 − 𝑦3
Resolver el sistema para con un h=0.1.Indicar los valores de x(t) y(t) z(t) para t=[0, 4, 5, 11.3,
11.8, 20], Las condiciones son: x(0)=10;y(0)=0;z(0)=1
Ecuaciones diferenciales de orden superior
Una ecuación diferencial ordinaria de orden es una ecuación que liga la variable independiente
, una función incógnita y sus derivadas sucesivas es decir, es una expresión, bien de la forma:
𝑦 𝑛 = 𝑓 ( 𝑥 , 𝑦 , 𝑦 ′ , 𝑦 ′ ′ , … … . , 𝑦 𝑛−1 )
Con las condiciones iniciales , en el intervalo , con un paso Indicar los valores de para
𝑦 ′ ′ ′ =−6 𝑦 4
Con las condiciones iniciales , en el intervalo , con un paso Indicar los valores de para
′′′ 36 60 5
𝑦 =− − +3( 𝑥 − 1)
𝑦 ( 𝑥 −1 )6
Con las condiciones iniciales , en el intervalo , con un paso Indicar los valores de para
Sistema de Ecuaciones Diferenciales de orden superior
𝑦 1𝑛 = 𝑓 (𝑥 , 𝑦 1 , 𝑦 2 … … 𝑦 𝑛 , 𝑦 ′ 1 , 𝑦 ′2 , … … . 𝑦 ′ 𝑛 , 𝑦 1𝑛− 1 , 𝑦 2𝑛 −1 … … .. 𝑦 𝑛𝑛− 1)
𝑦 2𝑛 = 𝑓 ( 𝑥 , 𝑦 1 , 𝑦 2 … … 𝑦 𝑛 , 𝑦 ′ 1 , 𝑦 ′ 2 , … … . 𝑦 ′𝑛 , 𝑦 1𝑛 −1 , 𝑦 2𝑛−1 … … .. 𝑦 𝑛𝑛− 1)
{
′′ 1 5 3
𝑢 + 𝑥𝑢 −𝑣 ´= 3
( 2− 3 𝑥 − 4 𝑥 )
𝑥
𝑣 ´ ´ − 𝑣=𝑥 (1− 𝑥2 )
Con las condiciones iniciales , en el intervalo , con un paso Indicar los valores de para
{
′′ ′ 2
𝑢 − 𝑢 +𝑣= 3 𝑡 +1
𝑣 ′ ′ +𝑢+𝑢 ′ = 4
Con las condiciones iniciales , en el intervalo , con un paso Indicar los valores de para
{𝑢′ ′ − 𝑢′ +𝑣 +4 =3 𝑡 2 +𝑒 𝑡 (2+ 4 𝑡 )
′′ ′
𝑣 +𝑢 − 𝑢 =6 t −2 t e
𝑡
Con las condiciones iniciales , en el intervalo , con un paso Indicar los valores de para
Encontrar la solución de la siguientes ecuaciones diferenciales:
𝑦 (9 𝑥 −15)
𝑦′ ′=
𝑥−1
Sea u(r) el potencial electrostático entre dos esferas metálicas concéntricas de radios R y 2R.
El potencial entre las dos esferas viene determinado por la expresión.
′′ 2 ′
𝑢 + 𝑢 =0
𝑟
Con las condiciones iniciales , en el intervalo , con un paso
Indicar los valores de para rTomar
Dos masas sujetas a resortes (sin masa) con constantes y , donde y es el desplazamiento. Determinar el
desplazamiento 1 y 2 sabiendo que:
¿
Con las condiciones iniciales , en el intervalo , con Indicar los valores de para t