Difference between revisions of "Applets:Linear Distortions of Periodic Signals"
Line 34: | Line 34: | ||
*the Distortion Power PD and | *the Distortion Power PD and | ||
*the Signal–to–Distortion Ratio (abbr.: SDR) | *the Signal–to–Distortion Ratio (abbr.: SDR) | ||
− | :$$\rho_{\rm D}=\frac{\rm Signal | + | :$$\rho_{\rm D}=\frac{\rm Signal Power}{\rm Distortion Power} = \frac{P_x}{P_{\rm D} }.$$ |
Line 44: | Line 44: | ||
[[File:LZI_T_2_2_S3_vers2.png|center|frame|For clarification of nonlinear distortions |class=fit]] | [[File:LZI_T_2_2_S3_vers2.png|center|frame|For clarification of nonlinear distortions |class=fit]] | ||
− | [[File:P_ID899__LZI_T_2_3_S1_neu.png|right |frame| | + | [[File:P_ID899__LZI_T_2_3_S1_neu.png|right |frame| Description of a linear system|class=fit]] |
− | *''' | + | *'''Linear distortions'''occur, if the transmission channel is characterized by a frequency response H(f)≠const.. Various frequencies are attenuated and delayed differently. Characteristic of this is that although frequencies can disappear (for example, through a low pass or a high pass), no new frequencies can arise. |
− | In | + | In this applet only linear distortions are considered. |
− | === | + | === Description forms for the frequency response === |
<br> | <br> | ||
− | + | The generally complex valued frequency response can be represented as follows: | |
:$$H(f) = |H(f)| \cdot {\rm e}^{-{\rm j} \hspace{0.05cm} \cdot | :$$H(f) = |H(f)| \cdot {\rm e}^{-{\rm j} \hspace{0.05cm} \cdot | ||
\hspace{0.05cm} b(f)} = {\rm e}^{-a(f)}\cdot {\rm e}^{-{\rm j} | \hspace{0.05cm} b(f)} = {\rm e}^{-a(f)}\cdot {\rm e}^{-{\rm j} | ||
\hspace{0.05cm} \cdot \hspace{0.05cm} b(f)}.$$ | \hspace{0.05cm} \cdot \hspace{0.05cm} b(f)}.$$ | ||
− | + | This results in the following description variables: | |
− | * | + | *The absolute value |H(f)| is called '''Amplitude response''' and in logarithmic form '''Attenuation curve''' bezeichnet: |
:$$a(f) = - \ln |H(f)|\hspace{0.2cm}{\rm in \hspace{0.1cm}Neper | :$$a(f) = - \ln |H(f)|\hspace{0.2cm}{\rm in \hspace{0.1cm}Neper | ||
\hspace{0.1cm}(Np) } = - 20 \cdot \lg |H(f)|\hspace{0.2cm}{\rm in | \hspace{0.1cm}(Np) } = - 20 \cdot \lg |H(f)|\hspace{0.2cm}{\rm in | ||
− | \hspace{0.1cm} | + | \hspace{0.1cm}Decibel \hspace{0.1cm}(dB) }.$$ |
− | * | + | *The '''Phase response''' b(f) indicates the negative frequency–dependent angle of H(f) in the complex plane based on the real axis: |
:$$b(f) = - {\rm arc} \hspace{0.1cm}H(f) \hspace{0.2cm}{\rm in | :$$b(f) = - {\rm arc} \hspace{0.1cm}H(f) \hspace{0.2cm}{\rm in | ||
\hspace{0.1cm}Radian \hspace{0.1cm}(rad)}.$$ | \hspace{0.1cm}Radian \hspace{0.1cm}(rad)}.$$ | ||
− | === | + | === <i>N^{th}</i> grade low pass === |
<br> | <br> | ||
− | + | The frequency response of a realizable <i>N^{th}</i> grade low pass is: | |
:H(f)=[11+j⋅f/f0]N. | :H(f)=[11+j⋅f/f0]N. | ||
− | + | An example is the RC–low pass with N=1. Consequently we can obtain | |
− | * | + | *the attenuation curve: |
:a(f)=N/2⋅ln[1+(f/f0)2], | :a(f)=N/2⋅ln[1+(f/f0)2], | ||
− | * | + | *the phase curve: |
:b(f)=N⋅arctan(f/f0), | :b(f)=N⋅arctan(f/f0), | ||
− | * | + | *the attenuation factor for the frequency f=fi: |
:αi=|H(f=fi)|=[1+(f/f0)2]N/2⇒x(t)=Ai⋅cos(2πfit)→y(t)=αi⋅Ai⋅cos(2πfit), | :αi=|H(f=fi)|=[1+(f/f0)2]N/2⇒x(t)=Ai⋅cos(2πfit)→y(t)=αi⋅Ai⋅cos(2πfit), | ||
− | * | + | *the phase delay for the frequency f=fi: |
:τi=b(fi)2πfi=N⋅arctan(fi/f0)2πfi⇒x(t)=Ai⋅cos(2πfit)→y(t)=Ai⋅cos(2πfi(t−τi)), | :τi=b(fi)2πfi=N⋅arctan(fi/f0)2πfi⇒x(t)=Ai⋅cos(2πfit)→y(t)=Ai⋅cos(2πfi(t−τi)), | ||
− | [[File:Tiefpass.png|center|frame| | + | [[File:Tiefpass.png|center|frame|Attenuation and phase curve of an <i>N^{th}</i> grade low pass]] |
=== Hochpass <i>N</i>–Ordnung === | === Hochpass <i>N</i>–Ordnung === |
Revision as of 23:36, 14 January 2018
Contents
Applet description
This applet showcases the effects of linear distortions(attenuation distortions and phase distortions) with
- the input signal x(t) ⇒ Power Px:
- x(t)=x1(t)+x2(t)=A1⋅cos(2πf1⋅t−φ1)+A2⋅cos(2πf2⋅t−φ2),
- the output signal y(t) ⇒ output signal Py:
- y(t)=α1⋅x1(t−τ1)+α2⋅x2(t−τ2),
- the matching output signal z(t) ⇒ Power Pz:
- z(t)=kM⋅y(t−τM)+α2⋅x2(t−τ2),
- the differential signal ε(t)=z(t)−x(t) ⇒ Power Pε.
The adjustment of the output signal's amplitude and phase y(t) ⇒ „Matching” allows for a differentiation between
- attenuation distortion and frequency–independant attenuation, as well as
- phase distortion and pure running time.
The Distortion Power PD is used to measure the strength of the linear distortion.
Theoretical background
Distortions refer to generally unwanted alterations of a message signal through a transmission system. Together with the strong stochastic effects (noise, crosstalk, etc.), they are a crucial limitation for the quality and rate of transmission.
Just as the „Stärke” of noise can be assessed through
- the Noise Power PN and
- the Signal–to–Noise Ratio (abbr.: SNR) ρN,
distortions can be quantified through
- the Distortion Power PD and
- the Signal–to–Distortion Ratio (abbr.: SDR)
- ρD=SignalPowerDistortionPower=PxPD.
Linear and nonlinear distortions
A distinction is made between linear and nonlinear distortions:
- Nonlinear distortions occur, if at all times t the nonlinear correlation y=g(x)≠const.⋅x exists between the signal values x=x(t) at the input and y=y(t) at the output, whereby y=g(x) is defined as the system's nonlinear characteristic. By creating a cosine signal at the input with frequency f0 the output signal value includes f0 as well as multiple harmonic waves. We conclude that new frequencies arise through nonlinear distortion.
- Linear distortionsoccur, if the transmission channel is characterized by a frequency response H(f)≠const.. Various frequencies are attenuated and delayed differently. Characteristic of this is that although frequencies can disappear (for example, through a low pass or a high pass), no new frequencies can arise.
In this applet only linear distortions are considered.
Description forms for the frequency response
The generally complex valued frequency response can be represented as follows:
- H(f)=|H(f)|⋅e−j⋅b(f)=e−a(f)⋅e−j⋅b(f).
This results in the following description variables:
- The absolute value |H(f)| is called Amplitude response and in logarithmic form Attenuation curve bezeichnet:
- a(f)=−ln|H(f)|inNeper(Np)=−20⋅lg|H(f)|inDecibel(dB).
- The Phase response b(f) indicates the negative frequency–dependent angle of H(f) in the complex plane based on the real axis:
- b(f)=−arcH(f)inRadian(rad).
N^{th} grade low pass
The frequency response of a realizable N^{th} grade low pass is:
- H(f)=[11+j⋅f/f0]N.
An example is the RC–low pass with N=1. Consequently we can obtain
- the attenuation curve:
- a(f)=N/2⋅ln[1+(f/f0)2],
- the phase curve:
- b(f)=N⋅arctan(f/f0),
- the attenuation factor for the frequency f=fi:
- αi=|H(f=fi)|=[1+(f/f0)2]N/2⇒x(t)=Ai⋅cos(2πfit)→y(t)=αi⋅Ai⋅cos(2πfit),
- the phase delay for the frequency f=fi:
- τi=b(fi)2πfi=N⋅arctan(fi/f0)2πfi⇒x(t)=Ai⋅cos(2πfit)→y(t)=Ai⋅cos(2πfi(t−τi)),
Hochpass N–Ordnung
Der Frequenzgang eines realisierbaren Hochpasses N–Ordnung lautet:
- H(f)=[j⋅f/f01+j⋅f/f0]N.
Ein einfacher LC–Tiefpass hat diesen Verlauf mit N=1. Damit erhält man
- den Dämpfungsverlauf:
- a(f)=N/2⋅ln[1+(f0/f)2],
- den Phasenverlauf:
- b(f)=−N⋅arctan(f0/f),
- den Dämpfungsfaktor für die Frequenz f=fi:
- αi=|H(f=fi)|=[1+(f0/f)2]N/2⇒x(t)=Ai⋅cos(2πfit)→y(t)=αi⋅Ai⋅cos(2πfit),
- die Phasenlaufzeit für die Frequenz f=fi:
- τi=b(fi)2\pifi=−N⋅arctan(f0/fi)2\pifi⇒x(t)=Ai⋅cos(2πfit)→y(t)=Ai⋅cos(2πfi(t−τi)),
Laufzeiten bei Tiefpass und Hochpass
Mache ich noch
Dämpfungsverzerrungen und Phasenverzerrungen
Überarbeite ich noch
Lineare Verzerrungen treten üblicherweise in Form von
- Dämpfungsverzerrungen αi und
- Phasenverzerrungen τi auf.
Ist α1≠α2 und τ1=τ2, so liegen ausschließlich Dämpfungsverzerrungen vor.
Dagegen führt α1=α2 und τ1≠τ2 zu reinen Phasenverzerrungen.
Ein Signal y(t) ist gegenüber x(t) unverzerrt, wenn α1=α2 und τ1undτ2 gilt.
Vorschlag für die Versuchsdurchführung
BlaBla
(1) Für das Sendesignal x(t) gelte A1=0.8 V, A2=0.6 V, f1=0.5 kHz, f2=1.5 kHz, φ1=90∘, φ2=0∘.
- Wie groß ist die Periodendauer T0? Welche Leistung Px weist dieses Signal auf? Wo können Sie diesen Wert im Programm ablesen?
⇒T0=[ größter gemeinsamer Teiler (0.5 kHz, 1.5 kHz)]−1=2.0 ms_;
Px=A21/2+A22/2=0.5 V2_=Pε, wenn kM=0_ ⇒ z(t)≡0.
(2) Variieren Sie bei sonst gleicher Einstellung wie unter (1) die Phase φ2 im gesamten möglichen Bereich ±180∘. Wie ändern sich T0 und Px?
⇒Keine Veränderungen:T0=2.0 ms;Px=0.5 V2_.
(3) Variieren Sie bei sonst gleicher Einstellung wie unter (1) die Frequenz f2 im Bereich 0≤f2≤5 kHz. Wie ändert sich die Signalleistung Px?
⇒Keine Veränderungen, falls f2≠0 oder f2≠f1:Px=0.5 V2_.T0 ändert sich, falls f2 kein Vielfaches von f1.
Falls f2=0:Px=A21/2+A22=0.68 V2_. Allgemeine Formel noch überprüfen
Falls f2=f1:Px=[A1cos(φ1)+A2cos(φ2)]2/2+[A1sin(φ1)+A2sin(φ2)]2/2. Mit φ1=90∘, φ2=0∘:Px=0.5 V2_.
(4) Ausgehend vom bisherigen Sendesignal x(t) gelte für den Kanal: α1=α2=0.5, τ1=τ2=0.5 ms. Zudem sei kM=1 und τM=0 .
- Gibt es lineare Verzerrungen? Wie groß ist die Empfangsleistung Py und die Leistung Pε des Differenzsignals ε(t)=z(t)−x(t)?
⇒y(t)=0.5⋅x(t−1 ms)_ ist unverzerrt, nur gedämpft und verzögert.
Empfangsleistung:Py=(A1/2)2/2+(A2/2)2/2=0.125 V2_. Pε ist deutlich größer:Pε=0.625 V2_.
(5) Variieren Sie bei sonst gleicher Einstellung wie unter (4) die Matchingparameter kM und τM. Wie groß ist die Verzerrungsleistung PD?
⇒PD ist gleich der Leistung Pε des Differenzsignals bei bestmöglicher Anpassung:kM=2 und τM=T0−0.5 ms=1.5 ms
⇒z(t)=x(t)⇒ε(t)=0⇒PD=Pε=0_⇒weder Dämpfungs- noch Phasenverzerrungen.
(6) Für den Kanal gelte nun α1=0.5,α2=0.2_, τ1=τ2=0.5 ms. Wie groß sind nun die Verzerrungsleistung PD und das SDR ρD?
⇒PD=Pε bei bestmöglicher Anpassung:kM=2.24_ und τM=1.5 ms_:PD=0.059 V2_.
Nur Dämpfungsverzerrungen.Signal-zu-Verzerrung-Leistungsverhältnis ρD=Px/Pε≈8.5_.
(7) Für den Kanal gelte nun α1=α2=0.5, τ1=2 ms_, τ2=0.5 ms. Wie groß sind nun PD und ρD?
⇒PD=Pε bei bestmöglicher Anpassung:kM=1.82_ und τM=0.15 ms_:PD=0.072 V2_.
Nur Phasenverzerrungen.Signal-zu-Verzerrung-Leistungsverhältnis ρD=Px/Pε≈7_.
(8) Für den Kanal gelte nun α1=0.5_,α2=0.2_, τ1=2 ms_, τ2=0.5 ms_. Wie groß sind nun PD und ρD? Wie lässt sich y(t) annähern?
⇒Dämpfungs- und Phasenverzerrungen. Bestmögliche Anpassung:kM=2.06_, τM=0.15 ms_:PD=0.136 V2_,ρD≈3.7_.
Zusammenfassen von φ- und τ-Parameter: y(t)=0.4 V⋅sin (2πf1t)−0.12 V⋅sin (2π⋅3f1⋅t)≈0.52 V⋅sin3(2πf1t)_.
(9) Nun gelte A1=A2=1 V, f1=1 kHz, f2=1 kHz, φ1=0∘, φ2=0∘_. Der Kanal sei ein Tiefpass erster Ordnung (f0=1 kHz)_.
- Gibt es Dämpfungsverzerrungen? Gibt es Phasenverzerrungen? Wie groß ist nun die Verzerrungsleistung PD?
⇒Dämpfungsverzerrungen, da α1=0.71≠α2=0.45_; geringere Phasenverzerrungen, da τ1=0.13 ms≈τ2=0.09 ms_.
Verzerrungsleistung PD=0.074 V2_ bei bestmöglicher Anpassung:kM=1.6_ und τM=0.9 ms_.
(10) Wie ändern sich die Kanalparameter durch einen Tiefpass zweiter Ordnung gegenüber einem Tiefpass erster Ordnung (f0=1 kHz).
- Wie groß ist nun die Verzerrungsleistung PD? Wie groß ist nun die Verzerrungsleistung PD?
⇒Es gilt α1=0.712≈0.5,α2=0.452≈0.5,τ1=2⋅0.13≈0.25 msτ2=2⋅0.09 ms≈0.18 ms_.
PD=0.228 V2 ist größer und der 2 kHz-Anteil wird im Vergleich zum 2 kHz-Anteil noch mehr unterdrückt.
(11) Welche Unterschiede ergeben sich bei einem Hochpass zweiter Ordnung gegenüber einem Tiefpass zweiter Ordnung (f0=1 kHz).
⇒???????????????
Zur Handhabung des Applets
(A) Parametereingabe per Slider
(B) Bereich der graphischen Darstellung
(C) Variationsmöglichkeit für die graphische Darstellung
(D) Abspeichern und Zurückholen von Parametersätzen
(E) Numerikausgabe des Hauptergebnisses T0; graphische Verdeutlichung durch rote Linie
(F) Ausgabe von xmax und der Signalwerte x(t∗)=x(t∗+T0)=x(t∗+2T0)
(G) Darstellung der Signalwerte x(t∗)=x(t∗+T0)=x(t∗+2T0) durch grüne Punkte
(H) Einstellung der Zeit t∗ für die Signalwerte x(t∗)=x(t∗+T0)=x(t∗+2T0)
Details zum obigen Punkt (C)
(*) Zoom–Funktionen „+” (Vergrößern), „−” (Verkleinern) und o (Zurücksetzen)
(*) Verschieben mit „←” (Ausschnitt nach links, Ordinate nach rechts), „↑” „↓” und „→”
Andere Möglichkeiten:
(*) Gedrückte Shifttaste und Scrollen: Zoomen im Koordinatensystem,
(*) Gedrückte Shifttaste und linke Maustaste: Verschieben des Koordinatensystems.
Über die Autoren
Dieses interaktive Berechnungstool wurde am Lehrstuhl für Nachrichtentechnik der Technischen Universität München konzipiert und realisiert.
- Die erste Version wurde 2005 von Bettina Hirner im Rahmen ihrer Diplomarbeit mit „FlashMX–Actionscript” erstellt (Betreuer: Günter Söder ).
- 2018 wurde dieses Programm von Jimmy He im Rahmen seiner Bachelorarbeit (Betreuer: Tasnád Kernetzky) auf „HTML5” umgesetzt und neu gestaltet.