Ecuaciones Diferenciales Guia 3
Ecuaciones Diferenciales Guia 3
Ecuaciones Diferenciales Guia 3
TUTOR
BOGOTÁ, D.C,
1
TABLA DE CONTENIDOS
CONTENIDO
1.Introducción…………………………………………………………..3
3.Conclusiones………………………………………………………… 16
4.Referente bibliográfico……………………………………………….17
2
INTRODUCCION
En esta guía de trabajo se estudia los diferentes métodos que tienen las
Ecuaciones diferenciales ordinarias de orden superior, su forma analítica y grafica
de su solución, también se profundizara con el concepto de ODE45 en Matlab la
cual mojara las capacidades para utilizar las técnicas, habilidades y herramientas
modernas.
3
GUÍA 3:
1. Actividad No. 2.
1.2 reproduzca las gráficas con el código de 0DE45 y obtenga otras graficas
variando los valores de parámetros
x=2 x+ y ; y =x−2 y
4
nos muestra recorrido como de silla y la solución de equilibrio en esta
ecuación y(t)=0 y x(t)=0.8, a continuación, algunas graficas modificando
unos parámetros del código, empezamos haciendo modificaciones a la
línea 8 del editor de Matlab.
Sys=inline (‘[4*x(1)+x(2)+x(2);x(1)-4*x(2)]´,´t´,´x´):
for x0=-1:1.5:1
for y0=-1:1.5:1
[ts,xs] = ode45(sys,[0 5],[x0 y0]);
plot(xs(:,1),xs(:,2))
end
end
for x0=-1:1.5:1
for y0=-1:1.5:1
[ts,xs] = ode45(sys,[0 -5],[x0 y0]);
plot(xs(:,1),xs(:,2))
end
5
figura: 3 modificación 2 ODE45
sys=inline('[8*x(1)+x(2);x(1)-7*x(2)]','t', 'x');
%warning off MATLAB: divideByZero
vectorfield(sys,-3:.25:3,-3:.25:3)
hold on
for x0=-4:1.5:4
for y0=-4:1.5:4
[ts,xs] = ode45(sys,[0 8],[x0 y0]);
plot(xs(:,1),xs(:,2))
end
end
for x0=-4:1.5:4
for y0=-4:1.5:4
[ts,xs] = ode45(sys,[0 -8],[x0 y0]);
plot(xs(:,1),xs(:,2))
end
6
end
hold off
axis([-8 8 -8 8])
fsize=15;
set(gca,'xtick',[-8:1:8],'FontSize',fsize)
set(gca,'ytick',[-8:1:8],'FontSize',fsize)
xlabel('x(t)','FontSize',fsize)
ylabel('y(t)','FontSize',fsize)
hold off
figura 4
Podemos concluir que al modificar los valores de una forma arbitraria la solución y
la grafica depende de los parámetros entre las líneas 12 y 18, donde todos los
datos de la Ecuación diferencial ingresados por medio de la ODE45 en el
programa de Matlab, los cuales al ingresar las ecuaciones salen las soluciones.
7
que contengan dos ecuaciones de primer orden con dos variables
dependientes que se apliquen a un área relacionada con su carrera y
desarrolle lo siguiente:
Obtenga la soluciones para cada una de las EDO a partir del método
de eliminación para sistemas con coeficientes constantes
2.1
ejercicio 1
Tenemos una ecuación de una red eléctrica básica, para las corrientes i,(t) e
i(t).
−R1 + R2 R2
d i1 L2 L2 i 1
()
dt i 2
=(
R2
)( )
−R2 i 2
L1 L1
Los parámetros para resolver son:
R1=8∩ , R 2=3 ∩, L1 =1h , L2=1 h ,
−8+ 3 3
i̇ 1 1 1 ) i̇1 −11 3 ) i̇ 1
()
i̇ 2
=(
3 ()
−3 i̇2
= ( i̇ )=(
()
3 −3 i̇ 2
1 1
ẏ=ay+ b
8
−14 + √ 196−96
λ 1= =−2
2
−14− √ 196−96
λ 2= =−12
2
3 9i −3 i2 9 −3 1
i 1=
| (−2 ) +11
3 | = 1
(−2 ) +3 −3 i 1 i2
= |
−3 1
=
3 |()
3 −1i 1 −3 i2 −1 −3
= 3
i 2=
| (−12 ) +11
3 |
=
(−12 ) +3 −3 i1 −9 i2
=
−3|−9 −1 |( )
AC =C1 1 e−2∗t +C 2 3 e−12∗t
() ( )
3 −1
9
figura 5 ejericcio1
Ejercicio 2
dx
=x+2 y
dt
dy
=4 x +3 y
dt
Primero determinaremos los valores propios y vectores propios de la matriz de
coeficientes en la ecuación característica
det ( A− AI ) = 1−λ 2
| 4 3−λ |
(1− λ) ( 3−λ )−8 = 0
10
λ 2−4 λ−5=0
( λ−5)(λ+1)=0
λ 1=5 ; λ2=−1
4 k 1−2 k 2=0
4 k 1=2 k 2
k 2=2 k 1
k 1= 1
()
2
1−(−1) 2 2 2
| 4
=
3−(−1) 4 4
=> || |
2 k 1 +2 k 2=0
4 k 1 +4 k 2=0
2 k 1=−2k 2
k 1=−k 2
k 2= −1
( )
1
x =( 1) . e 5t
1
2
x =( −1) . e −t
2
1
11
Concluimos que la solución general del sistema es:
x=c 1 1 .e 5 t +c 2 −1 . e−t
()
2 ( )
1
E igualamos los elementos, con los elementos correspondientes de los vectores:
y (t)=2 c1 e 5 t +c 2 e−t
Ejemplo 3
dx
=6 x− y
dx
dx
=6 x− y
dx
−b ± √ b2−4 ac
x=
2a
10 ± √ −16 10± 4 i
λ= = = 2¿¿
2 2
12
Ahora se sustituye λ 1 en la matriz
(1−2 i¿ K 1 – K 2 = 0
5 K 1−( 1+2 i ) K 2=0
Ejercicio 4
dx
=2 x +3 y
dt
dy
=2 x + y
dt
Primero hallamos el determinante
det ( A− AI ) = 2− λ 3
| 2 1−λ |
(2− λ) ( 1− λ )−6 = 0
λ 2−3 λ−4=0
( λ+ 1)( λ−4)=0
λ 1=−1; λ2=4
2−(−1) 3 3 3
| 2
=
1−(−1) 2 2
=> || |
3 k 1 +3 k 2 =0
13
2 k 1 +2 k 2=0
Por consiguiente:
3 k 1=−3 k 2
k 1=−k 2
k 1= 1
−1 ( )
Cuando λ 2=4 :
2−( 4) 3 −2 3
| 2 1−(4 )
=
2 −3||
=> |
−2 k 1 +3 k 2=0
2 k 1−3 k 2=0
2 k 1=3 k 2
3k2
k 1=
2
Cuando k 2=2, el vector propio relacionado es:
k 2= 3
() 2
x =( 1 ) . e −t
1
−1
x =( 3) . e 4t
2
2
Concluimos que la solución general del sistema es:
x=c 1 1 .e−t + c 2 3 . e 4 t
( )
−1 2 ()
E igualamos los elementos, con los elementos correspondientes del vector de la
izquierda, tenemos la solución general:
x (t)=c 1 e−t +3 c 2 e 4 t
y (t)=c 1 e−t + 2 c2 e 4 t
14
3. Conclusiones
15
BIBLIOGRAFIA
16