Wie kann man die Einschwingzeit eines überdämpften Systems abschätzen?

Ich habe versucht herauszufinden, wie ich die Einschwingzeit eines Systems zweiter Ordnung als Reaktion auf eine Schritteingabe der Größe 5 abschätzen kann. Die Übertragungsfunktion des Systems ist

G ( S ) = 1 ( S + 2 ) ( S + 4 )
und das Zeitverhalten habe ich bereits mit der Sprungeingabe R(s) ermittelt:
C ( S ) = R ( S ) G ( S ) C ( T ) = 5 8 + 5 8 e 4 T 5 4 e 2 T
Jetzt muss ich die 2% Einschwingzeit der Antwort anhand dieser Informationen schätzen, aber ich bin mir nicht sicher, wie. Ich weiß, dass das System als ζ> 1 überdämpft ist, daher kann ich die normale Einschwingzeitgleichung nicht verwenden
T S = 4 ζ ω N
Ich habe mir diesen Beitrag angesehen: ( Über- und kritisch gedämpfte Systemeinschwingzeit ), aber die Antworten erklären nur langwierige Wege, um ein genaues Ergebnis zu erhalten. Ich habe bereits MATLAB verwendet, um ein genaues Ergebnis von 2,3 Sekunden zu erhalten, aber ich muss es ohne MATLAB schätzen können.

Ich dachte, ich könnte Trial-and-Error mit verschiedenen Werten von t ausprobieren, bis c (t) innerhalb von 2% des Steady-State-Werts liegt (was 0,625 ist), aber obwohl dies funktionieren würde, bezweifle ich, dass es der richtige Weg ist Wer kann mir mit einer besseren Methode helfen?

Antworten (3)

Bei Systemen mit echten Polen in der linken Halbebene können Sie sie normalerweise abschätzen, indem Sie nur den dominanten Pol (den Pol mit der niedrigsten Frequenz) berücksichtigen. In Ihrem Fall wäre dies P D = 2 . Das Ergebnis wird genauer, wenn der nicht dominante Pol ( P N D ) entfernt sich weiter vom dominanten Pol.

Indem Sie nur den dominanten Pol betrachten, erhalten Sie eine ziemlich einfache Gleichung:

5 4 e 2 T = 0,02 5 8 T = 1 2 ln ( 0,02 5 8 4 5 ) 2.30258509299 S

Die Idee ist, dass der nicht dominante Pol bei P N D = 4 führt zu einem Begriff e 4 T das so schnell dämpft, dass es die Gesamteinschwingzeit nicht beeinflusst. Der Vorteil ist die Einfachheit der Gleichung und die Tatsache, dass es in elektronischen Schaltungen tatsächlich ziemlich häufig vorkommt, dass es einen sehr dominanten Pol und weit entfernte nicht dominante Pole gibt.

In Ihrem speziellen Fall ist es möglich, die Einschwingzeit analytisch zu berechnen. Die Zeit, die die zeitabhängigen Terme benötigen, um auf 2 % des Endwerts abzufallen, kann wie folgt berechnet werden (ähnlich wie bei Andys Antwort, aber unter Verwendung des Absolutwerts):

| e 4 T 2 e 2 T | = 0,02 ( j = e 2 T ) j 2 2 j = ± 0,02 ( Es gibt 4 verschiedene Lösungen, aber ich nehme nur die relevante ) j = e 2 T = 1 7 5 2 T = 1 2 ln ( 1 7 5 2 ) 2.30006613189 S

Also Faktor 2 für P N D / P D führt zu einem Fehler von etwa 0,1 % bei der berechneten Einschwingzeit, wenn die Dominant-Pol-Näherung verwendet wird. Ob dies ausreicht, überlasse ich Ihnen.

Das ist perfekt danke!

Ja, Ihre inverse Laplace-Berechnung ist korrekt.

Der endgültige stationäre Wert beträgt 5/8 – dies ist der DC-Wert nach langer Zeit. Sie suchen also wirklich nach dem Rest der Gleichung, der in der Größenordnung auf 2 % von 5/8 fällt: -

5 8 e 4 T 5 4 e 2 T = 5 8 0,02

= 8 8 e 4 T 8 4 e 2 T = 8 8 0,02

= e 4 T 2 e 2 T = 0,02

Hilft das?

Ich war damals ziemlich nah dran, da ich gerade etwas Ähnliches gemacht habe, aber ich habe 0,98 anstelle von 0,02 verwendet. Warum verwenden Sie 2 % und nicht 98 %? Wenn wir 2 % weniger als den Endwert wollen, wäre es doch 0,98 x 5/8?
@Sam, dein Ansatz und meiner sind gleichermaßen gültig, es ist nur so, dass ich so dachte und du den anderen.
Wir müssten verwenden 0,02 anstatt + 0,02 um die richtige Einschwingzeit zu finden (oder die Koeffizienten der Exponentiale umkehren oder den von @Sam vorgeschlagenen Wert von 0,98 verwenden und den DC-Term einbeziehen). Wahrscheinlich ein Versehen, aber vielleicht die Quelle der Verwirrung?

Nun, lassen Sie uns das in einem allgemeineren Fall lösen. Wir haben die folgende Übertragungsfunktion (unter der Annahme eines echten positiven Werts für ϵ ):

(1) H ( S ) := Y ( S ) X ( S ) = 1 ( S + ϵ ) ( S + 2 ϵ )

Wenn wir uns die Sprungantwort ansehen, die wir verwenden X ( S ) = L T [ θ ( T ) ] ( S ) = 1 S , also ist die Ausgabe gegeben durch:

(2) Y ( S ) = 1 S 1 ( S + ϵ ) ( S + 2 ϵ )

Unter Verwendung der inversen Laplace-Transformation finden wir:

(3) j ( T ) = L S 1 [ 1 S 1 ( S + ϵ ) ( S + 2 ϵ ) ] ( T ) = exp ( 2 ϵ T ) ( exp ( ϵ T ) 1 ) 2 2 ϵ 2

Es ist nicht schwer zu zeigen, wann T (unter der Annahme eines echten positiven Werts für ϵ ), wir bekommen:

(4) lim T j ( T ) = 1 2 ϵ 2

Nun, für die Eingewöhnungszeit wollen wir die Zeit finden T Wenn N % des Endwertes erreicht ist:

(5) j ( T N ) = N % 100 1 2 ϵ 2   T N =

Das Lösen ergibt:

(6) T N = 1 ϵ ln ( 10 10 N % )