Bedeutung des "Pols" einer Übertragungsfunktion mit Zeitverzögerung

Ich habe in der grundlegenden Steuerungstheorie gelernt, dass wir die Stabilität eines LTI-Systems durch die Vorzeichen des Realteils der Pole seiner Übertragungsfunktion bestimmen können H ( S ) = N ( S ) D ( S ) . Für ein rationales System, wo beides N ( S ) Und D ( S ) Polynome s sind, ist dies leicht verständlich, da ein positiver Realteil des Pols eine exponentielle Zunahme des Zeitbereichs bewirkt und somit instabil ist. Meine Frage ist, ob das System dann irgendwo eine Verzögerung von T Sekunden hat H ( S ) enthält einen Faktor exp ( S T ) , manchmal im Nenner. Warum können wir in diesem Fall die Stabilität immer noch bestimmen, indem wir den Realteil der "Pole" von betrachten H ( S ) (Punkte in der s-Ebene, wo D ( S ) = 0 )?

Ein Beispiel:

N ( S ) = S

D ( S ) = A S 2 + B S + C + D exp ( S T )

wobei a, b, c, d und T Konstanten sind.

In diesem Fall können wir keine Partialbruchentwicklung durchführen, wie wir es für rationale Übertragungsfunktionen tun. Dann, wie man das Zeitbereichsverhalten mit den Wurzeln von verknüpft D ( S ) ?

Antworten (2)

Ich glaube, die Eigenschaft, auf die Sie sich beziehen, ist

L { F ( T τ ) } = e S τ L { F ( T ) }

Dies bedeutet, dass Ihre Übertragungsfunktion dieser Regel nicht wirklich folgt und keine wirkliche Zeitverzögerung hat:

H ( S ) = S A S 2 + B S + C + D e S T

Das Lösen der inversen Laplace-Transformation für diese Funktion wird höchstwahrscheinlich numerische Methoden beinhalten. Was möglich ist , ist so etwas zu haben:

H ( S ) = S + D e S T A S 2 + B S + C

Weil es getrennt werden kann in:

H ( S ) = S A S 2 + B S + C + e S T D A S 2 + B S + C

Sie können sehen, dass der zweite Term auch in Teilbrüche zerlegt werden kann, wie der erste Term. Die Pole des Nenners sagen dann etwas über die Stabilität der Übertragungsfunktion aus.

Stabilität im Allgemeinen

Die Formel für die inverse Laplace-Transformation lautet

F ( T ) = 1 2 π J lim T σ J T σ + J T F ( S ) e S T D S

Wo σ größer als alle Singularitäten von gewählt wird F ( S ) auf der komplexen Ebene (in unserem Fall so, dass sie alle Pole umfasst).

Dieses Integral wird äquivalent durch den Cauchy-Residuensatz gelöst :

F ( T ) = L 1 { F ( S ) } = A l l   P Ö l e S   Ö F   F ( S ) R e S [ F ( S ) e S T ]
Denken Sie daran, dies ist allgemein , dh. es funktioniert immer für eine Übertragungsfunktion mit beliebigen Singularitäten!

Es spielt also keine Rolle, ob die Pole von einem Polynom oder von einer transzendenten Funktion stammen. Solange sich die Singularitäten alle in der Ebene der linken Hälfte befinden, enthält ihr Rest immer eine Exponentialfunktion, die auf 0 und nicht auf unendlich abfällt. Alle Singularitäten im RHP führen immer zu einer Exponentialfunktion, die explodiert.


Anhang

Das kann man anmerken e S T kann in seiner Taylor-Reihe erweitert werden:

e z = N = 0 + z N N !

Das bedeutet also, dass Ihre Beispielübertragungsfunktion geschrieben werden kann als

H ( S ) = S A S 2 + B S + C + D ( N = 0 + ( S T ) N N ! ) = S N = 0 + A N S N

Wo A 0 = C + D , A 1 = B D T , A 2 = A + D T 2 2 , A N = D ( T ) N N ! , N > 2 .

Diese Art von Übertragungsfunktion hat also überall "Pole". Dies soll veranschaulichen, dass Sie dies definitiv nicht wie eine normale Übertragungsfunktion zweiter Ordnung behandeln können.


Anhang

Der Rest eines (einfachen) Pols ist

R e S S = A [ H ( S ) ] = lim S A ( S A ) F ( S )

Für Pole mit Multiplizität N :

R e S S = A [ H ( S ) ] = 1 ( N 1 ) ! lim S A D N 1 D S N 1 ( ( S A ) N H ( S ) )

Die Taylor-Erweiterung in Ihrem Anhang ist meiner Meinung nach ein guter Punkt, da sie zeigt, dass es viele Pole über dem Raum geben wird, wie Sie vorgeschlagen haben. In meinem Fall ist die Exponentialfunktion tatsächlich im Nenner, da die Verzögerung in einer Rückkopplungsschleife auftritt. Numerische Methoden sind im Allgemeinen erforderlich, da stimme ich zu, aber ich würde gerne wissen, warum wir konzeptionell die Stabilität bestimmen können, indem wir uns die Wurzeln des Nenners ansehen.
@GeorgeC Ich habe meine Antwort auf das bearbeitet, was Sie meiner Meinung nach gemeint haben.
George, die Wurzeln des Nenners sind die Pole der Übertragungsfunktion. Und keiner dieser Pole darf nicht in der rechten Hälfte der S-Ebene (RHP) erscheinen. Daher beobachten wir die Positionen aller Pole, während wir die Verstärkung variieren. Mit dieser Methode können wir sehen, bei welcher Verstärkung sich die Pole zum RHP bewegen (Instabilität). Dies zeigt an, wo die Stabilitätsgrenze liegt.

Der Ausdruck für eine Zeitverzögerung (Totzeit) im Frequenzbereich exp(-sT) ist eine transzendente Funktion. Sie kann nicht mit konzentrierten Elementen modelliert werden und kann daher nicht direkt mit einer Übertragungsfunktion kombiniert werden, die aus Polynomen in "s" besteht. Zu diesem Zweck kann es jedoch durch den folgenden Ausdruck angenähert werden (Pade-Approximation zweiter Ordnung):

H D ( S ) = 1 S T 2 + S 2 T 2 12 1 + S T 2 + S 2 T 2 12

Dies ist eine rationale Funktion (Allpass) und erscheint als Faktor in der Schleifenverstärkung und erscheint natürlich in der Übertragungsfunktion der geschlossenen Schleife. Hier kann er wie jeder andere Block in der geschlossenen Schleife behandelt werden.

Ja, ich habe gelernt, dass Pade eine gängige Technik zur Annäherung von Verzögerungen ist. Aber selbst wenn wir Pade nicht machen, bestimmen wir trotzdem die Stabilität, indem wir die "Pole" der ursprünglichen Übertragungsfunktion finden. Beispielsweise ist die Übertragungsfunktion im offenen Regelkreis H(s)=KG(s), wobei G(s)=exp(-sT)/(s(s+1)). Das Lehrbuch zeichnet nur den Wurzelort des Pols mit geschlossener Schleife gegenüber der Verstärkung K (keine Pade-Näherung). Warum können wir das tun?
Okay - ich weiß was du meinst. In meiner Antwort habe ich die Pade-Näherung zweiter Ordnung angegeben. Genauer ist natürlich die dritte Ordnung oder sogar die 15. Ordnung. Ein Simulator kann höhere Ordnungen als n = 2 verarbeiten - aber ich weiß nicht genau, was der Simulator tut. Ich nehme an, die Näherungsordnung ist auf einen beliebigen großen Wert begrenzt, und schließlich verbindet der Simulator alle berechneten Pole und ist somit in der Lage, einen Wurzelort zu erzeugen.