martes, 26 de junio de 2012

4.2.1 METODO DE LOS OPERADORES


Un operador es un objeto matem atico que convierte una funci on en otra, por ejemplo, el operador

derivada convierte una funci on en una funci on diferente llamada la funci on derivada. Podemos

de nir el operador derivada

D que al actuar sobre una funci on diferenciable produce la derivada


de esta, esto es:


D

0f(x) = f(x) ; D1f(x) = f0(x) ; D2f(x) = f00(x) ; : : : ;Dnf(x) = f(n)(x) :


Es posible construir la siguiente combinaci on lineal con los operadores diferenciales:


P

(D) = a0 + a1D + a2D2 + + anDn ; an 6= 0 : (1)


donde

a2; a1; a2; : : : an son constantes. A este nuevo objeto lo podemos llamar el Operador Polinomial
de orden

n.


Por otro lado, recordemos que una ecuacion diferencial lineal de orden
n con coe cientes constantes
es una ecuacion de la forma A



ny(n) + an􀀀1y(n􀀀1) + + a2y00 + a1y0 + a0y = Q(x) ; (3)


por lo tanto, (

3) se puede escribir de una manera compacta como


P

(D)y = Q(x) : (4)


El operador polinomial es lineal, esto signi ca que tiene las siguientes propiedades

Si

f1(x) y f2(x) son dos funciones diferenciables de orden n, entones


P

(D) [
f
1(x) + f2(x)] =
P
(D)f1(x) + P(D)f2(x)


donde


y son constantes. Adem as:


Si

y1(x); y2(x); : : : ; yn(x) son n soluciones de la ecuaci on diferencial homog enea P(D)y = 0


entonces

yh(x) = C1y1(x) + C2y2(x) + + Cnyn(x) es tambi en una soluci on.


Si

yh(x) es una soluci on de P(D)y = 0 y yp(x) es una soluci on de P(D)y = Q(x) entonces


y

(x) = yh(x) + yp(x) es una soluci on de P(D)y = Q(x).

4.2.2 UTILIZACION DE LA TRANSFORMADA DE LAPLACE

La transformada de Laplace recibe su nombre en honor del matemático francés Pierre-Simon Laplace, que la presentó dentro de su teoría de la probabilidad. En 1744, Leonhard Euler había investigado un conjunto de integrales de la forma:

 z = \int X(x) e^{ax}\, dx
 z = \int X(x) x^A \, dx

— como soluciones de ecuaciones diferenciales, pero no profundizó en ellas y pronto abandonó su investigación. Joseph Louis Lagrange, admirador de Euler, también investigó ese tipo de integrales, y las ligó a la teoría de la probabilidad en un trabajo sobre funciones de densidad de probabilidad de la forma:

 \int X(x) e^{- a x } a^x\, dx,

— que algunos historiadores interpretan como auténticas transformadas de Laplace.

Este tipo de integrales atrajeron la atención de Laplace cuando, en 1782, y siguiendo la idea de Euler, trató de emplear estas integrales como soluciones de ecuaciones diferenciales. Parece ser que en 1785 dio un paso más allá, y reenfocó el problema para en vez de usar las integrales como soluciones, aplicarlas a las ecuaciones dando lugar a las transformadas de Laplace tal y como hoy en día se entienden. Usó una integral de la forma:

 \int x^s \phi (s)\, dx,

— análoga a la transformada de Mellin, con la que transformó una ecuación diferencial en una ecuación algebraica de la que buscó su solución. Planteó alguna de las principales propiedades de su transformada, y de alguna forma reconoció que el método de Joseph Fourier para resolver por medio de series de Fourier la ecuación de difusión podría relacionarse con su transformada integral para un espacio finito con soluciones periódicas.

Pese al logro, las transformadas de Laplace pronto cayeron en un relativo olvido, al haber sido presentadas en el campo de la probabilidad –ajeno a su moderna aplicación en la física y la ingeniería–, y ser tratadas sobre todo como objetos matemáticos meramente teóricos.

La moderna aplicación de las transformadas de Laplace y toda su teoría subyaciente surge en realidad en la segunda mitad del siglo XIX. Al tratar de resolver ecuaciones diferenciales relacionadas con la teoría de vibraciones, el ingeniero inglés Oliver Heaviside (1850-1925) descubrió que los operadores diferenciales podían tratarse analíticamente como variables algebraicas. De acuerdo con el "cálculo operacional", si se tiene una ecuación diferencial de la forma:

(D-a)y=f(t)\,

— donde D es el operador diferencial, esto es, D=d/dx, entonces la solución general a dicha ecuación es de la forma:

y=e^{at} \int e^{-at} f(t) dt +c_1 e^{at}.

Heaviside observó que si se trataba al operador D como una variable algebraica, era posible alcanzar igualmente la solución de toda ecuación pareja a la de arriba. En efecto, según la solución general, se cumple que:

y=\frac 1 {D-a} f(t)= e^{at} \int e^{-at} f(t) dt +c_1 e^{at}

Entonces, si se considera una ecuación diferencial de segundo orden como la siguiente:

y''-3y'+2y=e^t\,

— ésta puede reescribirse en para resaltar el operador D como:

(D^2-3D+2)y=e^t\,

Heaviside propuso despejar y y tratar a D algebraicamente, en cuyo caso se tendría que:

y= \frac {e^t} {D^2-3D+2} = \frac {e^t} {(D-1)(D-2)}=\frac 1 {D-2} e^t - \frac 1 {D-1} e^t

Sustituyendo las fracciones en D por la expresión integral de las mismas arriba presentada, se llega a la solución de la ecuación diferencial:

y=(e^{2t} \int e^{-2t} f(t) dt +c_1 e^{2t})-(e^{t} \int e^{-t} f(t) dt +c_2 e^{t})=e^{2t}(-e^{-t})+c_1 e^{2t})-(e^t (t) +c_2 e^t)
\Big (y=c_1 e^{2t}- (c_2+1) e^t -t e^t \Big)

Heaviside publicó sus resultados, cuya utilidad a la hora de resolver ecuaciones de la física y la ingeniería hizo que pronto se extendieran. Sin embargo, el trabajo de Heaviside, formal y poco riguroso, atrajo las críticas de algunos matemáticos puristas que los rechazaron argumentando que los resultados de Heaviside no podían surgir de tal forma. No obstante, el éxito del método hizo que pronto fuera adoptado por ingenieros y físicos de todo el mundo, de manera que al final atrajo la atención de cierto número de matemáticos tratando de justificar el método de manera rigurosa. Tras varias décadas de intentos, se descubrió que la Transformada descubierta por Laplace hacía un siglo no sólo ofrecía un fundamento teórico al método de cálculo operacional de Heaviside, sino que además ofrecía una alternativa mucho más sistemática a tales métodos.

Hacia principios del siglo XX, la transformada de Laplace se convirtió en una herramienta común de la teoría de vibraciones y de la teoría de circuitos, dos de los campos donde ha sido aplicada con más éxito. En general, la transformada es adecuada para resolver sistemas de ecuaciones diferenciales lineales con condiciones iniciales en el origen. Una de sus ventajas más significativas radica en que la integración y derivación se convierten en multiplicación y división. Esto transforma las ecuaciones diferenciales e integrales en ecuaciones polinómicas, mucho más fáciles de resolver.También se aplica en EDP y ecuaciones diferenciales en diferencias.

[editar]Propiedades


[editar]Linealidad


\mathcal{L}\left\{a f(t) + b g(t) \right\}
  = a \mathcal{L}\left\{ f(t) \right\} +
    b \mathcal{L}\left\{ g(t) \right\}

[editar]Derivación


\mathcal{L}\{f'(t)\}
  = s \mathcal{L}\{f(t)\} - f(0)
\mathcal{L}\{f''(t)\}
  = s^2 \mathcal{L}\{f(t)\} - s f(0) - f'(0)
\mathcal{L}\left\{ f^{(n)}(t) \right\}
  = s^n \mathcal{L}\{f(t)\} - s^{n - 1} f(0) - \cdots - f^{(n - 1)}(0) =
=\mathcal{L}\left\{ f^{(n)}(t) \right\}
  = s^n \mathcal{L}\{f(t)\} - \sum_{i=1}^{n} s^{n - i} f^{(i - 1)}(0)

[editar]Integración


\mathcal{L}\left\{ \int_0^t f(\tau)\, d\tau \right\}
  = {1 \over s} \mathcal{L}\{f\}

[editar]Dualidad


\mathcal{L}\{ t f(t)\}
  = -F'(s)

[editar]Desplazamiento de la frecuencia


\mathcal{L}\left\{ e^{at} f(t) \right\} =
  F(s-a)

Desplazamiento temporal


\mathcal{L}\left\{ f(t - a) u(t - a) \right\}
  = e^{-as} F(s)
\mathcal{L}^{-1} \left\{ e^{-as} F(s) \right\}
  = f(t - a) u(t - a)


[editar]Desplazamiento potencia n-ésima


\mathcal{L}\{\,t^nf(t)\} = (-1)^nD_s^n[F(s)]

[editar]Convolución


\mathcal{L}\{f*g\}
  = \mathcal{L}\{ f \} \mathcal{L}\{ g \}

[editar]Transformada de Laplace de una función con periodo p


\mathcal{L}\{ f \}
  = {1 \over 1 - e^{-ps}} \int_0^p e^{-st} f(t)\,dt



[editar]Condiciones de convergencia


\mathcal{L}\{(e^{t^2})\} (que crece más rápido que e^{-st}) no pueden ser obtenidas por Laplace, ya que e^{t^2}, es una función de orden exponencial de ángulos.

[editar]Teorema del valor inicial


Sea una función  f\in\varepsilon derivable a trozos y que f^{\prime}\in\varepsilon. Entonces :

f(0^{+})=\lim_{s\to\infty}{sF(s)}

 \varepsilon es el conjunto de funciones continuas a trozos con orden exponencial.

[editar]Teorema del valor fina


Seaf\in\varepsilon una función derivable a trozos tal que f^{\prime}\in\varepsilon.Entonces :


f(\infty)=\lim_{s\to0}{sF(s)}

 \varepsilon es el conjunto de funciones continuas a trozos con orden exponencial.

Tabla de las transformadas de Laplace más comunes

La siguiente tabla provee la mayoría de las transformaciones de Laplace para funciones de una sola variable.

Debido a que la transformada de Laplace es un operador lineal, la transformada de Laplace de una suma es la suma de la transformada de Laplace de cada término

\mathcal{L}\left\{f(t) + g(t) \right\}  = \mathcal{L}\left\{f(t)\right\} + \mathcal{L}\left\{ g(t) \right\}


\mathcal{L}\left\{a f(t)\right\}  = a \mathcal{L}\left\{ f(t)\right\}

4.3 APLICACIONES



Aplicaciones de los Sistemas de ecuaciones diferenciales lineales

Los sistemas de ecuaciones diferenciales lineales encuentran sus aplicaciones en varios problemas que surgen en el sistema del mundo real. Algunos de estos problemas se discuten a continuación.

1. Problema mecánico del acoplamiento de los resortes: Dos cuerpos con masa m1, m2, respectivamente, yacen sobre una mesa. La mesa está libre de fricción. Los dos cuerpos están conectados entre sí con la ayuda de un resorte. Esteresorte está en una posición no estirada. También cada uno de estos cuerpos está conectado a una superficie estática con la ayuda de los resortes. Una vez más, estos resortes no están estirados. La constante elástica de cada uno de los resortes es k1, k2, k3, respectivamente. La situación anterior puede ilustrarse como,



Aquí O1 es la posición inicial del primer cuerpo y O2 es la posición inicial del segundo cuerpo. Los cuerpos pueden ser cambiados de su posición de equilibrio mediante mover cualquiera delos cuerpos en cualquier dirección y luego soltarlos. Unejemplo de estoes,



En la figura anterior, x1 es la cantidad de distancia recorrida por el primer cuerpo cuando este se mueve desde la posición de equilibrio y x2 es la cantidad de distancia recorrida por el segundo cuerpo cuando este se mueve desde la posición de equilibrio. Esto implica que el primer resorte se alarga desde la posición estática por una distancia de x1 y el segundo resorte se alarga desde la posición estática por una distancia de x2 – x1.Esto implica que dos fuerzas restauradoras están actuando sobre el primer cuerpo, estas son:

• La fuerza del primer resorte la cual actúa en dirección izquierda. Esta fuerza por la ley de Hookes igual ak1×1.

• La fuerza del segundo resorte que actúa en dirección derecha. Esta fuerza es igual a k2(x2 – x1).
 Esto nos da la ecuación del movimiento,



 De manera similar, la ecuación del movimiento para el segundo cuerpo es,



Las dos ecuaciones anteriores forman un sistema de ecuaciones diferenciales lineales de segundo orden y pueden resolverse mediante el uso de las técnicas de solución de un sistema de ecuaciones diferenciales lineales de segundo orden.

2. Problemas eléctricos: Muchos de los circuitos eléctricos pueden ser reducidos para solucionar un sistema de ecuaciones diferenciales. Sea un circuito eléctrico dado como,



Ahora, mediante la aplicación de Kirchhoff se tiene la ecuación del flujo de corriente en un nodo como,



Esta ecuación puede ser reducida como,



De manera similar, la ecuación del flujo de corriente del nodo dos se da como,



Esta ecuación puede ser reducida como,



Ahora, aplicando la ley de Kirchoff a la parte izquierda del circuito dado. Por lo tanto tenemos,



Del mismo modo, mediante la aplicación de la ley de Kirchoff a la parte derecha del circuito dado obtenemos,



Ahora, diferencia lasdos últimas ecuaciones para obtener el sistema de ecuaciones como,





Las ecuaciones anteriores pueden ser resueltas para las variables i1, i2y el valor de la variable i puede determinarse con la ayuda de estas dos variables. Un punto importante a mencionar es que pueden existir más que ecuaciones para el ejemplo anterior. Por ejemplo, una de las ecuaciones puede ser,