Difference between revisions of "Digital Signal Transmission/Linear Nyquist Equalization"
Line 70: | Line 70: | ||
Für den Detektionsgrundimpuls <i>g<sub>d</sub></i>(<i>t</i>) am Filterausgang ergeben sich demzufolge zu den Zeitpunkten <i>ν</i><i>T</i> mit den Abkürzungen <i>g</i><sub>0</sub> = <i>g<sub>d</sub></i>(<i>t</i> = 0), <i>g</i><sub>1</sub> = <i>g<sub>d</sub></i>(<i>t</i> = ±<i>T</i>), <i>g</i><sub>2</sub> = <i>g<sub>d</sub></i>(<i>t</i> = ±2<i>T</i>) folgende Werte: | Für den Detektionsgrundimpuls <i>g<sub>d</sub></i>(<i>t</i>) am Filterausgang ergeben sich demzufolge zu den Zeitpunkten <i>ν</i><i>T</i> mit den Abkürzungen <i>g</i><sub>0</sub> = <i>g<sub>d</sub></i>(<i>t</i> = 0), <i>g</i><sub>1</sub> = <i>g<sub>d</sub></i>(<i>t</i> = ±<i>T</i>), <i>g</i><sub>2</sub> = <i>g<sub>d</sub></i>(<i>t</i> = ±2<i>T</i>) folgende Werte: | ||
− | :<math>t = 0\hspace{-0.1cm}:\hspace{0.2cm}g_0 = k_0 \cdot g_m(0) + k_1 \cdot 2 | + | :<math> t = 0\hspace{-0.1cm}:\hspace{0.2cm}g_0 = k_0 \cdot g_m(0) + k_1 \cdot 2 |
− | \cdot g_m(1) \hspace{1.23cm}+k_2 \cdot 2 \cdot g_m(2),\hspace{0.05cm}</math> | + | \cdot g_m(1) \hspace{1.23cm}+k_2 \cdot 2 \cdot g_m(2),\hspace{0.05cm} </math> |
− | :<math>t = \pm T\hspace{-0.1cm}:\hspace{0.2cm}g_1 = k_0 \cdot g_m(1) + k_1 | + | :<math> t = \pm T\hspace{-0.1cm}:\hspace{0.2cm}g_1 = k_0 \cdot g_m(1) + k_1 |
− | \cdot [g_m(0)+g_m(2)]+ k_2 \cdot [g_m(1)+g_m(3)],</math> | + | \cdot [g_m(0)+g_m(2)]+ k_2 \cdot [g_m(1)+g_m(3)], </math> |
− | :<math>t = \pm 2T\hspace{-0.1cm}:\hspace{0.2cm}g_2 = k_0 \cdot g_m(2) + k_1 | + | :<math> t = \pm 2T\hspace{-0.1cm}:\hspace{0.2cm}g_2 = k_0 \cdot g_m(2) + k_1 |
\cdot [g_m(1)+g_m(3)]+ k_2 \cdot [g_m(2)+g_m(4)] | \cdot [g_m(1)+g_m(3)]+ k_2 \cdot [g_m(2)+g_m(4)] | ||
− | \hspace{0.05cm}.</math> | + | \hspace{0.05cm}. </math> |
Aus diesem System mit drei linear unabhängigen Gleichungen kann man nun die Filterkoeffizienten <i>k</i><sub>0</sub>, <i>k</i><sub>1</sub> und <i>k</i><sub>2</sub> so bestimmen, dass der Detektionsgrundimpuls <i>g<sub>d</sub></i>(<i>t</i>) durch die normierten Stützstellen | Aus diesem System mit drei linear unabhängigen Gleichungen kann man nun die Filterkoeffizienten <i>k</i><sub>0</sub>, <i>k</i><sub>1</sub> und <i>k</i><sub>2</sub> so bestimmen, dass der Detektionsgrundimpuls <i>g<sub>d</sub></i>(<i>t</i>) durch die normierten Stützstellen | ||
Line 86: | Line 86: | ||
vollständig gegeben ist. Auf der nächsten Seite wird die Optimierung der Filterkoeffizienten an einem einfachen Beispiel verdeutlicht.<br> | vollständig gegeben ist. Auf der nächsten Seite wird die Optimierung der Filterkoeffizienten an einem einfachen Beispiel verdeutlicht.<br> | ||
+ | == Wirkungsweise des Transversalfilters (2) == | ||
+ | <br> | ||
+ | {{Beispiel}}''':''' Wir gehen von dem symmetrischen Eingangssignal entsprechend dem oberen Diagramm aus. Mit der Abkürzung <i>g<sub>m</sub></i>(<i>ν</i>) = <i>g<sub>m</sub></i>(± <i>ν</i> · <i>T</i>) gibt es folgende Abtastwerte im Abstand der Symboldauer <i>T</i>: | ||
+ | :<math>g_m(t) = {\rm exp }\left ( - \sqrt{2 \cdot |t|/T}\right )</math> | ||
+ | :<math>\Rightarrow \hspace{0.3cm} g_m(0) = 1 ,\hspace{0.15cm}g_m(1)= | ||
+ | 0.243,\hspace{0.15cm}g_m(2)= 0.135,\hspace{0.15cm}g_m(3)= 0.086, | ||
+ | \hspace{0.15cm}g_m(4)= 0.059 \hspace{0.05cm}.</math> | ||
+ | |||
+ | Für den Ausgangsimpuls soll <i>g<sub>d</sub></i>(0) = 1 und <i>g<sub>d</sub></i>(±<i>T</i>) = 0 gelten. Hierzu eignet sich ein Laufzeitfilter erster Ordnung mit den Koeffizienten <i>k</i><sub>0</sub> und <i>k</i><sub>1</sub>, die folgende Bedingungen erfüllen müssen: | ||
+ | |||
+ | :<math>t = \pm T\hspace{-0.1cm} : \hspace{0.2cm}g_1 = k_0 \cdot 0.243 + k_1 \cdot | ||
+ | [1.000 +0.135] = 0\hspace{0.3cm}\Rightarrow | ||
+ | \hspace{0.3cm}{k_1} = | ||
+ | -0.214 \cdot {k_0}\hspace{0.05cm},</math> | ||
+ | :<math> t = 0 \hspace{-0.1cm} : \hspace{0.2cm}g_0 = k_0 \cdot 1.000 + k_1 \cdot 2 \cdot | ||
+ | 0.243= 1\hspace{0.3cm}\Rightarrow \hspace{0.3cm}0.896 \cdot {k_0} | ||
+ | = 1 \hspace{0.05cm}.</math> | ||
+ | |||
+ | [[File:P ID1425 Dig T 3 5 S2b version1.png|rechts|Eingangs- und Ausgangsimpuls des optimalen Nyquistentzerrers]]<br> | ||
+ | |||
+ | Daraus erhält man die optimalen Filterkoeffizienten <nobr><i>k</i><sub>0</sub> = 1.116</nobr> und <i>k</i><sub>1</sub> = 0.239. Das mittlere Diagramm zeigt, dass damit der erste Vorläufer und der erste Nachläufer kompensiert werden können und zugleich <i>g<sub>d</sub></i>(0) = 1 gilt (gelbe Hinterlegung). Die weiteren Detektionsgrundimpulswerte (blaue Kreise) sind aber von 0 verschieden und bewirken Impulsinterferenzen.<br><br> | ||
+ | |||
+ | Das untere Diagramm zeigt, dass mit einem Filter zweiter Ordnung (<i>N</i> = 2) Nulldurchgänge bei ±<i>T</i> und bei ±2<i>T</i> erzwungen werden, wenn die Koeffizienten <i>k</i><sub>0</sub> = 1.127, <i>k</i><sub>1</sub> = 0.219 und <i>k</i><sub>2</sub> = 0.075 geeignet gewählt sind. Das Gleichungssystem zur Bestimmung der optimalen Koeffizienten lautet dabei: | ||
+ | |||
+ | :<math>t = 0\hspace{-0.1cm}:\hspace{0.2cm}g_0 = k_0 \cdot 1.000 + k_1 \cdot 2 | ||
+ | \cdot 0.243 + k_2 \cdot 2 \cdot 0.135 = 1\hspace{0.05cm},\\ | ||
+ | t= \pm T\hspace{-0.1cm}:\hspace{0.2cm}g_1 = k_0 \cdot 0.243 + k_1 \cdot | ||
+ | [1.000+0.135]+ k_2 \cdot [0.243+0.086] = 0\hspace{0.05cm}, | ||
+ | \\ | ||
+ | t = \pm 2 T\hspace{-0.1cm}:\hspace{0.2cm}g_2 = k_0 \cdot 0.135 + k_1 \cdot | ||
+ | [0.243+0.086]+ k_2 \cdot [1.000 + 0.059]= 0 \hspace{0.05cm}.</math>{{end}}<br> | ||
+ | |||
+ | Die Ergebnisse können wie folgt verallgemeinert werden: | ||
+ | *Mit einem Laufzeitfilter <i>N</i>–ter Ordnung können der Hauptwert <i>g<sub>d</sub></i>(0) zu 1 (normiert) sowie die ersten <i>N</i> Nachläufer und die ersten <i>N</i> Vorläufer zu Null gemacht werden.<br> | ||
+ | |||
+ | *Weitere Vor– und Nachläufer (|<i>ν</i>| > <i>N</i>) lassen sich so nicht kompensieren. Es ist auch möglich, dass diese außerhalb des Kompensationsbereichs vergrößert werden oder sogar neu entstehen.<br> | ||
+ | |||
+ | *Im Grenzübergang <i>N</i> → ∞ (in der Praxis heißt das: ein Filter mit sehr vielen Koeffizienten) ist eine vollständige Nyquistentzerrung und damit eine impulsinterferenzfreie Übertragung möglich.<br> | ||
{{Display}} | {{Display}} |
Revision as of 22:00, 26 December 2016
Struktur des optimalen Nyquistentzerrers
In diesem Abschnitt gehen wir von folgendem Blockschaltbild eines Binärsystems aus.
Hierzu ist anzumerken:
- Die Diracquelle liefert die zu übertragende Nachricht (Amplitudenkoeffizienten aν) in binärer bipolarer Form. Sie wird als redundanzfrei vorausgesetzt.
- Die Sendeimpulsform gs(t) wird durch den Senderfrequenzgang HS(f) berücksichtigt. Bei allen Beispielen ist HS(f) = si(π f T) zugrunde gelegt.
- Bei manchen Herleitungen werden Sender und Kanal – hierfür wird meist ein Koaxialkabel angenommen – durch den gemeinsamen Frequenzgang HSK(f) = HS(f) · HK(f) zusammengefasst.
- Das Empfangsfilter HE(f) setzt sich multiplikativ aus dem Matched–Filter HMF(f) = HSK∗(f) und dem Transversalfilter HTF(f) zusammen, zumindest kann es gedanklich so aufgespalten werden.
- Der Gesamtfrequenzgang zwischen der Diracquelle und dem Schwellenwertentscheider soll die erste Nyquistbedingung erfüllen. Es muss also gelten:
- \[H_{\rm S}(f) \cdot H_{\rm K}(f) \cdot H_{\rm MF}(f) \cdot H_{\rm TF}(f) = H_{\rm Nyq}(f) \hspace{0.05cm}.\]
- Mit dieser Bedingung ergibt sich die maximale Augenöffnung (keine Impulsinterferenzen). Deshalb gelten für das Detektions–SNR und den Systemwirkungsgrad bei binärer Signalisierung:
- \[\rho_d = \frac{2 \cdot s_0^2 \cdot T}{\sigma_d^2} = \frac{2 \cdot s_0^2 \cdot T}{N_0}\cdot \frac{1}{\sigma_{d,\hspace{0.05cm} {\rm norm}}^2} \hspace{0.3cm}\Rightarrow \hspace{0.3cm} \eta = \frac{\rho_d }{\rho_{d,\hspace{0.05cm} {\rm max}}} = \frac{\rho_d }{2 \cdot s_0^2 \cdot T/N_0} = \frac{1}{\sigma_{d,\hspace{0.05cm} {\rm norm}}^2} \hspace{0.05cm}.\]
- Die Optimierungsaufgabe beschränkt sich also darauf, das Empfangsfilter HE(f) so zu bestimmen, dass die normierte Rauschleistung vor dem Entscheider den kleinstmöglichen Wert annimmt:
- \[\sigma_{d,\hspace{0.05cm} {\rm norm}}^2 = \frac{\sigma_d^2}{N_0/ T} =T \cdot \int_{-\infty}^{+\infty} |H_{\rm E}(f)|^2 \,{\rm d} f \stackrel {!}{=} {\rm Minimum}\hspace{0.05cm}.\]
- Wir bezeichnen die Konfiguration als Optimale Nyquistentzerrung (ONE). Obwohl diese auch – und besonders effektiv – bei Mehrstufensystemen anwendbar ist, setzen wir zunächst M = 2.
Wirkungsweise des Transversalfilters (1)
Verdeutlichen wir uns zunächst die Aufgabe des symmetrischen Transversalfilters
\[H_{\rm TF}(f) \hspace{0.4cm}\bullet\!\!-\!\!\!-\!\!\!-\!\!\circ \hspace{0.4cm} h_{\rm TF}(t) = \sum_{\lambda = -N}^{+N} k_\lambda \cdot \delta(t - \lambda \cdot T) \hspace{0.05cm}.\]
N gibt die Ordnung des Filters an. Für die Filterkoeffizienten gilt k–λ = kλ. Dieses Filter ist somit durch die Koeffizienten k0, ... , kN vollständig bestimmt. Die Grafik zeigt ein Filter zweiter Ordnung (N = 2).
Für den Eingangsimpuls gm(t) setzen wir ohne Einschränkung der Allgemeingültigkeit voraus, dass dieser
- symmetrisch um t = 0 ist (Ausgang des Matched–Filters),
- zu den Zeiten νT und –νT den Wert gm(ν) besitzt.
Damit sind die Eingangsimpulswerte:
\[...\hspace{0.2cm} , g_m(3),\hspace{0.15cm}g_m(2),\hspace{0.15cm}g_m(1),\hspace{0.15cm}\hspace {0.15cm}g_m(0),\hspace{0.15cm}g_m(1),\hspace{0.15cm}g_m(2),\hspace{0.15cm}g_m(3),\hspace{0.1cm} ... \hspace{0.05cm}.\]
Für den Detektionsgrundimpuls gd(t) am Filterausgang ergeben sich demzufolge zu den Zeitpunkten νT mit den Abkürzungen g0 = gd(t = 0), g1 = gd(t = ±T), g2 = gd(t = ±2T) folgende Werte:
\[ t = 0\hspace{-0.1cm}:\hspace{0.2cm}g_0 = k_0 \cdot g_m(0) + k_1 \cdot 2 \cdot g_m(1) \hspace{1.23cm}+k_2 \cdot 2 \cdot g_m(2),\hspace{0.05cm} \] \[ t = \pm T\hspace{-0.1cm}:\hspace{0.2cm}g_1 = k_0 \cdot g_m(1) + k_1 \cdot [g_m(0)+g_m(2)]+ k_2 \cdot [g_m(1)+g_m(3)], \] \[ t = \pm 2T\hspace{-0.1cm}:\hspace{0.2cm}g_2 = k_0 \cdot g_m(2) + k_1 \cdot [g_m(1)+g_m(3)]+ k_2 \cdot [g_m(2)+g_m(4)] \hspace{0.05cm}. \]
Aus diesem System mit drei linear unabhängigen Gleichungen kann man nun die Filterkoeffizienten k0, k1 und k2 so bestimmen, dass der Detektionsgrundimpuls gd(t) durch die normierten Stützstellen
\[...\hspace{0.15cm} , g_3,\hspace{0.25cm}g_2 = 0 ,\hspace{0.15cm}g_1 = 0 ,\hspace{0.15cm}g_0 = 1,\hspace{0.15cm}g_1 = 0 ,\hspace{0.15cm}g_2 = 0 ,\hspace{0.25cm}g_3 ,\hspace{0.15cm} ...\]
vollständig gegeben ist. Auf der nächsten Seite wird die Optimierung der Filterkoeffizienten an einem einfachen Beispiel verdeutlicht.
Wirkungsweise des Transversalfilters (2)
\[g_m(t) = {\rm exp }\left ( - \sqrt{2 \cdot |t|/T}\right )\]
\[\Rightarrow \hspace{0.3cm} g_m(0) = 1 ,\hspace{0.15cm}g_m(1)= 0.243,\hspace{0.15cm}g_m(2)= 0.135,\hspace{0.15cm}g_m(3)= 0.086, \hspace{0.15cm}g_m(4)= 0.059 \hspace{0.05cm}.\]
Für den Ausgangsimpuls soll gd(0) = 1 und gd(±T) = 0 gelten. Hierzu eignet sich ein Laufzeitfilter erster Ordnung mit den Koeffizienten k0 und k1, die folgende Bedingungen erfüllen müssen:
\[t = \pm T\hspace{-0.1cm} : \hspace{0.2cm}g_1 = k_0 \cdot 0.243 + k_1 \cdot [1.000 +0.135] = 0\hspace{0.3cm}\Rightarrow \hspace{0.3cm}{k_1} = -0.214 \cdot {k_0}\hspace{0.05cm},\] \[ t = 0 \hspace{-0.1cm} : \hspace{0.2cm}g_0 = k_0 \cdot 1.000 + k_1 \cdot 2 \cdot 0.243= 1\hspace{0.3cm}\Rightarrow \hspace{0.3cm}0.896 \cdot {k_0} = 1 \hspace{0.05cm}.\]
Daraus erhält man die optimalen Filterkoeffizienten <nobr>k0 = 1.116</nobr> und k1 = 0.239. Das mittlere Diagramm zeigt, dass damit der erste Vorläufer und der erste Nachläufer kompensiert werden können und zugleich gd(0) = 1 gilt (gelbe Hinterlegung). Die weiteren Detektionsgrundimpulswerte (blaue Kreise) sind aber von 0 verschieden und bewirken Impulsinterferenzen.
Das untere Diagramm zeigt, dass mit einem Filter zweiter Ordnung (N = 2) Nulldurchgänge bei ±T und bei ±2T erzwungen werden, wenn die Koeffizienten k0 = 1.127, k1 = 0.219 und k2 = 0.075 geeignet gewählt sind. Das Gleichungssystem zur Bestimmung der optimalen Koeffizienten lautet dabei:
\[t = 0\hspace{-0.1cm}:\hspace{0.2cm}g_0 = k_0 \cdot 1.000 + k_1 \cdot 2 \cdot 0.243 + k_2 \cdot 2 \cdot 0.135 = 1\hspace{0.05cm},\\ t= \pm T\hspace{-0.1cm}:\hspace{0.2cm}g_1 = k_0 \cdot 0.243 + k_1 \cdot [1.000+0.135]+ k_2 \cdot [0.243+0.086] = 0\hspace{0.05cm}, \\ t = \pm 2 T\hspace{-0.1cm}:\hspace{0.2cm}g_2 = k_0 \cdot 0.135 + k_1 \cdot [0.243+0.086]+ k_2 \cdot [1.000 + 0.059]= 0 \hspace{0.05cm}.\]
Die Ergebnisse können wie folgt verallgemeinert werden:
- Mit einem Laufzeitfilter N–ter Ordnung können der Hauptwert gd(0) zu 1 (normiert) sowie die ersten N Nachläufer und die ersten N Vorläufer zu Null gemacht werden.
- Weitere Vor– und Nachläufer (|ν| > N) lassen sich so nicht kompensieren. Es ist auch möglich, dass diese außerhalb des Kompensationsbereichs vergrößert werden oder sogar neu entstehen.
- Im Grenzübergang N → ∞ (in der Praxis heißt das: ein Filter mit sehr vielen Koeffizienten) ist eine vollständige Nyquistentzerrung und damit eine impulsinterferenzfreie Übertragung möglich.