Laboratorio de Automatización y Control de Sistemas Dinámicos
Para esta semana consistió en elegir un problema y resolverlo:
Resolver la siguiente ecuación diferencial:
Con las siguientes características:
Primeramente explicar algunas expresiones que vienen en el problema:
ẍ = Segunda Derivada
Wn = Frecuencia Natural (Para mi problema sera tratado como constante)


Para resolverlo seria mediante las Transformadas de Laplace.

Teniendo que:
L{X''} = s²x(s) - sx(0) - x'(0)
L{X'} = sx(s) - s¹⁻¹x(0) = sx(s) - 1*x(0) = sx(s) - x(0)
L{X} = s¹⁻¹ x(s) = s⁰x(s) = 1*x(s) = x(s)
Nos da como resultado:
![s²x(s) - sx(0) - x'(0) + 2 \zeta Wn[sx(s) - x(0)] + Wn²[x(s)]](http://rinconmatematico.com/latexrender/pictures/fccaae18f4040fbf421da08d37e3996b.gif)
![s²x(s) - sx(0) - x'(0) + 2 \zeta Wn[sx(s)] - 2 \zeta Wnx(0) + Wn²[x(s)]](http://rinconmatematico.com/latexrender/pictures/09203b180e7983302f8aa43cd1299dac.gif)
Sustituyendo los valores de x' y de x'' en la ecuación anterior:
![s²x(s) - sa - b + 2 \zeta Wnsx(s) - 2 \zeta Wna + Wn²[x(s)]](http://rinconmatematico.com/latexrender/pictures/3c90369c41915ba12786260809d4d341.gif)
![s²[x(s)] - as - b + 2 \zeta Wn s[x(s)] - 2 \zeta Wna + Wn²[x(s)]](http://rinconmatematico.com/latexrender/pictures/3ff72ee031a372d392c35382361e5988.gif)
![x(s)[s²x + 2 \zeta Ws + Wn²] = as + b + 2 \zeta Wna](http://rinconmatematico.com/latexrender/pictures/5ee1aa7a487039298d043b6b50784a1b.gif)

Hay que factorizar el denominador la forma sencilla para mi es mediante la formula para resolver ecuaciones de 2do grado.

De ahí tenemos que:




Sustituyendo en la ecuación tenemos:
![s = \frac{-2 \delta Wn \pm \sqrt[2]{(2 \zeta Wn)² - 4(1) (Wn²)} }{2(1)}](http://rinconmatematico.com/latexrender/pictures/eca0a01cce187efcd3959338c5866487.gif)
Resolviéndolo:
![s = \frac{-2 \zeta Wn \pm \sqrt[2]{(2 \zeta Wn)²-4(1)(Wn²)}}{2(1)}](http://rinconmatematico.com/latexrender/pictures/69eefecb6bcde794b3d99e6cd2d4e0c8.gif)
![s = \frac{-2 \zeta Wn \pm \sqrt[2]{4 \zeta² Wn²-4Wn²}}{2}](http://rinconmatematico.com/latexrender/pictures/c46e2c30499bf4c1503c712993e1819f.gif)
![s = \frac{-2 \zeta Wn \pm \sqrt[2]{4 (\zeta² Wn²-Wn²)}}{2}](http://rinconmatematico.com/latexrender/pictures/f3da3f4bdd9801aac07b3c7cb9cccc78.gif)
![s = \frac{-2 \zeta Wn \pm \sqrt[2]{4}\sqrt[2]{\zeta² Wn²-Wn²}}{2}](http://rinconmatematico.com/latexrender/pictures/cf4265049468ed25f8a360cf3af7d84d.gif)
![s = \frac{-2 \zeta Wn \pm (2)\sqrt[2]{\zeta² Wn²-Wn²}}{2}](http://rinconmatematico.com/latexrender/pictures/8a369d15ff751539ecde75963de4a378.gif)
![s = -\zeta Wn \pm \sqrt[2]{\zeta² Wn²-Wn²}](http://rinconmatematico.com/latexrender/pictures/fd37f7e6b3ab3f3e32e923acf0ca8feb.gif)
Como ya no se puede simplificar más la factorización es igual a:
![(s + \zeta Wn - \sqrt[2]{\zeta² Wn²-Wn²})(s + \zeta Wn + \sqrt[2]{\zeta² Wn²-Wn²})](http://rinconmatematico.com/latexrender/pictures/efd13f81825c7aea06a62032b1f35158.gif)
Sustituyendo esto en la ecuación original:
![x(s) = \frac{as + b + 2 \zeta Wna}{(s + \zeta Wn - \sqrt[2]{\zeta² Wn²-Wn²})(s + \zeta Wn + \sqrt[2]{\zeta² Wn²-Wn²})}](http://rinconmatematico.com/latexrender/pictures/a6ce082d5b76c2a520785b284345aac9.gif)
Para que las operaciones se nos resulten mas fácil y los estos son constantes los sustituimos por variables:
![x = \zeta Wn - \sqrt[2]{\zeta² Wn²-Wn²})](http://rinconmatematico.com/latexrender/pictures/a73d25fb39e7d12c567af27f64c664e7.gif)
![y = \zeta Wn + \sqrt[2]{\zeta² Wn²-Wn²})](http://rinconmatematico.com/latexrender/pictures/460f8f17f46cefa9bac4dc91a70bc6ba.gif)

Por lo tanto nos queda:

Aplicando fracciones parciales:

Tenemos que:


Sustituyendo los valores de A y B en la ecuación original tenemos:

A continuación aplicamos la transformada inversa de laplace a la ecuación, ya que las variables x, y y z son constantes solo sera hará con los términos de s:

Resultado:

Sustituyendo los valores de x, y y z en la ecuación:
![x(t) = (\frac{ b + 2 \zeta Wna - a( \zeta Wn - \sqrt[2]{\zeta² Wn²-Wn²})}{(\zeta Wn + \sqrt[2]{\zeta² Wn²-Wn²})-( \zeta Wn - \sqrt[2]{\zeta² Wn²-Wn²})})e^-( \zeta Wn - \sqrt[2]{\zeta² Wn²-Wn²})t+(\frac{ b + 2 \zeta Wna - a(\zeta Wn + \sqrt[2]{\zeta² Wn²-Wn²})}{( \zeta Wn - \sqrt[2]{\zeta² Wn²-Wn²})-y})e^-(\zeta Wn + \sqrt[2]{\zeta² Wn²-Wn²})t](http://rinconmatematico.com/latexrender/pictures/974ecbe55b595bd3482d3110c7f13e54.gif)
Referencias:
http://es.wikipedia.org/wiki/Ecuaci%C3%B3n_de_segundo_grado
http://matematicas.udea.edu.co/~jescobar/docs/cap6.pdf
OK; 15 pts.
ResponderEliminar