Loading [MathJax]/jax/output/HTML-CSS/fonts/TeX/fontdata.js

Applets:Physical Signal & Analytic Signal

From LNTwww
Revision as of 18:50, 10 July 2018 by Xiaohan (talk | contribs)

Open Applet in a new tab

Applet Description


Dieses Applet zeigt den Zusammenhang zwischen dem physikalischen Bandpass–Signal x(t) und dem dazugehörigen analytischen Signal x+(t). Ausgegangen wird stets von einem Bandpass–Signal x(t) mit frequenzdiskretem Spektrum X(f):

x(t)=xU(t)+xT(t)+xO(t)=AUcos(2πfUtφU)+ATcos(2πfTtφT)+AOcos(2πfOtφO).

Das physikalische Signal x(t) setzt sich also aus drei harmonischen Schwingungen zusammen, einer Konstellation, die sich zum Beispiel bei der Zweiseitenband-Amplitudenmodulation des Nachrichtensignals xN(t)=ANcos(2πfNtφN) mit dem Trägersignal xT(t)=ATcos(2πfTtφT) ergibt. Die Nomenklatur ist ebenfalls an diesen Fall angepasst:

  • xO(t) bezeichnet das „Obere Seitenband” mit der Amplitude AO=AN/2, der Frequenz fO=fT+fN und der Phase φO=φT+φN.
  • Entsprechend gilt für das „Untere Seitenband” xU(t) mit fU=fT+fN, AU=AO und φU=φO.


Das dazugehörige analytische Signal lautet:

x+(t)=xU+(t)+xT+(t)+xO+(t)=AUej(2πfUtφU)+ATej(2πfTtφT)+AOej(2πfOtφO).
Analytische Signal zur Zeit t=0

Im Programm dargestellt wird x+(t) als vektorielle Summe dreier Drehzeiger (alle mit positiver Drehrichtung) als violetter Punkt (siehe beispielhafte Grafik für den Startzeitpunkt t=0):

  • Der (rote) Zeiger des Trägers xT+(t) mit der Länge AT und der Nullphasenlage φT=0 dreht mit konstanter Winkelgeschwindigkeit 2πfT (eine Umdrehung in der Zeit 1/fT.
  • Der (blaue) Zeiger des Oberen Seitenbandes xO+(t) mit der Länge AO und der Nullphasenlage φO dreht mit der Winkelgeschwindigkeit 2πfO, also etwas schneller als xT+(t).
  • Der (grüne) Zeiger des Unteren Seitenbandes xU+(t) mit der Länge AU und der Nullphasenlage φU dreht mit der Winkelgeschwindigkeit 2πfU, also etwas langsamer als xT+(t).


Den zeitlichen Verlauf von x+(t) bezeichnen wir im Folgenden auch als Zeigerdiagramm. Der Zusammenhang zwischen dem physikalischen Bandpass–Signal x(t) und dem dazugehörigen analytischen Signal x+(t) lautet:

x(t)=Re[x+(t)].

Hinweis:   Die Grafik gilt für φO=+30. Daraus folgt für den Startzeitpunkt t=0 der Winkel gegenüber dem Koordinatensystem:   ϕO=φO=30. Ebenso folgt aus der Nullphanlage φU=30 des unteren Seitenbandes für den in der komplexen Ebene zu berücksichtigenden Phasenwinkel:   ϕU=+30.


Englische Beschreibung (muss noch angepasst werden)


Theoretical Background


Beschreibungsmöglichkeiten von Bandpass-Signalen

Bandpass–Spektrum X(f)

Wir betrachten hier Bandpass-Signale x(t) mit der Eigenschaft, dass deren Spektren X(f) nicht im Bereich um die Frequenz f=0 liegen, sondern um eine Trägerfrequenz fT. Meist kann auch davon ausgegangen werden, dass die Bandbreite BfT ist.

Die Grafik zeigt ein solches Bandpass–Spektrum X(f). Unter der Annahme, dass das zugehörige x(t) ein physikalisches Signal und damit reell ist, ergibt sich für die Spektralfunktion X(f) eine Symmetrie bezüglich der Frequenz f=0. Ist x(t) eine gerade Funktion   ⇒   x(t)=x(t), so ist auch X(f) reell und gerade.


Neben dem physikalischen Signal x(t)  X(f) verwendet man zur Beschreibung von Bandpass-Signalen gleichermaßen:

  • das analytische Signal x+(t)  X+(f), wie im nächsten Unterabschnitt beschrieben,
  • das äquivalente Tiefpass–Signal xTP(t)  XTP(f), siehe Applet Physikalisches Signal und Äquivalentes Tiefpass–Signal.



Analytisches Signal – Spektralfunktion

Das zum physikalischen Signal x(t) gehörige analytische Signal x+(t) ist diejenige Zeitfunktion, deren Spektrum folgende Eigenschaft erfüllt:

Konstruktion der Spektralfunktion X+(f)
X+(f)=[1+sign(f)]X(f)={2X(f)f¨urf>0,0f¨urf<0.

Die so genannte Signumfunktion ist dabei für positive Werte von f gleich +1 und für negative f–Werte gleich 1.

  • Der (beidseitige) Grenzwert liefert sign(0)=0.
  • Der Index „+” soll deutlich machen, dass X+(f) nur Anteile bei positiven Frequenzen besitzt.


Aus der Grafik erkennt man die Berechnungsvorschrift für X+(f): Das tatsächliche BP–Spektrum X(f) wird

  • bei den positiven Frequenzen verdoppelt, und
  • bei den negativen Frequenzen zu Null gesetzt.


Aufgrund der Unsymmetrie von X+(f) bezüglich der Frequenz f=0 kann man bereits jetzt schon sagen, dass die Zeitfunktion x+(t) bis auf einen trivialen Sonderfall x+(t)=0  X+(f)=0 stets komplex ist.

Analytisches Signal – Zeitverlauf

An dieser Stelle ist es erforderlich, kurz auf eine weitere Spektraltransformation einzugehen.

Definition:  Für die Hilberttransformierte H{x(t)} einer Zeitfunktion x(t) gilt:

y(t)=H{x(t)}=1π+x(τ)tτdτ.

Dieses bestimmte Integral ist nicht auf einfache, herkömmliche Art lösbar, sondern muss mit Hilfe des Cauchy–Hauptwertsatzes ausgewertet werden.

Entsprechend gilt im Frequenzbereich:

Y(f)=jsign(f)X(f).


Das obige Ergebnis lässt sich mit dieser Definition wie folgt zusammenfassen:

  • Man erhält aus dem physikalischen BP–Signal x(t) das analytische Signal x+(t), indem man zu x(t) einen Imaginärteil gemäß der Hilberttransformierten hinzufügt:
x+(t)=x(t)+jH{x(t)}.
  • H{x(t)} verschwindet nur für den Fall x(t)=const.   ⇒   Gleichsignal. Bei allen anderen Signalformen ist somit das analytische Signal x+(t) komplex.


  • Aus dem analytischen Signal x+(t) kann das physikalische Bandpass–Signal in einfacher Weise durch Realteilbildung ermittelt werden:
x(t)=Re[x+(t)].

Beispiel 1:  Das Prinzip der Hilbert–Transformation wird durch die nachfolgende Grafik nochmals verdeutlicht:

  • Nach der linken Darstellung (A) kommt man vom physikalischen Signal x(t) zum analytischen Signal x+(t), indem man einen Imaginärteil jy(t) hinzufügt.
  • Hierbei ist y(t)=H{x(t)} eine reelle Zeitfunktion, die sich im Spektralbereich durch die Multiplikation des Spektrums X(f) mit jsign(f) angeben lässt.
Zur Verdeutlichung der Hilbert–Transformierten

Die rechte Darstellung (B) ist äquivalent zu (A). Nun gilt x+(t)=x(t)+z(t) mit der rein imaginären Funktion z(t). Ein Vergleich der beiden Bilder zeigt, dass tatsächlich z(t)=jy(t) ist.



Darstellung der harmonischen Schwingung als analytisches Signal

Die Spektralfunktion X(f) einer harmonischen Schwingung x(t)=Acos(2πfTtφ) besteht bekanntlich aus zwei Diracfunktionen bei den Frequenzen

  • +fT mit dem komplexen Gewicht A/2ejφ,
  • fT mit dem komplexen Gewicht A/2e+jφ.


Somit lautet das Spektrum des analytischen Signals (also ohne die Diracfunktion bei der Frequenz f=fT, aber Verdoppelung bei f=+fT):

X+(f)=Aejφδ(ffT).

Die dazugehörige Zeitfunktion erhält man durch Anwendung des Verschiebungssatzes:

x+(t)=Aej(2πfTtφ).

Diese Gleichung beschreibt einen mit konstanter Winkelgeschwindigkeit ωT=2πfT drehenden Zeiger.

Beispiel 2:  Aus Darstellungsgründen wird das Koordinatensystem entgegen der üblichen Darstellung um 90 gedreht (Realteil nach oben, Imaginärteil nach links).

Zeigerdiagramm einer harmonischen Schwingung

Anhand dieser Grafik sind folgende Aussagen möglich:

  • Zum Startzeitpunkt t=0 liegt der Zeiger der Länge A (Signalamplitude) mit dem Winkel φ in der komplexen Ebene. Im gezeichneten Beispiel gilt φ=45.
  • Für Zeiten t>0 dreht der Zeiger mit konstanter Winkelgeschwindigkeit (Kreisfrequenz) ωT in mathematisch positiver Richtung, das heißt entgegen dem Uhrzeigersinn.
  • Die Spitze des Zeigers liegt somit stets auf einem Kreis mit Radius A und benötigt für eine Umdrehung genau die Zeit T0, also die Periodendauer der harmonischen Schwingung x(t).
  • Die Projektion des analytischen Signals x+(t) auf die reelle Achse, durch rote Punkte markiert, liefert die Augenblickswerte von x(t).



x+(t)–Darstellung einer Summe aus drei harmonischen Schwingungen

In unserem Applet setzen wir stets einen Zeigerverbund aus drei Drehzeigern voraus. Das physikalische Signal lautet:

x(t)=xU(t)+xT(t)+xO(t)=AUcos(2πfUtφU)+ATcos(2πfTtφT)+AOcos(2πfOtφO).
  • Jede der drei harmonischen Schwingungen harmonischen Schwingungen xT(t), xU(t) und xO(t) wird durch eine Amplitude (A), eine Frequenz (f) und einen Phasenwert (φ) charakterisiert.
  • Die Indizes sind an das Modulationsverfahren Zweiseitenband–Amplitudenmodulation angelehnt. „T” steht für „Träger”, „U” für „Unteres Seitenband” und „O” für „Oberes Seitenband”. Entsprechend gilt stets fU<fT und fO>fT. Für die Amplituden und Phasen gibt es keine Einschränkungen.

Das dazugehörige analytische Signal lautet:

x+(t)=xU+(t)+xT+(t)+xO+(t)=AUej(2πfUtφU)+ATej(2πfTtφT)+AOej(2πfOtφO).

Beispiel 3:  Die hier angegebene Konstellation ergibt sich zum Beispiel bei der Zweiseitenband-Amplitudenmodulation des Nachrichtensignals xN(t)=ANcos(2πfNtφN) mit dem Trägersignal xT(t)=ATcos(2πfTtφT). Hierauf wird in der Versuchsdurchführung häufiger eingegangen.


Bei dieser Betrachtungsweise gibt es einige Einschränkungen bezüglich der Programmparameter:

  • Für die Frequenzen gelte stets fO=fT+fN und fU=fTfN.
  • Ohne Verzerrungen sind die Amplitude der Seitenbänder AO=AO=AN/2.
  • Die jeweiligen Phasenverhältnisse können der nachfolgenden Grafik entnommen werden.
Spektum X+(f) des analytischen Signals für verschiedene Phasenkonstellationen


Exercises

Exercises verzerrungen.png
  • First select the task number.
  • A task description is displayed.
  • Parameter values are adjusted.
  • Solution after pressing „Hide solition”.


The number „0” will reset to the same setting as the program start and will output a text with further explanation of the applet.


In the following, Green denotes the lower sideband   ⇒   (AU,fU,φU),   Red the carrier   ⇒   (AT,fT,φT) and Blue the upper sideband   ⇒   (AO,fO,φO).

(1)   Consider and interpret the analytic signal x+(t) for Red:AT=1.5 V, fT=50 kHz, φT=0. In addition, AU=AO=0.

Which Signal values x+(t) result for t=0, t = 5 \ \rm µ s and t = 20 \ \rm µ s? How lange are the corresponding signal values of x(t)?
  For a cosine signal x_+(t= 0) = A_{\rm T} = 1.5\ \text{V}. Then x_+(t) rotates in a mathematically positive direction (one revolution per period T_0 = 1/f_{\rm T}):
  x_+(t= 20 \ {\rm µ s}) = x_+(t= 0) = 1.5\ \text{V}\hspace{0.3cm}\Rightarrow\hspace{0.3cm}x(t= 20 \ {\rm µ s}) = 1.5\ \text{V,}\hspace{0.5cm} x_+(t= 5 \ {\rm µ s}) = {\rm j} \cdot 1.5\ \text{V}\hspace{0.3cm}\Rightarrow\hspace{0.3cm}x(t= 5 \ {\rm µ s}) = {\rm Re}[x_+(t= 5 \ {\rm µ s})] = 0.


(2)   How do the ratios change for \text{Red:} \hspace{0.15cm} A_{\rm T} = 1.0\ \text{V}, \ f_{\rm T} = 100 \ \text{kHz}, \ \varphi_{\rm T} = 90^\circ?

The signal x(t) is now a sine signal with a smaller amplitude. The analytic signal now starts because of \varphi_{\rm T} = 90^\circ   ⇒   \phi_{\rm T} = -90^\circ bei x_+(t= 0) = -{\rm j} \cdot A_{\rm T}. After that, x_+(t) rotates again in a mathematically positive direction, but twice as fast because of T_0 = 10 \ \rm µ s as in \rm (1).


(3)   Now applies   \text{Red:} \hspace{0.15cm} A_{\rm T} = 1\ \text{V}, \ f_{\rm T} = 100 \ \text{kHz}, \ \varphi_{\rm T} = 0^\circ,   \text{Green:} \hspace{0.15cm} A_{\rm U} = 0.4\ \text{V}, \ f_{\rm U} = 80 \ \text{kHz}, \ \varphi_{\rm U} = 0^\circ,   \text{Blue:} \hspace{0.15cm} A_{\rm O} = 0.4\ \text{V}, \ f_{\rm O} = 120 \ \text{kHz}, \ \varphi_{\rm O} = 0^\circ.

Consider and interpret the physical signal x(t) the analytic signal x_+(t).
The Signal x(t) results in the double sideband–Amplitude modulation (ZSB–AM) of the message signals A_{\rm N}\cdot \cos\left(2\pi f_{\rm N}\cdot t\right) with A_{\rm N} = 0.8\ \text{V}, f_{\rm N} = 20\ \text{kHz}. The carrier x_{\rm T}(t) with f_{\rm T} = 100\ \text{kHz} is also cosinusoidal. The degree of modulation is m = A_{\rm N}/A_{\rm T} = 0.8 and the period T_{\rm 0} = 50\ \text{µs}.
In the phasor diagram, the (red) carrier rotates faster than the (green) lower sideband and slower than the (blue) upper sideband. The analytic signal x_+(t) results as the geometric sum of the three rotating hands. It seems that the blue pointer is leading the wearer and the green pointer is following the wearer.


(4)   The settings of task (3)continue to apply. Which signal values are obtained at t=0, t=2.5 \ \rm µ s, t= 5 \ \rm µ s and t=10 \ \rm µ s?

At time t=0, all the pointers are in the direction of the real axis, so that x(t=0) = {\rm Re}\big [x+(t= 0)\big] = A_{\rm U} + A_{\rm T} + A_{\rm O} = 1.8\ \text{V}.
Until the time t=2.5 \ \rm µ s, the red carrier has rotated by 90^\circ, the blue pointer by 108^\circ and the green by 72^\circ. We have x(t=2.5 \ \rm µ s) = {\rm Re}\big [x_+(t= 2.5 \ \rm µ s)\big] = 0, because now the pointer group points in the direction of the imaginary axis. The other sought signal values are x(t=5 \ \rm µ s) = {\rm Re}\big [x_+(t= 5 \ \rm µ s)\big] = -1.647\ \text{V} and x(t=10 \ \rm µ s) = {\rm Re}\big [x_+(t= 10 \ \rm µ s)\big] = 1.247\ \text{V}.
For x_+(t) a spiral shape results, alternating with a smaller radius and then with a larger radius.


(5)   How should the phase parameters \varphi_{\rm T}, \varphi_{\rm U} and \varphi_{\rm O} be set if both the carrier x_{\rm T}(t) and the message signal x_{\rm N}(t) are sinusoidal?

The parameter selection \varphi_{\rm T} = \varphi_{\rm U} = \varphi_{\rm O}=90^\circ describes the signals x_{\rm T}(t) = A_{\rm T}\cdot \sin\left(2\pi f_{\rm T}\cdot t\right) and x_{\rm N}(t) = A_{\rm N}\cdot \cos\left(2\pi f_{\rm N}\cdot t\right). If, in addition, the message x_{\rm N}(t) is sinusoidal, then \varphi_{\rm O}=\varphi_{\rm T} - 90^\circ = 0 and \varphi_{\rm U}=\varphi_{\rm T} + 90^\circ = 180^\circ must be set.


(6)   The settings of task (3) apply except A_{\rm T} = 0.6\ \text{V}. Which modulation method is described here?

What are the consequences of this? What changes with A_{\rm T} = 0?
It is a ZSB–AM with carrier with the modulation degree m=0.8/0.6 = 1.333. For m > 1, however, Synchrondemodulation is required. Hüllkurvendemodulation no longer works.
With A_{\rm T} = 0   ⇒   m \to \infty results in a ZSB–AM without carrier. Also for this you absolutely need the synchronous demodulation.


(7)     Now applies   \text{Red:} \hspace{0.15cm} A_{\rm T} = 1\ \text{V}, \ f_{\rm T} = 100 \ \text{kHz}, \ \varphi_{\rm T} = 0^\circ,   \text{Green:} \hspace{0.15cm} A_{\rm U} = 0,   \text{Blue:} \hspace{0.15cm} A_{\rm O} = 0.8\ \text{V}, \ f_{\rm O} = 120 \ \text{kHz}, \ \varphi_{\rm O} = 90^\circ.

Which constellation is described here? What changes with A_{\rm U} = 0.8\ \text{V} and A_{\rm O} = 0?
In both cases, it is a Einseitenbandmodulation (ESB–AM) with the modulation degree \mu = 0.8 (in ESB we denote the degree of modulation with \mu instead of m). The carrier signal is cosinusoidal and the message signal is sinusoidal.
A_{\rm O} = 0.8\ \text{V}, A_{\rm U} = 0 is an OSB–Modulation. The green pointer is missing and the blue pointer rotates faster compared to the red carrier. ::A_{\rm U} = 0.8\ \text{V}, A_{\rm O} = 0 is a USB–Modulation. The blue pointer is missing and the green pointer rotates slower compared to the red carrier. <div class="bluebox"> '''(8)'''   Now applies   \text{Red:} \hspace{0.05cm} A_{\rm T} = 1\ \text{V}, \ f_{\rm T} = 100 \ \text{kHz}, \ \varphi_{\rm T} = 0^\circ,   \text{Green:} \hspace{0.05cm} A_{\rm U} = 0.4\ \text{V}, \ f_{\rm U} = 80 \ \text{kHz}, \ \varphi_{\rm U} = -90^\circ,   \text{Blue:} \hspace{0.05cm} A_{\rm O} = 0.2\ \text{V}, \ f_{\rm O} = 120 \ \text{kHz}, \ \varphi_{\rm O} = +90^\circ. :Which constellation could be described here? Which figure is given for the equivalent low-pass –signal x_{\rm TP}(t)?   ⇒   „locus”? <div style="clear:both;"> </div> </div> ::It could be a ZSB–AM of a sinusoidal signal with cosinusoidal carrier and modulation degree m=0.8, in which the upper sideband is attenuated by a factor of 2. The equivalent low-pass –signal x_{\rm TP}(t) has an elliptical course in the complex plane. =='"`UNIQ--h-8--QINU`"'Applet Manual== [[File:Handhabung_verzerrungen.png|center]] <br>     '''(A)'''     Parametereingabe für das Eingangssignal x(t) per Slider: Amplituden, Frequenzen, Phasenwerte     '''(B)'''     Vorauswahl für die Kanalparameter: per Slider, Tiefpass oder Hochpass     '''(C)'''     Eingabe der Kanalparameter per Slider: Dämpfungsfaktoren und Phasenlaufzeiten     '''(D)'''     Eingabe der Kanalparameter für Hoch– und Tiefpass: Ordnung n, Grenzfrequenz f_0     '''(E)'''     Eingabe der Matching–Parameter k_{\rm M} und \varphi_{\rm M}     '''(F)'''     Auswahl der darzustellenden Signale: x(t), y(t), z(t), \varepsilon(t), \varepsilon^2(t)     '''(G)'''     Graphische Darstellung der Signale     '''(H)'''     Eingabe der Zeit t_* für die Numerikausgabe     '''( I )'''     Numerikausgabe der Signalwerte x(t_*), y(t_*), z(t_*) und \varepsilon(t_*)     '''(J)'''     Numerikausgabe des Hauptergebnisses P_\varepsilon     '''(K)'''     Abspeichern und Zurückholen von Parametersätzen     '''(L)'''     Bereich für die Versuchsdurchführung: Aufgabenauswahl, Aufgabenstellung und Musterlösung     '''(M)'''     Variationsmöglichkeiten für die grafische Darstellung \hspace{1.5cm}Zoom–Funktionen „+” (Vergrößern), „-” (Verkleinern) und \rm o (Zurücksetzen) \hspace{1.5cm}Verschieben mit „\leftarrow” (Ausschnitt nach links, Ordinate nach rechts), „\uparrow” „\downarrow” und „\rightarrow\hspace{1.5cm}'''Andere Möglichkeiten''': \hspace{1.5cm}Gedrückte Shifttaste und Scrollen: Zoomen im Koordinatensystem, \hspace{1.5cm}Gedrückte Shifttaste und linke Maustaste: Verschieben des Koordinatensystems.

About the Authors

Dieses interaktive Berechnungstool wurde am Lehrstuhl für Nachrichtentechnik der Technischen Universität München konzipiert und realisiert.

  • Die erste Version wurde 2005 von Ji Li im Rahmen ihrer Diplomarbeit mit „FlashMX–Actionscript” erstellt (Betreuer: Günter Söder).
  • 2018 wurde dieses Programm von Xiaohan Liu im Rahmen ihrer Bachelorarbeit (Betreuer: Tasnád Kernetzky) neu gestaltet und erweitert.

Once again: Open Applet in new Tab

Open Applet in a new tab