Zero-Order-Hold und Z-Transformation

J. W. (arx) OP 26.11.2010 20:46 #1952719

Lesenswert?

•

Hallo,

ich versuche gerade eine Übertragungsfunktion zu diskretisieren und 
bekomme irgendwie nicht das Richtige Ergebnis raus:

Meine kontinuierliche Übertragungsfunktion ist ein IT1-Glied

$H(s) = \frac{70}{s^2 + 134 s}$


Ich erweitere um das Zero-Order-Glied und erhalte

$H_0(s) = \frac{70}{s^2 + 134 s} \cdot \frac{1 - e^{-Ts}}{s} = \frac{1}{s^2} \cdot \frac{70}{s + 134} \cdot (1 - e^{-Ts})$


Die Z-Transformation davon liefert

$H_0(z) = \frac{T z}{(z - 1)^2} \cdot \frac{70 z}{z - e^{-134 T}} \cdot (1 - z^{-1}) = \frac{T}{(z - 1)^2} \cdot (z - 1) \cdot \frac{70 z}{z - e^{-134 T}} = \frac{70 T z}{z^2 - (e^{-134 T} + 1) z + e^{-134 T}}$


Mit einer Abtastzeit von T = 0.001s erhalte ich dann

$H_0(z) = \frac{0.07 z}{z^2 - 1.875 z + 0.875}$


Matlab gibt aber nach dem Aufruf von c2d(H(s), 0.001, 'zoh') folgendes 
aus

$H_m(z) = \frac{3.353 \cdot 10^{-5} z + 3.206 \cdot 10^{-5}}{z^2 - 1.875 z + 0.875}$


Der Zähler stimmt also schon mal, aber der Nenner ist anders. Ein Test 
durch eine Sprungantwort zeigt schnell: Matlab hat recht, ich nicht. Ich 
weiß aber nicht, was ich falsch gemacht habe. Hat jemand einen Tipp?

Antwort Zitieren Bearbeiten Thread verschieben Thread sperren Anmeldepflicht aktivieren Anpinnen Thread löschen Thread mit anderem zusammenführen Beitrag melden

J. W. (arx) OP 27.11.2010 00:41 #1952947

Lesenswert?

•

▲
▼

Ok, die Lösung führt hier wohl über die Partialbruchzerlegung. Ist etwas 
mühsamer, bringt aber am Ende das gesuchte Ergebnis...

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Florian Werner Gast 17.05.2017 16:45 #5013671

Lesenswert?

•

▲
▼

Ich war eben mit der selben Problematik konfrontiert. Ich glaube die 
Lösung gefunden zu haben. Auch wenn der Thread schon alt ist poste ich 
das jetzt hier einfach mal, vielleicht bringt es ja jemanden weiter:

Ich habe die folgende Formel auf 
https://de.wikipedia.org/wiki/Sample-and-Hold-Schaltung gefunden, die 
auf den ersten Blick funktioniert:

$G_{h}(z)=(1-z^{{-1}})\cdot Z\left\{{\frac {G(s)}{s}}\right\}$


Hierbei muss man aufpassen, dass man den Term

$\frac{G(s)}{s}$

 nicht in Faktoren aufteilt, die einzeln Transformiert und das ganze in 
der Z-Welt wieder zusammenmultipliziert. Man darf höchstens konstante 
Parameter ausmultiplizieren, aber keine Terme mit s.
Das das nicht zulässig ist erkennt man schon daran, dass

$\mathcal{Z}\{ \frac{1}{s^2} \} = \mathcal{Z}\{ \frac{1}{s} \} \cdot \mathcal{Z}\{ \frac{1}{s} \} = (\frac{1}{1-z^{-1}})\cot(\frac{1}{1-z^{-1}}) = \frac{1}{(1-z^{-1})^2}$


nicht zum selben Ergebnis führt wie

$\mathcal{Z}\{ \frac{1}{s^2} \} = \frac{Tz^{-1}}{(1-z^{-1})^2}.$

Die erste Vorgehensweise ist falsch, die zweite richtig.

Für dein Beispiel bedeutet das also:

$H(s) = \frac{a}{s^2+bs} = \frac{a}{s(s+b)}; \; a=70; \; b=134$

$G_{h}(z)=(1-z^{{-1}}) Z\left\{{\frac{a}{s^2(s+b)}}\right\} = (1-z^{{-1}})\cdot Z\left\{{\frac{a}{b^2}\frac{b^2}{s^2\cdot(s+b)}}\right\}$


Aus Transformationstabelle:

$\mathcal{Z}\{\frac{a^2}{s^2(s+a)}\} = \frac{[(aT-1+e^{-aT})+(1-e^{-aT}-aTe^{-aT})z^{-1}]z^{-1}}{(1-z^{-1})^2(1-e^{-aT}z^{-1})}$


Also:

$G_{h}(z)=(1-z^{{-1}}) Z\left\{{\frac{a}{s^2(s+b)}}\right\} = (1-z^{{-1}})\frac{a}{b^2} \frac{[(bT-1+e^{-bT})+(1-e^{-bT}-bTe^{-bT})z^{-1}]z^{-1}}{(1-z^{-1})^2(1-e^{-bT}z^{-1})} =$

$\frac{a}{b^2} \frac{[(bT-1+e^{-bT})+(1-e^{-bT}-bTe^{-bT})z^{-1}]z^{-1}}{(1-z^{-1})(1-e^{-bT}z^{-1})} = \frac{[\frac{a}{b^2}(bT-1+e^{-bT})]z^{-1}+[\frac{a}{b^2}(1-e^{-bT}-bTe^{-bT})]z^{-2}}{1-(1+e^{-bT})z^{-1}+e^{-bT}z^{-2}} =$

$\frac{3.3488\cdot 10^{-5}\cdot z^{-1}+3.2025\cdot 10^{-5}\cdot z^{-2}}{1+1.8746\cdot z^{-1}+0.8746\cdot z^{-2}} = \frac{3.3488\cdot 10^{-5}z+3.2025 10^{-5}}{z^2+1.8746\cdot z+0.8746}$


Tadaaaa ;)

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Thread beobachten

Antwort schreiben

Bitte melde dich an, um einen Beitrag zu schreiben.

Bestehender Account

oder

Mit Google-Account einloggen

Noch kein Account?

Die Registrierung ist kostenlos und dauert nur eine Minute.

Jetzt registrieren