Daniel Rodríguez CI27462451
Sección 3 Matemática 4
Ejercicio 1 (Valor 10 puntos):
Para resolver ecuaciones diferenciales de orden superior seguimos utilizando el
comando dsolve('eq', 'cond', 'var'). Cabe resaltar que cuando no se le
indica a matlab, cuál es la variable independiente, toma a t por defecto.
Por ejemplo, si queremos resolver la ecuación diferencial: y ' ' −2 y '+ y=e ^x
sujeta a las condiciones iniciales y ' (0)=2 ; y (0)=0 y también deseamos
conocer la gráfica de la soluciónobtenida. Creamos un script en Matlab, con lo
siguiente:
Código en Matlab
% Solucionamos la ecuación diferencial:
y=dsolve('D2y-2*Dy+y=exp(x)','Dy(0)=2','y(0)=0','x')
% Simplificamos la expresión obtenida
simplify (y)
% Visualizamos en forma ordenada la solución:
pretty (y)
% Utilizamos el comando ezplot, para graficar en el intervalo [-2,2]
ezplot(y,[-2,2])
En la Figura 1 se muestra la solución de la ED y en la Figura 2 se muestra la
gráfica de la dicha solución.
a) Transcriba en un Script del Matlab y posteriormente ejecute el código,
insertando las instrucciones que se mostraron anteriormente. Una vez
compruebe el resultado, reutilice el código, ahora para que encuentre la
solución al problema de valor inicial:
4 y ' '+4 y '+17 y=x sujeta a y ' (0)=1; y (0)=0
b) Grafique la solución obtenida
Respuestas:
Código:
%% Daniel Rodriguez CI27462451
%EJERCICIO 1
%Solucion:
y=dsolve('4*D2y+4*Dy+17*y=x','Dy(0)=1','y(0)=0','x')
y=simplify(y)
ezplot(y,[-3,3])
>> ejercicio1
y = x/17 + (4*cos(2*x)*exp(-x/2))/289 + (137*sin(2*x)*exp(-x/2))/289 - 4/289
y = x/17 + (4*cos(2*x)*exp(-x/2))/289 + (137*sin(2*x)*exp(-x/2))/289 - 4/289
>>
Ejercicio 2 (Valor 10 puntos):
La ecuación del movimiento forzado de un resorte debido a una fuerza exterior
que actúa sobre la masa oscilante sujeta a dicho resorte, aplicando la segunda
ley de Newton, viene dada por la ecuación diferencial siguiente en el que se
relaciona la posición, x, de la masa con el tiempo, t: donde m es la masa, b es
la constante la constante del resorte y f la fuerza exterior.
c) Si m=1 ;b=6 ;k=10 y f (t )=25 cos( 4t ) determine la gráfica que modela la
posición del resorte considerando las condiciones iniciales y ' (0)=0; y
(0)= ½ y un intervalo de tiempo de [0,10]
Respuestas:
Código:
%% Daniel Rodriguez CI27462451
%EJERCICIO 2
%Solucion:
y=dsolve('D2y=-6*Dy-10*y+25*cos(4*t)','Dy(0)=0','y(0)=1/2','t')
y=simplify(y)
ezplot(y,[0,10])
>> ejercicio2
y = sin(t)*((25*cos(3*t))/12 + (75*cos(5*t))/68 + (25*sin(3*t))/12 +
(125*sin(5*t))/68) - cos(t)*((25*cos(3*t))/12 - (125*cos(5*t))/68 -
(25*sin(3*t))/12 + (75*sin(5*t))/68) + (38*exp(-3*t)*cos(t))/51 - (86*exp(-
3*t)*sin(t))/51
y = (50*sin(4*t))/51 - (25*cos(4*t))/102 + (38*exp(-3*t)*cos(t))/51 - (86*exp(-
3*t)*sin(t))/51
>>
Ejercicio 3 (Valor 10 puntos):
Para resolver un sistema de ecuaciones diferenciales, utilizando matlab, se
sigue el código que se muestra a continuación. Observe que las ecuaciones
diferenciales se han adaptado al ejemplo siguiente:
Código en Matlab
% Solucionamos el sistema de ecuaciones diferenciales:
sistema=dsolve('Dx=y','Dy=-x','x(0)=1','y(0)=1')
% Accedemos a la primera fila del vector que contiene la solución para y
y=sistema.y
% Accedemos a la primera fila del vector que contiene la solución para x
x=sistema.x
El resultado del código anterior se muestra en la figura 3.
a) Utilice el código anterior para encontrar la solución al sistema:
Resultados:
Código:
%% Daniel Rodriguez CI27462451
%EJERCICIO 3
%Solucion:
sistema=dsolve('Dx=(-2*x/25)+(1*y/50)','Dy=(2*x/25)-(2*y/
25)','x(0)=25','y(0)=0');
y=sistema.y
x=sistema.x
>> ejercicio3
y=
exp(-(3*t)/25)*(25*exp((2*t)/25) - 25)
x=
(exp(-(3*t)/25)*(25*exp((2*t)/25) + 25))/2
>>
