Wie schließt die Laplace-Transformation das Einschwingverhalten ein?

Wenn wir eine sinusförmige Eingabe liefern, wird die Bewertung der Lösung der Differentialgleichung frontal horrend. Wir geben also einen imaginären, komplexen Exponentialeingang und lösen die Differentialgleichung und nehmen schließlich den Imaginärteil als Lösung (Überlagerungssatz). Sobald wir also die komplexe Amplitude gefunden haben, können wir sie einfach damit multiplizieren e ich ω T , Imaginärteil nehmen und zur Lösung kommen.

Mit diesem netten Trick können wir die sinusförmige stationäre Antwort für das gegebene System finden. Wenn wir die Laplace-Transformation verwenden, erhalten wir, soweit ich das verstehen kann, auch eine komplexe Amplitude, nicht nur für den stationären Zustand, sondern auch für das Einschwingverhalten für jeden Eingang! . Die Idee hinter dem wertvollen Trick, einen sinusförmigen Eingang mit komplexen Zahlen zu modellieren, war vollkommen intuitiv. Aber ich habe keine Ahnung, wie die Laplace-Transformation das alles macht, ich verstehe nicht einmal, warum sie die Verwendung komplexer Zahlen erfordert? vielmehr wie findet es dieses alternative Modell für jeden gegebenen Eingang , wie wir es im vorherigen "Trick" nur für sinusförmigen Eingang getan haben.

@jDAQ danke, aber die Antwort erklärt alles, vorausgesetzt, wir verstehen, warum die Laplace-Transformation eine transiente Antwort enthält, oder?
Liegt das nicht daran, dass Sie Ihren Anfangs- und Endzustand in die Differentialgleichung einbeziehen müssen? Da Sie also das Signal von t=0bis haben t=infinity, ist der Übergang enthalten.

Antworten (3)

Meine Laplace-Transformationsmathematik ist ziemlich eingerostet, aber der nette Trick hinter einer Laplace-Transformation besteht darin, dass die Testwellenform, mit der die Transformation Ihre Eingangsfunktion vergleicht, die Amplitude im Laufe der Zeit ändert.

Untersuchen Sie die Fourier-Transformation:

F ^ ( ω ) = F ( T )   e ich ω T D X ,

Im Wesentlichen dauert es eine Test-Sinuswelle ( e ich ω ) mit Häufigkeit ω und bestimmt, wie ähnlich diese Sinuswelle Ihrer Eingabefunktion 1 ist . Es bestimmt diese Ähnlichkeit oder "Kreuzkorrelation", indem es diese Sinuswelle mit Ihrer Eingabefunktion über alle Zeiten hinweg multipliziert und integriert. Wenn es wenig Ähnlichkeit zwischen dieser Frequenz gibt ω und Ihrer Eingabe, dann ist diese Integration gleich Null.

Wie Sie beobachtet haben, ändert dieses Testsignal niemals die Amplitude und kann nur verwendet werden, um das stationäre Verhalten Ihres Eingangs zu untersuchen.

Vergleichen Sie nun mit der Laplace-Transformation:

F ( S ) = 0 F ( T ) e S T D T

Sehr ähnlich, außer jetzt S ist komplex und ersetzt das Reale (nicht komplex) ω . Zudem ist die Integration nun nicht für alle Zeiten, sondern erst die Zukunft.

Das Schöne daran, zu einer komplexen Potenz zu werden S ist, dass es sowohl eine reale als auch eine imaginäre Komponente gibt, also wenn S = A + ich B , Dann e S T = e A T e ich B T . Jetzt haben wir also eine Sinuswelle, die je nach Wert von entweder exponentiell anwächst, exponentiell schrumpft oder eine konstante Amplitude beibehält A .

Durch dieses zeitlich veränderliche Testsignal und die Tatsache, dass die Integration jetzt erst zum Zeitpunkt 0 beginnt, können wir nun das Einschwingverhalten aus dem Eingangssignal extrahieren.

A Und B bestimmen zusammen einen Punkt auf der komplexen Ebene (hier als s-Ebene bekannt), wohingegen ω wird immer auf der realen Linie sein. Dies ist der zusätzliche Freiheitsgrad, den @hotpaw2 erwähnt hat. Dieser Punkt definiert das von der Laplace-Transformation verwendete Testsignal, und Sie können die Auswirkungen der Bewegung des Punkts um die komplexe s-Ebene in diesem Bild sehen:

s-Ebene für Laplace-Transformationssignale

Diagramme, die die Fourier-Reihenzerlegung eines Signals darstellen, sind üblich (wie das hier: https://en.wikipedia.org/wiki/Fourier_series#/media/File:Fourier_series_and_transform.gif )

Weniger verbreitet sind Diagramme, die eine Laplace-Transformationszerlegung zeigen, aber ich habe dieses Beispiel gefunden, das eine RLC-Schaltung analysiert:Einschwingverhalten der RLC-Schaltung

Auf dem Diagramm oben links sehen Sie den Induktorstrom über die Zeit mit einer großen abklingenden Transiente von 0 bis 5 Sekunden. Auf dem Diagramm unten links können Sie das Ende des anfänglichen Übergangs und das endgültige stationäre Verhalten sehen (beachten Sie die große Änderung des Diagrammmaßstabs!).

Und die entsprechende Laplace-Zerlegung:Laplace-Zerlegung der RLC-Schaltung

Hier können Sie sehen, dass die Gesamtantwort im ersten Diagramm die Summe der beiden hier aufgetragenen Terme ist und dass die Auswirkung des anfänglichen Transienten mit der Zeit auf Null abfällt.

Auf dieser Seite finden Sie viele weitere Informationen:

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

1 : In Wirklichkeit e ich ω = cos X + ich Sünde X , also wird tatsächlich gleichzeitig eine Kosinus- und eine Sinuswelle (ein rotierender komplexer Punkt) als Testfunktion verwendet. Auf diese Weise erhält die Fourier-Transformation Phaseninformationen. Aus diesem Grund zeigt das obige Bild auch geschwungene 3D-Spiralen, da die Laplace-Transformation auch einen rotierenden komplexen Punkt hat. Aber es ist einfacher, es sich einfach als einzelne Sinuswelle vorzustellen.

Bildquellen: https://www.dsprelated.com/freebooks/mdft/Comparing_Analog_Digital_Complex.html https://www.mathworks.com/help/symbolic/solve-differential-equations-using-laplace-transform.html

Ihre Antwort ist die einzige, die die Frage des OP teilweise beantwortet, aber ich denke nicht, dass das hinzugefügte Sigma relevant ist (es lässt das Integral nur für eine große Klasse von Signalen konvergieren). Ich glaube, es ist die auf 0 gesetzte Integrationsgrenze, die den Unterschied ausmacht. Die gesamte Prägressionsgeschichte ist in den Anfangsbedingungen bei 0 eingeschlossen, und diese werden, wie Swedgin betonte, automatisch eingeschlossen, wenn Sie die Differentialgleichungen mit Hilfe der Differenzierungseigenschaft lösen. L{f'} = sF - f(0). Tatsächlich kann IIRC auch eine einseitige Fourier-Transformation definieren, aber die Konvergenz wird problematisch sein.
Wow, vielen Dank, ich habe das nicht ganz verstanden, bitte entblößen Sie mich. Ich kann nicht verstehen, warum die Integration von 0 plötzlich Transienten enthält und warum Fourier-Transformationen keine Transienten enthalten. Können Sie bitte ein Beispiel geben, vielleicht eine RC-Schaltung mit einem 5-V-DC-Eingang.
Wenn Sie eine Differentialgleichung in y transformieren, müssen Sie y und seine Ableitungen y', y'', ... transformieren. Was ist nun die (bilaterale) Fourier-Transformation von y'? Hinweis: math.stackexchange.com/questions/430858/… . Und was ist die (einseitige) Laplace-Transformation von y'? Hinweis: mathalino.com/reviewer/advance-engineering-mathematics/… . Also, in welches können Sie den Anfangswert y(0) einfügen und warum? Die (einseitige) Laplace-Transformation tut dies, und der Grund dafür ist, dass eines der Extreme der Integration 0 ist und f 2 B 0 nicht benötigt.
(Null sein). Wenn Sie also Laplace verwenden, können Sie die Anfangsbedingungen direkt angeben, wenn Sie y(t) in Y(S) verwandeln, oder im Falle einer Schaltung den Kondensator mit einer Anfangsspannung als Impedanz 1/sC mit einer Spannung modellieren Quelle, um den Anfangszustand darzustellen. So oder so, Sie können die Anfangsbedingung von Anfang an einsetzen und die Bestimmung von Konstanten im Ausdruck der allgemeinen vollständigen Lösung vermeiden, indem Sie am Ende Anfangsbedingungen angeben.
@AravindhVasu Wie ich zu Beginn meiner Antwort sagte, ist meine Laplace-Transformationsmathematik jetzt verrostet (viel intuitiver als mathematisch streng), daher werde ich es jetzt nicht wagen, mein eigenes Beispiel zu erstellen. Es gibt viele da draußen. Ich werde in meiner Antwort einen Link zu einem hinzufügen.
@SredniVashtar Ich denke, dass sowohl die Änderung der Integrationsgrenzen als auch die Hinzufügung des reellen Exponentialterms wichtig sind, um die transiente Reaktion zu modellieren. Aus Chus Antwort auf der anderen Seite: „Das Anlegen eines sinusförmigen Eingangs an z. B. ein stabiles SISO-, LTI-System erzeugt einen Ausgang mit einer transienten Komponente und einer stationären Komponente. Der Transient macht das, was alle Transienten tun – er fällt ab Null nach einer gewissen Zeit." Dieser Verfall ist das, was die e A T Begriff darstellt. electronic.stackexchange.com/a/208494/109442
Ich habe gerade eine Reihen-RL-Schaltungs-ODE gesehen, die unter Verwendung der (bilateralen) Fourier-Transformation mit einer Zwangsquelle gelöst wurde, die für t < 0 Null ist und bei = 0 auf den Anfangswert springt. Ich bekomme den Transienten, wie ich es bei Laplace tun würde. Das einzige Problem ist, dass es bei der bilateralen Transformation keine einfache Möglichkeit gibt, andere Anfangsbedingungen festzulegen, die nicht die „natürlichen“ sind. Tatsächlich war die Zwangsfunktion in diesem speziellen Beispiel ... exponentiell (ich brauchte etwas, das das Integral konvergieren lässt). Und es wurde perfekt modelliert, indem nicht abklingende Sinuskurven verwendet wurden. Laplace macht es einfach einfacher, Signale zu modellieren ...
Korrektur: Das Eingangssignal war Exp(-|t|), Symmetrie half beim Lösen. Der Ausgang war am Ende konstant auf dem Anfangswert für t < 0 und sprang dann zu Beginn des Übergangs. Ich werde versuchen, es mit einer einseitigen Fourier-Transformation zu lösen, um zu sehen, ob ich beliebige ICs einstellen kann. Oh, übrigens solltest du deine Variablen überarbeiten (x oder w?)

Wenn Sie eine integrale Transformation vornehmen, projizieren Sie die interessierende Funktion auf ein anderes Koordinatensystem, das immer noch alle Eigenschaften der ursprünglichen Funktion beibehält. Denken Sie daran, wie Sie einen Vektor sowohl in einem rechtwinkligen als auch in einem Polarkoordinatensystem schreiben können, während Sie immer noch derselbe Vektor sind. Das ist es, was Sie mit integralen Transformationen tun, aber anstatt in den zweidimensionalen Polarraum zu gehen, projizieren Sie in einen unendlich dimensionalen Hilbert-Raum. Also, anstatt zu haben X , j usw., Sie haben unendlich viele Basis der Form e J ω .

Um zu verstehen, was passiert, betrachten Sie eine einfache Vektorprojektion. So erhalten Sie die Komponenten desselben Vektors in einem anderen Koordinatensystem.

Projektion

Quelle

Um die Komponente des Vektors zu erhalten A im B Richtung nehmen Sie einfach das Skalarprodukt

A C = ich A ich C ich

mit dem Einheitsvektor C , neben B . Wenn der Vektor C komplex ist, ist das Skalarprodukt mit seinem komplexen Konjugat. Somit

A C = ich A ich C ich .

Wenn Sie eine Funktion projizieren, haben Sie unendlich viele Komponenten A ( T ) . Mit einer integralen Transformation, wie der Fourier-Transformation, haben Sie unendlich viele Basis e J ω . Jede Basis hat eine Komponente e J ω T . Damit ist die Summe vorbei ich wird das Integral über T

A ( T ) ( e J ω T )   D T = A ( T ) e J ω T   D T .

Wenn Sie also eine integrale Transformation vornehmen, haben Sie dieselbe Funktion, aber in einem anderen Koordinatensystem. Da es sich um dieselbe Funktion handelt, enthält sie das Einschwingverhalten in beiden Koordinatensystemen.

Wie beantwortet dies die Frage des OP? Das Eingangssignal enthält nicht unbedingt Informationen über den Anfangszustand des Systems (es sei denn, sie sind alle Null). Die Transferfunktion auch nicht. Wie passen diese Bedingungen in Ihr Bild? Zum Beispiel mit einer RLC-Schaltung, bei der in einem Fall C 5 V hat, L bei t = 0 1 A hat und in einem anderen Fall Vc (0) = -1 V und IL (0) = 0,2 mA ... Wo sind diese Bedingungen? notwendig, um die Transienten zu berechnen, gespeichert? Außerdem erwägen Sie die Integration von minus unendlich bis plus unendlich. Wo platzieren Sie Anfangsbedingungen?
@SredniVashtar Es gibt viele Möglichkeiten, das RLC-System in Ihrem Beispiel einzurichten. Wir müssen das System so aufstellen, dass es bei t = 0 nicht stabil ist. Andernfalls wäre es bereits in seinen stationären Zustand konvergiert. Eine Möglichkeit, dies zu erreichen, besteht darin, das System bei den Anfangsbedingungen für t < 0 zu halten. Dies bedeutet, dass es immer eine Art Fenstereffekt gibt, z. B. setzen wir den Strom durch den Kondensator für t < 0 auf Null, indem wir a einführen Schrittfenster. Meine Antwort besagt lediglich, dass die integrale Transformation alle Informationen beibehält, daher habe ich die spezifischen verwendeten Techniken nicht erklärt.
Mein Einwand ist, dass solche Informationen überhaupt nicht vorhanden sind. Wenn Sie die ODE im Zeitbereich lösen, erhalten Sie eine Lösung in Form beliebiger Konstanten, deren Wert durch Auferlegen der Anfangsbedingungen bestimmt wird. Wenn Sie es mit Laplace lösen, geben Sie die Anfangsbedingung ein, wenn Sie die ODE in AE transformieren. Die Anfangsbedingungen (entweder in Form des vorherigen Ladezustands/Stroms in den dynamischen Komponenten oder als durch die Anregungen auferlegte Anfangswerte) bestimmen den Einschwingvorgang, sind jedoch normalerweise nicht in der Übertragungsfunktion (oder im Eingangssignal) enthalten.
Andernfalls sollte man für jeden möglichen Satz von ICs eine Übertragungsfunktion angeben. Wenn wir die Differentialgleichung nicht direkt transformieren, werden die ICs stattdessen verwendet, um die Koeffizienten der generischen natürlichen Antwort zu bestimmen, die den Polen der Übertragungsfunktion zugeordnet sind. Trotzdem müssen sie hinzugefügt werden, wenn auch zu einem späteren Zeitpunkt der Lösung.
@SredniVashtar Abhängig von den Anfangsbedingungen erhalten Sie unterschiedliche Übertragungsfunktionen, die sich jedoch nur auf die Gleichstromfrequenz in Bezug auf Dirac-Deltas auswirken. Das Dirac-Delta wird eingeführt, wenn Sie eine Integration im Zeitbereich haben. Siehe hier: thefouriertransform.com/transform/fourierintegralproperty.gif .
@SredniVashtar Ich sollte auch hinzufügen, dass Sie je nach Anfangsbedingungen Phasenverschiebungen erhalten. Die Tatsache, dass Sie für unterschiedliche Anfangsbedingungen eine andere Fourier-Transformation erhalten, sollte erwartet werden. Denn wenn das Zeitbereichssignal anders ist, dann ist auch seine Fourier-Transformation anders.
Ich hätte den Begriff "Netzwerkfunktion" anstelle von "Übertragungsfunktion" verwenden sollen. Dies ist die Funktion, die ein Netzwerk charakterisiert. Wenn Sie also jedes Element durch sein Gegenstück im Frequenzbereich ersetzen und Spannungs- und Stromquellen verwenden, um die Anfangsbedingungen darzustellen, muss ich Ihnen zustimmen, dass Sie für jeden Satz von ICs eine andere Übertragungsfunktion erhalten können. Aber wenn man das tut, ist es normalerweise einfacher, nach den interessierenden Variablen ohne eine Übertragungsfunktion zu lösen. Oder, wie ich es bei Franco gesehen habe, kann man die Ausgabe in Bezug auf die Netzwerkfunktion und eine Funktion von ICs ausdrücken.

Die Steady-State-Fourier-Analyse ist praktisch, wenn Sie eine einzelne unendlich lange reine Sinuskurve mit einer einzigen Frequenz haben. Sie können die Reaktion eines LTI-Systems auf diese unendlich lange Sinuskurve leicht berechnen, da sie einem einzelnen Punkt im Fourier-Bereich zugeordnet ist.

Das Problem besteht darin, dass Signale in der realen Welt eine endliche Länge (eine Hüllkurve) haben, die einem unendlichen Stapel (Integral) von Frequenzen im Fourier-Bereich zugeordnet wird. Und die Berechnung der Antwort auf eine unendliche Reihe von Frequenzen könnte eine unendliche Menge an Berechnungen erfordern, was sehr lange dauern könnte.

Die Laplace-Transformation fügt einen weiteren Freiheitsgrad hinzu, indem sie die Zeit in ihren integralen Kern einbezieht. Es kann somit eine unendliche Reihe von Frequenzen in einem einzigen komplexen (oder 2D-) Punkt darstellen und ermöglicht so die Berechnung einer Systemantwort auf diese unendlichen Reihen von unendlichen Sinuskurven mit weniger als einer unendlichen Menge an Tafel. Transiente Signale können in Chargen dieser unendlichen Chargen zerlegt werden.

Betrachten Sie es (Laplace-Transformation) also als einen mathematischen Trick, um eine unendliche Menge an Einzelfrequenz-Steady-State-Analysen (Fourier-Transformation) in endlicher Zeit (und Tafel) durchzuführen, indem Sie einen weiteren Freiheitsgrad hinzufügen.

Hinzugefügt: Beachten Sie, dass die Fourier-Transformation einer abfallenden Exponentialfunktion unendlich lang ist. Unendlich im Fourier-Raum, aber ein einzelner (komplexer) Punkt im Laplace-Raum.

Vielen Dank für die Antwort, aber ich verstehe nicht: "Die Laplace-Transformation fügt einen weiteren Freiheitsgrad hinzu, indem sie die Zeit in ihren integralen Kern einbezieht. Sie kann daher eine unendliche Menge von " diesem Teil darstellen. Was meinst du damit, dass es eine unendliche Menge von Frequenzen darstellen kann?
Es stellt ein unendliches Integral dar und nicht nur über eine einzelne Sinusfrequenz. Während ein Punkt in der Fourier-Domäne nur für eine Frequenz steht.