Ecuaciones Diferenciales

Descargar como pptx, pdf o txt
Descargar como pptx, pdf o txt
Está en la página 1de 20

MÉTODOS

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

El método de Euler consiste en encontrar iterativamente la solución de una ecuación


diferencial de primer orden, utilizando como punto de inicio, un valor inicial, un paso h y un
intervalo, para lo cual se utiliza la siguiente formula.

𝑦 𝑖 +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:

Encontrar la solución de la siguientes ecuaciones diferencial con el método de Euler.

𝑦 ′ =2 𝑥 − 3 𝑦 +1 𝑥 𝜖 [ 1 ; 1.5 ] 𝑦 0=5 h= 0.1 𝑖𝑛𝑑𝑖𝑐𝑎𝑟 𝑡𝑜𝑑𝑜𝑠 𝑙𝑜𝑠 𝑣𝑎𝑙𝑜𝑟𝑒𝑠

𝑦 ′ = 𝑦𝐶𝑜𝑠 ( 𝑥 ) − 𝑥 𝑥 ∈ [ − 𝜋 ; 𝜋 ] 𝑦 0= 0 𝑛=10 𝑖𝑛𝑑𝑖𝑐𝑎𝑟 𝑙𝑜𝑠 𝑣𝑎𝑙𝑜𝑟𝑒𝑠 𝑝𝑎𝑟𝑎 {− 𝜋 , 0 , 𝜋 }

𝑦 ′ = ( 1− 2 𝑥 ) 𝑦 𝑥 ∈ [ 0 ; 3 ] 𝑦 ( 0 )=1 𝑛= 16 𝑖𝑛𝑑𝑖𝑐𝑎𝑟 𝑙𝑜𝑠 𝑣𝑎𝑙𝑜𝑟𝑒𝑠 𝑝𝑎𝑟𝑎 { 0,0.75,2 .25,3 }


Método de Runge-Kuta

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)

h=(b-a)/N; %tamaño de paso


t=a:h:b; % crea el vector como fila
t=t(:); % : 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=h*feval(f,t(k),y(k,:))';
k2=h*feval(f,t(k)+h/2,y(k,:)+k1/2)';
k3=h*feval(f,t(k)+h/2,y(k,:)+k2/2)';
k4=h*feval(f,t(k)+h,y(k,:)+k3)';
y(k+1,:)=y(k,:)+(k1+2*k2+2*k3+k4)/6;
end
Ejemplos:

Encontrar la solución de la siguientes ecuaciones diferencial con el Runge Kutta, .

𝑦 ′ =2 𝑥 − 3 𝑦 +1 𝑥 𝜖 [ 1 ; 1.5 ] 𝑦 0=5 h= 0.1 𝑖𝑛𝑑𝑖𝑐𝑎𝑟 𝑡𝑜𝑑𝑜𝑠 𝑙𝑜𝑠 𝑣𝑎𝑙𝑜𝑟𝑒𝑠

𝑦 ′ = 𝑦𝑆𝑖𝑛 ( 𝑥 ) − 𝑥 𝑥 ∈ [ − 𝜋 ; 𝜋 ] 𝑦 0= 0 𝑛=10 𝑖𝑛𝑑𝑖𝑐𝑎𝑟 𝑙𝑜𝑠 𝑣𝑎𝑙𝑜𝑟𝑒𝑠 𝑝𝑎𝑟𝑎 {− 𝜋 , 0 , 𝜋 }

−𝑦 𝑒𝑥
𝑥 ∈ [ 0 ;5 ] 𝑦 ( 0 )= −1𝑛=10 𝑖𝑛𝑑𝑖𝑐𝑎𝑟 𝑙𝑜𝑠 𝑣𝑎𝑙𝑜𝑟𝑒𝑠 𝑝𝑎𝑟𝑎 {0,1,4.5,5 }

𝑦 = + 3
3 3𝑦
Sistema de Ecuaciones Diferenciales de primer orden

Un sistema de ecuaciones diferenciales es un conjunto de varias ecuaciones diferenciales con


varias funciones incógnitas y un conjunto de condiciones de contorno. Una solución del mismo
es un conjunto de funciones diferenciables que satisfacen todas y cada una de las ecuaciones
del sistema, en esta caso todas las ecuaciones son de primer orden

𝑑 𝑦1
=𝐹 1 (𝑡 , 𝑦 1 , 𝑦 2 , … … , 𝑦 𝑛)
𝑑𝑡

𝑑 𝑦2
=𝐹 2 (𝑡 , 𝑦 1 , 𝑦 2 ,… … , 𝑦 𝑛 )
𝑑𝑡

𝑑 𝑦𝑛
=𝐹 𝑛 (𝑡 , 𝑦 1 , 𝑦 2 , … … , 𝑦 𝑛 )
𝑑𝑡
Ejemplos:

Del siguiente sistema de ecuaciones:

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

Del siguiente sistema de ecuaciones:

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

Encontrar la solución de la siguientes ecuación diferencial

𝑦′ ( 𝑥 ) =𝑦2 − 𝑦3

𝑥 𝜖 [ 0 ; 200 ] 𝑦 0= 0.01h=0.01 𝑖𝑛𝑑𝑖𝑐𝑎𝑟 𝑙𝑜𝑠 𝑣𝑎𝑙𝑜𝑟𝑒𝑠 𝑝𝑎𝑟𝑎 { 0 , 3 0 ,65,200 }

Del siguiente sistema de ecuaciones:

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 )

A la función se le llama función incógnita.


Ejemplos:

Encontrar la solución de la siguientes ecuación diferencial:


𝑦 ′ ′ =− 𝑦 −2 sin ⁡(𝑥 )

Con las condiciones iniciales , en el intervalo , con un paso


Indicar los valores de para

Encontrar la solución de la siguientes ecuación diferencial:


𝑦 ′ ′ =− 2 𝑦 ′ + 𝑥 2 𝑦

Con las condiciones iniciales , en el intervalo , con un paso


Indicar los valores de para
Encontrar la solución de la siguientes ecuaciones diferenciales:
𝑦 ′ ′ ′ =−2 𝑥𝑦 ´ ´ − 𝑦 ′ − 𝑦+ 𝑥𝑦

Con las condiciones iniciales , en el intervalo , con un paso Indicar los valores de para

𝑦 ′ ′ ′ =−12 𝑦 ´ + 𝑦 +1− 𝑥 [ cos ( 2 𝑥 )+16 sin ⁡(2 𝑥)]


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

Un sistema de ecuaciones diferenciales es un conjunto de varias ecuaciones diferenciales con


varias funciones incógnitas y un conjunto de condiciones de contorno. Una solución del mismo
es un conjunto de funciones diferenciables que satisfacen todas y cada una de las ecuaciones
del sistema, en esta caso todas las ecuaciones son de orden superior

𝑦 1𝑛 = 𝑓 (𝑥 , 𝑦 1 , 𝑦 2 … … 𝑦 𝑛 , 𝑦 ′ 1 , 𝑦 ′2 , … … . 𝑦 ′ 𝑛 , 𝑦 1𝑛− 1 , 𝑦 2𝑛 −1 … … .. 𝑦 𝑛𝑛− 1)

𝑦 2𝑛 = 𝑓 ( 𝑥 , 𝑦 1 , 𝑦 2 … … 𝑦 𝑛 , 𝑦 ′ 1 , 𝑦 ′ 2 , … … . 𝑦 ′𝑛 , 𝑦 1𝑛 −1 , 𝑦 2𝑛−1 … … .. 𝑦 𝑛𝑛− 1)

𝑦 3𝑛 = 𝑓 ( 𝑥 , 𝑦 1 , 𝑦 2 … … 𝑦 𝑛 , 𝑦 ′ 1 , 𝑦 ′ 2 ,… … . 𝑦 ′𝑛 , 𝑦 1𝑛−1 , 𝑦 2𝑛−1 … … .. 𝑦 𝑛𝑛 −1)




𝑦 𝑛 𝑛= 𝑓 (𝑥 , 𝑦 1 , 𝑦 2 … … 𝑦 𝑛 , 𝑦 ′ 1 , 𝑦 ′ 2 , … … . 𝑦 ′ 𝑛 , 𝑦 1𝑛 −1 , 𝑦 2𝑛−1 … … .. 𝑦 𝑛𝑛 −1)
Encontrar la solución de los siguientes sistemas de ecuaciones diferenciales:

{
′′ 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

Con las condiciones iniciales , en el intervalo , con un paso


Indicar los valores de para

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

También podría gustarte