Applets:Physical Signal & Analytic Signal
Contents
Applet Description
This applet shows the relationship between the physical bandpass–signal x(t) and the associated analytic signal x+(t). The starting point is always a bandpass signal–signal x(t) with frequency-discrete spectrum X(f):
- x(t)=xU(t)+xT(t)+xO(t)=AU⋅cos(2πfU⋅t−φU)+AT⋅cos(2πfT⋅t−φT)+AO⋅cos(2πfO⋅t−φO).
The physical signal x(t) is thus composed of three harmonic oscillations, a constellation that can be found, for example, in the Double-sideband amplitude modulation of the message signal xN(t)=AN⋅cos(2πfN⋅t−φN) with the carrier signal xT(t)=AT⋅cos(2πfT⋅t−φT) returns. The nomenclature is also adapted to this case:
- xO(t) denotes the „upper sideband” with the amplitude AO=AN/2, the frequency fO=fT+fN and the phase φO=φT+φN.
- Similarly, for the „lower sideband ” xU(t) with fU=fT+fN, AU=AO and φU=−φO.
The associated analytical signal is:
- x+(t)=xU+(t)+xT+(t)+xO+(t)=AU⋅e−j⋅(2π⋅fU⋅t−φU)+AT⋅e−j⋅(2π⋅fT⋅t−φT)+AO⋅e−j⋅(2π⋅fO⋅t−φO).
The program displays x+(t) as the vectorial sum of three rotation hands (all with positive rotation) as a purple dot (see example graphic for start time t=0):
- The (red) pointer of the carrier xT+(t) with the length AT and the zero phase position φT=0 rotates at constant angular velocity 2π⋅fT (one turn in time 1/fT.
- The (blue) pointer of the upper sideband xO+(t) with the length AO and the zero phase position φO rotates at the angular velocity 2π⋅fO, which is slightly faster than xT+(t).
- The (green) pointer of the lower sideband xU+(t) with the length AU and the zero phase position φU rotates at the angular velocity 2π⋅fU, which is slightly faster than xT+(t).
The time course of x+(t) is also referred to below as Pointer Diagram. The relationship between the physical bandpass signal x(t) and the associated analytic signal x+(t) is:
- x(t)=Re[x+(t)].
Note: The graphic applies to φO=+30∘. From this follows for the start time t=0 the angle with respect to the coordinate system: ϕO=−φO=−30∘. Similarly, from the null phantom φU=−30∘ of the lower sideband follows for the phase angle to be considered in the complex plane: ϕU=+30∘.
German description (muss noch angepasst werden)
Theoretical Background
Description possibilities of bandpass signals
We consider here bandpass signals x(t) with the property that their spectra X(f) are not in the range around the frequency f=0, but by a carrier frequency fT. In most cases it can also be assumed that the bandwidth is B≪fT.
The graph shows such a bandpass spectrum X(f). Assuming that the associated x(t) is a physical signal and thus real, the spectral function X(f) has a symmetry with respect to the frequency f=0. X(t) is an even function ⇒ x(−t)=x(t), so X(f) is real and even.
Beside the physical signal x(t) ∘−−−∙ X(f) we use for the description of Bandpass signals alike:
- the analytic signal x+(t) ∘−−−∙ X+(f), as in next subsection,
- the equivalent low-pass –signal xTP(t) ∘−−−∙ XTP(f), see Applet Physical signal and equivalent low pass –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:
- X+(f)=[1+sign(f)]⋅X(f)={2⋅X(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)=−j⋅sign(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)+j⋅H{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 j⋅y(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 −j⋅sign(f) angeben lässt.
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)=j⋅y(t) ist.
Darstellung der harmonischen Schwingung als analytisches Signal
Die Spektralfunktion X(f) einer harmonischen Schwingung x(t)=A⋅cos(2πfT⋅t−φ) besteht bekanntlich aus zwei Diracfunktionen bei den Frequenzen
- +fT mit dem komplexen Gewicht A/2⋅e−jφ,
- −fT mit dem komplexen Gewicht A/2⋅e+jφ.
Somit lautet das Spektrum des analytischen Signals (also ohne die Diracfunktion bei der Frequenz f=−fT, aber Verdoppelung bei f=+fT):
- X+(f)=A⋅e−jφ⋅δ(f−fT).
Die dazugehörige Zeitfunktion erhält man durch Anwendung des Verschiebungssatzes:
- x+(t)=A⋅ej⋅(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).
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)=AU⋅cos(2πfU⋅t−φU)+AT⋅cos(2πfT⋅t−φT)+AO⋅cos(2πfO⋅t−φ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)=AU⋅e−j⋅(2π⋅fU⋅t−φU)+AT⋅e−j⋅(2π⋅fT⋅t−φT)+AO⋅e−j⋅(2π⋅fO⋅t−φO).
Beispiel 3: Die hier angegebene Konstellation ergibt sich zum Beispiel bei der Zweiseitenband-Amplitudenmodulation des Nachrichtensignals xN(t)=AN⋅cos(2πfN⋅t−φN) mit dem Trägersignal xT(t)=AT⋅cos(2πfT⋅t−φ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=fT−fN.
- Ohne Verzerrungen sind die Amplitude der Seitenbänder AO=AO=AN/2.
- Die jeweiligen Phasenverhältnisse können der nachfolgenden Grafik entnommen werden.
Exercises
- 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, synchronous demodulation is required. envelope detection 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} und A_{\rm O} = 0?
- In both cases, it is a Single sideband (ESB–AM) with the modulation degree \mu = 0.8 (in ESB we denote the degree of modulation with \mu instead m). he 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.
(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 lowpass–signal x_{\rm TP}(t)? ⇒ „locus”?
- 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 lowpass–signal x_{\rm TP}(t) has an elliptical course in the complex plane.
Applet Manual
(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
This interactive calculation was designed and realized at the Lehrstuhl für Nachrichtentechnik of the Technischen Universität München .
- The original version was created in 2005 by Ji Li as part of her Diploma thesis using „FlashMX–Actionscript” (Supervisor: Günter Söder).
- In 2018 this Applet was redesigned and updated to „HTML5” by Xiaohan Liu as part of her Bachelor's thesis (Supervisor: Tasnád Kernetzky).