Solución Examen Final Matemática IV (ING)

1. Considere la ecuación diferencial y '' t y ' y  1; con y  0   1 , y '  0   2 . Determine la solución y  t  mediante:

a) Un método analítico visto a lo largo del curso, y que sea el más conveniente para su solución

b) El método de la transformada de Laplace, considerando que L 1

s   0 , k 
k 
0

Nota: Debe constatar inequívocamente que la solución debe ser la misma en ambos casos

Solución

Dado que la ecuación diferencial es de segundo orden y de coeficientes variables, se hará uso del método de series
de potencias.

Suponiendo una solución de la forma: y  C
n 0
n t n  C0  C1 t  C2 t 2  C3 t 3  C4 t 4  C5 t 5  ...

Determinando la primera y segunda derivada de la serie de potencias


y '   n Cn t n 1  C1  2C2 t  3C3 t 2  4C4 t 3  5C5 t 4  ...
n 1


y ''   n  n  1 Cn t n  2  2C2  6C3 t  12C4 t 2  20C5 t 3  ...
n2

Sustituyendo en la ecuación diferencial:

  

 n  n  1 Cn t n2  t  n Cn t n1   Cn t n  1
n2 n 1 n 0

  

 n  n  1 Cn t n2   n Cn t n   Cn t n  1
n2 n 1 n 0

Unificando parámetros:

  

 n  n  1 Cn t n2 
n2
 n Cn t n 
n 1
C
n 0
n tn  1
" k  n  2" "k n" "k n"

  

  k  2  k  1 Ck 2 t k   k Ck t k   Ck t k  1
k 0 k 1 k 0

Para sumar series de potencias, deben empezar con el mismo valor k-ésimo; por lo tanto hacemos que todas las series
de potencias empiecen en k=1, ya que una de las series de potencias empieza en dicho valor y es la serie que tiene
menor cantidad de términos.

Ciclo 01/2017
   
2C2    k  2  k  1 Ck  2 t k   k Ck t k  C0   Ck t k  1
k 1 k 1 k 1

2C2  C0    k  2  k  1 Ck  2  k Ck  Ck  t k  1
k 1

 2C2  C0     k  2  k  1 Ck 2  k Ck  Ck  t k  1  0 t  0 t 2  0 t 3  ...
k 1

Por lo que obtenemos los siguientes resultados:

1  C0
C0  2C2  1  C2 
2
Utilizando las condiciones iniciales para determinar las constantes y  0   1 , y '  0   2

 y  C0  C1 t  C2 t 2  C3 t 3  C4 t 4  C5 t 5  ...

 y  0   C0  C1  0   C2  0   C3  0   C4  0   C5  0   ...  C0  1
2 3 4 5

 y '  C1  2C2 t  3C3 t 2  4C4 t 3  5C5 t 4  ...

 y '  0   C1  2C2  0   3C3  0   4C4  0   5C5  0   ...  C1  2

2 3 4

Dado que C0  1  C2  0

Determinando los demás coeficientes, utilizando la fórmula de recurrencia

Fórmula de recurrencia
 k  1 C
 k  2  k  1 Ck 2  k Ck  Ck  0  Ck  2  k  1
 k  2  k  1 k
Con dicha expresión obtenemos los demás coeficientes de la serie de potencias.

Evaluando:

Si k  1  C3 
1  1 C  C  0
1  2 1  1 1 3

 2  1 C   2  1
Si k  2  C4   0   C4  0
 2  2  2  1 2  2  2  2  1
 3  1 C   3  1 0  C  0
Si k  3  C5   
 3  2  3  1 3  3  2  3  1 5

 4  1 C   4  1
Si k  4  C6   0   C6  0
 4  2  4  1 4  4  2  4  1
 5  1 C   5  1 0  C  0
Si k  5  C7   
 5  2  5  1 5  5  2  5  1 7

Ciclo 01/2017
Se observa que los coeficientes de la serie de potencias serán cero a partir de C2 .

Por lo que la solución de la ecuación diferencial tendrá la forma y  C0  C1 t , donde C0  1 y C1  2 .

 La solución de la ecuación diferencial es y  t   1  2t

Usando la Transformada de Laplace

L  y ''  L t y '  L  y  L 1

d 1
 s 2Y  s   s y  0   y '  0    sY  s   y  0   Y  s  
ds s

Sustituyendo las condiciones iniciales y  0   1 , y '  0   2

d 1
 s 2Y  s   s  2   sY  s   1  Y  s  
ds s

Determinando Y  s 

1
s 2Y  s   s  2  Y  s   sY '  s   Y  s  
s
sY '  s    s 2  2  Y  s   s  2 
1
s
 2 2 1
Y '  s    s   Y  s   1   2 (Ecuación Diferencial Lineal)
 s s s

Determinando un Factor Integrante para resolver la E.D.L

 2 s2
  s  s ds
 s  e s e
2 2

d  2 s2   2 1 2 s
2 2

 s e Y     1   2  s e 2
s
ds    s s 
2 2
1 s  2 1 2 2
s
Y s  2 e 2   s s 2 
1   s e ds
s
2
s2
1  s2
Y s  2 e  s  2s  1 e ds
2 2
s
 2 s 
2 2 2
1  s2 s
Y s  2 e 
 s e 2
 2 s e 2
 1  ds
s  

La integral anterior no puede resolverse por los métodos tradicionales, por lo que se hará uso de series de potencias
para re expresar la integral.

Ciclo 01/2017
 s2 s2
3 5 1
s e  2s e  1  1  2s  s 2  s 3  s 4  s 5  ...
2 2 2
2 8 4
 2 
2 2
 
s s
3 2 3 5 4 1 5
  s e 2  2s e 2  1 ds   1  2s  2 s  s  8 s  4 s  ...  ds
 
 2 s 
2 2
s
1 3 1 4 1 5 1 6
    ds  C  s  s  s  s  s  s  ...
2 2 2
s e 2 s e 1
  2 4 8 24

Sustituyendo en la expresión para Y  s 

2
1 s  1 3 1 4 1 5 1 6 
Y s  2 e 2 C  s  s  2 s  4 s  8 s  24 s  ...
2

s2 s2
  1 1 1
Re expresando e 2 como una serie de potencias  e 2
 1  s 2  s 4  s 6  ...
2 8 48

1  1 2 1 4 1 6  1 3 1 4 1 5 1 6 
Y s  1  2 s  8 s  48 s  ... C  s  s  2 s  4 s  8 s  24 s  ...
2

s2

 1 1 1 1  1 1 1 1 
Y  s      2  s 2  s 4  ... C  s  s 2  s 3  s 4  s 5  s 6  ...
 2 s 8 48  2 4 8 24 
C C 1 1 2 C 2 1 4 C 4 1 6 C 6
Y  s  1    s  s  s  s  s  s  ...
2 s2 s 4 8 24 48 192 384

Aplicando la transformada inversa de Laplace

 C 1 C 1 2 C 2 1 4 C 4 
L 1
Y  s   L 1
1    2  s  s  s  s 
1 6 C 6
s  s  ...
 2 s s 4 8 24 48 192 384 
C  1  1  1  C  1  1 2 
L 1
Y  s   L 1  L 1
 L
1
 L  2  L  s 
2 s s  4 
C  1  C   1 6  C 6
 L 1  s 2   L 1  s 4   L 1  s 4   L 1  s L 1
 s   ...
8   24   48  192   384 

Considerando el hecho que L

1
s   0 , k 
k 
0

1  C 
L 1
Y  s   L 1
 L
1
 2   y t   1  C t
s s 

Usando la condición y '  0   2  y '  t   C  y '  0   C  C2

 La solución de la ecuación diferencial es y  t   1  2t

Ciclo 01/2017
 2
  1 
 4 
 
2. Haciendo uso de Funciones Gamma o Beta demuestre que  4  7  x  x  3 dx  
7 2
3 3 
Demostración

Sea x  3  4 y  x  4 y  3  dx  4 dy
Cambio de límites de integración:

Si x  3  y  0 Si x  7  y  1
1

0  4 y  4  7   4 y  3  4 4 dy
1 1
Planteando la nueva integral

1 1 1

 4 4 y 4  7  4 y  3 4 4 dy
1

0

1 1 1

 4 4 y 4  4  4 y  4 4 dy
1

0

1 1 1 1 3 1 1

 4 y 4 1  y  4 dy y 1  y  4 dy
1

1
 4 4 4 4  4 2 4
0 0

La integral presenta la forma de la definición de la función Beta   m, n    x 1  x 

1
m1 n 1
dx , con m  1  1 4 y
0

n 1  1 4 , por lo que los valores respectivos serán: m  5 4 y n  5 4 .

5 5 5 5

   3   
5 5
3 3
 4  ,   4
2  4 4
4 2  4 4
2

4 4 5 5 5
    
4 4 2
5 1  1 1
        1    
4 4  4 4
5  3  3  1   3  1   1  3
        1     1         
2  2  2  2   2  2   2  4

2 2
 1  1   1  1   1   1    1 
3 
         4  2  
 4  4   4  4   2      4 
 42  
3

3
 3 
4 4
2
  1 
 
2   4  
   7  x  x  3 dx 
7
4 lqqd
3 3 

Ciclo 01/2017
   sin k t  k t cos k t
1  1     
3. Haciendo uso del teorema de convolución demuestre que L  2 
  s  k  
2 2 2k 3

Demostración

Usando el teorema de la convolución f  t   g  t   L

1
F  s  G  s 
  

1 1  1  1 1 
L  2   L   2 
  s  k     s  k   s  k  
2 2 2 2 2

  
 1  1 1  k k 
L 1    L  2  2 
  s  k   k   s  k   s  k  
2 2 2 2 2 2

 k 
Recordando que L 1
 2 2
 sin  k t 
s  k 
 
  1
  2  sin  k t   sin  k t  
1 1
L  2
  s  k   k
2 2

 
  1
 sin  k u  sin  k  t  u   du
1 1 t
L  2  2
  s  k   k
2 2 0

Para resolver la integral usamos la identidad sin  x  sin  y  

1
2
 cos  x  y   cos  x  y 
Siendo x  k u y y  k  t  u 

 
 
L 1
 2
1

1 t

 2  cos  k  2u  t    cos  k t  du 
  s  k   2k
2 2 0

 
1  cos  k  2u  t  
t

1  1  
L  2  2   u cos  kt  
  s  k   2k 
2 2 2k  0

  sin k t  k t cos k t
 L 1  1 

    lqqd
 2 
  s  k  
2 3
2 2k

Ciclo 01/2017
 t
4. Resuelva la siguiente ecuación integro-diferencial y '' 2 y ' y  2 y  u  du  5 , y  0   y '  0   0
0

Solución
Aplicando la transformada de Laplace

 t 
L  y ''  L 2 y '  L  y  L 2 y  u  du   L 5
 0 
2 5
 s 2Y  s   s y  0   y '  0   2  sY  s   y  0   Y  s   Y  s  
s s
Sustituyendo las condiciones iniciales y  0   0 , y '  0   0

2 5
s 2Y  s   2sY  s   Y  s   Y  s  
s s
 2 5
Y  s   s 2  2s  1   
 s s

Determinando Y  s 

5 5
Y s   Y s 
 2 s  2s 2  s  2
3
s  s 2  2s  1  
 s
s 2 1
Re expresando Y  s     2 
s 1 s 1 s  2
2

Aplicando la transformada inversa de Laplace

 s   2   1 
L 1
Y  s   L 1
 2 L
1
 2 L
1
 
 s  1  s  1 s  2

 y  t    cos  t   2sin  t   e 2t

3 y ''  t   3z ''  t   t e t  3cos  t 

5. Resolver el siguiente sistema de ecuaciones  , dado que:
t y ''  t   z '  t   sin  t 
y  0   1, y '  0   2, z  0   4, z ''  0   0
Nota: Debe especificar las expresiones para y  t  , z  t 

Solución
Aplicando la transformada de Laplace a la ecuación 1:

3L  y ''  t   3L  z ''  t   L t e   3L cos  t 

t

Ciclo 01/2017
 1 3s
3  s 2Y  s   s y  0   y '  0    3  s 2 Z  s   s z  0   z '  0   
 s  1 s 1
2 2

Sustituyendo las condiciones iniciales y  0   1 , y '  0   2 , z  0   4 y z '  0  A

1 3s
3  s 2Y  s   s  2  3  s 2 Z  s   4s  A  
 s  1 s 1
2 2

Desarrollando y despejando Z  s 

1 1 5 3  2  A
Z s      Y s
3s  s  1 s  s  1 s
2
2 2
s2

Aplicando la transformada de Laplace en la ecuación 2:

L t y ''  t   L z '  t   L sin  t 

d 2 1
 s Y  s   s y  0   y '  0    sZ  s   z  0   2
ds s 1
Sustituyendo las condiciones iniciales y  0   1 , y '  0   2 y z  0   4

d 2 1
 s Y  s   s  2   sZ  s   4  2
ds s 1
1
s 2Y '  s   2sY  s   1  sZ  s   4 
s 1
2

Desarrollando y despejando Y '  s 

2 1 3 Z s
Y ' s  Y  s   2 2  2
s s  s  1 s s

Sustituyendo la expresión de Z  s  se obtiene:

2 1 3 1 1 1 5 3  2  A 
Y ' s  Y  s   2 2  2  2     Y  s  
s  s  1 s s  3s  s  1 s  s  1 s 
2 2
s s2 
3 2 1 3  2  A
Y ' s  Y  s   2  3  (Ecuación Diferencial Lineal)
s s 3s  s  1 2
s3
3
 ds
Determinando un Factor Integrante para resolver la E.D.L   s   e s  s3

 2 3  2  A 
d 3
ds
 s Y  s    s3   2  3
1
 s 3s  s  12

s 3


 

1  3 2 1 3  2  A   1  1 
Y s  3 
 s       ds  Y  s   3   2 s   3  2  A   ds
s   s 3s3  s  1
2 2
s 3

 s 
 3  s  1
2




Ciclo 01/2017
 1 2 1  1 1 3  2  A C
Y s  3 
 s   3  2  A  s  C   Y s    3   3
s  3  s  1  s 3s  s  1 s2 s

Sustituyendo Y  s  en la expresión de Z  s  se obtiene:

1 1 5 3  2  A  1 1 3  2  A C 
Z s           3
3s 2  s  1 s  s 2  1 s  s 3s  s  1
2
s2 3
s2 s 

1 1 4 1 C
Z s     3  3
3s  s  1 s  s  1 s 3s  s  1 s
2 2 2

Aplicando la transformada inversa de Laplace para determinar Y  s  y Z  s 


 1 3  2  A C 

L 1
Y  s   L 1
  3
1
  3
 s 3s  s  1
2
 s s 


 1 1  1   1   1 
L 1
Y  s   L 1
 L
1
 L
1
 2 L
1
 3 L
1
 
 s  3s   3s   3s   3  s  1 
1  3  2  A  1 C 
L  L  3
s 
2
 s 
1 1 1 1 C 1 C  1 2 1
y  t   1   t  t 2  et  3  2  A t  t 2  y  t      t 2   3  2  A    t   e t
3 3 6 3 2 6 2   3 3 3


 C

L 1
Z  s   L 1
 2
1

1

4

1
 3
 3s  s  1 s  s  1

2 2
s 3s  s  1 s 
3


 s  
 1   
 
 8 C   1 
L 1
Z  s   L 1
 2 L
1
 L
1

1
 L 1
 L
1
 3 L
1
 3
 s  1  3  s  1   3  s  1   3s  s   3s 
2
   

1 1 8 C 1 1 C 8 1 1
z  t   cos  t   et  t et   t 2  t 2  z  t      t 2   t e t  e t  cos  t 
3 3 3 2 6 6 2  3 3 3

Para encontrar las constantes A y C se utilizan las condiciones iniciales y '  0   2 y z ''  0   0 .

4
Con lo que se obtiene A  y C 1
3

 2 2 5 1 t 2
 y  t   3 t  3 t  3 e  3
 
 z  t   2 t 2  1 t e t  1 e  t  cos  t   8
 3 3 3 3

Ciclo 01/2017
6. Como recordará, en Matemática III, usted se familiarizó con las series de potencias, las cuales permiten representar
una función a través de una serie polinomial infinita. En el actual curso de Matemática IV, usted se ha familiarizado
con la transformada de Laplace, la cual es una herramienta muy poderosa para resolver ecuaciones diferenciales
sumamente complejas de una forma algebraica transformando las funciones involucradas en el dominio del tiempo
al dominio de la frecuencia compleja. Haciendo uso de manipulaciones de la serie de Maclaurin de la función coseno
demuestre que:
 cos t 
   
  41s
L   e
 t  s

Demostración
 1 t 2k
k

Sabemos que la serie de Maclaurin para f  t   cos  t  es 
k  0  2k  !

Si sustituimos t , por t en la serie de Maclaurin del cos  t  , obtendremos la representación en series de potencias
de f  t   cos  t
 1 t 2k
k
 1 t k
k
 t   1 t k
k

 t
  cos 1 
 cos  t   
k 0  2k  !
 cos
k  0  2k  !

t t

k  0  2k  !

cos  t  
 1 t k 1 2
k


t

k  0  2k  !

Aplicando la transformada de Laplace

 cos t 
      1k k 1 2 
L    L  t 
 t   k  0  2k  ! 
 

L 
 cos t 
  
 L
 t 1 2 t1 2 t 3 2 t 5 2 t 7 2 t 9 2
    

 ...
 
 t   0! 2! 4! 6! 8! 10! 

 cos t 
    t 1 2   t1 2  t3 2  t5 2  t7 2  t9 2 
L  L    L    L    L    L    L    ...
 t   0!   2!   4!   6!   8!  10! 


t   s 
 n 1
Recordando que L n
n 1

 cos t 
  
 1  1 2  1   3 2  1   5 2  1   7 2  1   9 2  1  11 2 
L   12
      ...
 t  0! s 2! s 3 2 4! s 5 2 6! s 7 2 8! s 9 2 10! s11 2

Re expresando los valores de las funciones gammas a valores conocidos

  1 2   
1 1
   3 2    1 2  1   1 2   
2 2

Ciclo 01/2017
 3 31  3
   5 2     3 2  1    3 2      
2 22  4
5 53  15
   7 2     5 2  1    5 2      
2 24  8
7 7  15  105
   9 2     7 2  1    7 2      
2 2 8  16
9 9  105  945
  11 2     9 2  1    9 2      
2 2  16  32

 
1 3 15 105 945
 cos t      
  1  1 2 1 4 1 8 1 16 1 32
L    12
 32
 52
 72
 92
 11 2
 ...
 t  0! s 2! s 4! s 6 ! s 8 ! s 10! s


Re expresando y sacando factor común 12
s
  
 cos t 
   1 1 1 1 1 1 1 1 1 
L    1 2 1   2
 3
 4
 5
 ...
 t  s  4s 2!16 s 3! 64 s 4! 256 s 5!1024 s 


  
 cos t 
   1 1 1 1 1 1 1 1 1 
L   1  2 2
 3 3
 4 4
 5 5
 ...
 t  s  4s 2! 4 s 3! 4 s 4! 4 s 5! 4 s 


  
 cos t 
  1  1  1 1  1  1  1 1  1  1  1 1 
0 1 2 3 4 5

L                      ...
 t  s  0!  4s  1!  4s  2!  4s  3!  4 s  4!  4 s  5!  4s  

1 2 1 3 1 4 1 5 1
Recordando la serie de Maclaurin de la función e x  1  x  x  x  x  x  ... donde x  s
2! 3! 4! 5! 4

 cos t 

 L 


 
  41s
lqqd
 e
 t  s

t
Y u 
7. Dada la ecuación integral  t  u  13
du  t 1  t 
0

Demuestre que al resolver la ecuación integral por medio de la transformada de Laplace se obtiene como solución

3 3 13
y t   t  2  3t 
4
Nota: Este ejercicio no debe demostrarlo sustituyendo y  t  en la ecuación integral, sino que utilizando la transformada de Laplace

Demostración
t

Re expresando la ecuación integral Y  u  t  u 

1 3

0
du  t  t 2

Ciclo 01/2017
Se puede apreciar que se genera la integral de convolución, por lo tanto se aplica el teorema de convolución quedando
de la siguiente manera:
y  t   t 1 3  t  t 2

Aplicando la transformada de Laplace L  y  t   t   L t  t 

1 3 2

  1 3  1 1 2 s2 3  1 2  1  1 2 
 Y s   2 3  Y s      Y s   4 3  7 3 
s 1 31 s s   2 3  s 2 s 3    2 3  s s 

Aplicando la transformada inversa de Laplace


 1  1 2  
L 1
Y  s   L 1
   4 3  7 3 
   2 3  s
 s  

1    1 3  1     4 3  1  
L 1
Y  s    
1
L 1
 1 31  
2
L 1
  
 2 3   1 3  1  s    4 3  1  s
4 31


1  1 2 
 y t    t1 3  t 4 3 
  2 3    4 3    7 3 

Re expresando los valores de las funciones gammas a valores conocidos

1
   4 3   1 3  1   1 3
3
4 4
   7 3    4 3  1    4 3   1 3
3 9

1  1 2 
 y t    t1 3  t 4 3 
  2 3    4 3   7 3 

3 t1 3 18 t4 3
 y t   
  2 3  1 3 4   2 3  1 3


Pero   2 3  1 3   1 3  1  1 3 y recordando la propiedad   n   1  n  
sin  n  

 2
   2 3  1 3   1 3  1  1 3  
sin  3 3

3 3 13 9 3 43 3 3 13 3  3 3 13 1
 y t   t  t  y t   t 1  t   y  t   t  2  3t 
2 4 2  2  2 2

3 3 13
 y t   t  2  3t  lqqd
4

Ciclo 01/2017
8. Una viga uniforme de longitud L cuya sección transversal tiene un momento de inercia I , y el material de fabricación
posee un módulo de elasticidad E , se encuentra empotrada en sus dos extremos, en x  0 y x  L . Si se aplica una
carga distribuida por unidad de longitud w  x  , la cual viene dada por:

  2 
w 1  x , 0  x  L 2
w  x    0  L 
0 , L 2 xL


Si la deflexión estática " y " (medida hacia abajo) en términos de " x " obedece a la ecuación diferencial:

d4y
EI  w x , 0  x  L
dx 4
Determine una expresión para y  x  haciendo uso de la transformada de Laplace. Las condiciones frontera de la viga
doblemente empotrada son: y  0   y  L   0 , y '  0   y '  L   0

Solución

Re expresando la carga distribuida por unidad de longiutd en términos de la función escalón unitario

 2   2   L
w  x   w0 1  x   w0 1  x  u  x  
 L   L   2

2w0  L  L  L 
w x    x  1   u  x   
L 2  2  2 

Por lo que la ecuación diferencial que describe la deflexión estática de la viga será:

d 4 y 2w0  L  L  L 
EI    x  1   u  x   
 2  2 
4
dx L 2
Aplicando la transformada de Laplace

 d4y  2w  L  L  L  
L  EI 4 
 L  0   x  1   u  x    
 dx   L 2  2  2  

2w0  L 2 1 1  Ls 2 
EI  s 4Y  s   s 3 y  0   s 2 y '  0   sy ''  0   y '''  0     2 2e
L  s s s 
Ciclo 01/2017
Sustituyendo las condiciones iniciales y  0   0 , y '  0   0 , y ''  0   C1 y y '''  0   C2

2w0  L 2 1 1  Ls 2 
EI  s 4Y  s   C1 s  C2    s  s 2  s 2 e 
L

Desarrollando y despejando Y  s 

C1 C2 2w0  L 2 1 1  Ls 2 
Y s      e
s 3 s 4 EIL  s 5 s 6 s 6 

Aplicando la transformada inversa de Laplace

 C1 C2 2w0  L 2 1 1  Ls 2  
L 1
Y  s   L 1
 3 4    e
EIL  s 5 s 6 s 6  
s s 

 2!  C2  3!  2w0  L 2  4!  1  5!  1  1  Ls 2  
L 1
Y  s   C2! L
1 1
 3  L
1
 4   L 1
 5  L
1
 6  L
1
 6e 
 s  3!  s  EIL  4!  s  5!  s  5! s 

w0  5L 4 5  L 
5
C1 2 C2 3 L 
 y  x  x  x   x  x   x   u  x  
2 6 60 EIL  2  2  2  

Para encontrar las constantes C1 y C2 se utilizan las condiciones y  L   0 y y '  L   0

23w0 L2 9w L
Con lo que se obtiene C1  y C2   0
960 EI 40 EI

Por lo que la deflexión estática de la viga está dada por:

w0  5L 4 5  L 
5
23w0 L2 2 3w0 L 3 L 
 y  x  x  x   x  x   x   u x  
1920EI 80EI 60 EIL  2  2  2  

Ciclo 01/2017