¿Cómo incluye la transformada de Laplace la respuesta transitoria?

Cuando proporcionamos una entrada sinusoidal, la evaluación de la solución de la ecuación diferencial de frente se vuelve horrenda. Así que damos una entrada exponencial compleja e imaginaria y resolvemos la ecuación diferencial y finalmente tomamos la parte imaginaria como la solución (teorema de superposición). Entonces, una vez que encontramos la amplitud compleja, podemos simplemente multiplicarla por mi i ω t , tomar parte imaginaria y llegar a la solución.

Podemos encontrar la respuesta de estado estable sinusoidal para el sistema dado, con este ingenioso truco. Cuando usamos la transformada de Laplace, hasta donde puedo entender, también nos da una amplitud compleja, no solo para el estado estable, sino que también incluye la respuesta transitoria, ¡ para cualquier entrada! . La idea detrás del precioso truco de modelar una entrada sinusoidal con números complejos era perfectamente intuitiva. Pero no tengo idea de cómo la transformada de Laplace hace todo esto, ni siquiera entiendo, ¿por qué exige el uso de números complejos? más bien, ¿cómo encuentra este modelo alternativo para cualquier entrada dada , como hicimos en el "truco" anterior solo para la entrada sinusoidal?

@jDAQ gracias, pero la respuesta lo explica todo, suponiendo que entendamos por qué la transformada de Laplace incluye una respuesta transitoria, ¿verdad?
¿No es porque tienes que incluir tu estado inicial y final en la ecuación diferencial? Entonces, dado que tiene la señal de t=0a t=infinity, el transitorio está incluido.

Respuestas (3)

Mis matemáticas de transformada de Laplace están bastante oxidadas, pero el buen truco detrás de una transformada de Laplace es que la forma de onda de prueba con la que la transformada compara su función de entrada cambia la amplitud con el tiempo.

Examine la transformada de Fourier:

F ^ ( ω ) = F ( t )   mi i ω t d X ,

Esencialmente, se necesita una onda sinusoidal de prueba ( mi i ω ) con frecuencia ω y determina qué tan similar es esa onda sinusoidal a su función de entrada 1 . Determina esta similitud, o "correlación cruzada", al multiplicar esa onda sinusoidal con su función de entrada a lo largo de todo el tiempo e integra. Si hay poca similitud entre esta frecuencia ω y su entrada, entonces esta integración será igual a cero.

Como ha observado, esta señal de prueba nunca cambia de amplitud y solo puede usarse para examinar el comportamiento de estado estable de su entrada.

Ahora compare con la transformada de Laplace:

F ( s ) = 0 F ( t ) mi s t d t

Muy similar, excepto ahora. s es complejo y reemplaza lo real (no complejo) ω . Además, la integración ahora no es para siempre, sino solo el futuro.

Lo bueno de elevar a una potencia compleja s es que hay un componente tanto real como imaginario, de modo que si s = a + i b , entonces mi s t = mi a t mi i b t . Así que ahora tenemos una onda sinusoidal que crece exponencialmente, se reduce exponencialmente o mantiene una amplitud constante dependiendo del valor de a .

Esta señal de prueba que puede cambiar con el tiempo, más el hecho de que la integración ahora solo comienza en el tiempo 0 significa que ahora podemos extraer la respuesta transitoria de la señal de entrada.

a y b juntos determinan un punto en el plano complejo (conocido aquí como el plano s), mientras que ω siempre estará en la línea real. Este es el grado adicional de libertad que mencionó @hotpaw2. Este punto define la señal de prueba utilizada por la transformada de Laplace, y puede ver los efectos de mover el punto alrededor del plano s complejo en esta imagen:

Plano s para señales de transformada de Laplace

Los diagramas que representan la descomposición en series de Fourier de una señal son comunes (como el que se muestra aquí: https://en.wikipedia.org/wiki/Fourier_series#/media/File:Fourier_series_and_transform.gif )

Menos comunes son los diagramas que muestran una descomposición por transformada de Laplace, pero encontré este ejemplo analizando un circuito RLC:Respuesta transitoria del circuito RLC

En el gráfico superior izquierdo, puede ver la corriente del inductor a lo largo del tiempo, con un gran transitorio decreciente de 0 a 5 segundos. En el gráfico inferior izquierdo, puede ver el final del comportamiento transitorio inicial y el estado estable final (¡observe el gran cambio en la escala del gráfico!).

Y la correspondiente descomposición de Laplace:Descomposición de Laplace del circuito RLC

Aquí puede ver que la respuesta total en el primer gráfico es la suma de los dos términos trazados aquí, y que el impacto del transitorio inicial se reduce a cero con el tiempo.

Consulte esta página para obtener mucha más información:

https://en.wikipedia.org/wiki/Laplace_transform#Formal_definition

1 : En realidad, mi i ω = porque X + i pecado X , por lo que en realidad está usando una onda coseno y una sinusoidal simultáneamente (un punto complejo giratorio) como función de prueba. Así es como la transformada de Fourier obtiene información de fase. Esa es también la razón por la que la imagen de arriba muestra espirales rizadas en 3D, ya que la transformada de Laplace también tiene un punto complejo giratorio. Pero es más fácil pensar en ello como una sola onda sinusoidal.

Fuentes de imágenes: https://www.dsprelated.com/freebooks/mdft/Comparing_Analog_Digital_Complex.html https://www.mathworks.com/help/symbolic/solve-differential-equations-using-laplace-transform.html

La suya es la única respuesta para abordar parcialmente la pregunta del OP, pero no creo que la sigma agregada sea relevante (solo hace que la integral converja para una gran clase de señales). Creo que es el límite de integración colocado en 0 lo que marca la diferencia. Toda la historia previa está encerrada en las condiciones iniciales en 0 y estas, como señaló Swedgin, se incluyen automáticamente cuando resuelves las ecuaciones diferenciales por medio de la propiedad de diferenciación. L{f'} = sF - f(0). De hecho, IIRC también puede definir una transformada de Fourier unilateral, pero la convergencia será problemática.
Wow, muchas gracias, no he entendido completamente esto, por favor, sé sincero conmigo. No puedo entender por qué la integración desde 0 de repente incluye transitorios y por qué las transformadas de Fourier no incluyen transitorios. ¿ Puede dar un ejemplo, tal vez un circuito RC con una entrada de CC de 5V?
Cuando transformas una ecuación diferencial en y necesitas transformar y y sus derivadas y', y'',... Ahora, ¿cuál es la Transformada de Fourier (bilateral) de y'? Sugerencia: math.stackexchange.com/questions/430858/… . ¿Y cuál es la transformada de Laplace (unilateral) de y'? Sugerencia: mathalino.com/reviewer/advance-engineering-mathematics/… . Entonces, ¿cuál te permite poner el valor inicial y(0) y por qué? La transformada de Laplace (unilateral) sí lo hace, y la razón es que uno de los extremos de la integración es 0 y f no necesita 2 B 0.
(ser cero). Entonces, cuando usa Laplace, puede especificar las condiciones iniciales directamente cuando convierte y (t) en Y (S), o en el caso de un circuito, puede modelar el capacitor con un voltaje inicial como impedancia 1 / sC con un voltaje fuente para representar el estado inicial. De cualquier manera, puede poner la condición inicial desde el principio y evitar determinar constantes en la expresión de la solución general completa especificando las condiciones iniciales al final.
@AravindhVasu Como dije al comienzo de mi respuesta, mi matemática transformada de Laplace está oxidada ahora (mucho más intuitiva que matemáticamente rigurosa), por lo que no me aventuraré a crear mi propio ejemplo en este momento. Hay muchos por ahí. Agregaré un enlace a uno en mi respuesta.
@SredniVashtar Creo que tanto el cambio de los límites de integración como la adición del término exponencial real son importantes para modelar la respuesta transitoria. De la respuesta de Chu en la otra página: "Aplicar una entrada sinusoidal a, por ejemplo, un sistema estable, SISO, LTI producirá una salida con un componente transitorio y un componente de estado estable. El transitorio hace lo que hacen todos los transitorios: decae a cero después de cierto tiempo". Esa decadencia es lo que mi a t término representa. electronics.stackexchange.com/a/208494/109442
Acabo de ver una ODE de circuito RL en serie resuelta usando la transformada de Fourier (bilateral) con una fuente de forzamiento que es cero para t<0 y salta al valor inicial en = 0. Obtengo el transitorio, como lo haría con Laplace. El único problema es que al ser la transformada bilateral no hay manera fácil de establecer otras condiciones iniciales que no sean las 'naturales'. De hecho, la función forzada en este ejemplo específico era... exponencial (necesitaba algo que hiciera converger la integral). Y fue modelado perfectamente usando sinusoides que no decaen. Laplace simplemente hace que sea más fácil modelar señales...
corrección: la señal de entrada era Exp(-|t|), la simetría ayudó a resolver. La salida terminó siendo constante en el valor inicial para t<0 y luego saltó al comienzo del transitorio. Intentaré resolverlo con una transformada de Fourier unilateral para ver si puedo establecer circuitos integrados arbitrarios. Oh, por cierto, deberías revisar tus variables (¿x o w?)

Cuando toma una transformación integral, está proyectando la función de interés a un sistema de coordenadas diferente que aún conserva todas las propiedades de la función original. Piense en cómo puede escribir un vector en un sistema de coordenadas tanto rectangular como polar, sin dejar de ser el mismo vector. Esto es lo que está haciendo con las transformaciones integrales, pero, en lugar de ir al espacio polar bidimensional, está proyectando en un espacio de Hilbert de dimensión infinita. Entonces, en lugar de tener X , y etc., tienes infinitas bases de la forma mi j ω .

Para comprender lo que está sucediendo, considere una proyección vectorial simple. Así es como obtienes las componentes de un mismo vector en un sistema de coordenadas diferente.

proyección

fuente

Para obtener la componente del vector a en el b dirección simplemente toma el producto punto

a C = i a i C i

con el vector unitario C , Paralelo a b . Si el vector C es complejo, el producto escalar es con su complejo conjugado. Por eso

a C = i a i C i .

Cuando estás proyectando una función, tienes una cantidad infinita de componentes. a ( t ) . Con una transformada integral, como la transformada de Fourier, tienes infinitas bases mi j ω . Cada base tiene un componente. mi j ω t . Así la suma sobre i se convierte en la integral sobre t

a ( t ) ( mi j ω t )   d t = a ( t ) mi j ω t   d t .

Entonces, cuando toma una transformación integral, se queda con la misma función pero en un sistema de coordenadas diferente. Como es la misma función, contiene la respuesta transitoria en ambos sistemas de coordenadas.

¿Cómo responde esto a la pregunta del OP? La señal de entrada no necesariamente lleva información sobre el estado inicial del sistema (a menos que sean todos cero). La función de transferencia tampoco. ¿Dónde encajan estas condiciones en tu imagen? Por ejemplo con un circuito RLC donde en un caso C tiene 5V, L tiene 1A en t=0, y otro caso donde Vc(0)=-1V e IL(0)=.2 mA... ¿Dónde están esas condiciones? necesario para calcular el transitorio, almacenado? Además, está considerando la integración de menos infinito a más infinito. ¿Dónde pones las condiciones iniciales?
@SredniVashtar Hay muchas formas de configurar el sistema RLC en su ejemplo. Necesitamos configurar el sistema de tal manera que no sea estable en t = 0. De lo contrario, ya habría convergido a su estado estacionario. Una forma de lograr esto es mantener el sistema en las condiciones iniciales para t < 0. Esto significa que siempre hay algún tipo de efecto de ventana, por ejemplo, configuramos la corriente a través del capacitor para que sea cero para t < 0, introduciendo un ventana de paso Mi respuesta simplemente establece que la transformación integral conserva toda la información, por lo que no expliqué las técnicas específicas utilizadas.
Mi objeción es que dicha información no existe en absoluto, en primer lugar. Cuando resuelve la EDO en el dominio del tiempo, obtiene una solución en términos de constantes arbitrarias cuyo valor se determina imponiendo las condiciones iniciales. Cuando lo resuelves con Laplace, ingresas la condición inicial cuando transformas la ODE en AE. Las condiciones iniciales (ya sea en forma de estado previo de carga/corriente en los componentes dinámicos o como valores iniciales impuestos por las excitaciones) determinan el transitorio pero normalmente no se incluyen en la función de transferencia (ni en la señal de entrada).
De lo contrario, se debe especificar una función de transferencia para cada conjunto posible de circuitos integrados. En cambio, si no transformamos directamente la ecuación diferencial, los IC se utilizan para determinar los coeficientes de la respuesta natural genérica asociados con los polos de la función de transferencia. Aún así, deben agregarse, incluso en una etapa posterior de resolución.
@SredniVashtar Obtiene una función de transferencia diferente según las condiciones iniciales, pero solo afecta la frecuencia de CC en términos de deltas de Dirac. El delta de Dirac se introduce cuando tiene integración en el dominio del tiempo. Consulte aquí: thefouriertransform.com/transform/fourierintegralproperty.gif .
@SredniVashtar También debo agregar que obtienes cambios de fase según las condiciones iniciales. Se debe esperar el hecho de que obtenga una transformada de Fourier diferente para diferentes condiciones iniciales. Después de todo, si la señal en el dominio del tiempo es diferente, su transformada de Fourier también será diferente.
Debería haber usado el término "función de red" en lugar de "función de transferencia". Esta es la función que caracteriza una red. Entonces, si reemplaza cada elemento con su contraparte en el dominio de la frecuencia y coloca las fuentes de voltaje y corriente para representar las condiciones iniciales, estoy de acuerdo con usted, que para cada conjunto de circuitos integrados puede obtener una función de transferencia diferente. Pero, por lo general, cuando uno hace eso, es más fácil resolver las variables de interés sin una función de transferencia. O, como lo he visto en Franco, uno puede expresar la salida en términos de la función de red y una función de circuitos integrados.

El análisis de Fourier de estado estacionario es útil si tiene una sola sinusoide pura de frecuencia única infinitamente larga. Puede calcular fácilmente la respuesta de un sistema LTI a esta sinusoide infinitamente larga porque se asigna a un solo punto en el dominio de Fourier.

El problema es que en el mundo real, las señales tienen una longitud finita (una envolvente), que se asigna a un lote infinito (integral) de frecuencias en el dominio de Fourier. Y calcular la respuesta a un lote infinito de frecuencias podría requerir una cantidad infinita de cálculos, lo que podría llevar mucho tiempo.

La transformada de Laplace agrega otro grado de libertad al incluir el tiempo en su núcleo integral. Por lo tanto, puede representar un lote infinito de frecuencias en un solo punto complejo (o 2D), lo que permite calcular una respuesta del sistema a estos lotes infinitos de sinusoides infinitas con menos de una cantidad infinita de pizarra. Las señales transitorias se pueden descomponer en lotes de estos lotes infinitos.

Así que considérelo (transformada de Laplace) como un truco matemático para hacer una cantidad infinita de análisis de estado estacionario de frecuencia única (transformada de Fourier) en tiempo finito (y pizarra), agregando otro grado de libertad.

Agregado: tenga en cuenta que la transformada de Fourier de un exponencial en descomposición es infinitamente larga. Infinito en el espacio de Fourier, pero un solo punto (complejo) en el espacio de Laplace.

Gracias por responder, pero no entiendo, "La transformada de Laplace agrega otro grado de libertad, al incluir el tiempo en su núcleo integral. Por lo tanto, puede representar un lote infinito de" esta parte. ¿Qué quiere decir con que puede representar un lote infinito de frecuencias?
Representa una integral infinita, y no solo sobre una sola frecuencia de sinusoide. Mientras que un punto en el dominio de Fourier es para una sola frecuencia.