Harmonic Oscillation

Definition und Eigenschaften

Besondere Bedeutung für die Nachrichtentechnik – aber auch in vielen Naturwissenschaften – haben harmonische Schwingungen. Die folgende Grafik zeigt einen beispielhaften Signalverlauf.

Ihre Bedeutung hängt auch damit zusammen, dass die harmonische Schwingung die Lösung einer in vielen Disziplinen vorkommenden Differentialgleichung darstellt, die wie folgt lautet:

$$ x(t) + k \cdot\ddot{x} (t) = 0.$$

Hierbei kennzeichnen die beiden Punkte die zweite Ableitung der Funktion $x(t)$ nach der Zeit.

Eine jede harmonische Schwingung kann man in allgemeinster Form wie folgt darstellen:

$$x(t)= C \cdot \cos(2\pi f_0 t - \varphi).$$

Hierbei sind folgende Signalparameter verwendet:

die Amplitude $C$ - gleichzeitig der Maximalwert des Signals,
die Signalfrequenz $f_{0}$, und
der Nullphasenwinkel (oder kurz die Phase) $\phi$ der Schwingung.

Anmerkung:

In diesem Tutorial geht – wie auch in anderer Literatur üblich – bei der Beschreibung von harmonischen Schwingungen, Fourierreihe und Fourierintegral die Phase mit negativem Vorzeichen in die Gleichungen ein, während in Zusammenhang mit allen Modulationsverfahren die Phase stets mit einem Pluszeichen angesetzt wird.
Zur Unterscheidung der beiden Varianten benutzen wir in LNTwww $\varphi$ und $\phi$. Beide Symbole kennzeichnen das kleine griechische „phi”, wobei die Schreibweise $\varphi$ vorwiegend im deutschen und $\phi$ im anglo-amerikanischen Sprachraum angewandt wird.
Die Angaben $\varphi = 90^{\circ}$ und $\phi = -90^{\circ}$ sind somit äquivalent und stehen beide für die Sinusfunktion:

$$\cos(2 \pi f_0 t - 90^{\circ}) = \cos(2 \pi f_0 t - \varphi) = \cos(2 \pi f_0 t + \phi) = \sin(2 \pi f_0 t ).$$

Das folgende Lernvideo verdeutlicht die Eigenschaften harmonischer Schwingungen anhand von Tonleitern:

Harmonische Schwingungen (Dauer Teil 1: 4:33 – Teil 2: 6:15)

Zeitsignaldarstellung

Signalparameter einer harmonischen Schwingung

Die Amplitude $C$ kann aus der nebenstehenden Grafik direkt abgelesen werden. Die Signalfrequenz $f_0$ ist gleich dem Kehrwert der Periodendauer $T_0$. Schreibt man die obige Gleichung in der Form

$$x(t) = C \cdot \cos(2\pi f_0 t - \varphi) = C \cdot \cos(2\pi f_0 (t - \tau)), $$

so wird klar, dass der Nullphasenwinkel $\phi$ und die Verschiebung $\tau$ gegenüber einem cosinusförmigen Signal wie folgt zusammenhängen:

$$\varphi = \frac{\tau}{T_0} \cdot 2{\pi}. $$

Bei einem Cosinussignal sind die Kenngrößen $\tau$ und $\varphi$ jeweils 0. Demgegenüber ist ein Sinussignal um $\tau = T_0/4$ verschoben und entsprechend gilt für den Nullphasenwinkel $\varphi = \pi/2$ (im Bogenmaß) bzw. $90^{\circ}$.

Es ist also festzustellen, dass – wie für das obige Beispiel vorausgesetzt – bei einem positiven Wert von $\tau$ bzw. $\varphi$ das (bezüglich $t = 0$) nächstgelegene Signalmaximum später kommt als beim Cosinussignal und bei negativen Werten früher. Liegt am Systemeingang ein Cosinussignal an und ist das Ausgangssignal demgegenüber um einen Wert $\tau$ verzögert, so bezeichnet man $\tau$ auch als die Laufzeit des Systems.

Da eine harmonische Schwingung durch lediglich drei Signalparameter eindeutig festliegt, kann der gesamte Zeitverlauf von $-\infty$ bis $+\infty$ aus nur drei Signalwerten $x_1=x(t_1)$, $x_2=x(t_2)$, $x_3=x(t_3)$ analytisch beschrieben werden, wenn die Zeiten $t_1$, $t_2$ und $t_3$ geeignet festgelegt wurden.

Harmonische Schwingung, festgelegt durch drei Abtastwerte

Aus den drei Abtastwerten,

$x_1 = x(t_1 = 3.808 \;{\rm ms}) = +1.609$,
$x_2 = x(t_2 = 16.696 \;{\rm ms})= –0.469$,
$x_3 = x(t_3 = 33.84 \;{\rm ms}) = +1.227$.

erhält man folgendes Gleichungssystem:

$C \cdot \cos(2\pi \hspace{0.05 cm} f_0 \hspace{0.05 cm} t_1 - \varphi) \hspace{-0.15 cm} = \hspace{-0.15 cm}+1.609\hspace{0.05 cm},\\ C \cdot \cos(2\pi \hspace{0.05 cm} f_0 \hspace{0.05 cm} t_2 - \varphi) \hspace{-0.15 cm} = \hspace{-0.15 cm}-0.469\hspace{0.05 cm},\\ C \cdot \cos(2\pi \hspace{0.05 cm} f_0 \hspace{0.05 cm} t_3 - \varphi) \hspace{-0.15 cm} = \hspace{-0.15 cm}+1.227\hspace{0.05 cm}.$

Nach Lösen dieses nichtlinearen Gleichungssystem ergeben sich folgende Signalparameter:

Signalamplitude $C = 2$,
Periodendauer $T_0 = \;{\rm ms}$ ⇒ Signalfrequenz $f_0 = 125 \;{\rm Hz}$,
Verschiebung gegenüber einem Cosinus $\tau = 3 \;{\rm ms}$ ⇒ Nullphasenwinkel $\varphi = 3\pi /4 = 135^\circ$.

Hinweis: Legt man alle Abtastzeitpunkte $t_1$,$t_2$, $t_3$ in Maxima, Minima und/oder Nullstellen, so gibt es für das nichtlineare Gleichungssystem keine eindeutige Lösung.

Darstellung mit Cosinus- und Sinusanteil

Eine weitere Darstellungsform der harmonischen Schwingung lautet wie folgt:

$$x(t)=A\cdot\cos(2\pi f_0 t)+ B\cdot\sin(2\pi f_0 t).$$

Die Bezeichnungen $A$ und $B$ für die Amplituden von Cosinus– und Sinusanteil sind so gewählt, dass sie mit der Nomenklatur des folgenden Kapitels Fourierreihe übereinstimmen.

Durch Anwendung trigonometrischer Umformungen erhalten wir aus der Darstellung auf der letzten Seite:

$$x(t)=C\cdot \cos(2\pi f_0 t-\varphi)=C\cdot\cos(\varphi)\cdot\cos(2\pi f_0 t)+C\cdot\sin(\varphi)\cdot\sin(2\pi f_0 t).$$

Daraus folgt durch Koeffizientenvergleich direkt:

$$A=C\cdot\cos(\varphi),\hspace{0.5cm}B=C\cdot\sin(\varphi).$$

Der Betrag und der Nullphasenwinkel der harmonischen Schwingung können aus den Parametern $A$ und $B$ ebenfalls nach einfachen trigonometrischen Überlegungen berechnet werden:

$$C=\sqrt{A^2+B^2}, \hspace{0.5 cm}\varphi = \arctan\left({-B}/{A}\right).$$

Das Minuszeichen bei der Berechnung des Nullphasenwinkels $\varphi$ hängt damit zusammen, dass $\varphi$ in das Argument der Cosinusfunktion mit negativem Vorzeichen eingeht. Würde man statt „$\cos(2\pi f_0 t - \varphi)$” die Schreibweise „$\cos(2\pi f_0 t + \phi)$” verwenden, so würde $\phi= \arctan(B/A)$ gelten. Bitte beachten Sie:

Bei Fourierreihe und Fourierintegral ist in der Literatur die $\varphi$–Darstellung üblich.
Zur Beschreibung der Modulationsverfahren verwendet man fast immer die $\phi$–Darstellung.

Die in der linken Grafik dargestellte Schwingung ist durch die Parameter $f_0 = 125 \;{\rm Hz}$ und $\varphi = +135^{\circ}$ ⇒ $\phi = -135^{\circ}$ gekennzeichnet.

Die Schwingung wird durch jede der beiden Gleichungen vollständig beschrieben:

$$x(t) = 2\cdot \cos(2\pi f_0 t-135^{\circ})\hspace{0.05cm},$$ $$x(t) = -\sqrt{2}\cdot \cos(2\pi f_0 t) +\sqrt{2}\cdot \sin(2\pi f_0 t)\hspace{0.05cm}.$$

Die rechte Skizze verdeutlicht die trigonometrische Umformung:

$$A = 2\cdot \cos(-135^{\circ}) = -\sqrt{2}\hspace{0.05cm},$$ $$B = 2\cdot \sin(-135^{\circ}) = +\sqrt{2}\hspace{0.05cm}.$$

Spektraldarstellung eines Cosinussignals

Zur Herleitung der Spektralfunktion beschränken wir uns zunächst auf ein Cosinussignal, das mit der komplexen Exponentialfunktion und dem Satz von Euler auch in folgender Weise geschrieben werden kann:

$$x(t)=A \cdot \cos(2\pi f_0 t)={A}/{2}\cdot({\rm e}^{\rm -j2 \pi \it f_{\rm 0} t} + {\rm e}^{\rm j2\pi \it f_{\rm 0} t}).$$

Bereits aus dieser Zeitbereichsdarstellung ist ersichtlich, dass das Cosinussignal – spektral gesehen – nur eine einzige (physikalische) Frequenz beinhaltet, nämlich die Frequenz $f_0$.

Zur mathematischen Herleitung der Spektralfunktion benutzen wir folgende Beziehungen:

den auf der Seite Diracfunktion im Frequenzbereich hergeleiteten Funktionalzusammenhang:

$x(t)=A \circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\, X(f)=A \cdot \rm \delta (\it f).$

den Verschiebungssatz (für den Frequenzbereich) im Vorgriff auf ein späteres Kapitel:

$x(t) \cdot {\rm e}^{\rm j2\pi\it f_{\rm 0} t}\circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\,X(f-f_0). $

Daraus ergibt sich die folgende Fourierkorrespondenz:

$$x(t)=A\cdot \cos(2\pi f_0t)\circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\,X(f)={A}/{\rm 2}\cdot \rm \delta (\it f+f_{\rm 0})+{A}/{\rm 2}\cdot \rm \delta(\it f-f_{\rm 0}).$$

Das bedeutet: Die Spektralfunktion $X(f)$ eines Cosinussignals mit der Frequenz $f_0$ setzt sich aus zwei Diracfunktionen bei $\pmf_0$ zusammen. Die Impulsgewichte sind jeweils gleich der halben Signalamplitude.

Das Bild zeigt das Spektrum einer Cosinusschwingung mit Amplitude $A$ = 4 V und Frequenz $f_0$ = 5 kHz ⇒ $T_0$ = 200 μs. Die Diracfunktion bei $-f_0$ gehört zum ersten Term obiger Gleichung (ableitbar aus der Bedingung $f + f_0 = 0$), die bei $+f_0$ zum Term $\delta (f_0)$. Die Impulsgewichte sind jeweils 2 V.

Es wird darauf hingewiesen, dass die Spektralfunktion einer jeden reellen Zeitfunktion mit Ausnahme des Gleichsignals sowohl Anteile bei positiven als auch bei negativen Frequenzen aufweist. Diese Tatsache, die Studienanfängern oft Probleme bereitet, ergibt sich ganz formal aus dem Satz von Euler (siehe oben). Durch die Erweiterung des Frequenzwertebereichs von $f \ge 0$ auf die Menge der reellen Zahlen kommt man von der physikalischen zur mathematischen Frequenz. Allerdings ist für eine negative Frequenz die vorne angegebene Definition nicht mehr anwendbar: Man kann –5 kHz nicht als „minus 5000 Schwingungen pro Sekunde” interpretieren. Im Verlauf dieses Kurses werden Sie feststellen, dass durch die Verkomplizierung des einfachen Sachverhaltes später kompliziertere Sachverhalte sehr elegant und einfach beschrieben werden können.

Allgemeine Spektraldarstellung

Für ein sinusförmiges Signal gilt mit dem Satz von Euler in ähnlicher Weise:

$$x(t)=B\cdot \sin(2 \pi f_0 t)= \frac{\it B}{2 \cdot \rm j}({\rm e}^{{\rm j} 2 \pi \it f_{\rm 0} t}-{\rm e}^{-\rm j2 \pi \it f_{\rm 0} t})=\rm j\cdot {\it B}/{2}({\rm e}^{-j2 \pi \it f_{\rm 0} t}-{\rm e}^{\rm j2 \pi \it f_{\rm 0} t}) .$$

Daraus folgt für die Spektralfunktion, die jetzt rein imaginär ist:

$$x(t)=B\cdot \sin(2\pi f_0 t)\circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\,X(f)={\rm j} \cdot \left [ {B}/{2} \cdot \delta (f+f_0)- {B}/{2} \cdot \delta (f-f_0) \right ].$$

Das Bild zeigt die rein imaginäre Spektralfunktion einer Sinusschwingung $x(t)$ mit Amplitude $B$ = 3 V und Frequenz $f_0$ = 5 kHz. Die Phase ist $\phi=90^circ$ ⇒ $\phi = 90^circ$. Beachten Sie: Bei der positiven Frequenz (+f_0) ist der Imaginärteil negativ und bei der negativen Frequenz (-f_0) positiv.

Bei Überlagerung von Cosinus– und Sinusanteil entsprechend der Beziehung

$$x(t)=A \cdot \cos(2\pi f_0 t) +B \cdot \sin(2 \pi f_0 t)$$

überlagern sich auch die einzelnen Spektralfunktionen und man erhält:

$$X(f)=\frac{A+{\rm j} \cdot B}{2}\, {\rm \delta} (f+f_0)+\frac{A-{\rm j} \cdot B}{2}\, \delta (f-f_0).$$

Mit dem Betrag $C$ und der Phase $\phi$ lautet diese Fourierkorrespondenz:

$$x(t)=C\cdot \cos(2\pi f_0 t-\varphi) \circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\, X(f)=\frac{C}{2}\cdot {\rm e}^{{\rm j} \varphi} \cdot \delta(f+f_0) + \frac{C}{2} \cdot {\rm e}^{\rm{-j} \varphi} \cdot \delta (f-f_0).$$

Man erkennt, dass die Spektralfunktion $X(f)$ nicht nur für positive und negative Frequenzen definiert ist, sondern im Allgemeinen auch noch komplexwertig ist. Beispiel: Mit den Parametern $C$ = 5 V, $f_0$ = 5 kHz und $\phi$ = 30 Grad (im Bogenmaß $\pi$/6) ergibt sich wegen 2.5 · cos(30°) = 2.165 und 2.5 · sin(30°) = 1.25 für der Real- bzw. der Imaginärteil von $X(f)$ gemäß folgender Grafik: