Z-Transformation

Die z-Transformation ist eine mathematische Methode, um zyklisch abgetastete kontinuierliche Signale im Zusammenhang mit dynamischen Systemen im z-Bereich berechenbar zu machen .

Die z-Transformation ist aus der Laplace-Transformation entstanden und hat auch ähnliche Eigenschaften und Berechnungsregeln. Sie erlaubt für die digitale Signalverarbeitung mit Computern (Mikrocomputern) die Berechnung von impulsförmig abgetasteten und digitalisierten Signalfolgen (Wertefolgen) in Verbindung mit kontinuierlichen dynamischen Systemen G(s) oder den zugehörigen systembeschreibenden linearen Differenzengleichungen y_{(k)}.

Ein Vorteil der Anwendung der z-Transformation ergibt sich, wenn eine Wertefolge und eine systembeschreibende Differenzengleichung in eine algebraisch zusammengefasste z-Übertragungsfunktion überführt wird. Die Übertragungsfunktion des z-Bereiches G(z) dient der Systemanalyse und dem Systemverhalten. Der Verlauf der Systemausgangsgröße y_{(k)} kann bei gegebener Eingangsgröße u_{(k)} durch verschiedene Methoden der inversen z-Transformation in den zeitdiskreten Bereich f_{{(k)}} und dann im Zeitbereich f(t) dargestellt werden.

Die Anwendung der z-Transformation bezieht sich ausschließlich auf Digitalrechner-geführte Anlagen, die Steuer- und Regelungstechnik und digitale Filter.

Einführung in die z-Transformation

Geschichtliche Entwicklung

Die grundsätzlichen Ideen zur z-Transformation gehen auf Pierre-Simon Laplace zurück und wurde 1947 durch Witold Hurewicz zur Lösung von linearen Differentialgleichungen mit konstanten Koeffizienten verwendet.

Ursprünglich wurde sie als „Laplace-Transformation von Abtastfunktionen“ eingeführt, im Jahr 1952 erfolgte die heute übliche Begriffsfestlegung z-Transformation durch John R. Ragazzini und Lotfi A. Zadeh bei Arbeiten mit zeitdiskreten Daten im Rahmen der Regelungstechnik an der Columbia University.

Die Modifizierte z-Transformation geht auf Arbeiten von Eliahu Ibrahim Jury aus dem Jahr 1958 zurück.

Definition der Systemgrößen

Folgende Systemgrößen werden für die verschiedenen f(t)-, F(s)-, f_{{(k)}}- und F(z)-Bereiche verwendet:

Werden in der Regelungstechnik Regler und Regelstrecke gleichzeitig betrachtet, wird die Eingangsgröße als e(t) (Regelabweichung) bezeichnet, die Ausgangsgröße u(t) ist gleichzeitig die Eingangsgröße der Regelstrecke. Die Ausgangsgröße der Regelstrecke ist y(t). Im Zeitbereich werden diese Größen als Kleinbuchstaben, im Bildbereich als Großbuchstaben geschrieben.
Bei der rekursiven Simulation mit Differenzengleichungen 1. Ordnung eines dynamischen Systems mit dem Digitalrechner erfolgt keine Abtastung. Jedes berechnete Glied von k der endlichen Ausgangsfolge bezieht sich auf das um k-1 zurückliegende Glied. Auf diese Weise entstehen Eingangsfolgen und Ausgangsfolgen aller benutzten Differenzengleichungen von Teilsystemen für eine bestimmte Folge von k. Die gleichen Differenzengleichungen werden bis k_{\mathrm {max} } wiederholt berechnet und gespeichert.
Die Abtastfolge k=(0,1,2,3,\dots ) bedeutet eine Nummerierung der Folgeglieder der Wertefolge des Eingangssignals (Eingangsfolge) oder des Ausgangssignals (Ausgangsfolge) eines Systems.

Grundlagen z-Transformation

Die Laplace-Transformation ist ein mathematisches Verfahren der Systemtheorie zur Behandlung und Berechnung von kontinuierlichen Signalen und linearen zeitinvarianten dynamischen Systemen.

Die z-Transformation ist ein mathematisches Verfahren der Systemtheorie zur Behandlung und Berechnung von kontinuierlich abgetasteten Signalen und linearen zeitinvarianten zeitdiskreten dynamischen Systemen. Ein zeitdiskretes dynamisches System wird durch Differenzengleichungen oder als z-Transformierte beschrieben.

Zeitinvariante lineare dynamische Systeme g(t) haben mindestens einen Systemeingang und einen Systemausgang für das Eingangssignal u(t) und Ausgangssignal y(t). Das Systemverhalten wird im Zeitbereich durch gewöhnliche Differentialgleichungen beschrieben. Mit der Laplace-Transformation kann eine gewöhnliche Differentialgleichung in den sogenannten Bildbereich G(s) (s-Bereich, komplexer Frequenzbereich) überführt werden. Das Systemverhalten im Bildbereich wird durch die s-Übertragungsfunktion bestimmt und kann algebraisch behandelt werden.

Diskrete Zeit: Die Zeitdiskretisierung eines kontinuierlichen Signals oder eines dynamischen zeitinvarianten Übertragungssystems bedeutet der Übergang der Berechnung eines kontinuierlichen Signals oder Systems f(t) mit unendlicher hoher Auflösung zu einem Signal oder Systems f(kT_{A}) mit einer endlichen Auflösung eines fortlaufenden konstanten Zeitintervalls T_A. Das Zeitintervall T_A muss genügend klein sein, damit dominante Systembewegungen auch erfasst werden können, bzw. der Approximationsfehler gegenüber dem Verlauf der analytischen Funktion gering ist. Es werden hier zur Kennzeichnung der physikalischen Unterschiede der Zeitdiskretisierung folgende Definitionen festgelegt:

Ideale Abtastung eines kontinuierlichen Signals f(t) im Abstand T_A durch eine Sample-and-Hold-Schaltung nullter Ordnung.

Abtastfolge und Wertefolge: In der Mathematik wird eine Auflistung von endlich und unendlich vielen fortlaufend nummerierten Objekten als Folge bezeichnet.

Die Abtastfolge k=(0,1,2,3,\dots ) bedeutet eine Nummerierung der Folgeglieder der Wertefolge des Eingangssignals (Eingangsfolge) und des Ausgangssignals (Ausgangsfolge) eines Systems.

Eine Wertefolge besteht aus k_{\mathrm {max} } oder k=\infty vielen Folgegliedern. Das Objekt mit der Nummer i wird i-tes Folgeglied oder i-te Komponente der Folge genannt. Die abgetasteten und digitalisierten Signale entsprechen einer Folge von modulierten impulsförmigen Signalen u(kT_{A}) im Abstand T_A, die erst nach der A/D-Wandlung mit einer Haltestufe zu einem gestuften quasi kontinuierlichem Signal u(t) werden.

Die z-Transformierten und die Laplace-Transformierten verhalten sich ähnlich bei der algebraische Behandlung der Einzelsysteme des z- und s-Bereichs und in den Transformations-Rechenregeln. Sowohl im z- wie auch im s-Bereich lassen sich die Übertragungsfunktionen G(z) oder G(s) der dynamischen Systeme mittels Nullstellenzerlegung der Systemanalyse und Systemsynthese unterziehen.

Trotz der ähnlichen Eigenschaften der z-Transformierten im Vergleich zur s-Transformierten besteht bei der Behandlung von Signalen und dynamischen Systemen mit der z-Transformation ein höherer mathematischer und signaltechnischer Schwierigkeitsgrad.

Während die Laplace-Transformation sich mit Signalen und dynamischen Systemen f(t) des Zeitbereichs zur Wandlung in den Bildbereich F(s) und umgekehrt, mit der inversen Laplace-Transformation in den Zeitbereich befasst, hat die Behandlung mit der z-Transformation zum z-Bereich weitere mathematische Zwischenformen. Es handelt sich bei der z-Transformation um die Beziehungen der Signale und dynamische Systeme vom Zeitbereich in den:

Bezeichnung der Transformationen im f(t)-, F(s)-, f_{{(k)}}- und F(z)-Bereich:

F(s)={\mathcal  L}\{f(t)\}
f(t)={\mathcal  L}^{{-1}}\{F(s)\}
F(z)={\mathcal  Z}\{f_{{(k)}}\}
f_{{(k)}}={\mathcal  Z}^{{-1}}\{F(z)\}

Z-Transformation von Abtastsignalen

Die z-Transformation eines abgetasteten Signals F^{*}(s)={\mathcal  L}\{f^{*}{(t)}\} entspricht dem Austausch der komplexen Variable s=\sigma +j\omega der s-Transformierten durch die komplexe z-Variable z=e^{{sT_{A}}}. Dadurch wird eine unendliche Summe der Exponentialterme in eine Potenzreihe in z^{{-k}} überführt (spezielle Laurent-Reihe).

{\displaystyle F(z)={\mathcal {Z}}\{f^{*}{(t)}\}=\left.{\mathcal {L}}\{f^{*}{(t)}\}\ \right|_{z=e^{sT_{A}}}=\sum _{k=0}^{\infty }f(kT_{A})\cdot e^{-kT_{A}s}=\sum _{k=0}^{\infty }f(kT_{A})\cdot z^{-k}=f_{0}+f_{1}\cdot z^{-1}+f_{2}\cdot z^{-2}+f_{3}\cdot z^{-3}+\dots }

Für die z-Transformation eines zeitdiskreten Eingangssignals u(kT_{A}) in vereinfachter Schreibweise u_{(k)} gilt:

U(z)=\sum_{k=0} ^\infty u_{(k)} \cdot z^{-k}

Grundlagen Differenzengleichungen für lineare zeitinvariante Systeme

Für die numerische Berechnung des Systemverhaltens eines dynamischen Systems g(t) oder im Zusammenhang mit der z-Transformation werden Differenzengleichungen {\displaystyle y_{(k)}=f(\mathrm {System} ,k,k-1,\Delta t,u_{(k)})} benötigt. Mit ihrer Hilfe lässt sich das Systemverhalten, die Systemausgangsgröße y(t), dynamischer Systeme g(t) für ein gegebenes Eingangssignal u(t) im zeitdiskreten Bereich f(k\Delta t) berechnen.

Differenzengleichungen entstehen meist aus systembeschreibenden gewöhnlichen Differenzialgleichungen, deren Differentialquotienten durch Differenzenquotienten ersetzt werden. Die kontinuierlichen mathematischen Operationen der Integration und Differentiation werden zeitdiskret durch Summen- und Differenzenbildung angenähert.

Die Differenzengleichungen beschreiben mit dem Approximationsalgorithmus für ein kleines Zeitintervall \Delta t die Signaländerungen nach jedem Zeitintervall als Funktion des betreffenden Teilsystems (Linearfaktoren im s-Bereich) und des Eingangssignals. Mit der fortlaufenden Wiederholung der Berechnung mit dem Zeitintervall \Delta t und Addition der Änderungsergebnisse zum vorherigen Ergebnis ergibt sich der Signalverlauf eines Systems über die Zeit k_{{\mathrm  {max}}}\cdot \Delta t.

Es bestehen verschiedene mathematische Verfahren, zeitkontinuierliche Systeme in zeitdiskrete Systeme zu beschreiben und umzuwandeln.

Differenzengleichungen der einfachsten Art beziehen sich auf die den Linearfaktoren der Übertragungsfunktion G(s) zugehörigen Differenzialgleichungen erster Ordnung, deren Differentialquotienten durch Differenzenquotienten ersetzt werden. Diese Beziehung ist von großer Bedeutung, weil nur 4 verschiedene Typen von Differenzengleichungen erster Ordnung existieren, mit denen alle Formen von linearen Übertragungssystemen gebildet werden können, auch solche mit Schwingungsanteilen mit konjugiert komplexen Polen oder Nullstellen. Diese Teilsysteme können beliebig multiplikativ, additiv, zurück gekoppelt oder strukturell vermascht sein und gelten sowohl für den s-Bereich wie auch im diskreten Zeitbereich.

Meistens wird zur Aufstellung der Differenzengleichungen das explizite Euler-Rückwärtsverfahren der Rechteckapproximation als einfachstes Verfahren verwendet. Nach diesem Verfahren können aus den zugehörigen Differenzialgleichungen der 4 Elementarsysteme G(s) erster Ordnung der Übertragungsfunktionen Differenzengleichungen gebildet werden, indem an Stelle des Differenzialquotienten mit dy/dt der Differenzenquotient {\frac  {\Delta y}{\Delta t}}={\frac  {y_{{(k)}}-y_{{(k-1)}}}{\Delta t}} näherungsweise eingeführt wird.

In der Regel wird davon ausgegangen, dass die inneren Systemspeicher des Übertragungssystems sich im Ruhezustand befinden und die Anfangswerte bei t = 0 für y_{{(0)}}=0 und alle Ableitungen von y_{{(0)}} Null sind.

Beispiel der Entwicklung der Differenzengleichung der Integration (I-Glied) aus der Differenzialgleichung:

Die Übertragungsfunktion des I-Gliedes lautet:
\ G(s)={\frac  {Y(s)}{U(s)}}={\frac  1{T_{I}\cdot s}}

Die zugehörige Differenzialgleichung lautet:

T_{I}\cdot {\dot  y}(t)=u(t)

Der Differenzenquotient wird an Stelle des Differenzialquotienten {\dot  y}(t) eingesetzt.
Anstelle der kontinuierlichen Systemgrößen y(t) und u(t) treten die aktuellen zeitdiskreten Werte y(k\Delta t)=y_{{(k)}} und u(k\Delta t)=u_{{(k)}}:

T_{I}\cdot {\frac  {y_{{(k)}}-y_{{(k-1)}}}{\Delta t}}=u_{{(k)}}.

Damit lautet die nach y_{(k)} umgestellte Differenzengleichung des I-Gliedes:

y_{{(k)}}=y_{{(k-1)}}+u_{{(k)}}\cdot {\frac  {\Delta t}{T_{I}}}

In gleicher Weise können die Differenzengleichungen von Systemen erster Ordnung aus den zugehörigen Differenzialgleichungen abgeleitet werden.

Tabelle der Differenzengleichungen (Euler-Rückwärts) der Elementarsysteme G(s) erster Ordnung
Elementarsysteme P-Glied I-Glied D-Glied PD1-Glied PT1-Glied
Übertragungsfunktion {\frac {Y}{U}}(s)=K {\frac {Y}{U}}(s)={\frac {1}{T\cdot s}} {\frac {Y}{U}}(s)=T\cdot s {\frac {Y}{U}}(s)=K\cdot (T\cdot s+1) {\frac {Y}{U}}(s)={\frac {K}{T\cdot s+1}}
Differenzengleichungen
y_{(k)}=\,
 K \cdot u_{(k)} y_{(k-1)}+u_{(k)}\cdot {\frac {\Delta t}{T}} [u_{(k)}-u_{(k-1)}]\cdot {\frac {T}{\Delta t}} K_{{PD1}}\cdot \left[u_{{(k)}}+(u_{{(k)}}-u_{{(k-1)}})\cdot {\frac  {T}{\Delta t}}\right] y_{(k-1)}+[K_{PT1}\cdot u_{(k)}-y_{(k-1)}]\cdot {\frac {\Delta t}{T+\Delta t}}

(Mit K = Verstärkungsfaktor, y_{(k)} = aktuelle Ausgangsgröße, y_{(k-1)} = vorherige Ausgangsgröße, T = Zeitkonstante, u_{(k)} = aktuelle Eingangsgröße)

Die einmalige Anwendung einer Differenzengleichung y_{{(k)}}=f(T,\Delta t,k,k-1,u_{{(k)}}) zum Zeitpunkt k\cdot \Delta t ergibt für eine gegebene Eingangsfolge u_{(k)} ein Folgeglied der Ausgangsfolge y_{(k)}. Jedes Folgeglied y_{(k)} bezieht sich auf eine zurückliegende Folge y_{(k-1)}. Deshalb wird eine solche Differenzengleichung als Rekursionsgleichung bezeichnet, weil jedes Folgeglied eine Funktion des vorherigen Folgegliedes ist.

Die rekursive Anwendung von Differenzengleichungen zur Berechnung von Eingangs-Wertefolgen zu Ausgangs-Wertefolgen bedeutet die angenäherte Lösung der systembeschreibenden Differentialgleichung des Systemausgangssignals y_{(k)} von Wertefolgen (Berechnungspunkten) y_{{(k)}}=y(0),y(1),y(2),y(3),\dots .

Mit Hilfe eines Personal Computers kann das Systemverhalten eines dynamisches Systems oder eines Regelkreises mit Differenzengleichungen vollständig simuliert werden. Dabei wird eine endliche Anzahl k_{\mathrm {max} } von Berechnungsfolgen (Wertefolgen) festgelegt und die Rechenergebnisse der Teilsysteme - das Systemverhalten - tabellarisch und grafisch als Berechnungspunkte im Abstand \Delta t dargestellt. Die Differenzengleichung enthält bereits die Lösungsvorschrift der Systemausgangsgröße in Annäherung an die systembeschreibende Differentialgleichung.

Handelt es sich bei dem dynamischen System um eine Hardware mit einer im zeitlichen Abstand T_A abgetasteten Eingangsfolge u_{(k)}, die über einen Mikrocomputer mit Differenzengleichungen zu einer Ausgangsfolge y_{(k)} berechnet wird, kann mit Hilfe eines Haltegliedes eine treppenförmig gestufte quasi kontinuierliche Ausgangsgröße y(t) als Beispiel der prinzipiellen Funktion eines digitalen Reglers erreicht werden. Regelstrecken liegen in der Praxis meist als kontinuierliche Systeme vor, die eine kontinuierliche Stellgröße benötigen.

Differenzengleichungen höherer Ordnung

Differenzengleichungen können auch aus gewöhnlichen Differenzialgleichungen höherer Ordnung entwickelt werden, wenn ab dem Zeitpunkt t=kT_{A}\geq 0 die letzten vergangenen Ausgangs-Wertefolgen mit k=[-1\to -n] und die Eingangs-Wertefolgen mit k=[-1\to -m] bekannt sind.

Zeitkontinuierliche lineare Systeme werden im Zeitbereich durch die gewöhnlichen Differenzialgleichungen n-ter Ordnung mit konstanten Koeffizienten für n\geq m beschrieben. Dabei sind n und m die höchsten Ableitungen der Ausgangssignale y(t) und Eingangssignale u(t).

Eine gegebene gewöhnliche Differentialgleichung wird durch den Koeffizienten a_{0} dividiert, um y(t) freistellen zu können. Diese Form der Differentialgleichung wird entsprechend der dargestellten Koeffizienten wie folgt neu geordnet.

y(t)+a_{{1}}{\dot  y}(t)+a_{{2}}{\ddot  y}(t)+\dots +a_{n}y^{{(n)}}(t)=b_{0}u(t)+b_{1}{\dot  u}(t)+b_{2}{\ddot  u}(t)+\dots +b_{m}u^{{(m)}}(t).

Diese Differentialgleichung kann in eine Differenzengleichung überführt werden:

Jede Ableitung der Systemgrößen wird im zeitdiskreten Bereich entsprechend der Ordnung als zurückliegende Folgeglieder der Eingangs- und Ausgangsfolgen k-1 bis k-n oder k-m berücksichtigt.

Daraus folgt die Differenzengleichung:

{y_{{(k)}}+a_{1}\cdot y_{{(k-1)}}+a_{2}\cdot y_{{(k-2)}}+\dots +a_{n}\cdot y_{{(k-n)}}=b_{0}\cdot u_{{(k)}}+b_{{1}}\cdot u_{{(k-1)}}+b_{2}\cdot u_{{(k-2)}}+\dots +b_{m}\cdot u_{{(k-m)}}}.

Damit kann die allgemeine Form der Differenzengleichung nach y_{(k)} aufgelöst werden:

{y_{{(k)}}=b_{0}\cdot u_{{(k)}}+b_{{1}}\cdot u_{{(k-1)}}+b_{2}\cdot u_{{(k-2)}}+\dots +b_{m}\cdot u_{{(k-m)}}-a_{1}\cdot y_{{(k-1)}}-a_{2}\cdot y_{{(k-2)}}-\dots -a_{n}\cdot y_{{(k-n)}}}.

Für die numerische Berechnung eines dynamischen Systems wird die s-Übertragungsfunktion oder die zugehörige Differentialgleichung benötigt. Die Umsetzung einer systembeschreibenden Differentialgleichung in eine angenäherte Differenzengleichung zur Beschreibung von Eingangsfolgen und Ausgangsfolgen eines dynamischen Systems wird ermöglicht, wenn die Differentiale der Differentialgleichung durch Rückwärts-Differenzenquotienten über die Abtastperiode ersetzt werden.

Die folgenden Ableitungen der Differentialquotienten in Differenzenquotienten der 1. 2. und 3. Ordnung sind gegeben:

Differenzenquotient 1. Ordnung:

{\dot  y}(t)={\frac  {dy}{dt}}\approx {\frac  {\Delta y}{\Delta t}}={\frac  {\Delta y_{k}}{T_{A}}}={\frac  {y_{{(k)}}-y_{{(k-1)}}}{T_{A}}}

Der Differenzenquotient 2. Ordnung entsteht aus Differenzen der Differenz:

\ddot y(t)=\frac {d^2y}{dt^2} \approx \frac {\Delta^2y_k}{T_A^2} = \frac {y_{(k)}-2y_{(k-1)} +y_{(k-2)} }{T_A^2}

Der Differenzenquotient 3. Ordnung lautet:

y^{(3)}(t)=\frac {d^3y}{dt^3} \approx \frac {\Delta^3y_k}{T_A^3} = \frac {y_{(k)}-3y_{(k-1)}+3y_{(k-2)}-y_{(k-3)}}{T_A^3}

Nach erfolgtem Einsetzen der Differenzenquotienten in die Differenzengleichung eines dynamischen Systems lassen sich die neuen Koeffizienten aus den Koeffizienten der Differentialgleichung berechnen.

Beispiel der Entwicklung einer Differenzengleichung zur Berechnung der Sprungantwort eines PT_{{2kk}}-Gliedes mit konjugiert komplexen Polen:

Sprungantwort eines PT2-Gliedes mit konjugiert komplexen Polen. TA = 0,01 s
0{,}25{\ddot  y}(t)+0{,}125{\dot  y}(t)+y(t)=1(t)
Gegeben: Übertragungsfunktion im s-Bereich:
G(s)={\frac  {Y(s)}{U(s)}}={\frac  1{0{,}25s^{2}+0{,}125s+1}}
Sprungfunktion: u(t)=1

Gesucht: Differenzengleichung y_{(k)} zur numerischen Bestimmung des System-Zeitverhaltens.

Zugehörige Differentialgleichung nach dem Differentiationssatz der Laplace-Transformation:

0{,}25{\ddot  y}(t)+0{,}125{\dot  y}(t)+y(t)=u(t)

Die Differenzenquotienten für {\ddot  y}(t) und {\dot  y}(t) werden in die nachfolgende Differenzengleichung eingesetzt:

0{,}25\cdot {\frac  {y_{{(k)}}-2y_{{(k-1)}}+y_{{(k-2)}}}{T_{A}^{2}}}+0{,}125\cdot {\frac  {y_{{(k)}}-y_{{(k-1)}}}{T_{A}}}+y_{{k}}=u_{{k}}\quad |\qquad y(t)\to y_{{(k)}};\ u(t)\to u_{{(k)}}

Die Brüche werden in einzelne additive Terme aufgelöst, um y_{(k)} freistellen zu können:

 y_{(k)} \bigg(\frac {0{,}25}{T_A^2} + \frac{0{,}125}{T_A} + 1 \bigg) = u_{(k)} + \frac {0{,}5y_{(k-1)}}{T_A^2}- \frac {0{,}25y_{(k-2)}}{T_A^2} + \frac {0{,}125y_{(k-1)}}{T_A}
 y_{(k)}=\frac {u_{(k)}}{\lambda} + \frac {0{,}5y_{(k-1)}}{\lambda T_A^2} - \frac {0{,}25y_{(k-2)}}{\lambda T_A^2}+ \frac {0{,}125y_{(k-1)}}{\lambda T_A}\quad \bigg|\qquad \lambda =\frac {0{,}25}{T_A^2}+\frac{0{,}125}{T_A}+1

Berechnungsbeispiel für einige Werte der Ausgangsfolge y_{(k)} mit T_A = 0{,}01; \quad \lambda = 2513{,}5; \quad u_{(k)}=1:

y_{(k=0)} = 0{,}0004 + 0 - 0 + 0 = 0{,}0004.
y_{(k=1)} = 0{,}0004 + 0{,}0008 - 0 + 0 = 0{,}0012.
y_{(k=2)} = 0{,}0004 + 0{,}0024 - 0{,}0004 + 0 = 0{,}0024.
y_{(k=157)} = 0{,}0004 + 3{,}286 - 1{,}643 + 0{,}008 = 1{,}652.

Diese Differenzengleichung entspricht einem Rekursionsalgorithmus eines dynamischen Systems, der schrittweise mit einem digitalen Rechner gelöst werden kann.

Die rekursive Berechnung der Differenzengleichung 2. Ordnung bezieht sich für die aktuelle Ausgangsfolge y_{(k)} durch Einsetzen der zurückliegenden Werte der Ausgangsfolge u_{{(k-1)}} und u_{(k-2)} in die Gleichung. Für das 1. Folgeglied der Berechnungsfolge k=0 sind die zurückliegenden Werte der Ausgangsfolge noch nicht verfügbar und damit Null. Die Anzahl k=(0, 1, 2, 3, \dots k_\mathrm{max}) der Glieder der Ausgangsfolge wird durch die diskrete Zeit T_A und durch die gewünschte zu beobachtende Gesamtzeit des Einschwingvorgangs bestimmt.

Z-Transformation zeitdiskreter Signale und dynamischer Systeme

Die z-Transformation wird auf zeitdiskrete Signale u_{(k)}, auf die systembeschreibenden Differenzengleichungen y_{(k)} und auf Übertragungsfunktionen des s- und z-Bereiches meist mit Hilfe der Korrespondenztabellen angewendet.

Funktionsblöcke eines digitalen Reglers.

Das nebenstehende Bild ist ein Beispiel der Darstellung der Signalarten und Systeme an einem aufgeschnittenen digitalen Regelkreis mit einer kontinuierlichen Regelstrecke.

Abtastsysteme wandeln in Verbindung mit A/D-Wandlern ein kontinuierliches Signal in ein zeitdiskretes Signal als Wertefolge um. D/A-Wandler in Verbindung mit Haltesystemen nullter Ordnung wandeln eine Wertefolge in ein gestuftes zeitkontinuierliches Signal um.

Rechenregeln der z-Transformation

Wie die Laplace-Transformation wird die z-Transformation durch Sätze und Rechenregeln definiert, dennoch bestehen in einigen Funktionen große Unterschiede.

Die Transformationen und Rücktransformationen der z-Transformation erfolgen meist mit Hilfe von Transformations-Tabellen. In der Fachliteratur werden in Tabellen die Zeitfunktionen f(t), die Laplace-Transformierten f(s) und die z-Transformierten f(z) dargestellt. Nicht vorhandene Zeitfunktionen für die inverse z-Transformation können wie bei der Laplace-Transformation durch die Partialbruchzerlegung bestimmt werden.

Die Verfahren der Anwendung der Rechenregeln der z-Transformationen sind umfangreich und können in diesem Abschnitt nur angedeutet werden. Die nachstehenden meist tabellarisch aufgeführten Verfahren sind in jedem guten Fachbuch der Regelungstechnik oder in Vorlesungsmanuskripten der z-Transformation zu finden.

Die nachfolgenden mathematischen Beziehungen gelten für Systeme mit einem Eingangssignal u_{(k)} und einem Ausgangssignal y_{(k)}. Für einen digitalen Regler müssen die zugehörigen Größen des Eingangssignals e_{{(k)}} und Ausgangssignals u_{(k)} in die Gleichungen eingesetzt werden.

Vergleich der z-Transformation und der Laplace-Transformation.
Mitz=e^{{s\cdot T_{A}}} wird die s-Ebene auf die z-Ebene abgebildet.

Die wichtigsten Eigenschaften der z-Transformation:

z-Variable:
z=e^{{T_{A}\cdot s}}
z-Transformierte einer Folge f_{{(k)}}:
F(z)={\mathcal  Z}\{f_{{(k)}}\}
Inverse z-Transformation von F(z):
f_{{(k)}}={\mathcal  Z}^{{-1}}\{F(z)\}
Für den s-Bildbereich gilt für die s-Variable.
{\displaystyle \lim _{t\to +0}f(t)=\lim _{s\to \infty }s\cdot F(s)}
{\displaystyle \lim _{t\to \infty }f(t)=\lim _{s\to 0}s\cdot F(s)}
Für den z-Bildbereich gilt für die z-Variable:
{\displaystyle f(k=0)=\lim _{z\to \infty }F(z)}
{\displaystyle \lim _{k\to \infty }f(k)=\lim _{z\to 1}(z-1)\cdot F(z)}
u(k=0)=\lim _{{z\to \infty }}U(z)
Totzeiten innerhalb eines Digitalreglers entstehen durch die Signalabtastung über A/D-Wandler, D/A-Wandler und durch die Rechenzeit des Mikrocomputers. Die Berechnung einer Totzeit T_{t}=d\cdot T_{A} einer Abtastfolge entspricht einer Rückwärtsverschiebung (Rechtsverschiebung) der Abtastfolge (Zahlenwerte) um d Abtastschritte von T_A. Diese mathematische Operation bedeutet im s-Bereich für die Totzeit e^{{-s\cdot T_{t}}} und im z-Bereich:
u_{{(k-d)}}\ \circ \!\!-\!\!\bullet \ z^{{-d}}\cdot U(z)
Eine Vorhersage um die Zeit T_{P}=d\cdot T_{A} einer Wertefolge entspricht einer Vorwärtsverschiebung nach links um d-Abtastschritte. Diese Operation entspricht der Laplace-Transformierten e^{{s\cdot T_{A}}} und bedeutet im z-Bereich:
u_{{(k+d)}}\ \circ \!\!-\!\!\bullet \ z^{{d}}\cdot U(z)
Die Differenz zweier aufeinander folgender Abtastwerte dividiert durch die Abtastzeit entspricht der Annäherung an einen Differenzialquotienten im Zeitbereich.
u_{{(k)}}-u_{{(k-1)}}\ \circ \!\!-\!\!\bullet \ {\frac  {z-1}z}\cdot U(z)
Wird die Summe aller Abtastwerte mit der Abtastzeit TA multipliziert, entsteht in Annäherung an den Zeitbereich die numerische Integration. Im s-Bereich entspricht die Integration 1 / s. Im z-Bereich gilt die Integration:
\sum _{{i=0}}^{k}u_{{(i)}}\ \circ \!\!-\!\!\bullet \ {\frac  z{z-1}}\cdot U(z)
Die z-Übertragungsfunktion G(z) ist die z-Transformierte der Gewichtsfolge g(k).
G(z)\ \bullet \!\!-\!\!\circ \ g(k)
In G(z) und g(k) sind alle Eigenschaften eines linearen, zeitdiskreten dynamischen Systems enthalten.
Ausgangssignal als Funktion des Eingangssignals und der z-Übertragungsfunktion
Y(z)=U(z)\cdot G(z)
{\displaystyle y_{(k)}=\sum _{i=0}^{k}g_{(i)}\cdot u_{(k-i)}}

Im Bildbereich steht die Multiplikation an Stelle der Faltungssumme, wie bei zeitkontinuierlichen Systemen an Stelle des Faltungsintegrals.

Tabelle der Korrespondenzen des Zeitbereichs f(t), des Laplace- und z-Bereichs (Auszüge)

Funktion im
Zeitbereich f(t)
Laplace-Transformierte
im Bildbereich F(s)
Diskrete Laplace-Transformierte
nach der z-Transformation F(z)={\mathcal  Z}\{F(s)\}
δ-Impuls 1 1
Einheits-
sprung 1
{\frac  1s} {\frac  z{z-1}}
t {\frac  1{s^{2}}} {\frac  {T_{A}\cdot z}{(z-1)^{2}}}
t^{2} {\frac  2{s^{3}}} {\frac  {T_{A}^{2}\cdot z\cdot (z+1)}{(z-1)^{3}}}
e^{{-a\cdot t}} {\frac  1{s+a}} {\frac  z{z-e^{{-a\cdot T_{A}}}}}
t\cdot e^{{-a\cdot t}} {\frac  1{(s+a)^{2}}} >{\frac  {e^{{-a\cdot T_{A}}}\cdot T_{A}\cdot z}{(z-e^{{-a\cdot T_{A}}})^{2}}}
1-e^{{-a\cdot t}} {\frac  a{s\cdot (s+a)}} {\frac  {(1-e^{{-a\cdot T_{A}}})\cdot z}{(z-1)\cdot (z-e^{{-a\cdot T_{A}}})}}
1-e^{{-t/T_{1}}} {\frac  1{s\cdot (1+T_{1}\cdot s)}} {\frac  {(1-e^{{-T_{A}/T_{1}}})\cdot z}{(z-1)\cdot (z-e^{{-T_{A}/T_{1}}})}}
\sin \omega_0t {\frac  {\omega _{0}}{s^{2}+\omega _{0}^{2}}} {\frac  {z\cdot \sin(\omega _{0}T_{A})}{z^{2}-2z\cdot \cos(\omega _{0}T_{A})+1}}
\cos \omega_0t {\frac  {s}{s^{2}+\omega _{0}^{2}}} {\frac  {z^{2}-z\cdot \cos(\omega _{0}T_{A})}{z^{2}-2z\cdot \cos(\omega _{0}T_{A})+1}}
1-(1+at)e^{{-at}} {\frac  {a^{2}}{s\cdot (s+a)^{2}}} {\frac  z{z-1}}-{\frac  z{z-c}}-{\frac  {acT_{A}z}{(z-c)^{2}}}\ ;\quad c=e^{{-aT_{A}}}
1+{\frac  {be^{{-at}}-ae^{{-bt}}}{a-b}} {\frac  {ab}{s\cdot (s+a)(s+b)}} {\frac  z{z-1}}+{\frac  {bz}{(a-b)(z-c)}}-{\frac  {az}{(a-b)(z-d)}}
c=e^{{-aT_{A}}}\ \quad d=e^{{-bT_{A}}}
e^{-at}\sin \omega_0t {\frac  {\omega _{0}}{(s+a)^{2}+\omega _{0}^{2}}} {\frac  {cz\sin(\omega _{0}T_{A})}{z^{2}-2cz\cos(\omega _{0}T_{A})+c^{2}}}\ ;\quad c=e^{{-aT_{A}}}
e^{-at}\cos \omega_0t {\frac  {s+a}{(s+a)^{2}+\omega _{0}^{2}}} {\frac  {z^{2}-cz\cos(\omega _{0}T_{A})}{z^{2}-2cz\cos(\omega _{0}T_{A})+c^{2}}}\ ;\quad c=e^{{-aT_{A}}}
a^{{t/T_{A}}} {\frac  1{s-(1/T_{A})\ln a}} {\frac  z{z-a}}
TA = Abtastzeit, a oder b = Zahlenwert der Nullstelle (Pol) im s-Bereich, T1 = Zeitkonstante der s-Übertragungsfunktion.

Tabelle der Korrespondenzen der zeitdiskreten Funktionen f_{{(k)}} zum z-Bereich F(z) (Auszüge)

Name der
Zeitfunktion
Wertefolge f_{{(k)}} z-Transformierte
\delta -Impulsfolge f_{{(k)}}=[0(-T_{A}),1(0T_{A}),0(T_{A}),0(2T_{A}),0(3T_{A}),\dots ] 1
Sprungfolge f_{{(k)}}=[0(-T_{A}),1(0T_{A}),1(T_{A}),1(2T_{A}),1(3T_{A}),\dots ] {\frac  z{z-1}}
Anstiegsfolge u_{{(k)}}=[-u_{1}(-T_{A}),u_{0}(0T_{A}),u_{1}(T_{A}),u_{2}(2T_{A}),u_{3}(3T_{A}),\dots ]
k=[-1,0,1,2,3,\dots ]
{\frac  z{(z-1)^{2}}}
Potenzfolge a^{k} {\frac  z{z-a}}
e-Funktionsfolge e^{{-ak}} {\frac  z{z-e^{{-a}}}}
Sinusfunktionsfolge \sin(ak) {\frac  {z\cdot \sin a}{z^{2}-2z\cdot \cos a+1}}
Kosinusfunktionsfolge \cos(ak) {\frac  {z\cdot (z-\cos a)}{z^{2}-2z\cdot \cos a+1}}
Abklingende Sinus-
funktionsfolge
e^{{-ak}}\cdot \sin(bk) \frac {z \cdot e^{-a}\cdot \sin b} {z^2-2z \cdot e^{-a} \cdot \cos \ b+e^{-2a}}
Abklingende Cosinus-
funktionsfolge
e^{{-ak}}\cdot \cos(bk) \frac {z \cdot (z-e^{-a}\cdot \cos b)} {z^2-2z \cdot e^{-a} \cdot \cos \ b+e^{-2a}}
Linearitätssatz a_{1}f_{1}(k)+a_{2}f_{2}(k)+\dots +a_{n}f_{n}(k) a_{1}F_{1}(z)+a_{2}F_{2}(z)+\dots +a_{n}F_{n}(z)
Rechts-Verschiebesatz f(k-d) z^{{-d}}\cdot F(z)
Faltungssumme \sum _{{i=0}}^{k}f_{1}(i)\cdot f_{2}(k-i) F_{1}(z)\cdot F_{2}(z)
Anfangswertsatz \lim _{{k\to 0}}f(k) \lim _{{z\to \infty }}F(z)
Endwertsatz \lim _{{k\to \infty }}f(k) \lim _{{z\to 1}}\,(z-1)\cdot F(z)

Z-Transformation einer Wertefolge (Impulsfolge)

Die Laplace-Transformation bezieht sich auf die Ableitungen einer gewöhnlichen Differentialgleichung und ersetzt diese nach dem Laplace-Differentiationssatz durch die komplexe Variable s=\sigma +j\omega . Ein Exponent von s kennzeichnet die Ordnung der Ableitung.

Die z-Transformation transformiert eine Impulsfolgefunktion f^{*}(t) oder eine Zahlenfolge f_{{(k)}} in eine Funktion F(z) mit der z-Variable z=e^{{T_{A}\cdot s}}.

Da die Berechnungen mit Impulsfunktionen oder Folgen aufwendig sind, ist es sinnvoll, diese durch einfachere Berechnungen im z-Bildbereich auszuführen. Die z-Transformation von Impulsfolgen kann als diskrete Laplace-Transformation aufgefasst werden.

Die zu diskreten Zeitpunkten abgetasteten kontinuierlichen Eingangssignale entsprechen mit der Eingangsgröße modulierten (gewichteten) \delta -Impulsfolgen, die mit der z-Transformation berechnet werden.

Beim Übergang von kontinuierlichen Systemen f(t) zu Abtastsystemen f_{{(k)}} mit der Abtastfolge k=(0,1,2,3,\dots ) gehen lineare zeitinvariante Differenzialgleichungen in zeitinvariante Differenzengleichungen der Abtastzeit TA beziehungsweise der Abtastfrequenz f_{{TA}}=1/T_{A} der Funktion f(kT_{A}) über.

Kontinuierliche Signale z.B. die Regelabweichung e(t) eines Regelkreises mit einer analog arbeitenden Regelstrecke werden zu gleichen Zeitabständen T_A abgetastet. Ein Abtaster mit A/D-Wandler erzeugt aus einem zeitkontinuierlichen Regler-Eingangssignal e(t) ein zeitdiskretes Signal e_{{(k)}}. Es können sowohl abgetastete Regler-Eingangssignale e_{{(k)}} als auch Differenzengleichungen f(k), die im diskreten Zeitbereich den Regelalgorithmus eines Reglers beschreiben, als z-Übertragungsfunktionen in den z-Bereich transformiert und als algebraische Gleichungen behandelt werden.

Wird eine inverse z-Transformation der z-Übertragungsfunktionen durchgeführt, entsteht die Lösung der zeitdiskreten Differenzengleichung f_{{(k)}}={\mathcal  Z}^{{-1}}\{F(z)\} im f(k)-Bereich. Mit Hilfe verschiedener Verfahren der Rücktransformation vom z-Bereich in den k-Bereich ergeben sich dann als Lösung die Differenzengleichungen des Regelalgorithmus für den diskreten Bereich f(k).

Für die häufig vorkommenden Transformationen des k- und z-Bereiches stehen in vielen Fachbüchern der Regelungstechnik Transformationstafeln zur Verfügung.

Eine einzelne Abtastung eines kontinuierlichen Signals u(t) an einer beliebigen Stelle des Signalverlaufs z.B. zum Zeitpunkt T_{1} wird als Modulation von u(t) mit einem Dirac-Impuls mit \delta \cdot (t-T_{1}) beschrieben.

Damit ergibt sich für die Multiplikation des Eingangssignals u(t) mit dem Dirac-Impuls:

u_{{T_{1}}}(t)=u(t)\cdot \delta (t-T1)=u(T_{1})\cdot \delta (t-T_{1})

Wird nun periodisch mit zu den Zeitpunkten kT_{A} das Signal u(t) abgetastet, kann das abgetastete Signal u(kT_{A})=u_{k} als Multiplikation einer \delta -Impulsfolge mit u(t) betrachtet werden.

Die Impulsfolgefunktion bezieht sich auf die Summengrenzen von k=-\infty zu k=\infty

u(t)=\sum _{{k=-\infty }}^{\infty }u(kT_{A})\cdot \delta (t-kT_{A})

Mit dem Übergang der Summengrenzen von k=0 zu k=\infty wird die Impulsfolgefunktion:

{\displaystyle u(t)=\left.\sum _{k=0}^{\infty }u_{k}\cdot \delta (t-kT_{A})\right|_{{\text{mit }}u_{k}=0{\text{ für }}k<0}}

Die Impulsfolgefunktion wird der Laplace-Transformation unterzogen:

U(s)={\mathcal  L}{\Bigl \{}\sum _{{k=0}}^{\infty }u_{k}\cdot \delta (t-kT_{A}){\Bigr \}}=\int \limits _{{0}}^{\infty }\sum _{{k=0}}^{\infty }u_{k}\cdot \delta (t-kT_{A})\cdot e^{{-st}}dt=\sum _{{k=0}}^{\infty }u_{k}\cdot e^{{-k\cdot s\cdot T_{A}}}

Mit der z-Variablen z=e^{{s\cdot T_{A}}} und z^{{-k}}=e^{{-ksT_{A}}} wird die Gleichung vereinfacht zu der z-Transformierten des kontinuierlichen Signals u(t).

Mit z=e^{{s\cdot T_{A}}} wird die s-Ebene auf die z-Ebene abgebildet.

Damit lautet die z-Transformation eines abgetasteten kontinuierlichen Signals u(t) mit der Folge k=(0,1,2,3,\dots ) und der Wertefolge u(kT_{A}):=u_{{k}}=[u_{{(0)}};u_{{(T_{A})}};u_{{(2T_{A})}};u_{{(3T_{A})}};\cdots ]:

\mathcal Z \{u_{(k)}\}=U(z)=\sum_{k=0} ^\infty u_k \cdot z^{-k}

oder allgemein als Funktion F(z):

F(z)=\sum_{k=0} ^\infty f_k \cdot z^{-k}

Z-Übertragungsfunktion (Impulsübertragungsfunktion) von zeitdiskreten Elementen

Übersicht:

Wendet man auf die einzelnen Terme der Differenzengleichung den Linearitäts- und den rechts-Verschiebungssatz an und bringt alle Terme mit y_{(k)} auf die linke Seite und alle Terme mit u_{(k)} auf die rechte Seite der Differenzengleichung, so lässt sich das Verhältnis G(z)=Y(z)/U(z) mit den verbleibenden Elementen als gebrochen-rationale Funktion darstellen.
Die z-Übertragungsfunktion lautet:
G(z)={\frac  {Y(z)}{U(z)}}={\frac  {{\text{Zählerpolynom}}\ (z)}{{\text{Nennerpolynom}}\ (z)}}

Mit der so errechneten Differenzengleichung des Übertragungssystems ist man nun in der Lage, für eine gegebene Eingangserregung u_{(k)} des Systems die Ausgangsfolgen des Systems y_{(k)} zu berechnen. Dies ist auf verschiedenen Arten möglich.

Entstehung der z-Übertragungsfunktion:

Ein zeitdiskretes dynamisches System wird durch eine rekursive Differenzengleichung beschrieben.

Wie bei zeitkontinuierlichen Systemen g(t) und der Übertragungsfunktion G(s) besteht eine vergleichbare Beziehung bei den zeitdiskreten Systemen zwischen der Gewichtsfolge g(k) und z-Übertragungsfunktion G(z). Die Übertragungsfunktion G(z) ist die z-Transformierte der Gewichtsfolge g(k).

Beziehung der Gewichtsfolge g(k) zur z-Übertragungsfunktion G(z).
G(z)\ \bullet \!\!-\!\!\circ \ g(k)

In G(z) und g(k) sind alle Eigenschaften eines dynamischen Systems enthalten. Zur Systemberechnung wird der Verlauf des Eingangssignals U(z) und das Systemverhalten G(z) oder g(k) benötigt.

Für die Ermittlung der z-Übertragungsfunktion ist die Gewichtsfunktion g(t) die Systemantwort (Impulsantwort) auf die Vorgabe eines Eingangs-DIRAC-Impulses \delta (t).

Bei der Laplace-Transformation eines Systems lautet die Übertragungsfunktion:

{\displaystyle G(s)={\frac {Y(s)}{U(s)}}={\frac {{\text{Zählerpolynom}}\ (s)}{{\text{Nennerpolynom}}\ (s)}}\ |\qquad {\text{mit }}U(s)=\delta (t)}

Aus der Gewichtsfunktion wird die Impulsfolgefunktion g^{*}(t) oder die Gewichtsfolge g(kTA) gebildet und in den z-Bereich transformiert. Die z-Übertragungsfunktion von zeitkontinuierlichen Systemen ergibt sich mit folgender Vorschrift:

G(z)={\mathcal  Z}\{g(t)|_{{t=kT_{A}}}\}={\mathcal  Z}\{{\mathcal  L}^{{-1}}\{G(s)\}|_{{t=kT_{A}}}\}

Die Berechnung der Impulsübertragung von linearen zeitdiskreten Systemen vereinfacht sich, wenn aus der Differenzengleichung mit der z-Transformation die z-Übertragungsfunktion bestimmt wird. So können z.B. die Regelalgorithmen eines digitalen Reglers mit Differenzengleichungen formuliert werden.

Ist die z-Übertragungsfunktion eines Systems G(z) gegeben, lässt sich die zugehörige z-Differenzengleichung bestimmen:

Es wird eine lineare Differenzengleichung eines Regelalgorithmus gegeben:

Die Eingangssignale der Differenzengleichung sind wieder u_{(k)} und die Ausgangssignale y_{(k)}.

a_{n}\cdot y_{{(k+n)}}+a_{{n-1}}\cdot y_{{(k+n-1)}}+\cdots +a_{1}\cdot y_{{(k+1)}}+a_{0}\cdot y_{{(k)}}=b_{m}\cdot u_{{(k+m)}}+b_{{m-1}}\cdot u_{{(k+m-1)}}+\cdots +b_{1}\cdot u_{{(k+1)}}+b_{0}\cdot u_{{(k)}}|\qquad n\geq m

Unter der Voraussetzung, dass die Anfangswerte des zeitdiskreten Systems gleich Null sind, wird die Differenzengleichung transformiert. Mit den Rechenregeln des rechts-Verschiebungssatzes ergibt sich die Transformationsvorschrift:

y_{{(k)}}\to y(z),\qquad y_{{(k+i)}}\to z^{i}\cdot y(z),\qquad y_{{(k-i)}}\to z^{{-i}}\cdot y(z),
u_{{(k)}}\to u(z),\qquad u_{{(k+i)}}\to z^{i}\cdot u(z),\qquad u_{{(k-i)}}\to z^{{-i}}\cdot u(z)

Die z-transformierte Differenzengleichung ergibt sich zu:

a_{n}\cdot z^{n}\cdot Y(z)+a_{{n-1}}\cdot z^{{n-1}}\cdot Y(z)+\cdots +a_{1}\cdot z\cdot Y(z)+a_{0}\cdot Y(z)=b_{m}\cdot z^{m}\cdot U(z)+b_{{m-1}}\cdot z^{{m-1}}\cdot U(z)+\cdots +b_{1}\cdot z\cdot U(z)+b_{0}\cdot U(z)

Durch Ausklammern des Quotienten Y(z) / U(z) kann die Impulsübertragungsfunktion G(z) gebildet werden.

Damit lautet die z-Übertragungsfunktion:

G(z)={\frac  {Y(z)}{U(z)}}={\frac  {b_{m}\cdot z^{m}+b_{{m-1}}\cdot z^{{m-1}}+\cdots +b_{1}\cdot z+b_{0}}{a_{n}\cdot z^{n}+a_{{n-1}}\cdot z^{{n-1}}+\cdots +a_{1}\cdot z+a_{0}}}={\frac  {{\text{Zähler}}(z)}{{\text{Nenner}}(z)}}

Bei der Berechnung von Differenzengleichungen mit dem Digitalrechner ist die von n-Schritten nach rechts verschobene Gleichung besser geeignet.

 a_n \cdot y_{(k)} + a_{n-1} \cdot y_{(k-1)} + \cdots + a_1 \cdot y_{(k-n+1)} + a_0 \cdot y_{(k-n)} = b_m \cdot u_{(k-n+m)}+ b_{m-1} \cdot u_{(k-n+m-1)}+ \cdots + b_1 \cdot u_{(k-n+1)}+b_0 \cdot u_{(k-n)}

Damit lautet die z-Übertragungsfunktion mit negativen Exponenten von z:

G(z)={\frac  {Y(z)}{U(z)}}={\frac  {b_{m}\cdot z^{{m-n}}+b_{{m-1}}\cdot z^{{m-1-n}}+\cdots +b_{1}\cdot z^{{-n+1}}+b_{0}\cdot z^{{-n}}}{a_{n}+a_{{n-1}}\cdot z^{{-1}}+\cdots +a_{1}\cdot z^{{-n+1}}+a_{0}\cdot z^{{-n}}}}

Beispiel: Für eine gegebenen Differenzengleichung 3. Ordnung lautet die z-Übertragungsfunktion für Ordnungen m = 2 und n = 3:

G(z)={\frac  {Y(z)}{U(z)}}={\frac  {b_{2}\cdot z^{2}+b_{1}\cdot z+b_{0}}{a_{3}\cdot z^{3}+a_{2}\cdot z^{2}+a_{1}\cdot z+a_{0}}}

Mit der rechts-Verschiebung der zugehörigen Differenzengleichung um n = 3 Abtastschritte T_A lautet die z-Übertragungsfunktion:

G(z)={\frac  {Y(z)}{U(z)}}={\frac  {b_{2}\cdot z^{{-1}}+b_{1}\cdot z^{{-2}}+b_{0}\cdot z^{{-3}}}{a_{3}+a_{2}\cdot z^{{-1}}+a_{1}\cdot z^{{-2}}+a_{0}\cdot z^{{-3}}}}

Es fällt auf, dass die z-Transformation einer Differenzengleichung eine große Ähnlichkeit mit der Polynomform einer s-Übertragungsfunktion für kontinuierliche Systeme hat. Ebenso fällt auf, dass die z-Übertragungsfunktion negative Potenzen hat. Allgemein werden z-Übertragungsfunktionen mit negativen Potenzen von z dargestellt.

Jede z-Übertragungsfunktion mit negativen Potenzen von z kann in eine Form mit positiven Potenzen von z gebracht werden, wenn Zähler und Nenner der z-Übertragungsfunktion mit der höchsten vorkommenden inversen z-Potenz multipliziert werden. Für die z-Übertragungsfunktion der Form mit positiven Exponenten muss gelten: Nennerpolynom > Zählerpolynom. Ist diese Bedingung erfüllt, ist G(z) kausal.

Es können sowohl abgetastete Eingangssignale u(kTA) als auch Differenzengleichungen f(kTA), die im diskreten Zeitbereich das Verhalten eines Systems (z.B. den Regelalgorithmus eines Reglers) beschreiben, als z-Übertragungsfunktionen in den z-Bereich transformiert und als algebraische Gleichungen behandelt werden.

Wird eine inverse z-Transformation der z-Übertragungsfunktionen durchgeführt, entsteht die Lösung der zeitdiskreten Differenzengleichung f_{{(k)}}={\mathcal  Z}^{{-1}}\{F(z)\} im f(kT_{A})-Bereich. Mit Hilfe verschiedener Verfahren der Rücktransformation vom z-Bereich in den k-Bereich ergeben sich dann als Lösung die Differenzengleichungen des Regelalgorithmus für den diskreten Bereich f(kTA).

Die typische Anwendung der z-Transformation eines digitalen Systems, eines digitalen Reglers oder eines digitalen Filters, für den Regelalgorithmus lautet wie folgt:

Die Analyse und die Synthese diskreter Signale und Systeme lässt sich mit der z-Transformation erleichtern, setzt aber auch umfangreiches mathematisches Spezialwissen voraus.

Für die häufig vorkommenden Transformationen des k- und z-Bereich stehen in vielen Fachbüchern der Regelungstechnik Transformationstafeln zur Verfügung.

Z-Übertragungsfunktion einer Totzeit

Eine Totzeit Tt hat im Zeitbereich das Verhalten \displaystyle y(t)={u}(t-T_{t}).

Im s-Bereich lautet die Übertragungsfunktion eines Totzeitsystems:

G(s)={\frac  {y(s)}{u(s)}}=e^{{-s\cdot T_{t}}}

Es handelt sich hier bei der Verbindung einer s-Übertragungsfunktion als gebrochen-rationalen Funktion mit einem Totzeitglied um eine transzendente Funktion, die als Anhang einer Übertragungsfunktion - beispielsweise einer Regelstrecke - multiplikativ zugeordnet wird, mit der Einschränkung, dass keine algebraische Behandlung erlaubt ist.

Im z-Bereich entspricht eine Totzeit T_{t}=d\cdot T_{A} einer Rückwärtsverschiebung (Verschiebung nach rechts) der Abtastfolgen um d Abtastschritte. Die z-Übertragungsfunktion des Totzeitgliedes in Verbindung mit weiteren z-Übertragungssystemen bleibt gebrochen-rational (Bruch mit Zähler- und Nennerpolynom).

Für die z-Transformation der Abtastfolge für das Signal u_{{(k-d)}} gilt:

u_{{(k-d)}}\ \circ \!\!-\!\!\bullet \ z^{{-d}}\cdot U(z)

Anders als bei der Laplace-Transformation bedeutet ein Totzeitglied in Verbindung mit einer z-Übertragungsfunktion keine Einschränkung der bestehenden Form der gebrochen-rationalen Funktion. Eine z-Übertragungsfunktion mit einer mit d Abtastschritten T_A definierten Totzeit T_{t}=d\cdot T_{A} wird durch hinzufügen von eines Terms z^{{-d}} berücksichtigt. Es können beliebige algebraische Operationen durchgeführt werden.

Beispiel einer z-Übertragungsfunktion mit Totzeit:

{\displaystyle G(z)={\frac {Y(z)}{U(z)}}={\frac {b_{m}\cdot z^{m}+b_{m-1}\cdot z^{m-1}+\cdots +b_{1}\cdot z+b_{0}}{a_{n}\cdot z^{n}+a_{n-1}\cdot z^{n-1}+\cdots +a_{1}\cdot z+a_{0}}}\cdot z^{-d}}

Z-Übertragungsfunktion von zeitkontinuierlichen Elementen G(s) mit Abtaster und Halteglied

Definition der Signale: Allgemein gilt für den Zeitbereich f(t), zeitdiskreten Bereich f(k) und Bildbereich f(s), f(z) die Eingangsgröße u und die Ausgangsgröße y.

Abtastsysteme erlauben eine Umwandlung kontinuierlicher Signale in zeitdiskrete Abtastfolgen (Impulsfolgen). Halteglieder erlauben eine Umwandlung zeitdiskreter Ausgangsfolgen (Impulsfolgen) in kontinuierliche gestufte Signale im Zeitbereich. Ein Abtaster setzt ein kontinuierliches Signal u(t) in Verbindung mit einem nachgeschalteten A/D-Wandler in eine Zahlenfolge u(k) mit digitalen Werten um.

Ein Halteglied nullter Ordnung setzt eine Zahlenfolge u_{(k)} mit einem vorgeschalteten D/A-Wandler in ein gestuftes kontinuierliches Signal u(t) um. Bei der Reihenschaltung von einem Abtastsystem (gewichteter δ-Abtaster) und einem Halteglied handelt es sich um die Umwandlung einer aus einem kontinuierlichen Eingangssignal modulierten Impulsfolge in eine gestufte Treppenfunktion f(t), die in ein kontinuierliches dynamisches System eingeleitet werden kann.

Liegt ein kontinuierliches dynamisches System vor, erlaubt die Abtastung eines kontinuierlichen Signals, z.B. eine Regelabweichung, die Verarbeitung einer zeitdiskreten Eingangsfolge in einem Mikrocomputer als Regelalgorithmus zu einer zeitdiskreten Ausgangsfolge. Die über ein Halteglied geleitete Ausgangsfolge erzeugt ein quasi stetiges gestuftes Ausgangssignal zu einer analogen zeitabhängigen Stellgröße.

Je nach Anwendung des kontinuierlichen Systems G(s), beispielsweise in der Regelungstechnik, werden unterschiedlich Funktionsblöcke zusammengefasst, die auch mehrere Abtaster und Halteglieder enthalten können.

Um ein kontinuierliches System G(s) in eine z-Übertragungsfunktion (Impulsübertragungsfunktion) zu definieren, muss es im Eingang durch eine Abtastfolge e_{{(k)}} als Eingangsfolge und am Ausgang des Systems zeitsynchron durch einen Abtaster die Ausgangsfolge u_{(k)} zur Verfügung stellen.

Berechnung der Varianten kontinuierlicher Systeme G(s) mit Abtaster und Halteglied.

Je nachdem wie die Reihenfolge der Funktionsböcke eines Haltegliedes G_{H}(s), eines Systems G(s) und eines Abtasters festgelegt sind, können diese Funktionsblöcke als ein Teil einer offenen Regelstrecke betrachtet werden. Die wären z.B.:

Das obere Grafikbild zeigt die Kombination der Funktionsblöcke Halteglied, kontinuierliches System G(s) und Abtaster.

Damit ein kontinuierliches System G(s) in eine (Impuls-)Übertragungsfunktion G(z) definiert werden kann, muss das System durch eine Impulsfolge (Wertefolge) gespeist und am Ausgang ein Abtaster eingesetzt sein, der synchron zur Eingangsfolge eine Ausgangsfolge realisiert. Das Halteglied vor dem dynamischen System wandelt die Impulsfolge in ein quasikontinuierliches Signal um. Das kontinuierliche System-Ausgangssignal wird über einen Abtaster als Ausgangsfolge definiert.

Die s-Übertragungsfunktion des Haltegliedes G_{H}(s) nullter Ordnung lautet:

G_{H}(s)={\frac  {{\bar  U}(s)}{{\hat  U}(s)}}={\frac  {1-e^{{T_{A}\cdot s}}}s}

Die z-Übertragungsfunktion (Impulsübertragungsfunktion) des Haltegliedes G_{H}(s) und der Regelstrecke G_{S}(s) wird durch die z-Transformation von {\mathcal  Z}\{G_{H}(s)\cdot G_{S}(s)\} berechnet.

Die z-Übertragungsfunktion von Halteglied und Regelstrecke lautet:

G_{{HS}}(z)=\left.{\mathcal  Z}\{G_{H}(s)\cdot G_{S}(s)\right|_{{t=kT_{A}}}\}={\mathcal  Z}\{G_{H}(s)\cdot G_{S}(s)\}/DD>

In diese Gleichung wird G_{H}(s) eingesetzt:

G_{{HS}}(z)={\mathcal  Z}{\biggl \{}{\frac  {1-e^{{-T_{A}s}}}{s}}\cdot G_{S}(s){\biggr \}}={\mathcal  Z}{\biggl \{}{\frac  {G_{S}(s)}s}{\biggr \}}-{\mathcal  Z}{\biggl \{}{\frac  {e^{{-T_{A}s}}\cdot G_{S}(s)}s}{\biggr \}}

Eine Multiplikation mit e^{{-sT_{A}}} entspricht im z-Bereich eine Rechtsverschiebung um einen Abtastzyklus T_A und damit einer Multiplikation mit z^{{-1}}.

G_{{HS}}(z)={\mathcal  Z}{\biggl \{}{\frac  {1-e^{{-T_{A}s}}}{s}}\cdot G_{S}(s){\biggr \}}=(1-z^{{-1}})\cdot {\mathcal  Z}{\biggl \{}{\frac  {G_{S}(s)}s}{\biggr \}}

Die z-Übertragungsfunktion der Reihenschaltung eines Haltegliedes nullter Ordnung mit einem kontinuierlichen System ergibt sich aus der Laplace-Transformierten der Sprungantwort G_{S}(s), multipliziert mit {\frac  {z-1}z}. (Siehe Tabelle {\displaystyle {\frac {1}{s}}\ \circ \!\!-\!\!\bullet \ {\frac {z}{z-1}}}):

G_{{HS}}(z)={\frac  {Y{(z)}}{U{(z)}}}={\frac  {z-1}z}\cdot {\mathcal  Z}{\biggl \{}{\frac  {G_{S}(s)}s}{\biggr \}}

Analog wie bei der Regelstrecke G_{S}(s) kann für die im 2. Teil des Bildes dargestellte Anordnung der Funktionsblöcke Abtaster, Regelalgorithmus des Reglers und Halteglied die gleiche Form der Gleichung der z-Transformation benutzt werden. Ein- und Ausgangsgröße und die Systemgröße des Reglers G_{R}(s) werden getauscht:

G_{{HR}}(z)={\frac  {U{(z)}}{E{(z)}}}={\frac  {z-1}z}\cdot {\mathcal  Z}{\biggl \{}{\frac  {G_{R}(s)}s}{\biggr \}}

Um die z-Transformation des kontinuierlichen Systems G(s) bilden zu können, werden die Korrespondenztabellen angewendet. Steht für die gegebene s-Übertragungsfunktion G(s) kein Eintrag für die z-Transformation zur Verfügung, muss eine Partialbruchzerlegung von G(s)/s vorgenommen werden.

Beispiel: z-Transformation eines IT_{1}-Gliedes.

G(s)={\frac  1{s\cdot (s+1)}}

Gesucht z-Übertragungsfunktion G_{{HS}}(z) des kombinierten Systems:

G_{{HS}}(z)={\frac  {Y{(z)}}{U{(z)}}}={\frac  {z-1}z}\cdot {\mathcal  Z}{\biggl \{}{\frac  1{s^{2}\cdot (s+1)}}{\biggr \}}

Partialbruchzerlegung von

{\frac  1{s^{2}\cdot (1+s)}}={\frac  1{s^{2}}}-{\frac  1s}+{\frac  1{(s+1)}}

z-Transformation der Partialbrüche in den z-Bereich:

G_{{HS}}(z)={\frac  {Y{(z)}}{U{(z)}}}={\frac  {z-1}z}\cdot {\mathcal  Z}{\biggl \{}{\frac  1{s^{2}}}-{\frac  1s}+{\frac  1{(s+1)}}{\biggr \}}={\frac  {z-1}z}\cdot {\biggl \{}{\frac  {T_{A}\cdot z}{(z-1)^{2}}}-{\frac  z{z-1}}+{\frac  z{z-e^{{-T_{A}}}}}{\biggr \}}

Durch Ausmultiplizieren der Brüche ergibt sich die z-Transformation:

G_{{HS}}(z)={\frac  {Y(z)}{U(z)}}={\frac  {z-1}z}-1+{\frac  {z-1}{z-e^{{-T_{A}}}}}

Z-Übertragungsfunktionen der Standardregler

Es gelten hier allgemein die formalen Begriffe der Ein- und Ausgangsgrößen der Systemtheorie: Eingangssignal = U(z), Ausgangssignal Y(z).
Für die Belange der Regler gelten die Eingangs- und Ausgangsgrößen des Reglers zu einem Regelkreis lauten: E(z) und U(z).

Die Regelalgorithmen mit den allgemeinen Koeffizienten lauten nach der rekursiven Rechenvorschrift:

y_{{(k)}}=a_{1}\cdot y_{{(k-1)}}+b_{0}\cdot u_{{(k)}}+b_{1}\cdot u_{{(k-1)}}+b_{2}\cdot u_{{(k-2)}}

Die z-Transformation der Differenzengleichung liefert die z-Übertragungsfunktion GR(z) des Reglers:

G_{R}(z)={\frac  {Y(z)}{U(z)}}={\frac  {b_{0}+b_{1}\cdot z^{{-1}}+b_{2}\cdot z^{{-2}}}{1-a_{1}\cdot z^{{-1}}}}={\frac  {b_{0}\cdot z^{2}+b_{1}\cdot z+b_{2}}{z^{2}-a_{1}\cdot z}}

Tabelle der z-Übertragungsfunktionen der Standardregler (Typ II, Obersumme)

(ai = Koeffizienten des Nennerpolynoms, bi = Koeffizienten der Zählerpolynoms, TI = TN = 1 / KI = Zeitkonstante = Nachstellzeit, TV = Vorhaltezeit, TA = Abtastzeit)

Reglerart z-Übertragungsfunktion GR(z) der Standardregler
P-Regler G_{R}(z)=b_{0},\ b_{0}=K_{R}
I-Regler G_{R}(z)={\frac  {b_{0}\cdot z}{z-a_{1}}}={\frac  {b_{0}}{1-a_{1}\cdot z^{{-1}}}}

a_{1}=1;\ b_{0}=K_{I}\cdot T_{A}
PD-Regler G_{R}(z)={\frac  {b_{0}\cdot z+b_{1}}z}=b_{0}+b_{1}\cdot z^{{-1}}

b_{0}=K_{R}\cdot [1+{\frac  {T_{V}}{T_{A}}}];\ b_{1}=-K_{R}\cdot {\frac  {T_{V}}{T_{A}}}
PI-Regler G_{R}(z)={\frac  {b_{0}\cdot z+b_{1}}{z-a_{1}}}={\frac  {b_{0}+b_{1}\cdot z^{{-1}}}{1-a_{1}\cdot z^{{-1}}}}

a_{1}=1;\ b_{0}=K_{R}\cdot \left[1+{\frac  {T_{A}}{T_{N}}}\right];\ b_{1}=-K_{R}
PID-Regler G_{R}(z)={\frac  {b_{0}\cdot z^{2}+b_{1}\cdot z+b_{2}}{z^{2}-a_{1}\cdot z}}={\frac  {b_{0}+b_{1}\cdot z^{{-1}}+b_{2}\cdot z^{{-2}}}{1-a_{1}\cdot z^{{-1}}}}

a_{1}=1;\ b_{0}=K_{R}\cdot \left[1+{\frac  {T_{A}}{T_{N}}}+{\frac  {T_{V}}{T_{A}}}\right];\ b_{1}=-K_{R}\cdot \left[1+2\cdot {\frac  {T_{V}}{T_{A}}}\right];\ b_{2}=K_{R}\cdot {\frac  {T_{V}}{T_{A}}}

Umwandlung einer z-Übertragungsfunktion in eine Differenzengleichung

Eine gegebene z-Übertragungsfunktion lässt sich in eine Differenzengleichung umwandeln. Dazu sind folgende Schritte erforderlich.

Die Übertragungsfunktion wurde durch den Koeffizienten a_{n} dividiert, um z^{n} freistellen zu können.

Diese Form der Übertragungsfunktion wurde entsprechend der dargestellten Koeffizienten wie folgt neu geordnet.

G(z)={\frac  {Y(z)}{U(z)}}={\frac  {b_{m}\cdot z^{m}+b_{{m-1}}\cdot z^{{m-1}}+b_{{m-2}}\cdot z^{{m-2}}\cdots +b_{1}\cdot z+b_{0}}{z^{n}+a_{{n-1}}\cdot z^{{n-1}}+a_{{n-2}}\cdot z^{{n-2}}+\cdots +a_{1}\cdot z+a_{0}}}={\frac  {{\text{Zähler}}(z)}{{\text{Nenner}}(z)}}
{b_m \cdot z^m \cdot U(z) + b_{m-1} \cdot z^{m-1}\cdot U(z) + b_{m-2}\cdot z^{m-2}\cdot U(z)+ \cdots + b_1 \cdot z \cdot U(z)+ b_0 \cdot U(z)=
z^n\cdot Y(z) + a_{n-1} \cdot z^{n-1}\cdot Y(z) + a_{n-2} \cdot z^{n-2}\cdot Y(z) + \cdots + a_1 \cdot z\cdot Y(z) + a_0\cdot Y(z)}
{b_m \cdot z^{-(n-m)} \cdot U(z) + b_{m-1} \cdot z^{-(n-m+1)}\cdot U(z) + b_{m-2}\cdot z^{-(n-m+2)}\cdot U(z)+ \cdots + b_0 \cdot z^{-n} \cdot U(z)= Y(z) + a_1 \cdot z^{-(n-1)}\cdot Y(z) + a_2 \cdot z^{-(n-2)}\cdot Y(z) + a_0 \cdot z^{-n}\cdot Y(z)}
{b_m \cdot u[k-n+m]+b_{m-1} \cdot u[k-n+m-1]+\dots +b_1 \cdot u[k-n+1]+b_0 \cdot u[k-n]=y[k]+a_{n-1} \cdot y[k-1]+ \dots +a_1 \cdot y[k-n+1]+a_0 \cdot y[k-n]}
y_{k}= -a_{n-1} \cdot y_{(k-1)}- \dots - a_1 \cdot y_{(k-n+1)} -a_0 \cdot y_{(k-n)} + b_m \cdot u_{(k-n+m)} + \dots + b_1 \cdot u_{(k-n+1)}+ b_0 \cdot u_{(k-n)}

Rücktransformation in den diskreten Zeitbereich f(kTA)

Für die Programmierung des Mikrocomputers des Digitalreglers werden Differenzengleichungen benötigt, die durch Überführung der z-Übertragungsfunktion des z-Bildbereiches in den diskreten Zeitbereich f(k) gewonnen wird.

Die inverse z-Transformation liefert für F(z) wieder die Werte der Zahlenfolge:

{\mathcal  Z}^{{-1}}\{F(z)\}=f(kT_{A}) für k=(0;1;2;3;\cdots ).

Dazu sind drei Verfahren gegeben:

Durch die Polynomdivision wird die z-Transformierte F(z) in eine konvergierende Potenzreihe nach z^{{-i}} entwickelt. Dividiert man das Zählerpolynom von F(z) durch das Nennerpolynom nach den Regeln der Polynomdivision, so ergibt sich die gewünschte Potenzreihe. Die Vorfaktoren der Potenzen von z sind die gesuchten Werte der Zahlenfolge f(kT_{A}).
Die Rücktransformation von F(z) erfolgt durch die Partialbruchzerlegung mit Hilfe der zugehörigen Korrespondenztabelle. Die Partialbrüche sollen die Form {\frac  z{z-z_{i}}} haben. Die inverse z-Transformation ergibt f(kT_{A}) als Summe der rücktransformierten Partialbrüche.
Die inverse z-Transformation ist über das komplexe Umlaufintegral bestimmt. Die Lösung des Integrals lässt sich auf die Bestimmung der Residuen  \operatorname{Res} \{z_i \} zurückführen.
f(kT_A) = \sum_i \operatorname{Res} \{F(z) \cdot z^{k-1} \} \quad |_{Z=Z_i}

Berechnungsbeispiele zur z-Transformation

Beispiel 1: z-Transformierte der Sprungfunktion der Eingangsgröße u_{{\sigma }}(t):

Die z-Transformierte der normierten Sprungfunktion u_{{\sigma }}(t) = 1 für t > 0 ist zu ermitteln.

Lösung:

Die kontinuierliche Eingangsgröße u(t) wird als Sprungfunktion (Einheitssprung) zum Zeitpunkt t = 0 mit dem Wert 1 abgetastet.
Durch Abtastung entsteht eine Wertefolge mit konstanten Amplituden = 1.

u(kT_{A})=u_{{(k)}}=0(-1),1(0),1(T_{A}),1(2T_{A}),1(3T_{A})\dots

oder die Impulsfolgefunktion:

u^{*}(t)=\sum _{{k=0}}^{\infty }u_{{(k)}}\cdot \delta (t-kT_{A})

Die z-Transformation lautet damit:

{\mathcal  Z}\{u_{k}\}=u(z)=\sum _{{k=0}}^{\infty }u_{k}\cdot z^{{-k}}=\sum _{{k=0}}^{\infty }z^{{-k}}\quad \quad {\text{mit}}\ u_{k}=1

Der Grenzwert für z^{{-n}} für n\to \infty und |z|>1 (Konvergenz) ist Null.
Die z-Transformierte der Einheitssprungfunktion ist damit:

{\mathcal  Z}\{u_{k}\}=u(z)=\sum _{{k=0}}^{\infty }u_{k}\cdot z^{{-k}}={\frac  z{z-1}}

Beispiel 2: z-Transformierte des Produktes einer Exponentialfunktion f(t) mit dem Einheitssprung u(t)=1:

Gegeben: Exponentialfunktion y(t)=e^{{-t/T_{1}}}

Einheitssprung u_{{\sigma }} = 1 für t > 0.

Gesucht: Ergebnis der beiden Beziehungen als Produkt im z-Bereich:

Die Ausgangswertefolge y(kT_{A}) lautet:

y(kT_{A})=y_{{(k)}}=e^{{-{\frac  {kT_{A}}{T_{1}}}}}

Daraus ergibt sich die z-Transformierte:

y(z)=\sum _{{k=0}}^{\infty }y_{k}\cdot z^{{-k}}=\sum _{{k=0}}^{\infty }e^{{-{\frac  {kT_{A}}{T_{1}}}}}\cdot z^{{-k}}

Substituiert man

z^{'}=z\cdot e^{{T_{A}/T_{1}}},\quad (z^{'})^{{-k}}=z^{{-k}}\cdot e^{{-kT_{A}/T_{1}}},

so lautet die z-Transformierte

y(z)={\mathcal  Z}\{y(kT_{A})\}={\mathcal  Z}\{e^{{-kT_{A}/T_{1}}}\}={\frac  z{z-e^{{-T_{A}/T_{1}}}}}

Beispiel 3: Ermittlung der z-Übertragungsfunktion aus einem Verzögerungsglied 1. Ordnung

Gegeben: Übertragungsfunktion G(s) des PT1-Gliedes mit der Verstärkung K und der Zeitkonstante T1:

{\displaystyle G(s)={\frac {Y(s)}{U(s)}}={\frac {K}{1+s\cdot {T_{1}}}}={\frac {K\cdot a}{s+a}};\qquad a=1/T_{1}=} Polstelle des Nennerpolynoms.
Anmerkung: Die Zeitkonstantendarstellung der Übertragungsfunktion ist der Pol- Nullstellendarstellung mathematisch identisch.

Gesucht: z-Übertragungsfunktion G(z) aus G(s).

Dazu wird zuerst die Gewichtsfunktion g(t) des PT1-Gliedes ermittelt.
Aus einer Laplace-Korrespondenztabelle findet man für das Zeitverhalten f(t) der Übertragungsfunktion G(s) für die Gewichtsfunktion (Impulsantwort):

{\displaystyle g(t)={\mathcal {L}}^{-1}\{G(s)\}={\mathcal {L}}^{-1}{\bigg \{}\underbrace {\frac {1}{1+T_{1}\cdot s}} _{\text{s-Bereich}}{\bigg \}}=\underbrace {{\frac {1}{T_{1}}}\cdot e^{-{\frac {t}{T_{1}}}}} _{\text{Zeitbereich}}\qquad {oder}\qquad {\mathcal {L}}^{-1}{\bigg \{}\underbrace {\frac {1}{s+a}} _{\text{s-Bereich}}{\bigg \}}=\underbrace {e^{-a\cdot t}} _{\text{Zeitbereich}}}.

Die Verstärkung K wurde für den Transfer in den Zeitbereich nicht berücksichtigt und wird wieder eingesetzt.

{\displaystyle g(t)=K\cdot {\mathcal {L}}^{-1}\{G(s)\}=K\cdot {\mathcal {L}}^{-1}{\bigg \{}{\frac {1}{1+T_{1}\cdot s}}{\bigg \}}={\frac {K}{T_{1}}}\cdot e^{-{\frac {t}{T_{1}}}}={\frac {K}{T_{1}}}\cdot e^{-a\cdot t}}

Für die Eingangsimpulsfolge g(kT_{A}) wird die z-Transformation angewendet.

Aus einer in Fachbüchern der Regelungstechnik dargestellten Korrespondenztabelle wird für {\displaystyle f(t)=e^{-a\cdot t}} als Gewichtsfunktion (Impulsantwort)
die zugehörige F(z)-Funktion entnommen. Der Faktor K/T_{1} hat keine Auswirkungen auf die Transformation und wird später berücksichtigt:

{\displaystyle {\mathcal {Z}}\{e^{-at}\}{\bigg |}_{t=kT_{A}}={\frac {z}{z-e^{-a\cdot T_{A}}}};\qquad T_{A}={\text{Abtastzeit}};\quad k=Abtastfolge}

Daraus ergibt sich unter Berücksichtigung des Faktors K/T_{1} die gesuchte z-Übertragungsfunktion, die gerne vereinfacht dargestellt wird:

{\displaystyle G(z)={\frac {K}{T_{1}}}\cdot {\mathcal {Z}}\{g(kT_{A})\}={\frac {K}{T_{1}}}\cdot {\frac {z}{z-a_{1}}};\qquad {\text{mit }}a_{1}=e^{-{\frac {T_{A}}{T_{1}}}};\quad T_{1}=1/a}

Beispiel 4: z-Übertragungsfunktion einer Differenzengleichung für den Integralalgorithmus:

Definition der Signale:
Analoge Regeldifferenz = e(t), Abgetastete Regeldifferenz = e_{{(k)}},
Reglerausgangsgröße = u_{(k)}, Zeitkonstante K_{I}=1/T_{N}.

Gegeben: Differenzengleichung:

u_{{(k)}}=u_{{(k-1)}}+K_{I}\cdot T_{A}\cdot e_{{(k-1)}}
(Differenzengleichung nach Euler-Vorwärts, wegen e_{{(k-1)}}) entsteht "Untersumme".

Gesucht: z-Übertragungsfunktion.

Nach der Transformationsvorschrift ändern sich die Größen:

u_{{(k)}}\longrightarrow U(z),\qquad u_{{(k-1)}}\longrightarrow z^{{-1}}\cdot U(z),\qquad e_{{(k-1)}}\longrightarrow z^{{-1}}\cdot E(z).

Daraus folgt für die z-Transformation:

U(z)=z^{{-1}}\cdot U(z)+K_{I}\cdot T_{A}\cdot z^{{-1}}\cdot E(z)

Wird U(z) / E(z) ausgeklammert, entsteht die z-Übertragungsfunktion des Integralalgorithmus:

G_{R}(z)={\frac  {U(z)}{E(z)}}={\frac  {K_{I}\cdot T_{A}\cdot z{-1}}{1-z^{{-1}}}}={\frac  {K_{I}\cdot T_{A}}{z-1}}

Für die gleiche Differenzengleichung des Integrationsalgorithmus nach Euler-Rückwärts (Obersumme):

u_{{(k)}}=u_{{(k-1)}}+K_{I}\cdot T_{A}\cdot e_{{(k)}}

lautet die z-Übertragungsfunktion mit e_{{(k)}}\to E(z):

G_{R}(z)={\frac  {U{(z)}}{E{(z)}}}={\frac  {K_{I}\cdot T_{A}}{1-z^{{-1}}}}

Aus der Übertragungsfunktion GR(z) wird die Stellgröße des Reglers U(z) abgeleitet mit der Regelabweichung E(z) = W(z) - Y(z):

Anwendung der z-Übertragungsfunktion für einen digitalen PI-Regler mit Abtastfunktion und Halteglied

Das nachfolgende Berechnungsbeispiel zeigt die über die z-Transformation ermittelte Differenzengleichung eines digitalen PI-Reglers im Vergleich mit den Differenzengleichungen nach Euler-Rückwärts.

Digitale Regler in industriellen Erzeugnissen werden gegenüber analogen Reglern zunehmend in großen Stückzahlen produziert. Dafür sind verschiedene Gründe gegeben:

Mit der Anwendung der z-Transformation auf abgetastete Signale und des gewünschten Regelalgorithmus entsteht die z-Übertragungsfunktion.

Die Hauptaufgabe der Realisierung eines digitalen Reglers ist die Aufstellung der z-Übertragungsfunktion des Reglers bzw. das Finden der zur Programmierung des Mikrocomputers benötigten Rekursionsgleichung. Das nachstehende Berechnungsbeispiel aus einem Fachbuch zeigt einen digitalen PI-Regler, der durch eine Reihenschaltung von Halteglied, Reglerfunktion GR(s) und Abtaster approximiert wird.

Die gesuchte z-Übertragungsfunktion des digitalen Reglers D(z) besteht aus dem Produkt von Halteglied H_{0}(z)={\frac  {z-1}z} und Reglerfunktion GR(s):

D(z)={\frac  {U{(z)}}{E{(z)}}}={\frac  {z-1}z}\cdot {\mathcal  Z}{\biggl \{}{\frac  {G_{R}(s)}s}{\biggr \}}

Berechnungsbeispiel eines PI-Reglers als Reihenschaltung von idealem Abtaster, PI-Regleralgorithmus und Halteglied.

Gegeben: s-Übertragungsfunktion eines PI-Reglers GR:
G_{R}(s)={\frac  {U{(s)}}{E{(s)}}}={\frac  {K_{P}\cdot (1+T_{N}\cdot s)}{T_{N}\cdot s}}=K_{P}+{\frac  {K_{P}}{T_{N}\cdot s}}

Parameter des PI-Reglers mit folgenden Zahlenwerten:

Verstärkung: K_{P}=2; Zeitkonstante: T_N = 1{,}5[s]; Abtastzeit: T_A = 0{,}1[s];

Gesucht: z-Übertragungsfunktion und daraus die Differenzengleichung mit dem Regelalgorithmus des PI-Reglers:

D(z)={\frac  {U(z)}{E(z)}}=H_{0}(z)\cdot {\mathcal  Z}{\biggl \{}{\frac  {G_{R}(s)}s}{\biggr \}}

Nach Anwendung der oben stehenden abgeleiteten Gleichung mit Eintrag der PI-Komponenten für G_{R}(s) müssen die Komponenten des s-Bereichs in den z-Bereich transformiert werden:

D(z)={\frac  {z-1}z}\cdot {\mathcal  Z}{\biggl \{}{\frac  {K_{P}}s}+{\frac  {K_{P}}{T_{N}\cdot s^{2}}}{\biggr \}}={\frac  {z-1}z}\cdot {\mathcal  Z}{\biggl \{}K_{P}\cdot {\frac  1s}+{\frac  {K_{P}}{T_{N}}}\cdot {\frac  1{s^{2}}}{\biggr \}}

Die Terme des s-Bereichs werden z-transformiert:

D(z)={\frac  {z-1}z}\cdot {\biggl \{}K_{P}\cdot {\frac  z{z-1}}+{\frac  {K_{P}}{T_{N}}}\cdot {\frac  {T_{A}\cdot z}{(z-1)^{2}}}{\biggr \}}={\frac  {K_{P}\cdot T_{N}+K_{P}\cdot (T_{A}-T_{N})\cdot z^{{-1}}}{T_{N}-T_{N}\cdot z^{{-1}}}}

Damit lautet die z-Übertragungsfunktion des Reglers mit Zahlenwerten:

D(z)={\frac  {U{(z)}}{E{(z)}}}={\frac  {3-2,8\cdot z^{{-1}}}{1,5-1,5\cdot z^{{-1}}}}={\frac  {2-1{,}8667\cdot z^{{-1}}}{1-z^{{-1}}}}

Diese z-Übertragungsfunktion wird in Operatorendarstellung nach U(z) geordnet und als Differenzengleichung u_{(k)} im Mikroprozessor programmiert.

U(z)\cdot (1-z^{{-1}})=E(z)\cdot (2-1{,}8667\cdot z^{{-1}})
{\displaystyle U(z)-z^{-1}\cdot U(z)=2\cdot E(z)-1{,}8667\cdot z^{-1}\cdot E(z)}

Durch Anwendung der Rechtsverschiebung um einen Abtastschritt "-1" ergibt sich für Rücktransformation in den diskreten Zeitbereich die Differenzengleichung u_{(k)}.

u_{{(kT_{A})}}-u_{{(kT_{A}-1T_{A})}}=2\cdot e_{{(kT_{A})}}-1{,}8667\cdot e{(kT_{A}-1T_{A})}

Diese Gleichung nach u_{{(kT_{A})}} aufgelöst ergibt:

u_{{(kT_{A})}}=u_{{(kT_{A}-T_{A})}}+2\cdot e_{{(kT_{A})}}-1{,}8667\cdot e_{{(kT_{A}-T_{A})}}

In vereinfachter Schreibweise lautet die Differenzengleichung eines PI-Reglers für ein abgetastetes Regelsystem mit Halteglied bei TA = 0,1 [s]:

u_{{(k)}}=u_{{(k-1)}}+2\cdot e_{{(k)}}-1{,}8667\cdot e_{{(k-1)}}

Zum Vergleich die errechnete Differenzengleichung für TA = 0,01 [s] zur besseren Annäherung an die analytische Funktion:

u_{{(k)}}=u_{{(k-1)}}+2\cdot e_{{(k)}}-1{,}9867\cdot e_{{(k-1)}}
Sprungantwort eines digitalen PI-Reglers mit A/D-Wandler und Abtast-Halteglied am Reglerausgang.
Tabelle zur grafischen Darstellung der Sprungantwort des PI-Reglers mit TA = 0,1 [s] und e(k) = 1,
im Vergleich mit einer direkten Anwendung von Differenzengleichungen nach Euler-Rückwärts.
Folge k Zeit TA
[sec]
Differenzengleichung aus z-Transformation
Ausgangssignal u_{(k)}
Differenzengleichung
Euler-Rückwärts
u_{(k)} Untersumme
Differenzengleichung
Euler-Rückwärts
u_{(k)} Obersumme
0 0 u_{{(k)}}=0+2\cdot 1-1{,}8667\cdot 0=2{,}0000 2,0000 2,1333
1 0,1 u_{{(k)}}=2{,}0000+2\cdot 1-1{,}8667\cdot 1=2{,}1333 2,1333 2,2667
2 0,2 u_{{(k)}}=2{,}1333+2\cdot 1-1{,}8667\cdot 1=2{,}2666 2,2667 2,4000
3 0,3 u_{{(k)}}=2{,}2666+2\cdot 1-1{,}8667\cdot 1=2{,}3999 2,4000 2,5333
4 0,4 u_{{(k)}}=2{,}3999+2\cdot 1-1{,}8667\cdot 1=2{,}5332 2,5333 2,6667
5 0,5 u_{{(k)}}=2{,}5332+2\cdot 1-1{,}8667\cdot 1=2{,}6665 2,6667 2,8000
6 0,6 u_{{(k)}}=2{,}6665+2\cdot 1-1{,}8667\cdot 1=2{,}7998 2,8000 2,9333
       
10 1,0 3,3330 3,3330 3,4667
15 1,5 3,9995 4,0000 4,1333

Qualitativer Vergleich zwischen PI-Reglern mit Differenzengleichungen über die z-Transformation und Differenzengleichungen nach Euler-Rückwärts

Das Fachbuch-Rechenbeispiel der Bestimmung der Differenzengleichung eines PI-Reglers über die z-Transformation ist offensichtlich als Untersumme definiert. Dieses Beispiel setzt ideale Schnittstellen mit einem Halteglied nullter Ordnung (entspricht einem Zeitschritt TA) und einen idealen Mikrorechner voraus.
Die Grafik zeigt den Unterschied der Stellgröße für die Berechnungsmethoden der Differenzengleichungen nach der Ober- und Untersumme.
Die grafische Darstellung der Sprungantwort des PI-Reglers laut der Tabellen-Berechnungspunkte y_{(k)} entsprechen einer Rechteckfunktion, weil am Reglerausgang für die Wertefolge ein Halteglied für den Zeitraum TA wirkt, das nach jedem Rechenschritt gelöscht wird. Die errechneten Daten der Differenzengleichung über die z-Transformation und die Daten nach den Differenzengleichungen Euler-Rückwärts-Untersumme sind identisch.
Sowohl die errechneten Ausgangsgrößen des Reglers der Differenzengleichung als Wertefolgen y_{(k)} über die z-Übertragungsfunktion wie auch die errechneten Wertefolgen der Ausgangsgrößen y_{(k)} nach den Differenzengleichungen der Methode Euler-Rückwärts (Untersumme und Obersumme) entsprechen nur dann der analytischen Funktion der Sprungantwort des PI-Reglers, wenn die Abtastzeit T_{A}=\Delta t gegen Null geht.
Der Vergleich der Rechenergebnisse führt zu der Auffassung, dass für die Bestimmung der Differenzengleichungen der Standardregler die im Abschnitt Tabelle der Differenzengleichungen der Standardregler (Euler Rückwärts) aufgeführten Differenzengleichungen unmittelbar ohne z-Transformation verwendet werden können.
Aus der Übertragungsfunktion des s-Bereiches wird deutlich, dass der PI-Regler in der Reihenstruktur aus einem I-Glied und einem PD1-Glied besteht.
Übertragungsfunktion des PI-Reglers in der Parallelstruktur:
{\frac  U{E}}{(s)}=K_{P}+{\frac  {K_{P}}{T_{N}\cdot s}}=K_{P}\left(1+{\frac  1{T_{N}\cdot s}}\right)
Wird der Klammerausdruck der Gleichung auf einen gemeinsamen Nenner gebracht, entsteht die Produktdarstellung in der Reihenstruktur:
{\frac  U{E}}{(s)}=K_{{P}}{\frac  {T_{N}\cdot s+1}{T_{N}\cdot s}}=G_{I}(s)\cdot G_{{PD}}(s)={\frac  {K_{P}}{T_{N}\cdot s}}\cdot (T_{N}\cdot s+1)
KPI = KP / TN ist die Verstärkung des PI-Reglers.
Für die Berechnung des PI-Reglers werden zwei Differenzengleichungen 1. Ordnung benötigt. Dabei ist wegen der empfohlenen Reihenschaltung die Ausgangsgröße des I-Gliedes die Eingangsgröße des PD1-Gliedes. Die Ausgangsgröße des PD-Gliedes im zeitdiskreten Bereich ist das Folgeglied zum Zeitpunkt k\cdot \Delta t. Die Eingangsgröße des PI-Reglers ist die normierte zeitdiskrete Sprungfunktion 1(t)=1(0);1(\Delta t);1(2\Delta t);1(3\Delta t);\dots .
Die Differenzengleichung eines PI-Reglers nach der Untersumme als Reihenschaltung von I-Glied und PD1-Glied lautet:
{\text{I-Glied:}}\qquad e_{{1(k)}}=e_{{1(k-1)}}+{e_{{(k)}}\cdot K_{P}\cdot \Delta t}/{T_{N}}
{\text{PD-Glied:}}\quad u_{{(k)}}=e_{{1(k-1)}}+[e_{{1(k)}}-e_{{1(k-1)}}]\cdot {T_{N}}/{\Delta t}
Anmerkung zum PD-Glied: e_{{1(k-1)}} bedeutet wegen Untersumme: aktuelles e_{{1(k)}} des I-Gliedes um einen Schritt (k-1) zurückgesetzt.
Die Güte der Approximation ist durch das Verhältnis der Abtastzeit T_{A}=\Delta t zur dominanten Zeitkonstante T des Reglers und der Regelstrecke bestimmt. Soll ein PD-Glied des PI-Reglers ein PT1-Glied der Regelstrecke kompensieren, sollte dieses Verhältnis 0,1 bis 0,01 betragen.
Welche Genauigkeit des Regelvorgangs erreicht werden muss, hängt von der Vorgabe des Lastenheftes und von der Art der Regelstrecke ab.
Bei schnellen Regelstrecken können die realen Schnittstellen des Reglers durch die signaltechnischen Verzögerungen als Ersatztotzeit berücksichtigt werden. Im Gegensatz zur Laplace-Transformation kann eine z-transformierte Übertragungsfunktion mit einer Totzeit im Rahmen der Rechenregeln der z-Transformation durch Rückwärtsverschiebung "Verschiebung nach rechts" um d = Tt / TA-Abtastschritte als komplette z-Übertragungsfunktion berechnet werden.

Einsatz von z-Transformation zur Bestimmung expliziter Formeln von Rekursionsvorschriften

Motivation

Die z-Transformation kann zur Bestimmung expliziter Formeln für Zahlenfolgen eingesetzt werden, die rekursiv definiert sind. Ein Paradebeispiel hierfür ist die Fibonacci-Zahlenfolge

{\displaystyle 0,1,1,2,3,5,8,13,\dots }

Die Zahlenfolge ist rekursiv definiert und zwar mit der Rekursionsvorschrift

{\displaystyle y_{k}=y_{k-1}+y_{k-2},\quad k>1}

und dem Rekursionsanfang

{\displaystyle y_{0}=0\quad {\text{und}}\quad y_{1}=1.}

Die rekursive Definition ist einfach, hat jedoch zum Nachteil, dass sie es nicht erlaubt, das k-te Element der Zahlenfolge direkt zu bestimmen, ohne seine k-1 Vorgänger zu berechnen, was bei größer werdendem n immer zeitintensiver wird. Man fragt sich also, ob nicht auch eine explizite Formel existiert, die eine direkte Bestimmung eines jeden Elements der Zahlenfolge erlaubt. Für Fibonacci-Zahlen existiert eine solche explizite Formel und sie lautet

{\displaystyle y_{k}={\frac {1}{\sqrt {5}}}\cdot \left(\left({\frac {1+{\sqrt {5}}}{2}}\right)^{k}-\left({\frac {1-{\sqrt {5}}}{2}}\right)^{k}\right),\quad k\geq 0.}.

Auf den ersten Blick scheint es unmöglich zu sein, dass diese Formel mit {\sqrt {5}} Termen überhaupt ganze Zahlen erzeugt, doch durch Einsetzen von {\displaystyle k=0,1,2,3,\dots } kann sich jeder selbst davon überzeugen, dass die Formel in der Tat die entsprechenden Fibonacci-Zahlen liefert.

Doch wie kommt man nur auf eine solche Formel? Dazu gibt es verschiedenen Methoden und die z-Transformation ist eine davon. Im Folgenden wird eine Methodik gezeigt, wie man mit Hilfe von z-Transformation die explizite Formel einer rekursiven Berechnungsvorschrift bestimmen kann. Die allgemeinen Ausführungen werden am Beispiel der Fibonacci-Zahlen angewendet und verständlich gemacht.

Allgemeine Form der Rekursionsgleichung

Differenzengleichung und Rekursion

Eine Rekursionsformel vom Typ

{\displaystyle y_{k}=a_{1}\cdot y_{k-1}+a_{2}\cdot y_{k-2}+\dots +a_{n}\cdot y_{k-n}\quad {\text{mit }}y_{k}=0{\text{ für }}k<0}

kann als Beschreibung eines abgetasteten dynamischen Systems ohne Eingangssignal interpretiert werden, denn die weiter oben beschriebene Differenzengleichung

{\displaystyle {y_{(k)}=b_{0}\cdot u_{(k)}+b_{1}\cdot u_{(k-1)}+b_{2}\cdot u_{(k-2)}+\dots +b_{m}\cdot u_{(k-m)}\quad -a_{1}\cdot y_{(k-1)}-a_{2}\cdot y_{(k-2)}-\dots -a_{n}\cdot y_{(k-n)}}}.

geht direkt in die Rekursionsformel über, wenn man das Eingangssignal u zu Null setzt, d.h. {\displaystyle u_{k}=0} für alle k:

{\displaystyle {y_{(k)}=-a_{1}\cdot y_{(k-1)}-a_{2}\cdot y_{(k-2)}-\dots -a_{n}\cdot y_{(k-n)}}.}

Das bedeutet, dass eine Rekursionsformel als eine Art Eigendynamik des Systems betrachtet werden kann, wenn auf dieses keine äußeren Störungen einwirken. Die Eigendynamik wird allerdings erst sichtbar, wenn das System aus der Ruhelage ausgelenkt wird, und der genaue Verlauf der „Bewegung“ hängt von den Anfangsbedingungen ab. Man setzt also voraus, dass das System bis zum Zeitpunkt Null keinerlei Dynamik aufwies (d.h. {\displaystyle y_{k}=0,\ \forall k<0}), doch die Anfangsbedingungen (Rekursionsanfang) müssen unbedingt berücksichtigt werden.

Anfangsbedingungen und Rekursionsanfang

Um die Rekursionsformel mit Hilfe der z-Transformation lösen zu können, muss der Rekursionsanfang (also die Anfangsbedingung) in die Differenzengleichung eingefügt werden. Dies kann man mit Hilfe des zeitdiskreten Dirac-Impulses mit der Definition

{\displaystyle \delta _{k}={\begin{cases}1\ &{\text{falls }}k=0\\0\ &{\text{sonst}}\end{cases}}}

umsetzen. Man platziert eine Reihe von n diskreten \delta -Impulsen an den ersten n Zeitpunkten

{\displaystyle {\begin{aligned}y_{k}&=a_{1}\cdot y_{k-1}+a_{2}\cdot y_{k-2}+\dots +a_{n}\cdot y_{k-n}\quad +d_{0}\cdot \delta _{k}+d_{1}\cdot \delta _{k-1}+d_{2}\cdot \delta _{k-2}+\dots +d_{n-1}\cdot \delta _{k-(n-1)}\\&=\sum _{i=1}^{n}a_{i}\cdot y_{k-i}\quad +\sum _{j=0}^{n-1}d_{j}\cdot \delta _{k-j}\end{aligned}}}

und wählt deren Amplituden d_{i} so, dass sich die gewünschten n Anfangswerte der Zahlenfolge ergeben. Die Anzahl der Anfangswerte muss n sein, damit die Rekursionsformel starten kann, denn n ist die Anzahl der Rekursionselemente auf der rechten Seite der Rekursionsgleichung. Angabe weiterer Anfangswerte (mehr als n) ist überflüssig, weil diese Werte direkt aus der Rekursionsformel berechnet werden.

Unter Berücksichtigung dass {\displaystyle y_{k}=0} für alle negativen k angesetzt ist und der diskrete Dirac-Impuls {\displaystyle \delta _{k-j}} nur an der Stelle {\displaystyle k=j} den Wert 1 hat und ansonsten 0 erhält man etwa für die ersten drei Elemente die Beziehungen:

{\displaystyle {\begin{aligned}y_{0}&=\sum _{i=1}^{n}a_{i}\cdot \underbrace {y_{0-i}} _{=0}+\sum _{j=0}^{n-1}d_{j}\cdot \delta _{0-j}=d_{0}\cdot \delta _{0}=d_{0}\\y_{1}&=\sum _{i=1}^{n}a_{i}\cdot y_{1-i}+\sum _{j=0}^{n-1}d_{j}\cdot \delta _{1-j}=a_{1}\cdot \underbrace {x_{0}} _{=d_{0}}+d_{1}=a_{1}\cdot d_{0}+d_{1}\\y_{2}&=\sum _{i=1}^{n}a_{i}\cdot y_{2-i}+\sum _{j=0}^{n-1}d_{j}\cdot \delta _{2-j}=a_{1}\cdot x_{1}+a_{2}\cdot x_{0}+d_{2}=a_{1}\cdot (a_{1}\cdot d_{0}+d_{1})+a_{2}\cdot d_{0}+d_{2}\end{aligned}}}

Sind diese drei Anfangswerte {\displaystyle y_{0},\,y_{1}} und y_2 als Rekursionsanfang bekannt, so kann man aus diesen Gleichungen direkt die Amplituden {\displaystyle d_{0},\,d_{1}} und d_{2} bestimmen:

{\displaystyle {\begin{aligned}d_{0}&=y_{0}\\d_{1}&=y_{1}-a_{1}\cdot d_{0}\\d_{2}&=y_{2}-a_{1}\cdot (a_{1}\cdot d_{0}+d_{1})-a_{2}\cdot d_{0}=y_{2}-a_{1}\cdot (a_{1}\cdot d_{0}+\underbrace {y_{1}-a_{1}\cdot d_{0}} _{d_{1}})-a_{2}\cdot d_{0}=y_{2}-a_{1}\cdot y_{1}-a_{2}\cdot d_{0}\end{aligned}}}
Im Beispiel der Fibonacci-Zahlen gibt es zwei Anfangsbedingungen
y_{0}=0 und {\displaystyle y_{1}=1.}

Die Differenzengleichung wird also mit

{\displaystyle y_{k}=y_{k-1}+y_{k-2}+d_{0}\cdot \delta _{k}+d_{1}\cdot \delta _{k-1},\quad k\geq 0}

angesetzt. Das heißt, die beiden Koeffizienten a_{1} und a_2 sind beide 1, so dass man für die Amplituden der beiden Dirac-Impulse die Werte

{\displaystyle {\begin{aligned}d_{0}&=y_{0}=0\\d_{1}&=y_{1}-a_{1}\cdot d_{0}=1\end{aligned}}}

und damit die Differenzengleichung

{\displaystyle y_{k}=y_{k-1}+y_{k-2}+\delta _{k-1}}

erhält.

Lösung der Rekursionsgleichung mittels z-Transformation

Nun kann die gewonnene Differenzengleichung

{\displaystyle y_{k}=\sum _{i=1}^{n}a_{i}\cdot y_{k-i}\quad +\sum _{j=0}^{n-1}d_{j}\cdot \delta _{k-j}}

mittels z-Transformation gelöst werden. Im ersten Schritt erhält man für die z-Transformierte den Ausdruck

{\displaystyle Y(z)={\mathcal {Z}}\{y_{k}\}=\sum _{i=1}^{n}a_{i}\cdot Y(z)\cdot z^{-i}\quad +\sum _{j=0}^{n-1}d_{j}\cdot 1\cdot z^{-j},}

in dem man die Y(z)-Terme auf die linke Seite bringen

{\displaystyle Y(z)\cdot \left(1-\sum _{i=1}^{n}a_{i}\cdot Y(z)\cdot z^{-i}\right)=\sum _{j=0}^{n-1}d_{j}\cdot z^{-j}}

und anschließend nach Y(z)auflösen kann:

{\displaystyle Y(z)={\frac {\sum _{j=0}^{n-1}d_{j}\cdot z^{-j}}{1-\sum _{i=1}^{n}a_{i}\cdot Y(z)\cdot z^{-i}}}}

bzw. durch Multiplikation mit z^{n} im Zähler und Nenner:

{\displaystyle Y(z)={\frac {\sum _{j=0}^{n-1}d_{j}\cdot z^{n-j}}{z^{n}-\sum _{i=1}^{n}a_{i}\cdot Y(z)\cdot z^{n-i}}}}

Im Allgemeinen erhält man also Polynome gleicher Ordnung im Nenner und im Zähler. Das Zählerpolynom kann nie eine höhere Ordnung als das Nennerpolynom besitzen, weil wie bereits erwähnt, lediglich n Anfangswerte vorgegeben werden so dass die Rekursion beginnen kann. Es kann jedoch vorkommen, dass die Anfangsbedingungen zu Nullwerten bei manchen d_j Koeffizienten führen, so dass das Zählerpolynom eine kleinere Ordnung als das Nennerpolynom besitzt. Mit diesen Überlegungen kann man also davon ausgehen, dass eine Partialbruchzerlegung möglich ist und die z-Transformierte in Form

{\displaystyle {\begin{aligned}Y(z)&=\sum _{i=1}^{n}{\frac {\alpha _{i}\cdot z+\beta _{i}}{z-p_{i}}}\\&=\sum _{i=1}^{n}\left({\frac {\alpha _{i}\cdot z}{z-p_{i}}}+{\frac {\beta _{i}}{z-p_{i}}}\right)\\&=\sum _{i=1}^{n}\left({\frac {\alpha _{i}\cdot z}{z-p_{i}}}+{\frac {\beta _{i}\cdot z}{z-p_{i}}}\cdot z^{-1}\right)\end{aligned}}}

geschrieben werden kann. Wegen der Beziehung

{\displaystyle {\mathcal {Z}}\{a^{k}\}={\frac {z}{z-a}}}

und der Rückwärtsverschiebung (nach rechts) der z-Transformation

{\displaystyle {\mathcal {Z}}\{y_{k-1}\}=z^{-1}\cdot Y(z)}

erhält man also für die Rücktransformierte

{\displaystyle {\begin{aligned}y_{k}=\sum _{i=1}^{n}\left(\alpha _{i}\cdot p_{i}^{k}+\beta _{i}\cdot p_{i}^{k-1}\right)=\sum _{i=1}^{n}p_{i}^{k-1}\cdot \left(\alpha _{i}\cdot p_{i}+\beta _{i}\right)\end{aligned}}}

und somit die explizite Formel für die einzelnen Elemente y_{k} der rekursiv definierten Zahlenfolge.

Im Beispiel der Fibonacci-Zahlen mit der Differenzengleichung
{\displaystyle y_{k}=y_{k-1}+y_{k-2}+\delta _{k-1}}

erhält man zunächst durch z-Transformation

{\displaystyle Y(z)=Y(z)\cdot z^{-1}+Y(z)\cdot z^{-2}+1\cdot z^{-1}}

und daraus

{\displaystyle Y(z)={\frac {z}{z^{2}-z-1}}}

bzw. mit der Partialbruchzerlegung

{\displaystyle Y(z)={\frac {1}{\sqrt {5}}}\cdot \left({\frac {z}{z-{\frac {1+{\sqrt {5}}}{2}}}}-{\frac {z}{z-{\frac {1-{\sqrt {5}}}{2}}}}\right).}

Die Rücktransformation ergibt nun direkt die gesuchte explizite Form für die Fibonacci-Zahlenfolge:

{\displaystyle y_{k}={\frac {1}{\sqrt {5}}}\cdot \left(\left({\frac {1+{\sqrt {5}}}{2}}\right)^{k}-\left({\frac {1-{\sqrt {5}}}{2}}\right)^{k}\right).}

Dass diese „krumme“ Formel tatsächlich ganzzahlige Ergebnisse liefert wird übrigens sofort klar, wenn man die Binomialausdrücke ausschreibt. Wegen des Minuszeichens heben sich die Terme mit geraden Exponenten von {\sqrt {5}} auf während die Terme mit ungeraden Exponenten sich addieren und mit dem {\sqrt {5}}-Term im Nenner kürzen. Damit verbleiben lediglich gerade Exponenten von {\sqrt {5}}, d.h. Exponenten von 5. Man kann also die obige Gleichung auf mit Fünfer-Potenzen gewichtete Binomialkoeffizienten und das Pascalsche Dreieck zurückführen (siehe detaillierte Herleitung im Artikel Verwandtschaft der Fibonacci-Zahlen mit dem Pascalschen Dreieck).

Siehe auch

Trenner
Basierend auf einem Artikel in: Wikipedia.de
Seitenende
Seite zurück
©  biancahoegel.de
Datum der letzten Änderung:  Jena, den: 04.07. 2021