Presentación TEMA 6
Álgebra y Cálculo Matricial
Sistema de ecuaciones diferenciales y en
diferencias lineales
Alida Iralis Ortíz Pupo
Buscar un ejemplo de aplicación de sistemas de ecuaciones
diferenciales lineales no homogéneos, resolver y discretizar
con los métodos aprendidos. Comparar ambos sistemas.
2
Sistema masa-resorte-amortiguador con una fuerza externa
Para describir el sistema masa-resorte-amortiguador con una fuerza externa dependiente del tiempo, planteamos la
siguiente ecuación diferencial:
De manera abreviada, se puede escribir como:
Donde:
• : masa del objeto
• : constante del resorte
• : constante de amortiguamiento
• : fuerza externa aplicada al sistema que varía en el tiempo
Una posible definición de la fuerza externa es:
Donde
es la frecuencia angular en [
: amplitud de la fuerza externa
Despejando la aceleración :
3
Convertimos esta ecuación a un sistema de primer orden. Definimos:
Entonces:
Es decir que:
La solución de este sistema no homogéneo es:
Hallamos la descomposición espectral de :
Primero obtenemos los autovalores y autovectores de
4
En el caso del sistema donde
De aquí en adelante para simplificar los cálculos haremos la siguiente sustitución:
(√ ) √( )
2 2
Quedando los autovalores: 𝑐 𝑘 12 12
−4 = (𝑐 −4𝑘𝑚)= √𝑐 −4𝑘𝑚
𝑚 𝑚𝑚 𝑚
5
Para se tiene: Para se tiene:
La matriz es diagonalizable, por lo tanto:
6
Hallando :
Entonces:
7
Ya tenemos la descomposición de A:
Luego calculamos , y :
Entonces la solución del sistema es:
8
[[ ] ][ ][ [ ][ [ [] ] [ ][ ][] [ ][ ][ ]
𝑐− 𝑅 𝑐 − 𝑅 𝑐− 𝑅 𝑐 +𝑅 𝑐 − 𝑅 𝑐+ 𝑅 𝑐−𝑅
𝑐 −𝑅
− (𝑡 − 0 ) 𝑘 𝑐 − 𝑅 𝑘 𝑐 − 𝑅 𝑐 + 𝑅 𝑐 − 𝑅 𝑘 𝑐 − 𝑅 𝑐+ 𝑅 𝑐 − 𝑅 𝑘 𝑐 − 𝑅 𝑐 + 𝑅 𝑐 − 𝑅 𝑚 𝑚 𝑐+ 𝑅 𝑐 − 𝑅 𝑐 − 𝑅 ( )
𝑐− 𝑅
−𝑡
𝑐− 𝑅
𝜏
𝑐 −𝑅
− (𝑡 − 0 )
𝑐 +𝑅
− (𝑡 − 0)
𝑐 − 𝑅 𝑐 + 𝑅 𝜏 𝜏 − ( 𝑡 − 0)
− 𝑡 − 𝑡 2𝑚 2𝑚 2𝑚
− ( 𝑡 − 0)
2𝑚
− (𝑡− 0) − (𝑡− 0) 𝑐 − 𝑅
2𝑚 2𝑚 − 𝑡
𝑐+ 𝑅 𝜏
− 𝑡 2𝑚
𝑐+𝑅 𝑐−𝑅 𝑒 0 − − 𝑐+𝑅 𝑐−𝑅 𝑒 0 𝑒 0 − − 0 − 𝑒 − 𝑒 − − − 𝑒 − 𝑒 −𝑅𝑒 − 2𝑅 𝑒 0 2𝑅 𝑒 − 2𝑅 𝑒 𝑅𝑒 −𝑅𝑒 − 𝑒 − 𝑒 −2𝑚𝑅𝑓𝜏𝑒
2𝑚 2𝑚 𝑡 2𝑚 2𝑚 2𝑚 2𝑚 2𝑚 𝑡 2𝑚 2𝑚 𝑡
− − 𝑐+ 𝑅 𝑅 2 𝑅 − − 𝑐+ 𝑅 𝑐+ 𝑅 𝑅 2 𝑅 2 𝑘 2 𝑘 𝑅 2 𝑅 2 𝑘 2 𝑘 𝑘 𝑘
𝑢(𝑡)= 2𝑘 2𝑘 − (𝑡− ) 𝑢0+ 2𝑘 2𝑘 − 𝑡∫ 𝜏 𝑓(𝜏) 𝑑𝜏= 𝑢0+ ∫ 𝑐+𝑅 𝑐+𝑅 𝑓(𝜏) 𝑑𝜏= 𝑐−𝑅 𝑐+𝑅 𝑐−𝑅 𝑐+𝑅 𝑢0+ ∫ 𝑐+𝑅 𝑑𝜏 2 2
2𝑚 𝑘 𝑐 + 𝑅
0
2 𝑚 𝑡 2𝑚 𝑘 𝑐 + 𝑅 − 𝑐 − 𝑅 𝑡 − − 𝑐 + 𝑅 𝑡 − 𝑘 𝑐 + 𝑅 − 𝑐 − 𝑅 𝑡 − 𝑐 + 𝑅 𝑡 𝑡 𝑘 𝜏 𝑐 + 𝑅 𝜏 𝑘 − (𝑡 − ) 𝑘 − ( 𝑡 − ) 𝑐 − 𝑅 − ( 𝑡 − ) 𝑐 + 𝑅 − (𝑡 − ) − 𝑐 − 𝑅 𝑡 − 𝑐 + 𝑅 𝑡 𝑡 𝑐 + 𝑅 𝜏
0 𝑒 0 𝑒 0 𝑒 0 (2𝑚 ) 2𝑚( ) 2𝑚 2𝑚 2𝑚 2𝑚 2𝑚 2𝑚 2𝑚 2𝑚 2𝑚 2𝑚 2𝑚
1 1 1 1 𝑚 𝑒 𝑒 𝑒 𝑒 𝑒 𝑒 𝑚 − 𝑒 + 𝑒 − 𝑒 + 𝑒 𝑒 𝑒 𝑓 (𝜏 ) 𝑒 0 0 0
0 0 0 0
0
𝑅 2𝑅 𝑅 2 𝑅 𝑅 2 𝑅 𝑅 2 𝑅 𝑅 𝑅 2 𝑅 2 𝑅 2𝑚 𝑅 9
Resolviendo la integral considerando
Luego sustituyendo en la expresión de
10
Dado que solo interesa la primera componente de la solución:
Considerando las condiciones iniciales y sustituyendo :
11
Sustituyendo en la ecuación inicial:
12
Y dado que definimos anteriormente , se ha comprobado que la hallada es solución de la ecuación.
13
Ahora, hallamos la solución de por los métodos vistos en el grado, para comprobar que son equivalentes. Primeramente,
encontramos la solución de la ecuación homogénea:
Para hallar la solución hacemos las sustituciones:
La solución homogénea es de la forma:
La solución particular de la ecuación no homogénea entonces es:
, dado que
Luego,
14
Entonces la solución completa sería:
Para hallar las constantes sustituimos por los valores iniciales:
Luego, para hallar las constantes resolvemos el sistema:
15
+¿
Ahora sustituimos para hallar :
16
La solución entonces es:
La cual coincide con la solución hallada por el método visto en clases.
17
Considerando y y las condiciones iniciales:
[ ] [ ][ ]
𝑐−√𝑐−4𝑘𝑚 𝑐+√ −4𝑘𝑚 𝑐−√ −4𝑘𝑚 𝑐+√ −4𝑘𝑚 𝑐−√ −4𝑘𝑚 𝑐+√ −4𝑘𝑚
𝑚 1 𝑐−√ −4𝑘𝑚
2 𝑡 𝑐−√ −4𝑘𝑚𝜏
2 2 2 2 2 2 2
1 − (𝑡 − 0 ) − (𝑡− 0) − (𝑡− 0) − (𝑡− 0) − 𝑡 − 𝑡
𝑥(𝑡)= 2 (𝑐+√ −4𝑘𝑚)𝑒 −(𝑐 √ −4𝑘𝑚)𝑒 𝑢 + 𝑒 −𝑒 𝑢 − (𝑐+√ −4𝑘𝑚)𝑒 +(𝑐−√ −4𝑘𝑚)𝑒 − ∫𝑒 𝑓(𝜏)𝑑
2 2𝑚 2 2𝑚 2 𝑚 2 𝑚 2 2𝑚 2 2 𝑚 2𝑚
2√𝑐−4𝑘𝑚 √𝑐−4𝑘𝑚 2𝑚√𝑐−4𝑘𝑚
01 2 02
2 𝑘 2
𝑡0
18
La discretización del sistema viene dada por:
Donde:
y
Resolviendo la integral componente a componente:
19
De esta manera conformamos la matriz G:
20
Y el sistema discretizado sería:
Y sustituyendo :
21
GRACIAS
22
