Applets:Abtastung analoger Signale und Signalrekonstruktion: Unterschied zwischen den Versionen

Aus LNTwww
Wechseln zu:Navigation, Suche
Zeile 18: Zeile 18:
 
==Theoretischer Hintergrund==
 
==Theoretischer Hintergrund==
  
=== Signaldarstellung mit orthonormalen Basisfunktionen ===
 
  
Wir gehen von einer Menge  $\{s_i(t)\}$  möglicher Sendesignale aus, die den möglichen Nachrichten  $m_i$  eineindeutig zugeordnet sind. Mit  $i = 1$, ... , $M$  gelte:
 
:$$m \in \{m_i \}, \hspace{0.2cm} s(t) \in \{s_i(t) \}\hspace{-0.1cm}: \hspace{0.3cm} m = m_i  \hspace{0.1cm} \Leftrightarrow \hspace{0.1cm} s(t) = s_i(t) \hspace{0.05cm}.$$
 
  
Für das Folgende setzen wir weiter voraus, dass die  $M$ Signale  $s_i(t)$  [[Signaldarstellung/Klassifizierung_von_Signalen#Energiebegrenzte_und_leistungsbegrenzte_Signale| energiebegrenzt]]  sind, was meist gleichzeitig bedeutet, dass sie nur von endlicher Dauer sind.
+
===Beschreibung der Abtastung im Zeitbereich===
  
{{BlaueBox|TEXT= 
+
[[Datei:P_ID1120__Sig_T_5_1_S1_neu.png|center|frame|Zur Zeitdiskretisierung des zeitkontinuierlichen Signals  $x(t)$]]
$\text{Satz:}$  Eine jede Menge  $\{s_1(t), \hspace{0.05cm}  \text{...} \hspace{0.05cm} , s_M(t)\}$  energiebegrenzter Signale lässt sich in  $N \le M$   '''orthonormale Basisfunktionen'''  $\varphi_1(t),  \hspace{0.05cm} \text{...} \hspace{0.05cm} , \varphi_N(t)$  entwickeln.  Es gilt:
 
  
:$$s_i(t) = \sum\limits_{j = 1}^{N}s_{ij} \cdot \varphi_j(t) ,
+
Im Folgenden verwenden wir für die Beschreibung der Abtastung folgende Nomenklatur:
\hspace{0.3cm}i = 1,\hspace{0.05cm} \text{...}\hspace{0.1cm} , M, \hspace{0.3cm}j = 1,\hspace{0.05cm} \text{...} \hspace{0.1cm}, N
+
*Das zeitkontinuierliche Signal sei  $x(t)$.
\hspace{0.05cm}.$$
+
*Das in äquidistanten Abständen  $T_{\rm A}$  abgetastete zeitdiskretisierte Signal sei  $x_{\rm A}(t)$.
 +
*Außerhalb der Abtastzeitpunkte  $\nu \cdot T_{\rm A}$  gilt stets  $x_{\rm A}(t) = 0$.
 +
*Die Laufvariable  $\nu$  sei  [[Signaldarstellung/Zum_Rechnen_mit_komplexen_Zahlen#Reelle_Zahlenmengen|ganzzahlig]]:      $\nu \in \mathbb{Z} = \{\hspace{0.05cm} \text{...}\hspace{0.05cm} , –3, –2, –1, \hspace{0.2cm}0, +1, +2, +3, \text{...} \hspace{0.05cm}\} $.
 +
*Dagegen ergibt sich zu den äquidistanten Abtastzeitpunkten mit der Konstanten  $K$:
 +
 +
:$$x_{\rm A}(\nu \cdot T_{\rm A}) = K \cdot x(\nu \cdot T_{\rm A})\hspace{0.05cm}.$$
  
Jeweils zwei Basisfunktionen  $\varphi_j(t)$  und  $\varphi_k(t)$  müssen orthonormal zueinander sein, das heißt, dass gelten muss  $(\delta_{jk}$  nennt man das [https://de.wikipedia.org/wiki/Kronecker-Delta Kronecker–Symbol]$)$:
+
Die Konstante hängt von der Art der Zeitdiskretisierung ab. Für die obige Skizze gilt  $K = 1$.
  
:$$<\hspace{-0.1cm}\varphi_j(t), \hspace{0.05cm}\varphi_k(t) \hspace{-0.1cm}> = \int_{-\infty}^{+\infty}\varphi_j(t) \cdot \varphi_k(t)\,d \it t = {\rm \delta}_{jk} =
+
{{BlaueBox|TEXT=
\left\{ \begin{array}{c} 1 \\
+
$\text{Definitionen:}$&nbsp;
0  \end{array} \right.\quad
 
\begin{array}{*{1}c} {\rm falls}\hspace{0.4cm}j = k\hspace{0.1cm}
 
\\ {\rm falls}\hspace{0.4cm} j \ne k \hspace{0.1cm}\\ \end{array}
 
\hspace{0.05cm}.$$}}<br>
 
  
Der Parameter&nbsp; $N$&nbsp; gibt dabei an, wieviele Basisfunktionen&nbsp; $\varphi_j(t)$&nbsp; benötigt werden, um die&nbsp; $M$&nbsp; möglichen Sendesignale darzustellen.&nbsp; Mit anderen Worten: &nbsp; $N$&nbsp; ist die ''Dimension des Vektorraums'', der von den&nbsp; $M$&nbsp; Signalen aufgespannt wird.&nbsp; Dabei gilt:
+
Der&nbsp; '''Diracpuls (im Zeitbereich)'''&nbsp; besteht aus unendlich vielen Diracimpulsen, jeweils im gleichen Abstand&nbsp; $T_{\rm A}$&nbsp; und alle mit gleichem Impulsgewicht&nbsp; $T_{\rm A}$:
*Ist&nbsp; $N = M$, so sind alle Sendesignale zueinander orthogonal.&nbsp; Sie sind nicht notwendigerweise orthonormal, das heißt, die Energien&nbsp; $E_i = \ <\hspace{-0.01cm}s_i(t), \hspace{0.05cm}s_i(t) \hspace{-0.01cm}>$&nbsp; können durchaus ungleich Eins sein.<br>
+
*Der Fall&nbsp; $N < M$&nbsp; ergibt sich, wenn mindestens ein Signal&nbsp; $s_i(t)$&nbsp; als Linearkombination von Basisfunktionen&nbsp; $\varphi_j(t)$&nbsp; dargestellt werden kann, die sich bereits aus anderen Signalen&nbsp; $s_j(t) \ne s_i(t)$&nbsp; ergeben haben.<br>
+
:$$p_{\delta}(t) = \sum_{\nu = - \infty }^{+\infty} T_{\rm A} \cdot
 +
\delta(t- \nu \cdot T_{\rm A}
 +
)\hspace{0.05cm}.$$
 +
 +
* Unter&nbsp; '''Abtastung'''&nbsp; verstehen wir im Folgenden die Multiplikation des zeitkontinuierlichen Signals&nbsp; $x(t)$&nbsp; mit einem solchen&nbsp; '''Diracpuls''':
 +
 +
:$$x_{\rm A}(t) = x(t) \cdot p_{\delta}(t)\hspace{0.05cm}.$$}}
  
 +
Aufgrund dieser Definition ergeben sich für das abgetastete Signal folgende Eigenschaften:
 +
*Das abgetastete Signal zum betrachteten Zeitpunkt&nbsp; $(\nu \cdot T_{\rm A})$&nbsp; ist gleich&nbsp; $T_{\rm A} \cdot x(\nu \cdot T_{\rm A}) · \delta (0)$.
 +
*Da die Diracfunktion&nbsp; $\delta (t)$&nbsp; zur Zeit&nbsp; $t = 0$&nbsp; unendlich ist, sind eigentlich alle Signalwerte&nbsp; $x_{\rm A}(\nu \cdot T_{\rm A})$&nbsp; ebenfalls unendlich groß.
 +
*Somit ist auch der auf der letzten Seite eingeführte Faktor&nbsp; $K$&nbsp; eigentlich unendlich groß.
 +
*Zwei Abtastwerte&nbsp; $x_{\rm A}(\nu_1 \cdot T_{\rm A})$&nbsp; und&nbsp; $x_{\rm A}(\nu_2 \cdot T_{\rm A})$&nbsp; unterscheiden sich jedoch  im gleichen Verhältnis wie die Signalwerte&nbsp; $x(\nu_1 \cdot T_{\rm A})$&nbsp; und&nbsp; $x(\nu_2 \cdot T_{\rm A})$.
 +
*Die Abtastwerte von&nbsp; $x(t)$&nbsp; erscheinen in den Impulsgewichten der Diracfunktionen:
 +
 +
:$$x_{\rm A}(t) =  \sum_{\nu = - \infty }^{+\infty} T_{\rm A} \cdot x(\nu \cdot T_{\rm A})\cdot
 +
\delta (t- \nu \cdot T_{\rm A}
 +
)\hspace{0.05cm}.$$
  
[[Datei:P ID1993 Dig T 4 1 S2 version1.png|right|frame|Darstellung der drei Sendesignale durch zwei Basisfunktionen|class=fit]]
+
*Die zusätzliche Multiplikation mit&nbsp; $T_{\rm A}$&nbsp; ist erforderlich, damit&nbsp; $x(t)$&nbsp; und&nbsp; $x_{\rm A}(t)$&nbsp; gleiche Einheit besitzen. Beachten Sie hierbei, dass&nbsp; $\delta (t)$&nbsp; selbst die Einheit „1/s” aufweist.
{{GraueBox|TEXT= 
 
$\text{Beispiel 1:}$&nbsp; Wir betrachten&nbsp; $M = 3$&nbsp; energiebegrenzte Signale gemäß der Grafik.  
 
  
Man erkennt sofort:
 
*Die Signale&nbsp; $s_1(t)$&nbsp;  und &nbsp;$s_2(t)$&nbsp; sind zueinander orthogonal.<br>
 
  
*Die Energien sind&nbsp; $E_1 = A^2 \cdot T = E$&nbsp; und &nbsp;$E_2 = (A/2)^2 \cdot T = E/4$.<br>
+
Die folgenden Seiten werden zeigen, dass diese gewöhnungsbedürftigen Gleichungen durchaus zu sinnvollen Ergebnissen führen, wenn man sie konsequent  anwendet.
  
*Die Basisfunktionen&nbsp; $\varphi_1(t)$&nbsp; und &nbsp;$\varphi_2(t)$&nbsp; sind jeweils formgleich mit&nbsp; $s_1(t)$&nbsp;  bzw.&nbsp;  $s_2(t)$.
 
*Beide Signale besitzen jeweils die Energie &bdquo;Eins&rdquo;:
 
  
:$$\varphi_1(t)=\frac{s_1(t)}{\sqrt{E_1} } = \frac{s_1(t)}{\sqrt{A^2 \cdot T} } = \frac{1}{\sqrt{ T} \cdot \frac{s_1(t)}{A}$$
+
==Diracpuls im Zeit- und im Frequenzbereich==
:$$\hspace{0.5cm}\Rightarrow \hspace{0.1cm}s_1(t) = s_{11} \cdot \varphi_1(t)\hspace{0.05cm},\hspace{0.1cm}s_{11} = \sqrt{E}\hspace{0.05cm},$$
+
<br>
:$$\varphi_2(t) =\frac{s_2(t)}{\sqrt{E_2} } = \frac{s_2(t)}{\sqrt{(A/2)^2 \cdot T} } = \frac{1}{\sqrt{ T} \cdot \frac{s_2(t)}{A/2}\hspace{0.05cm}$$
+
{{BlaueBox|TEXT=
:$$\hspace{0.5cm}\Rightarrow \hspace{0.1cm}s_2(t) = s_{21} \cdot \varphi_2(t)\hspace{0.05cm},\hspace{0.1cm}s_{21} = {\sqrt{E} }/{2}\hspace{0.05cm}.$$
+
$\text{Satz:}$&nbsp; Entwickelt man den&nbsp; '''Diracpuls'''&nbsp; in eine&nbsp; [[Signaldarstellung/Fourierreihe|Fourierreihe]]&nbsp; und transformiert diese unter Anwendung des&nbsp; [[Signaldarstellung/Gesetzmäßigkeiten_der_Fouriertransformation#Verschiebungssatz|Verschiebungssatzes]]&nbsp; in den Frequenzbereich, so ergibt sich folgende Korrespondenz:
 +
   
 +
:$$p_{\delta}(t) = \sum_{\nu = - \infty }^{+\infty} T_{\rm A} \cdot
 +
\delta(t- \nu \cdot T_{\rm A}
 +
  )\hspace{0.2cm}\circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\, \hspace{0.2cm} P_{\delta}(f) = \sum_{\mu = - \infty }^{+\infty} \delta
 +
(f- \mu \cdot f_{\rm A} ).$$
  
*Das Signal&nbsp; $s_3(t)$&nbsp; kann durch die vorher bestimmten Basisfunktionen&nbsp; $\varphi_1(t)$&nbsp; und &nbsp;$\varphi_2(t)$&nbsp; ausgedrückt werden:
+
Hierbei gibt&nbsp; $f_{\rm A} = 1/T_{\rm A}$&nbsp; den Abstand zweier benachbarter Diraclinien im Frequenzbereich an. }}
:$$s_3(t) =s_{31} \cdot \varphi_1(t) + s_{32} \cdot \varphi_2(t)\hspace{0.05cm},$$
 
:$$\hspace{0.5cm}\Rightarrow \hspace{0.1cm}
 
s_{31} = {A}/{2} \cdot \sqrt {T}= {\sqrt{E} }/{2}\hspace{0.05cm}, \hspace{0.2cm}s_{32} = - A \cdot \sqrt {T} = -\sqrt{E}  \hspace{0.05cm}.$$
 
  
Trotz&nbsp; $M=3$&nbsp; gilt also im vorliegenen Fall nur&nbsp; $N=2$.
+
 +
{{BlaueBox|TEXT=
 +
$\text{Beweis:}$&nbsp; Die Herleitung der hier angegebenen Spektralfunktion&nbsp; $P_{\delta}(f)$&nbsp; geschieht in mehreren Schritten:
  
Im rechten unteren Bild sind die Signale in einer 2D&ndash;Darstellung mit den Basisfunktionen&nbsp; $\varphi_1(t)$&nbsp; und &nbsp;$\varphi_2(t)$&nbsp; als Achsen dargestellt, wobei&nbsp; $E = A^2 \cdot T$&nbsp; gilt und der Zusammenhang zu den anderen Grafiken durch die Farbgebung zu erkennen ist.
+
'''(1)'''&nbsp;&nbsp; Da&nbsp; $p_{\delta}(t)$&nbsp; periodisch mit dem konstanten Abstand&nbsp; $T_{\rm A}$&nbsp; zwischen zwei Diraclinien ist, kann die&nbsp; [[Signaldarstellung/Fourierreihe#Komplexe_Fourierreihe|(komplexe) Fourierreihendarstellung]]&nbsp; angewendet werden:
 +
 +
:$$p_{\delta}(t) =  \sum_{\mu = - \infty }^{+\infty} D_{\mu} \cdot
 +
{\rm e}^{\hspace{0.05cm}{\rm j} \hspace{0.05cm} \cdot 2 \hspace{0.05cm} \pi \hspace{0.05cm}\cdot \hspace{0.05cm}\mu \hspace{0.05cm}\cdot \hspace{0.05cm}t/T_{\rm A} }
 +
\hspace{0.3cm}{\rm mit}\hspace{0.3cm}
 +
D_{\mu} = \frac{1}{T_{\rm A} } \cdot \int_{-T_{\rm A}/2
 +
}^{+T_{\rm A}/2}p_{\delta}(t) \cdot {\rm e}^{- {\rm j} \hspace{0.05cm}\cdot \hspace{0.05cm} 2 \pi \hspace{0.05cm}\cdot \hspace{0.05cm}\mu \hspace{0.05cm}
 +
\cdot \hspace{0.05cm}t/T_{\rm A} }\hspace{0.1cm} {\rm d}t\hspace{0.05cm}.$$
  
Die vektoriellen Repräsentanten der Signale&nbsp; $s_1(t)$,&nbsp; $s_2(t)$&nbsp; und&nbsp; $s_3(t)$&nbsp; in diesem zweidimensionellen Vektorraum lassen sich daraus wie folgt ablesen:
+
'''(2)'''&nbsp;&nbsp; Im Bereich von&nbsp; $–T_{\rm A}/2$&nbsp; bis&nbsp; $+T_{\rm A}/2$&nbsp; gilt für den Diracpuls im Zeitbereich: &nbsp; $p_{\delta}(t) = T_{\rm A} \cdot \delta(t)$. Damit kann man für die komplexen Fourierkoeffizienten schreiben: &nbsp;  
:$$\mathbf{s}_1 = (\sqrt{ E}, \hspace{0.1cm}0), \hspace{0.5cm}
+
:$$D_{\mu} = \int_{-T_{\rm A}/2
\mathbf{s}_2 = (0, \hspace{0.1cm}\sqrt{ E}/2), \hspace{0.5cm} \mathbf{s}_3 = (\sqrt{ E}/2,\hspace{0.1cm}-\sqrt{ E} )   \hspace{0.05cm}.$$}}
+
}^{+T_{\rm A}/2}{\delta}(t) \cdot {\rm e}^{- {\rm j} \hspace{0.05cm}\cdot \hspace{0.05cm} 2 \pi \hspace{0.05cm}\cdot \hspace{0.05cm}\mu \hspace{0.05cm}
<br clear= all>
+
\cdot \hspace{0.05cm}t/T_{\rm A} }\hspace{0.1cm} {\rm d}t\hspace{0.05cm}.$$
 +
'''(3)'''&nbsp;&nbsp; Unter Berücksichtigung der Tatsache, dass für&nbsp; $t \neq 0$&nbsp; der Diracimpuls Null ist und für&nbsp; $t = 0$&nbsp; der komplexe Drehfaktor gleich&nbsp; $1$, gilt weiter:
 +
:$$D_{\mu} = \int_{- T_{\rm A}/2
 +
}^{+T_{\rm A}/2}{\delta}(t) \hspace{0.1cm} {\rm d}t = 1\hspace{0.5cm}{\Rightarrow}\hspace{0.5cm}
 +
p_{\delta}(t) =  \sum_{\mu = - \infty }^{+\infty} {\rm e}^{ {\rm j} \hspace{0.05cm}
 +
\cdot 2 \hspace{0.05cm} \pi \cdot \hspace{0.05cm}\mu \hspace{0.05cm}\cdot \hspace{0.05cm}t/T_{\rm A} }\hspace{0.05cm}.
 +
$$
 +
'''(4)'''&nbsp;&nbsp; Der&nbsp; [[Signaldarstellung/Gesetzmäßigkeiten_der_Fouriertransformation#Verschiebungssatz|Verschiebungssatz im Frequenzbereich]]&nbsp; lautet mit&nbsp; $f_{\rm A} = 1/T_{\rm A}$:
 +
:$${\rm e}^{ {\rm j} \hspace{0.05cm}
 +
\hspace{0.05cm} \cdot 2 \hspace{0.05cm} \pi \hspace{0.05cm}\cdot \hspace{0.05cm}\mu \hspace{0.05cm}\cdot \hspace{0.05cm}
 +
f_{\rm A}\hspace{0.05cm}\cdot \hspace{0.05cm}t}\hspace{0.2cm}\circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\, \hspace{0.2cm}
 +
\delta
 +
(f- \mu \cdot f_{\rm A}
 +
)\hspace{0.05cm}.$$  
 +
'''(5)'''&nbsp;&nbsp; Wendet man das Ergebnis auf jeden einzelnen Summanden an, so erhält man schließlich:
 +
 +
:$$P_{\delta}(f) =  \sum_{\mu = - \infty }^{+\infty} \delta
 +
(f- \mu \cdot f_{\rm A}
 +
)\hspace{0.05cm}.$$
 +
<div align="right">q.e.d.</div>}}
  
=== Das Verfahren nach Gram-Schmidt===
 
  
Im letzten &nbsp;$\text{Beispiel}$&nbsp; war die Bestimmung der beiden orthonormalen Basisfunktionen&nbsp; $\varphi_1(t)$&nbsp; und&nbsp; $\varphi_2(t)$&nbsp; sehr einfach, da diese formgleich mit&nbsp; $s_1(t)$&nbsp;  bzw.&nbsp;  $s_2(t)$&nbsp; waren. Das&nbsp; [https://de.wikipedia.org/wiki/Gram-Schmidtsches_Orthogonalisierungsverfahren Gram&ndash;Schmidt&ndash;Verfahren]&nbsp; findet die Basisfunktionen&nbsp; $\varphi_1(t)$, ... , $\varphi_N(t)$&nbsp; für beliebig vorgebbare Signale&nbsp; $s_1(t)$, ... , $s_M(t)$, und zwar wie folgt:
+
Das Ergebnis besagt:
 +
*Der Diracpuls&nbsp; $p_{\delta}(t)$&nbsp; im Zeitbereich besteht aus unendlich vielen Diracimpulsen, jeweils im gleichen Abstand&nbsp; $T_{\rm A}$&nbsp; und alle mit gleichem Impulsgewicht&nbsp; $T_{\rm A}$.
 +
*Die Fouriertransformierte von&nbsp; $p_{\delta}(t)$&nbsp;  ergibt wiederum einen Diracpuls, aber nun im Frequenzbereich  &nbsp; &nbsp; $P_{\delta}(f)$.
 +
*$P_{\delta}(f)$&nbsp; besteht ebenfalls aus unendlich vielen Diracimpulsen, nun aber im jeweiligen Abstand&nbsp; $f_{\rm A} = 1/T_{\rm A}$&nbsp; und alle mit dem Impulsgewicht&nbsp; $1$.
 +
*Die Abstände der Diraclinien in der Zeit– und Frequenzbereichsdarstellung folgen demnach dem&nbsp; [[Signaldarstellung/Gesetzmäßigkeiten_der_Fouriertransformation#Reziprozit.C3.A4tsgesetz_von_Zeitdauer_und_Bandbreite|Reziprozitätsgesetz]]: &nbsp;  
 +
:$$T_{\rm A} \cdot f_{\rm A} = 1 \hspace{0.05cm}.$$
  
*Die erste Basisfunktion&nbsp; $\varphi_1(t)$&nbsp; ist stets formgleich mit&nbsp; $s_1(t)$. Es gilt:
 
:$$\varphi_1(t) = \frac{s_1(t)}{\sqrt{E_1}} = \frac{s_1(t)}{|| s_1(t)||}
 
\hspace{0.3cm}\Rightarrow \hspace{0.3cm} || \varphi_1(t) || = 1, \hspace{0.2cm}s_{11} =|| s_1(t)||,\hspace{0.2cm}s_{1j} = 0 \hspace{0.2cm}{\rm f{\rm \ddot{u}r }}\hspace{0.2cm} j \ge 2
 
\hspace{0.05cm}.$$
 
  
{{BlaueBox|TEXT=
+
[[Datei:P_ID1121__Sig_T_5_1_S3_NEU.png|right|frame|Diracpuls im Zeit- und Frequenzbereich]]
$\text{Hinweise zur Nomenklatur:}$&nbsp;  
+
{{GraueBox|TEXT=
 +
$\text{Beispiel 1:}$&nbsp; Die Grafik verdeutlicht die obigen Aussagen für
 +
*$T_{\rm A} = 50\,{\rm &micro;s}$,
 +
*$f_{\rm A} = 1/T_{\rm A} = 20\,\text{kHz}$ .
  
'''(1)'''&nbsp; Ausgehend von zwei reellen und energiebegrenzten Zeitfunktionen &nbsp;$x(t)$&nbsp; und &nbsp;$y(t)$&nbsp; erhält man für das &nbsp;[https://de.wikipedia.org/wiki/Inneres_Produkt innere Produkt] allgemein:
 
:$$<\hspace{-0.01cm}x(t), \hspace{0.05cm}y(t) \hspace{-0.01cm}> \hspace{0.15cm}= \int_{-\infty}^{+\infty}x(t) \cdot y(t)\,d \it t
 
\hspace{0.05cm}.$$
 
  
'''(2)'''&nbsp; Daraus ergibt sich die&nbsp; [https://de.wikipedia.org/wiki/Euklidische_Norm Euklidische Norm]&nbsp; der Zeitfunktion $s_1(t)$:
+
Man erkennt aus dieser Skizze auch die unterschiedlichen Impulsgewichte von&nbsp; $p_{\delta}(t)$&nbsp; und&nbsp; $P_{\delta}(f)$.}}
:$$\vert \vert s_1(t) \vert \vert = \sqrt{<\hspace{-0.01cm}s_1(t), \hspace{0.15cm}s_1(t) \hspace{-0.01cm}>} $$}}
+
===Prinzip und Motivation===
 +
<br>
 +
Viele Nachrichtensignale sind analog und damit  gleichzeitig&nbsp; [[Signaldarstellung/Klassifizierung_von_Signalen#Zeitkontinuierliche_und_zeitdiskrete_Signale|zeitkontinuierlich]]&nbsp; und&nbsp; [[Signaldarstellung/Klassifizierung_von_Signalen#Wertkontinuierliche_und_wertdiskrete_Signale|wertkontinuierlich]]. Soll ein solches Analogsignal mittels eines Digitalsystems übertragen werden, so sind folgende Vorverarbeitungsschritte erforderlich:
 +
*die&nbsp; '''Abtastung'''&nbsp; des Nachrichtensignals&nbsp; $x(t)$, die zweckmäßigerweise – aber nicht notwendigerweise – zu äquidistanten Zeitpunkten erfolgt &nbsp; &rArr; &nbsp; '''Zeitdiskretisierung''',
 +
*die&nbsp; '''Quantisierung'''&nbsp; der Abtastwerte, um so die Anzahl&nbsp; $M$&nbsp; der möglichen Werte auf einen endlichen Wert zu begrenzen  &nbsp; &rArr; &nbsp; '''Wertdiskretisierung'''.
  
  
Es wird nun angenommen, dass aus den Signalen&nbsp; $s_1(t)$, ... , $s_{k-1}(t)$&nbsp; bereits die Basisfunktionen&nbsp; $\varphi_1(t)$, ... , $\varphi_{n-1}(t)$&nbsp; berechnet wurden &nbsp;$(n \le k)$.  
+
Die Quantisierung wird erst im Kapitel&nbsp; [[Modulationsverfahren/Pulscodemodulation|Pulscodemodulation]]&nbsp; des Buches „Modulationsverfahren” im Detail behandelt.
  
*Dann berechnen wir mittels der nächsten Funktion&nbsp; $s_k(t)$&nbsp; die Hilfsfunktion
+
[[Datei:P_ID1120__Sig_T_5_1_S1_neu.png|center|frame|Zur Zeitdiskretisierung des zeitkontinuierlichen Signals&nbsp; $x(t)$]]
:$$\theta_k(t) = s_k(t) - \sum\limits_{j = 1}^{n-1}s_{kj} \cdot \varphi_j(t) \hspace{0.4cm}{\rm mit}\hspace{0.4cm}
 
s_{kj} = \hspace{0.01cm} < \hspace{-0.1cm} s_k(t), \hspace{0.05cm}\varphi_j(t) \hspace{-0.01cm} >, \hspace{0.2cm} j = 1, \hspace{0.05cm} \text{...}\hspace{0.05cm}, n-1\hspace{0.05cm}.$$
 
  
*Hat diese Hilfsfunktion die Norm  &nbsp; $||\theta_k(t)|| = 0$, so liefert&nbsp; $s_k(t)$&nbsp; keine neue Basisfunktion.&nbsp; Vielmehr lässt sich dann&nbsp; $s_k(t)$&nbsp; durch die&nbsp; $n-1$&nbsp; bereits vorher gefundenen Basisfunktionen &nbsp;$\varphi_1(t)$, ... , $\varphi_{n-1}(t)$&nbsp;  ausdrücken:
+
Im Folgenden verwenden wir für die Beschreibung der Abtastung folgende Nomenklatur:
:$$s_k(t) = \sum\limits_{j = 1}^{n-1}s_{kj}\cdot \varphi_j(t) \hspace{0.05cm}.$$
+
*Das zeitkontinuierliche Signal sei&nbsp; $x(t)$.
 +
*Das in äquidistanten Abständen&nbsp; $T_{\rm A}$&nbsp; abgetastete zeitdiskretisierte Signal sei&nbsp; $x_{\rm A}(t)$.
 +
*Außerhalb der Abtastzeitpunkte&nbsp; $\nu \cdot T_{\rm A}$&nbsp; gilt stets&nbsp; $x_{\rm A}(t) = 0$.
 +
*Die Laufvariable&nbsp; $\nu$&nbsp; sei&nbsp; [[Signaldarstellung/Zum_Rechnen_mit_komplexen_Zahlen#Reelle_Zahlenmengen|ganzzahlig]]:  &nbsp; &nbsp; $\nu \in \mathbb{Z} =  \{\hspace{0.05cm} \text{...}\hspace{0.05cm} , –3, –2, –1, \hspace{0.2cm}0, +1, +2, +3, \text{...} \hspace{0.05cm}\} $.
 +
*Dagegen ergibt sich zu den äquidistanten Abtastzeitpunkten mit der Konstanten&nbsp; $K$:
 +
   
 +
:$$x_{\rm A}(\nu \cdot T_{\rm A}) = K \cdot x(\nu \cdot T_{\rm A})\hspace{0.05cm}.$$
  
*Eine neue Basisfunktion&nbsp; (nämlich die &nbsp;$n$&ndash;te)&nbsp; ergibt sich nur für den Fall&nbsp; $||\theta_k(t)|| \ne 0$:
+
Die Konstante hängt von der Art der Zeitdiskretisierung ab. Für die obige Skizze gilt&nbsp; $K = 1$.
  
:$$\varphi_n(t) =  \frac{\theta_k(t)}{|| \theta_k(t)||}
 
\hspace{0.3cm}\Rightarrow \hspace{0.3cm} || \varphi_n(t) || = 1\hspace{0.05cm}.$$
 
  
Diese Prozedur wird solange fortgesetzt, bis alle&nbsp; $M$&nbsp; Signale berücksichtigt wurden.
+
==Frequenzbereichsdarstellung==
*Danach hat man alle&nbsp; $N \le M$&nbsp; orthonormalen Basisfunktionen&nbsp; $\varphi_j(t)$&nbsp; gefunden.  
+
<br>
*Der Sonderfall&nbsp; $N = M$&nbsp; ergibt sich nur dann, wenn alle&nbsp; $M$&nbsp; Signale linear voneinander unabhängig sind.<br>
+
Zum Spektrum des abgetasteten Signals&nbsp; $x_{\rm A}(t)$&nbsp; kommt man durch Anwendung des&nbsp; [[Signaldarstellung/Faltungssatz_und_Faltungsoperation#Faltung_im_Frequenzbereich|Faltungssatzes]]. Dieser besagt, dass der Multiplikation im Zeitbereich die Faltung im Spektralbereich entspricht:
 +
 +
:$$x_{\rm A}(t) = x(t) \cdot p_{\delta}(t)\hspace{0.2cm}\circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\, \hspace{0.2cm}
 +
X_{\rm A}(f) = X(f) \star P_{\delta}(f)\hspace{0.05cm}.$$
  
 +
Aus dem Spektrum&nbsp; $X(f)$&nbsp; wird durch Faltung mit der um&nbsp; $\mu \cdot f_{\rm A}$&nbsp; verschobenen Diraclinie:
 +
 +
:$$X(f) \star \delta
 +
(f- \mu \cdot f_{\rm A}
 +
)= X (f- \mu \cdot f_{\rm A}
 +
)\hspace{0.05cm}.$$
  
{{GraueBox|TEXT=   
+
Wendet man dieses Ergebnis auf alle Diraclinien des Diracpulses an, so erhält man schließlich:
$\text{Beispiel 2:}$&nbsp; Wir betrachten die &nbsp;$M = 4$&nbsp; energiebegrenzten Signale &nbsp;$s_1(t)$, ... , $s_4(t)$&nbsp; entsprechend der Grafik. Zur Vereinfachung der Berechnungen sind hier sowohl die Amplituden als auch die Zeit normiert.
+
 +
:$$X_{\rm A}(f) = X(f) \star \sum_{\mu = - \infty }^{+\infty} \delta
 +
  (f- \mu \cdot f_{\rm A}
 +
) = \sum_{\mu = - \infty }^{+\infty} X (f- \mu \cdot f_{\rm A}
 +
)\hspace{0.05cm}.$$
  
[[Datei:Dig_T_4_1_S3_neu.png|center|frame|Zum Gram-Schmidt-Verfahren|class=fit]]
+
{{BlaueBox|TEXT=
 +
$\text{Fazit:}$&nbsp; Die Abtastung des analogen Zeitsignals&nbsp; $x(t)$&nbsp; in äquidistanten Abständen&nbsp; $T_{\rm A}$&nbsp; führt im Spektralbereich zu einer&nbsp; '''periodischen Fortsetzung'''&nbsp; von&nbsp; $X(f)$&nbsp; mit dem Frequenzabstand&nbsp; $f_{\rm A} = 1/T_{\rm A}$.}}
  
Man erkennt aus diesen Skizzen:
 
*Die Basisfunktion&nbsp; $\varphi_1(t)$&nbsp; ist formgleich mit&nbsp; $s_1(t)$.&nbsp; Wegen&nbsp; $E_1 = \vert \vert s_1(t) \vert \vert ^2 = 3 \cdot 0.5^2 = 0.75$&nbsp; ergibt sich&nbsp; $s_{11} = \vert \vert s_1(t) \vert \vert = 0.866$.&nbsp; $\varphi_1(t)$&nbsp; selbst besitzt abschnittsweise die Werte&nbsp; $\pm 0.5/0.866 = \pm0.577$.
 
  
*Zur Berechnung der Hilfsfunktion&nbsp; $\theta_2(t)$&nbsp; berechnen wir
+
{{GraueBox|TEXT=
 +
$\text{Beispiel 2:}$&nbsp;
 +
Die obere Grafik zeigt&nbsp; (schematisch!)&nbsp; das Spektrum&nbsp; $X(f)$&nbsp; eines Analogsignals&nbsp; $x(t)$, das Frequenzen bis&nbsp; $5 \text{ kHz}$&nbsp; beinhaltet.
  
:$$s_{21}  = \hspace{0.1cm} < \hspace{-0.01cm} s_2(t), \hspace{0.05cm}\varphi_1(t) \hspace{-0.1cm} > \hspace{0.01cm} = 0 \cdot (+0.577) + 1 \cdot (-0.577)+ 0 \cdot (-0.577)= -0.577$$
+
[[Datei:P_ID1122__Sig_T_5_1_S4_neu.png|center|frame|Spektrum des abgetasteten Signals]]
:$$ \Rightarrow  \hspace{0.3cm}\theta_2(t) = s_2(t) - s_{21} \cdot \varphi_1(t) = (0.333, \hspace{0.15cm} 0.667, \hspace{0.15cm} -0.333)
 
\hspace{0.3cm}\Rightarrow  \hspace{0.3cm}\vert \vert \theta_2(t) \vert \vert^2 = (1/3)^2 + (2/3)^2 + (-1/3)^2 = 0.667$$
 
:$$ \Rightarrow  \hspace{0.3cm} s_{22} = \sqrt{0.667} = 0.816,\hspace{0.3cm}
 
\varphi_2(t) = \theta_2(t)/s_{22} = (0.408, \hspace{0.15cm} 0.816, \hspace{0.15cm} -0.408)\hspace{0.05cm}. $$
 
  
*Die inneren Produkte zwischen&nbsp; $s_1(t)$&nbsp; mit&nbsp; $\varphi_1(t)$&nbsp; bzw. &nbsp;$\varphi_2(t)$&nbsp; liefern folgende Ergebnisse:
+
Tastet man das Signal mit der Abtastrate&nbsp; $f_{\rm A}\,\text{ = 20 kHz}$, also im jeweiligen Abstand&nbsp; $T_{\rm A}\, = {\rm 50 \, &micro;s}$&nbsp; ab, so erhält man das unten skizzierte periodische Spektrum&nbsp; $X_{\rm A}(f)$.
:$$s_{31}  \hspace{0.01cm} =  \hspace{0.1cm} < \hspace{-0.1cm} s_3(t), \hspace{0.07cm}\varphi_1(t) \hspace{-0.01cm} > \hspace{0.1cm} = 0.5 \cdot (+0.577) + 0.5 \cdot (-0.577)- 0.5 \cdot (-0.577)= 0.289,$$
+
*Da die Diracfunktionen unendlich schmal sind, beinhaltet das abgetastete Signal&nbsp; $x_{\rm A}(t)$&nbsp; auch beliebig hochfrequente Anteile.  
:$$s_{32}  \hspace{0.1cm} =  \hspace{0.01cm} < \hspace{-0.1cm} s_3(t), \hspace{0.07cm}\varphi_2(t) \hspace{-0.01cm} > \hspace{0.1cm} = 0.5 \cdot (+0.408) + 0.5 \cdot (+0.816)- 0.5 \cdot (-0.408)= 0.816$$
+
*Dementsprechend ist die Spektralfunktion&nbsp; $X_{\rm A}(f)$&nbsp; des abgetasteten Signals bis ins Unendliche ausgedehnt.}}
:$$\Rightarrow  \hspace{0.3cm}\theta_3(t) = s_3(t) - 0.289 \cdot \varphi_1(t)- 0.816 \cdot \varphi_2(t) = 0\hspace{0.05cm}.$$
 
  
Das bedeutet: &nbsp; Die grüne Funktion&nbsp; $s_3(t)$&nbsp; liefert keine neue Basisfunktion&nbsp; $\varphi_3(t)$, im Gegensatz zur Funktion&nbsp; $s_4(t)$. Die numerischen Ergebnisse hierfür können der Grafik entnommen werden.}}
 
  
 +
==Signalrekonstruktion==
 +
<br>
 +
Die Signalabtastung ist bei einem digitalen Übertragungssystem kein Selbstzweck, sondern sie muss irgendwann wieder rückgängig gemacht werden. Betrachten wir zum Beispiel das folgende System:
  
===Die verschiedenen Rubriken bei der Auswahl der Programmparameter===
+
[[Datei:P_ID1123__Sig_T_5_1_S5a_neu.png|center|frame|Signalabtastung und Signalrekonstruktion]]
Das Programm bietet insgesamt&nbsp; $4 \cdot 6 = 24$&nbsp; Möglichkeiten zur Einstellung der jeweiligen Menge &nbsp;$\{s_i(t)\}$&nbsp; möglicher Sendesignale.&nbsp; Diese&nbsp; $24$&nbsp; Parametersätze sind in vier Rubriken eingeteilt. Die vier Rubriküberschriften treffen den Sachverhalt nicht hundertprozentig und sind deshalb in Hochkommata gesetzt: 
 
  
'''(1)'''&nbsp; Rubrik&nbsp; <u>&bdquo;Basisband&rdquo;</u> &nbsp; &rArr; &nbsp; gültig für die Einstellungen &nbsp;$\rm (A)$&nbsp; ... &nbsp;$\rm (F)$:
+
*Das Analogsignal&nbsp; $x(t)$&nbsp; mit der  Bandbreite&nbsp; $B_{\rm NF}$&nbsp; wird wie oben beschrieben abgetastet.
[[Datei:Gram_1_version2.png|right|frame|Signalform bei &bdquo;Basisband&rdquo;]]
+
*Am Ausgang eines idealen Übertragungssystems liegt das ebenfalls zeitdiskrete Signal&nbsp; $y_{\rm A}(t) = x_{\rm A}(t)$&nbsp; vor.  
*Jedes Mustersignal&nbsp; $s_i(t)$&nbsp; besteht aus drei Rechteckfunktionen unterschiedlicher Höhen und jeweiliger Dauer&nbsp; $T$.&nbsp;
+
*Die Frage ist nun, wie der Block&nbsp; '''Signalrekonstruktion'''&nbsp; zu gestalten ist, damit auch&nbsp; $y(t) = x(t)$&nbsp; gilt.
*Die einzelnen Rechteckhöhen sind Vielfache von&nbsp; $\pm 0.25$&nbsp; und die gesamte Signaldauer ergibt&nbsp; $3T$.
+
 
*Mit dem seitlichen Slider kann man das Signal&nbsp; $s_i(t)$&nbsp; um Vielfache von&nbsp; $\pm 0.25$&nbsp; nach oben und unten verschieben.
+
 
*Solche Signale treten zum Beispiel bei der binären oder mehrstufigen&nbsp; [[Digitalsignalübertragung/Redundanzfreie_Codierung#Quatern.C3.A4rsignal_mit_rc_.3D_0_und_Tern.C3.A4rsignal_mit_rc_.E2.89.88_0|Basisbandübertragung]]&nbsp; auf.
+
Die Lösung ist relativ einfach, wenn man die Spektralfunktionen betrachtet: &nbsp; Man erhält aus&nbsp; $Y_{\rm A}(f)$&nbsp; das Spektrum&nbsp; $Y(f) = X(f)$&nbsp; durch einen Tiefpass mit dem&nbsp; [[Lineare_zeitinvariante_Systeme/Systembeschreibung_im_Frequenzbereich#.C3.9Cbertragungsfunktion_-_Frequenzgang|Frequenzgang]]&nbsp; $H(f)$, der&nbsp;
*Im&nbsp; $\text{Beispiel 2}$&nbsp; des hier angegebenen Links erkennt man zum Beispiel die grafischen Darstellungen
+
[[Datei:P_ID1124__Sig_T_5_1_S5b_neu.png|right|frame|Frequenzbereichsdarstellung der Signalrekonstruktion]]
:* eines binären Signals&nbsp; $q(t)$,
+
*die tiefen Frequenzen unverfälscht durchlässt:
:* eines ternären Signals&nbsp; $s_3(t)$,
+
:$$H(f) = 1 \hspace{0.3cm}{\rm{f\ddot{u}r}} \hspace{0.3cm} |f| \le B_{\rm
:* eines quaternären Signals&nbsp; $s_4(t)$.  
+
  NF}\hspace{0.05cm},$$
 +
*die hohen Frequenzen vollständig unterdrückt:
 +
:$$H(f) = 0 \hspace{0.3cm}{\rm{f\ddot{u}r}} \hspace{0.3cm} |f| \ge f_{\rm A} - B_{\rm
 +
  NF}\hspace{0.05cm}.$$
 +
 +
Weiter ist aus der Grafik zu erkennen, dass der Frequenzgang&nbsp; $H(f)$&nbsp; im Bereich von&nbsp; $B_{\rm NF}$&nbsp; bis&nbsp; $f_{\rm A}–B_{\rm NF}$&nbsp; beliebig geformt sein kann,  
 +
*beispielsweise linear abfallend (gestrichelter Verlauf)  
 +
*oder auch rechteckförmig,
 +
 
 +
 
 +
solange die beiden oben genannten Bedingungen erfüllt sind.
 
<br clear=all>
 
<br clear=all>
'''(2)'''&nbsp; Rubrik&nbsp; <u>&bdquo;''M''&ndash;ASK / BPSK&rdquo;</u>&nbsp; &rArr; &nbsp; gültig für die Einstellungen &nbsp;$\rm (G)$&nbsp; ... &nbsp;$\rm (L)$:
+
==Das Abtasttheorem==
[[Datei:Gram_2_version2.png|right|frame|Signalform bei &bdquo;''M''&ndash;ASK / BPSK&rdquo;]]
+
<br>
*Die Mustersignale&nbsp; $s_i(t)$&nbsp; haben ebenfalls die Dauer&nbsp; $3T$&nbsp; und sind ähnlich aufgebaut wie bei der Rubrik&nbsp; '''(1)'''.
+
Die vollständige Rekonstruktion des Analogsignals&nbsp; $y(t)$&nbsp; aus dem abgetasteten Signal&nbsp; $y_{\rm A}(t) = x_{\rm A}(t)$&nbsp; ist nur möglich, wenn die Abtastrate&nbsp; $f_{\rm A}$&nbsp; entsprechend der Bandbreite&nbsp; $B_{\rm NF}$&nbsp; des Nachrichtensignals richtig gewählt wurde.  
* Im Unterschied zu&nbsp; '''(1)'''&nbsp; wird jede Rechteckfunktion&nbsp; $($Dauer $T)$&nbsp; durch eine Periode einer Sinusfunktionen ersetzt.
+
 
*Der angegebene Zahlenwert gibt hier die Amplitude des sinusförmigen Teilstücks an.
+
Aus der Grafik der&nbsp; [[Signaldarstellung/Zeitdiskrete_Signaldarstellung#Signalrekonstruktion|letzten Seite]]&nbsp; erkennt man, dass folgende Bedingung erfüllt sein muss:
*Bei negativem Vorzeichen wird aus dem &bdquo;Sinus&rdquo; die Funktion &bdquo;Minus&ndash;Sinus&rdquo;.
+
 
*Mit dem seitlichen Slider kann man die Amplitude von&nbsp; $s_i(t)$&nbsp; um Vielfache von&nbsp; $\pm 0.25$&nbsp; vergrößern oder verkleinern.
+
:$$f_{\rm A} - B_{\rm  NF} > B_{\rm  NF} \hspace{0.3cm}\Rightarrow \hspace{0.3cm}f_{\rm A} > 2 \cdot  B_{\rm   NF}\hspace{0.05cm}.$$
*Solche Signale können zum Beispiel bei der&nbsp;[[Digitalsignalübertragung/Trägerfrequenzsysteme_mit_kohärenter_Demodulation#M.E2.80.93stufiges_Amplitude_Shift_Keying_.28M.E2.80.93ASK.29|''M''&ndash;ASK]]&nbsp; (mehrstufiges ''Amplitude Shift Keying'')&nbsp; auftreten, ebenso bei&nbsp; [[Digitalsignalübertragung/Trägerfrequenzsysteme_mit_kohärenter_Demodulation#Binary_Phase_Shift_Keying_.28BPSK.29|BPSK]]&nbsp;(''Binary Phase Shift Keying''). 
+
<br clear=all>
+
{{BlaueBox|TEXT=
'''(3)'''&nbsp; Rubrik&nbsp; <u>&bdquo;Nur eine Frequenz&rdquo;</u>&nbsp; &rArr; &nbsp; gültig für die Einstellungen &nbsp;$\rm (M)$&nbsp; ... &nbsp;$\rm (R)$:
+
$\text{Abtasttheorem:}$&nbsp; Besitzt ein Analogsignal&nbsp; $x(t)$&nbsp; Spektralanteile im Bereich&nbsp; $\vert f \vert < B_{\rm NF}$, so kann dieses aus seinem abgetasteten Signal nur dann vollständig rekonstruiert werden, wenn die Abtastrate hinreichend groß ist:
[[Datei:Gram_3_version2.png|right|frame|Signalform bei &bdquo;Nur eine Frequenz&rdquo;]]
+
:$$f_{\rm A} ≥ 2 \cdot B_{\rm NF}.$$  
*Alle Mustersignale&nbsp; $s_i(t)$&nbsp; haben die Dauer&nbsp; $T$&nbsp; und sind jeweils Harmonische Schwingungen der Form
+
 
:$$s_i(t) = A \cdot \cos(2\pi \cdot f \cdot t - \varphi).$$
+
Für den Abstand zweier Abtastwerte muss demnach gelten:
*In der Grafik dargestellt ist der Fall:&nbsp; $A=0.75, \hspace{0.3cm}f= 2 \cdot f_0 =2/T, \hspace{0.3cm}\varphi= 90^\circ$ &nbsp; &rArr; &nbsp;   '''sinusförmiger Verlauf'''.
+
*Die Eigenschaft &bdquo;Nur eine Frequenz&rdquo; bezieht sich auf die einzelnen Mustersignale&nbsp; $s_i(t)$.
+
:$$T_{\rm A} \le \frac{1}{ 2 \cdot B_{\rm  NF} }\hspace{0.05cm}.$$}}
*Die&nbsp; $M$&nbsp; Signalformen eines Sets können durchaus unterschiedliche Frequenzen&nbsp; $f$&nbsp; aufweisen.
+
 
*Mit dem Slider lässt sich die Phase von&nbsp; $s_i(t)$&nbsp; um Vielfache von&nbsp; $\pm 22.5^\circ$&nbsp; in beide Richtungen variieren.
+
 
*Solche Harmonische haben für alle (analogen und digitalen) Nachrichtensysteme große Bedeutung.
+
Wird bei der Abtastung der größtmögliche Wert &nbsp; &nbsp; $T_{\rm A} = 1/(2B_{\rm NF})$&nbsp; herangezogen,
<br clear=all>
+
*so muss zur Signalrekonstruktion des Analogsignals aus seinen Abtastwerten
'''(4)'''&nbsp; Rubrik&nbsp; <u>&bdquo;Mehrere Frequenzen&rdquo;</u>&nbsp; &rArr; &nbsp; gültig für die Einstellungen &nbsp;$\rm (S)$&nbsp; ... &nbsp;$\rm (X)$:
+
*ein idealer, rechteckförmiger Tiefpass mit der Grenzfrequenz&nbsp; $f_{\rm G} = f_{\rm A}/2 = 1/(2T_{\rm A})$&nbsp; verwendet werden.
[[Datei:Gram_4_version2.png|right|frame|Signalform bei &bdquo;Mehrere Frequenzen&rdquo;]]
+
 
*Alle Mustersignale&nbsp; $s_i(t)$&nbsp; haben die Dauer&nbsp; $T$&nbsp; und sind meist Harmonische Schwingungen der Form
+
 
:$$s_i(t) = A \cdot \cos(2\pi \cdot f_k \cdot t)\hspace{0.6cm}\text{bzw.}\hspace{0.6cm}A \cdot \sin(2\pi \cdot f_k \cdot t).$$
+
{{GraueBox|TEXT=
*Der Parameter&nbsp; $k$&nbsp; gibt die Anzahl der Schwingungen innerhalb der Periodendauer&nbsp; $T$&nbsp; an. <br>In der Grafik dargestellt ist der Fall:&nbsp; $k=4$ &nbsp; &rArr; &nbsp;  vier Schwingungen innerhalb von&nbsp; $T$.
+
$\text{Beispiel 3:}$&nbsp; Die Grafik zeigt oben das auf&nbsp; $\pm\text{ 5 kHz}$&nbsp; begrenzte Spektrum&nbsp; $X(f)$&nbsp; eines Analogsignals, unten das Spektrum&nbsp; $X_{\rm A}(f)$&nbsp; des im Abstand&nbsp; $T_{\rm A} =\,\text{ 100 &micro;s}$&nbsp; abgetasteten Signals &nbsp; &nbsp; $f_{\rm A}=\,\text{ 10 kHz}$.  
*Möglich sind auch Mustersignale mit mehreren Frequenzen, zum Beispiel&nbsp; $($mit&nbsp; $k=0$:&nbsp; Gleichsignal$)$:
+
[[Datei:P_ID1125__Sig_T_5_1_S6_neu.png|right|frame|Abtasttheorem im Frequenzbereich]]
:$$s_i(t) = 1 \cdot \cos(2\pi \cdot f_0 \cdot t) - 0.5 \cdot \cos(2\pi \cdot f_2 \cdot t)-0.5 \cdot \cos(2\pi \cdot f_3 \cdot t).$$ 
+
<br>Zusätzlich eingezeichnet ist der Frequenzgang&nbsp; $H(f)$&nbsp; des Tiefpasses zur Signalrekonstruktion, dessen Grenzfrequenz&nbsp; $f_{\rm G} = f_{\rm A}/2 = 5\,\text{ kHz}$&nbsp; betragen muss.
*Die Eigenschaft &bdquo;Mehrere Frequenzen&rdquo; bezieht sich auf einzelne Mustersignale&nbsp; $s_i(t)$&nbsp; oder auch auf den gesamten Set&nbsp; $\{s_i(t)\}$.
+
 
*Mit dem Slider können die Frequenzen .... .
+
 
*Auch Signale dieser Rubrik haben für viele (analoge und digitale) Nachrichtensysteme eine große Bedeutung.  
+
*Mit jedem anderen&nbsp; $f_{\rm G}$–Wert ergäbe sich&nbsp; $Y(f) \neq X(f)$.  
 +
*Bei&nbsp; $f_{\rm G} < 5\,\text{ kHz}$&nbsp; fehlen die oberen&nbsp; $X(f)$–Anteile.
 +
* Bei&nbsp; $f_{\rm G} > 5\,\text{ kHz}$&nbsp; kommt es aufgrund von Faltungsprodukten zu unerwünschten Spektralanteilen in&nbsp; $Y(f)$.
 
<br clear=all>
 
<br clear=all>
 +
Wäre die Abtastung beim Sender mit einer Abtastrate&nbsp; $f_{\rm A} < 10\,\text{ kHz}$&nbsp;  erfolgt  &nbsp; ⇒  &nbsp;  $T_{\rm A} >100 \,{\rm &micro;  s}$, so wäre das Analogsignal&nbsp; $y(t) = x(t)$&nbsp; aus den Abtastwerten&nbsp; $y_{\rm A}(t)$&nbsp; auf keinen Fall rekonstruierbar.}}
  
  
 +
''Hinweis'': &nbsp; Zu der hier behandelten Thematik gibt es ein interaktives Applet: &nbsp;
 +
[[Applets:Abtastung_periodischer_Signale_und_Signalrekonstruktion_(Applet)|Abtastung analoger Signale und Signalrekonstruktion]]
 
<br clear=all>
 
<br clear=all>
 
==Versuchsdurchführung==
 
==Versuchsdurchführung==

Version vom 3. Februar 2020, 16:07 Uhr

Applet in neuem Tab öffnen

Programmbeschreibung


Das Applet verdeutlicht das Gram–Schmidt–Verfahren. Dieses ermöglicht, eine Menge  $\{s_1(t), \hspace{0.05cm} \text{...} \hspace{0.05cm} , s_M(t)\}$  energiebegrenzter Signale mit Hilfe von   $N \le M$  orthonormalen Basisfunktionen   $\varphi_1(t), \hspace{0.05cm} \text{...} \hspace{0.05cm} , \varphi_N(t)$  in folgender Form darzustellen:

$$s_i(t) = \sum\limits_{j = 1}^{N}s_{ij} \cdot \varphi_j(t) , \hspace{0.3cm}i = 1,\hspace{0.05cm} \text{...}\hspace{0.1cm} , M, \hspace{0.3cm}j = 1,\hspace{0.05cm} \text{...} \hspace{0.1cm}, N \hspace{0.05cm}.$$

Der vektorielle Repräsentant der Musterfunktion  $s_1(t)$  lautet dann: $$\mathbf{s}_i = \big( s_{i1}\hspace{0.05cm}, \hspace{0.3cm}s_{i2}\hspace{0.05cm},\hspace{0.05cm} \text{...}\hspace{0.05cm},\hspace{0.05cm} s_{iN} \big ).$$

Das Applet zeigt alle Grafiken, die zum Verständnis des Gram–Schmidt–Verfahrens erforderlich sind, und als jeweiliges Ergebnis

  • die 2D–Darstellung der  $M$  vektoriellen Repräsentanten, falls  $N=2$,
  • die 3D–Darstellung der  $M$  vektoriellen Repräsentanten, falls  $N=3$.

Theoretischer Hintergrund

Beschreibung der Abtastung im Zeitbereich

Zur Zeitdiskretisierung des zeitkontinuierlichen Signals  $x(t)$

Im Folgenden verwenden wir für die Beschreibung der Abtastung folgende Nomenklatur:

  • Das zeitkontinuierliche Signal sei  $x(t)$.
  • Das in äquidistanten Abständen  $T_{\rm A}$  abgetastete zeitdiskretisierte Signal sei  $x_{\rm A}(t)$.
  • Außerhalb der Abtastzeitpunkte  $\nu \cdot T_{\rm A}$  gilt stets  $x_{\rm A}(t) = 0$.
  • Die Laufvariable  $\nu$  sei  ganzzahlig:     $\nu \in \mathbb{Z} = \{\hspace{0.05cm} \text{...}\hspace{0.05cm} , –3, –2, –1, \hspace{0.2cm}0, +1, +2, +3, \text{...} \hspace{0.05cm}\} $.
  • Dagegen ergibt sich zu den äquidistanten Abtastzeitpunkten mit der Konstanten  $K$:
$$x_{\rm A}(\nu \cdot T_{\rm A}) = K \cdot x(\nu \cdot T_{\rm A})\hspace{0.05cm}.$$

Die Konstante hängt von der Art der Zeitdiskretisierung ab. Für die obige Skizze gilt  $K = 1$.

$\text{Definitionen:}$ 

Der  Diracpuls (im Zeitbereich)  besteht aus unendlich vielen Diracimpulsen, jeweils im gleichen Abstand  $T_{\rm A}$  und alle mit gleichem Impulsgewicht  $T_{\rm A}$:

$$p_{\delta}(t) = \sum_{\nu = - \infty }^{+\infty} T_{\rm A} \cdot \delta(t- \nu \cdot T_{\rm A} )\hspace{0.05cm}.$$
  • Unter  Abtastung  verstehen wir im Folgenden die Multiplikation des zeitkontinuierlichen Signals  $x(t)$  mit einem solchen  Diracpuls:
$$x_{\rm A}(t) = x(t) \cdot p_{\delta}(t)\hspace{0.05cm}.$$

Aufgrund dieser Definition ergeben sich für das abgetastete Signal folgende Eigenschaften:

  • Das abgetastete Signal zum betrachteten Zeitpunkt  $(\nu \cdot T_{\rm A})$  ist gleich  $T_{\rm A} \cdot x(\nu \cdot T_{\rm A}) · \delta (0)$.
  • Da die Diracfunktion  $\delta (t)$  zur Zeit  $t = 0$  unendlich ist, sind eigentlich alle Signalwerte  $x_{\rm A}(\nu \cdot T_{\rm A})$  ebenfalls unendlich groß.
  • Somit ist auch der auf der letzten Seite eingeführte Faktor  $K$  eigentlich unendlich groß.
  • Zwei Abtastwerte  $x_{\rm A}(\nu_1 \cdot T_{\rm A})$  und  $x_{\rm A}(\nu_2 \cdot T_{\rm A})$  unterscheiden sich jedoch im gleichen Verhältnis wie die Signalwerte  $x(\nu_1 \cdot T_{\rm A})$  und  $x(\nu_2 \cdot T_{\rm A})$.
  • Die Abtastwerte von  $x(t)$  erscheinen in den Impulsgewichten der Diracfunktionen:
$$x_{\rm A}(t) = \sum_{\nu = - \infty }^{+\infty} T_{\rm A} \cdot x(\nu \cdot T_{\rm A})\cdot \delta (t- \nu \cdot T_{\rm A} )\hspace{0.05cm}.$$
  • Die zusätzliche Multiplikation mit  $T_{\rm A}$  ist erforderlich, damit  $x(t)$  und  $x_{\rm A}(t)$  gleiche Einheit besitzen. Beachten Sie hierbei, dass  $\delta (t)$  selbst die Einheit „1/s” aufweist.


Die folgenden Seiten werden zeigen, dass diese gewöhnungsbedürftigen Gleichungen durchaus zu sinnvollen Ergebnissen führen, wenn man sie konsequent anwendet.


Diracpuls im Zeit- und im Frequenzbereich


$\text{Satz:}$  Entwickelt man den  Diracpuls  in eine  Fourierreihe  und transformiert diese unter Anwendung des  Verschiebungssatzes  in den Frequenzbereich, so ergibt sich folgende Korrespondenz:

$$p_{\delta}(t) = \sum_{\nu = - \infty }^{+\infty} T_{\rm A} \cdot \delta(t- \nu \cdot T_{\rm A} )\hspace{0.2cm}\circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\, \hspace{0.2cm} P_{\delta}(f) = \sum_{\mu = - \infty }^{+\infty} \delta (f- \mu \cdot f_{\rm A} ).$$

Hierbei gibt  $f_{\rm A} = 1/T_{\rm A}$  den Abstand zweier benachbarter Diraclinien im Frequenzbereich an.


$\text{Beweis:}$  Die Herleitung der hier angegebenen Spektralfunktion  $P_{\delta}(f)$  geschieht in mehreren Schritten:

(1)   Da  $p_{\delta}(t)$  periodisch mit dem konstanten Abstand  $T_{\rm A}$  zwischen zwei Diraclinien ist, kann die  (komplexe) Fourierreihendarstellung  angewendet werden:

$$p_{\delta}(t) = \sum_{\mu = - \infty }^{+\infty} D_{\mu} \cdot {\rm e}^{\hspace{0.05cm}{\rm j} \hspace{0.05cm} \cdot 2 \hspace{0.05cm} \pi \hspace{0.05cm}\cdot \hspace{0.05cm}\mu \hspace{0.05cm}\cdot \hspace{0.05cm}t/T_{\rm A} } \hspace{0.3cm}{\rm mit}\hspace{0.3cm} D_{\mu} = \frac{1}{T_{\rm A} } \cdot \int_{-T_{\rm A}/2 }^{+T_{\rm A}/2}p_{\delta}(t) \cdot {\rm e}^{- {\rm j} \hspace{0.05cm}\cdot \hspace{0.05cm} 2 \pi \hspace{0.05cm}\cdot \hspace{0.05cm}\mu \hspace{0.05cm} \cdot \hspace{0.05cm}t/T_{\rm A} }\hspace{0.1cm} {\rm d}t\hspace{0.05cm}.$$

(2)   Im Bereich von  $–T_{\rm A}/2$  bis  $+T_{\rm A}/2$  gilt für den Diracpuls im Zeitbereich:   $p_{\delta}(t) = T_{\rm A} \cdot \delta(t)$. Damit kann man für die komplexen Fourierkoeffizienten schreiben:  

$$D_{\mu} = \int_{-T_{\rm A}/2 }^{+T_{\rm A}/2}{\delta}(t) \cdot {\rm e}^{- {\rm j} \hspace{0.05cm}\cdot \hspace{0.05cm} 2 \pi \hspace{0.05cm}\cdot \hspace{0.05cm}\mu \hspace{0.05cm} \cdot \hspace{0.05cm}t/T_{\rm A} }\hspace{0.1cm} {\rm d}t\hspace{0.05cm}.$$

(3)   Unter Berücksichtigung der Tatsache, dass für  $t \neq 0$  der Diracimpuls Null ist und für  $t = 0$  der komplexe Drehfaktor gleich  $1$, gilt weiter:

$$D_{\mu} = \int_{- T_{\rm A}/2 }^{+T_{\rm A}/2}{\delta}(t) \hspace{0.1cm} {\rm d}t = 1\hspace{0.5cm}{\Rightarrow}\hspace{0.5cm} p_{\delta}(t) = \sum_{\mu = - \infty }^{+\infty} {\rm e}^{ {\rm j} \hspace{0.05cm} \cdot 2 \hspace{0.05cm} \pi \cdot \hspace{0.05cm}\mu \hspace{0.05cm}\cdot \hspace{0.05cm}t/T_{\rm A} }\hspace{0.05cm}. $$

(4)   Der  Verschiebungssatz im Frequenzbereich  lautet mit  $f_{\rm A} = 1/T_{\rm A}$:

$${\rm e}^{ {\rm j} \hspace{0.05cm} \hspace{0.05cm} \cdot 2 \hspace{0.05cm} \pi \hspace{0.05cm}\cdot \hspace{0.05cm}\mu \hspace{0.05cm}\cdot \hspace{0.05cm} f_{\rm A}\hspace{0.05cm}\cdot \hspace{0.05cm}t}\hspace{0.2cm}\circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\, \hspace{0.2cm} \delta (f- \mu \cdot f_{\rm A} )\hspace{0.05cm}.$$

(5)   Wendet man das Ergebnis auf jeden einzelnen Summanden an, so erhält man schließlich:

$$P_{\delta}(f) = \sum_{\mu = - \infty }^{+\infty} \delta (f- \mu \cdot f_{\rm A} )\hspace{0.05cm}.$$
q.e.d.


Das Ergebnis besagt:

  • Der Diracpuls  $p_{\delta}(t)$  im Zeitbereich besteht aus unendlich vielen Diracimpulsen, jeweils im gleichen Abstand  $T_{\rm A}$  und alle mit gleichem Impulsgewicht  $T_{\rm A}$.
  • Die Fouriertransformierte von  $p_{\delta}(t)$  ergibt wiederum einen Diracpuls, aber nun im Frequenzbereich   ⇒   $P_{\delta}(f)$.
  • $P_{\delta}(f)$  besteht ebenfalls aus unendlich vielen Diracimpulsen, nun aber im jeweiligen Abstand  $f_{\rm A} = 1/T_{\rm A}$  und alle mit dem Impulsgewicht  $1$.
  • Die Abstände der Diraclinien in der Zeit– und Frequenzbereichsdarstellung folgen demnach dem  Reziprozitätsgesetz:  
$$T_{\rm A} \cdot f_{\rm A} = 1 \hspace{0.05cm}.$$


Diracpuls im Zeit- und Frequenzbereich

$\text{Beispiel 1:}$  Die Grafik verdeutlicht die obigen Aussagen für

  • $T_{\rm A} = 50\,{\rm µs}$,
  • $f_{\rm A} = 1/T_{\rm A} = 20\,\text{kHz}$ .


Man erkennt aus dieser Skizze auch die unterschiedlichen Impulsgewichte von  $p_{\delta}(t)$  und  $P_{\delta}(f)$.

Prinzip und Motivation


Viele Nachrichtensignale sind analog und damit gleichzeitig  zeitkontinuierlich  und  wertkontinuierlich. Soll ein solches Analogsignal mittels eines Digitalsystems übertragen werden, so sind folgende Vorverarbeitungsschritte erforderlich:

  • die  Abtastung  des Nachrichtensignals  $x(t)$, die zweckmäßigerweise – aber nicht notwendigerweise – zu äquidistanten Zeitpunkten erfolgt   ⇒   Zeitdiskretisierung,
  • die  Quantisierung  der Abtastwerte, um so die Anzahl  $M$  der möglichen Werte auf einen endlichen Wert zu begrenzen   ⇒   Wertdiskretisierung.


Die Quantisierung wird erst im Kapitel  Pulscodemodulation  des Buches „Modulationsverfahren” im Detail behandelt.

Zur Zeitdiskretisierung des zeitkontinuierlichen Signals  $x(t)$

Im Folgenden verwenden wir für die Beschreibung der Abtastung folgende Nomenklatur:

  • Das zeitkontinuierliche Signal sei  $x(t)$.
  • Das in äquidistanten Abständen  $T_{\rm A}$  abgetastete zeitdiskretisierte Signal sei  $x_{\rm A}(t)$.
  • Außerhalb der Abtastzeitpunkte  $\nu \cdot T_{\rm A}$  gilt stets  $x_{\rm A}(t) = 0$.
  • Die Laufvariable  $\nu$  sei  ganzzahlig:     $\nu \in \mathbb{Z} = \{\hspace{0.05cm} \text{...}\hspace{0.05cm} , –3, –2, –1, \hspace{0.2cm}0, +1, +2, +3, \text{...} \hspace{0.05cm}\} $.
  • Dagegen ergibt sich zu den äquidistanten Abtastzeitpunkten mit der Konstanten  $K$:
$$x_{\rm A}(\nu \cdot T_{\rm A}) = K \cdot x(\nu \cdot T_{\rm A})\hspace{0.05cm}.$$

Die Konstante hängt von der Art der Zeitdiskretisierung ab. Für die obige Skizze gilt  $K = 1$.


Frequenzbereichsdarstellung


Zum Spektrum des abgetasteten Signals  $x_{\rm A}(t)$  kommt man durch Anwendung des  Faltungssatzes. Dieser besagt, dass der Multiplikation im Zeitbereich die Faltung im Spektralbereich entspricht:

$$x_{\rm A}(t) = x(t) \cdot p_{\delta}(t)\hspace{0.2cm}\circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\, \hspace{0.2cm} X_{\rm A}(f) = X(f) \star P_{\delta}(f)\hspace{0.05cm}.$$

Aus dem Spektrum  $X(f)$  wird durch Faltung mit der um  $\mu \cdot f_{\rm A}$  verschobenen Diraclinie:

$$X(f) \star \delta (f- \mu \cdot f_{\rm A} )= X (f- \mu \cdot f_{\rm A} )\hspace{0.05cm}.$$

Wendet man dieses Ergebnis auf alle Diraclinien des Diracpulses an, so erhält man schließlich:

$$X_{\rm A}(f) = X(f) \star \sum_{\mu = - \infty }^{+\infty} \delta (f- \mu \cdot f_{\rm A} ) = \sum_{\mu = - \infty }^{+\infty} X (f- \mu \cdot f_{\rm A} )\hspace{0.05cm}.$$

$\text{Fazit:}$  Die Abtastung des analogen Zeitsignals  $x(t)$  in äquidistanten Abständen  $T_{\rm A}$  führt im Spektralbereich zu einer  periodischen Fortsetzung  von  $X(f)$  mit dem Frequenzabstand  $f_{\rm A} = 1/T_{\rm A}$.


$\text{Beispiel 2:}$  Die obere Grafik zeigt  (schematisch!)  das Spektrum  $X(f)$  eines Analogsignals  $x(t)$, das Frequenzen bis  $5 \text{ kHz}$  beinhaltet.

Spektrum des abgetasteten Signals

Tastet man das Signal mit der Abtastrate  $f_{\rm A}\,\text{ = 20 kHz}$, also im jeweiligen Abstand  $T_{\rm A}\, = {\rm 50 \, µs}$  ab, so erhält man das unten skizzierte periodische Spektrum  $X_{\rm A}(f)$.

  • Da die Diracfunktionen unendlich schmal sind, beinhaltet das abgetastete Signal  $x_{\rm A}(t)$  auch beliebig hochfrequente Anteile.
  • Dementsprechend ist die Spektralfunktion  $X_{\rm A}(f)$  des abgetasteten Signals bis ins Unendliche ausgedehnt.


Signalrekonstruktion


Die Signalabtastung ist bei einem digitalen Übertragungssystem kein Selbstzweck, sondern sie muss irgendwann wieder rückgängig gemacht werden. Betrachten wir zum Beispiel das folgende System:

Signalabtastung und Signalrekonstruktion
  • Das Analogsignal  $x(t)$  mit der Bandbreite  $B_{\rm NF}$  wird wie oben beschrieben abgetastet.
  • Am Ausgang eines idealen Übertragungssystems liegt das ebenfalls zeitdiskrete Signal  $y_{\rm A}(t) = x_{\rm A}(t)$  vor.
  • Die Frage ist nun, wie der Block  Signalrekonstruktion  zu gestalten ist, damit auch  $y(t) = x(t)$  gilt.


Die Lösung ist relativ einfach, wenn man die Spektralfunktionen betrachtet:   Man erhält aus  $Y_{\rm A}(f)$  das Spektrum  $Y(f) = X(f)$  durch einen Tiefpass mit dem  Frequenzgang  $H(f)$, der 

Frequenzbereichsdarstellung der Signalrekonstruktion
  • die tiefen Frequenzen unverfälscht durchlässt:
$$H(f) = 1 \hspace{0.3cm}{\rm{f\ddot{u}r}} \hspace{0.3cm} |f| \le B_{\rm NF}\hspace{0.05cm},$$
  • die hohen Frequenzen vollständig unterdrückt:
$$H(f) = 0 \hspace{0.3cm}{\rm{f\ddot{u}r}} \hspace{0.3cm} |f| \ge f_{\rm A} - B_{\rm NF}\hspace{0.05cm}.$$

Weiter ist aus der Grafik zu erkennen, dass der Frequenzgang  $H(f)$  im Bereich von  $B_{\rm NF}$  bis  $f_{\rm A}–B_{\rm NF}$  beliebig geformt sein kann,

  • beispielsweise linear abfallend (gestrichelter Verlauf)
  • oder auch rechteckförmig,


solange die beiden oben genannten Bedingungen erfüllt sind.

Das Abtasttheorem


Die vollständige Rekonstruktion des Analogsignals  $y(t)$  aus dem abgetasteten Signal  $y_{\rm A}(t) = x_{\rm A}(t)$  ist nur möglich, wenn die Abtastrate  $f_{\rm A}$  entsprechend der Bandbreite  $B_{\rm NF}$  des Nachrichtensignals richtig gewählt wurde.

Aus der Grafik der  letzten Seite  erkennt man, dass folgende Bedingung erfüllt sein muss:

$$f_{\rm A} - B_{\rm NF} > B_{\rm NF} \hspace{0.3cm}\Rightarrow \hspace{0.3cm}f_{\rm A} > 2 \cdot B_{\rm NF}\hspace{0.05cm}.$$

$\text{Abtasttheorem:}$  Besitzt ein Analogsignal  $x(t)$  Spektralanteile im Bereich  $\vert f \vert < B_{\rm NF}$, so kann dieses aus seinem abgetasteten Signal nur dann vollständig rekonstruiert werden, wenn die Abtastrate hinreichend groß ist:

$$f_{\rm A} ≥ 2 \cdot B_{\rm NF}.$$

Für den Abstand zweier Abtastwerte muss demnach gelten:

$$T_{\rm A} \le \frac{1}{ 2 \cdot B_{\rm NF} }\hspace{0.05cm}.$$


Wird bei der Abtastung der größtmögliche Wert   ⇒   $T_{\rm A} = 1/(2B_{\rm NF})$  herangezogen,

  • so muss zur Signalrekonstruktion des Analogsignals aus seinen Abtastwerten
  • ein idealer, rechteckförmiger Tiefpass mit der Grenzfrequenz  $f_{\rm G} = f_{\rm A}/2 = 1/(2T_{\rm A})$  verwendet werden.


$\text{Beispiel 3:}$  Die Grafik zeigt oben das auf  $\pm\text{ 5 kHz}$  begrenzte Spektrum  $X(f)$  eines Analogsignals, unten das Spektrum  $X_{\rm A}(f)$  des im Abstand  $T_{\rm A} =\,\text{ 100 µs}$  abgetasteten Signals   ⇒   $f_{\rm A}=\,\text{ 10 kHz}$.

Abtasttheorem im Frequenzbereich


Zusätzlich eingezeichnet ist der Frequenzgang  $H(f)$  des Tiefpasses zur Signalrekonstruktion, dessen Grenzfrequenz  $f_{\rm G} = f_{\rm A}/2 = 5\,\text{ kHz}$  betragen muss.


  • Mit jedem anderen  $f_{\rm G}$–Wert ergäbe sich  $Y(f) \neq X(f)$.
  • Bei  $f_{\rm G} < 5\,\text{ kHz}$  fehlen die oberen  $X(f)$–Anteile.
  • Bei  $f_{\rm G} > 5\,\text{ kHz}$  kommt es aufgrund von Faltungsprodukten zu unerwünschten Spektralanteilen in  $Y(f)$.


Wäre die Abtastung beim Sender mit einer Abtastrate  $f_{\rm A} < 10\,\text{ kHz}$  erfolgt   ⇒   $T_{\rm A} >100 \,{\rm µ s}$, so wäre das Analogsignal  $y(t) = x(t)$  aus den Abtastwerten  $y_{\rm A}(t)$  auf keinen Fall rekonstruierbar.


Hinweis:   Zu der hier behandelten Thematik gibt es ein interaktives Applet:   Abtastung analoger Signale und Signalrekonstruktion

Versuchsdurchführung


Aufgaben 2D-Gauss.png

Noch anpassen

  • Wählen Sie zunächst die Nummer  (1, ...)  der zu bearbeitenden Aufgabe.
  • Eine Aufgabenbeschreibung wird angezeigt. Die Parameterwerte sind angepasst.
  • Lösung nach Drücken von „Musterlösung”.


Die Nummer 0 entspricht einem „Reset”:

  • Gleiche Einstellung wie beim Programmstart.
  • Ausgabe eines „Reset–Textes” mit weiteren Erläuterungen zum Applet.

Bis hierher

(1)  Es gilt die Einstellung  $\rm A$.  Interpretieren Sie die ausgegebenen Grafiken.  Wählen Sie hierfür „Einzelschritt”.

  •  Einstellung  $\rm A$  beschreibt das $\text{Beispiel 2}$  im Theorieteil. Die Basisfunktion  $\varphi_1(t)$  ist identisch mit dem Signal  $s_1(t)$,  aber mit Signalenergie  $E=1$.
  •  Es gibt hier nur  $N=3$  Basisfunktionen, da die Hilfsfunktion  $\theta_3(t)$  identisch Null ist.
  •  Die vektoriellen Repräsentanten der Signale  $s_1(t)$,  ... , $s_4(t)$  können im 3D–Vektorraum abgelesen werden;  Beispiel:  $\mathbf{s}_4 = (-1.444, \hspace{0.15cm} -0.408, \hspace{0.15cm} +0.707)$.

(2)  Interpretieren Sie die ausgegebenen Grafiken für die Einstellung  $\rm B$.  Wählen Sie hierfür und bei den weiteren Aufgaben „Gesamtdarstellung”.

  •  Auch hier gibt es  $N=3$  Basisfunktionen.  Bei Änderung auf  $s_4 = (-1, \hspace{0.15cm} -1, \hspace{0.25cm} 0)$  nur mehr  $N=2$.

(3)  Bei der Einstellung  $\rm C$  ist die Reihenfolge der Signale gegenüber  $\rm B$  vertauscht.  Wie wirkt sich das auf die Basisfunktionen aus?

  •  Auch hier gibt es  $N=3$  Basisfunktionen, aber nun andere:  Nämlich  $\varphi_1(t) = s_1(t)$,  $\varphi_2(t) = s_2(t)$,  $\varphi_3(t) = s_3(t)$.

(4)  Die  $M=4$  Signale der Einstellung  $\rm D$  lassen sich durch nur  $N=2$  Basisfunktionen ausdrücken?  Begründen Sie dieses Ergebnis.

  •  Es gilt  $s_3(t) = s_1(t)/4 - s_2(t)/2$  und  $s_4(t) = -s_1(t) - s_2(t)$.  Das heißt:  $s_3(t)$  und  $s_4(t)$  liefern keine neuen Basisfunktionen.

(5)  Interpretieren Sie die ausgegebenen Grafiken für die Einstellung  $\rm E$  im Vergleich zur Einstellung  $\rm D$.

  •  Bei der Einstellung  $\rm E$  ist die Reihenfolge der Signale gegenüber der Einstellung   $\rm D$  vertauscht. Ähnlich wie zwischen  $\rm B$  und  $\rm C$.
  •  Auch diese  $M=4$  Signale lassen sich somit durch nur  $N=2$  Basisfunktionen ausdrücken, aber durch andere als in der Aufgabe  (4).

(6)  Welches Ergebnis liefern die vier Signale gemäß der Einstellung  $\rm F$?

  •  Die die Signale  $s_1(t)$, ... , $s_4(t)$  basieren alle auf einer einzigen Basisfunktion   $\varphi_1(t)$, die formgleich mit  $s_1(t)$  ist.  Es gilt  $N=1$.
  •  Die vektoriellen Repräsentanten der Signale  $s_1(t)$,  ... , $s_4(t)$  sind  $\pm 0.866$  und  $\pm 1.732$.  Sie liegen alle auf einer Linie.

(7)  Es gilt nun die „M–ASK / BPSK”–Einstellung  $\rm G$.  Interpretieren Sie das Ergebnis und versuchen Sie, einen Zusammenhang zu einer früheren Aufgabe herzustellen.

  •  Vergleicht man die angegebenen Zahlenwerte, so erkennt man, dass eine ähnliche Konstellation betrachtet wird wie bei der „Basisband”–Einstellung  $\rm A$.
  •  Der einzige Unterschied ist, dass nun alle Energien nur halb so groß sind wie vorher.  Bezüglich der Amplituden wirkt sich das um den Faktor  $\sqrt{2}$  aus.
  •  Somit ist nun der vektorielle Repräsentant des unteren Signals    $\mathbf{s}_4 = (-1.021, \hspace{0.15cm} -0.289, \hspace{0.15cm} +0.500)$  anstelle von  $\mathbf{s}_4 = (-1.444, \hspace{0.15cm} -0.408, \hspace{0.15cm} +0.707)$.
  •  Bei der Einstellung  $\rm H$  sind gegenüber  $\rm G$  alle Amplituden verdoppelt. Somit ergibt sich hier wieder  $\mathbf{s}_4 = (-1.444, \hspace{0.15cm} -0.408, \hspace{0.15cm} +0.707)$.

(8)  Es gelte die „M–ASK / BPSK”–Einstellung  $\rm I$.  Interpretieren Sie das Ergebnis.  Versuchen Sie wieder, einen Zusammenhang zu einer früheren Aufgabe herzustellen.

  •  Hier wird eine ähnliche Konstellation betrachtet wird wie bei der „Basisband”–Einstellung  $\rm C$, aber nun mit nur halb so großen Energien.
  •  Somit ist nun der vektorielle Repräsentant des unteren Signals    $\mathbf{s}_4 = (+0.707, \hspace{0.15cm} -0.707, \hspace{0.15cm} 0.000)$  anstelle von  $\mathbf{s}_4 = (+1.000, \hspace{0.15cm} -1.000, \hspace{0.15cm} 0.000)$.

(9)  Wählen Sie die Einstellungen  $M=4 \text{, nach Spalt–TP, }T_{\rm E}/T = 1$,  $10 \cdot \lg \ E_{\rm B}/N_0 = 10 \ {\rm dB}$  und  $12 \ {\rm dB}$.  Interpretieren Sie die Ergebnisse.

  •  Es gibt nun drei Augenöffnungen. Gegenüber  (5)  ist also  $ö_{\rm norm}$  um den Faktor  $3$  kleiner,  $\sigma_{\rm norm}$  dagegen nur um etwa den Faktor  $\sqrt{5/9)} \approx 0.75$.
  •  Für  $10 \cdot \lg \ E_{\rm B}/N_0 = 10 \ {\rm dB}$  ergibt sich nun die Fehlerwahrscheinlichkeit  $p_{\rm U} \approx 2.27\%$  und für  $10 \cdot \lg \ E_{\rm B}/N_0 = 12 \ {\rm dB}$  nur mehr  $0.59\%$.

(10)  Für die restlichen Aufgaben gelte stets  $10 \cdot \lg \ E_{\rm B}/N_0 = 12 \ {\rm dB}$. Betrachten Sie das Augendiagramm für  $M=4 \text{, CRO–Nyquist, }r_f = 0.5$.

  •  In  $d_{\rm S}(t)$  müssen alle „Fünf–Symbol–Kombinationen” enthalten sein   ⇒   mindestens  $4^5 = 1024$  Teilstücke   ⇒   maximal  $1024$  unterscheidbare Linien.
  •  Alle  $1024$  Augenlinien gehen bei  $t=0$  durch nur vier Punkte:  $ö_{\rm norm}= 0.333$. $\sigma_{\rm norm} = 0.143$  ist etwas größer als in  (9)  ⇒   ebenso  $p_{\rm U} \approx 1\%$.

(11)  Wählen Sie die Einstellungen  $M=4 \text{, nach Gauß–TP, }f_{\rm G}/R_{\rm B} = 0.48$  und variieren Sie  $f_{\rm G}/R_{\rm B}$.   Interpretieren Sie die Ergebnisse.

  •  $f_{\rm G}/R_{\rm B}=0.48$  führt zur minimalen Fehlerwahrscheinlichkeit  $p_{\rm U} \approx 0.21\%$.  Kompromiss zwischen  $ö_{\rm norm}= 0.312$  und  $\sigma_{\rm norm}= 0.109$.
  •  Bei zu kleiner Grenzfrequenz dominieren die Impulsinterferenzen.  Beispiel:  $f_{\rm G}/R_{\rm B}= 0.3$:  $ö_{\rm norm}= 0.157; $ $\sigma_{\rm norm}= 0.086$  ⇒    $p_{\rm U} \approx 3.5\%$.
  •  Bei zu großer Grenzfrequenz dominiert das Rauschen.  Beispiel:  $f_{\rm G}/R_{\rm B}= 1.0$:  $ö_{\rm norm}= 0.333; $ $\sigma_{\rm norm}= 0.157$  ⇒    $p_{\rm U} \approx 1.7\%$.
  •  Aus dem Vergleich mit  (9)  erkennt man:  Bei Quaternärcodierung ist es günstiger, Impulsinterferenzen zuzulassen.

(12)  Welche Unterschiede zeigt das Auge für  $M=3 \text{ (AMI-Code), nach Gauß–TP, }f_{\rm G}/R_{\rm B} = 0.48$  gegenüber dem vergleichbaren Binärsystem? Interpretation.

  •  Der Detektionsgrundimpuls  $g_d(t)$  ist in beiden Fällen gleich. Die Abtastwerte sind jeweils  $g_0 = 0.771, \ g_1 = 0.114$.
  •  Beim AMI–Code gibt es zwei Augenöffnungen mit je  $ö_{\rm norm}= 1/2 \cdot (g_0 -3 \cdot g_1) = 0.214$.  Beim Binärcode:  $ö_{\rm norm}= g_0 -2 \cdot g_1 = 0.543$.
  •  Die AMI–Folge besteht zu 50% aus Nullen. Die Symbole  $+1$  und  $-1$  wechseln sich ab   ⇒   es gibt keine lange  $+1$–Folge und keine lange  $-1$–Folge.
  •  Darin liegt der einzige Vorteil des AMI–Codes:  Dieser kann auch bei einem gleichsignalfreien Kanal   ⇒   $H_{\rm K}(f= 0)=0$  angewendet werden.

(13)  Gleiche Einstellung wie in  (12), zudem  $10 \cdot \lg \ E_{\rm B}/N_0 = 12 \ {\rm dB}$. Analysieren Sie die Fehlerwahrscheinlichkeit des AMI–Codes.

  •  Trotz kleinerem  $\sigma_{\rm norm} = 0.103$  hat der AMI–Code eine höhere Fehlerwahrscheinlichkeit  $p_{\rm U} \approx 2\%$  als der Binärcode:  $\sigma_{\rm norm} = 0.146, \ p_{\rm U} \approx \cdot 10^{-4}.$
  •  Für  $f_{\rm G}/R_{\rm B}<0.34$  ergibt sich ein geschlossenes Auge  $(ö_{\rm norm}= 0)$  ⇒    $p_{\rm U} =50\%$. Beim Binärcode:  Für  $f_{\rm G}/R_{\rm B}>0.34$  ist das Auge geöffnet.

(14)  Welche Unterschiede zeigt das Auge für  $M=3 \text{ (Duobinärcode), nach Gauß–TP, }f_{\rm G}/R_{\rm B} = 0.30$  gegenüber dem vergleichbaren Binärsystem?

  •  Redundanzfreier Binärcode:  $ö_{\rm norm}= 0.096, \ \sigma_{\rm norm} = 0.116 \ p_{\rm U} \approx 20\% $       Duobinärcode:  $ö_{\rm norm}= 0.167, \ \sigma_{\rm norm} = 0.082 \ p_{\rm U} \approx 2\% $.
  • Insbesondere bei kleinem  $f_{\rm G}/R_{\rm B}$  liefert der Duobinärcode gute Ergebnisse, da die Übergänge von  $+1$  nach  $-1$  (und umgekehrt) im Auge fehlen.
  • Selbst mit  $f_{\rm G}/R_{\rm B}=0.2$  ist das Auge noch geöffnet. Im Gegensatz zum AMI–Code  ist aber „Duobinär” bei gleichsignalfreiem Kanal nicht anwendbar.

Zur Handhabung des Applets


Anleitung Auge.png

    (A)     Auswahl:   Codierung
                   (binär,  quaternär,  AMI–Code,  Duobinärcode)

    (B)     Auswahl:   Detektionsgrundimpuls
                    (nach Gauß–TP,  CRO–Nyquist,  nach Spalt–TP}

    (C)     Prametereingabe zu  (B)
                   (Grenzfrequenz,  Rolloff–Faktor,  Rechteckdauer)

    (D)     Steuerung der Augendiagrammdarstellung
                   (Start,  Pause/Weiter,  Einzelschritt,  Gesamt,  Reset)

    (E)     Geschwindigkeit der Augendiagrammdarstellung

    (F)     Darstellung:  Detektionsgrundimpuls  $g_d(t)$

    (G)     Darstellung:  Detektionsnutzsignal  $d_{\rm S}(t - \nu \cdot T)$

    (H)     Darstellung:  Augendiagramm im Bereich  $\pm T$

    ( I )     Numerikausgabe:  $ö_{\rm norm}$  (normierte Augenöffnung)

    (J)     Prametereingabe  $10 \cdot \lg \ E_{\rm B}/N_0$  für  (K)

    (K)     Numerikausgabe:  $\sigma_{\rm norm}$  (normierter Rauscheffektivwert)

    (L)     Numerikausgabe:  $p_{\rm U}$  (ungünstigste Fehlerwahrscheinlichkeit)

    (M)     Bereich für die Versuchsdurchführung:   Aufgabenauswahl

    (N)     Bereich für die Versuchsdurchführung:   Aufgabenstellung

    (O)     Bereich für die Versuchsdurchführung:   Musterlösung einblenden

Ü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 2008 von  Thomas Großer  im Rahmen einer Werkstudententätigkeit mit „FlashMX–Actionscript” erstellt (Betreuer:  Günter Söder).
  • 2019 wurde das Programm von  Carolin Mirschina  im Rahmen einer Werkstudententätigkeit auf „HTML5” umgesetzt und neu gestaltet (Betreuer:  Tasnád Kernetzky).


Die Umsetzung dieses Applets auf HTML 5 wurde durch  Studienzuschüsse  der Fakultät EI der TU München finanziell unterstützt. Wir bedanken uns.


Nochmalige Aufrufmöglichkeit des Applets in neuem Fenster

Applet in neuem Tab öffnen