From eb2c2220bb7a6f17a99cd0cdd083bc94bee9520c Mon Sep 17 00:00:00 2001 From: wieerwill Date: Thu, 24 Feb 2022 10:58:25 +0100 Subject: [PATCH] Abtastung und Rekonstruktion --- Grundlagen der Biosignalverarbeitung.pdf | 4 +- Grundlagen der Biosignalverarbeitung.tex | 90 ++++++++++++------------ 2 files changed, 46 insertions(+), 48 deletions(-) diff --git a/Grundlagen der Biosignalverarbeitung.pdf b/Grundlagen der Biosignalverarbeitung.pdf index 02db5cb..0f0fa85 100644 --- a/Grundlagen der Biosignalverarbeitung.pdf +++ b/Grundlagen der Biosignalverarbeitung.pdf @@ -1,3 +1,3 @@ version https://git-lfs.github.com/spec/v1 -oid sha256:b319ebc60762163b65d2543f9a1d1cba949ce695f5f1f2972836b3cae7ffeb42 -size 3836653 +oid sha256:e3396c37476e7ad8eea1a8da387748f35026d588935b7be2d21990962a8858e3 +size 3812553 diff --git a/Grundlagen der Biosignalverarbeitung.tex b/Grundlagen der Biosignalverarbeitung.tex index 2985333..8493577 100644 --- a/Grundlagen der Biosignalverarbeitung.tex +++ b/Grundlagen der Biosignalverarbeitung.tex @@ -926,8 +926,7 @@ \item Spektrogramm des gefilterten Signals zeigt typischen Effekt nichtkonstanter Laufzeit: Die Frequenzen(-gruppen), die höher als ca. 20Hz liegen, werden deutlich verzögert und erzeugen das schnelle Nachschwingen. \end{itemize*} - \section{Signalkonditionierung, Abtastung und Digitalisierung}\label{signalkonditionierung-abtastung-und-digitalisierung} - + \section{Signalkonditionierung} \subsection{Pegelanpassung}\label{pegelanpassung} Pegelanpassung notwendig für \begin{itemize*} @@ -955,9 +954,9 @@ \begin{itemize*} \item Abtastung: Erfassung des momentanen Wertes eines Signals zu definierten Zeitpunkten \item Üblicherweise ist Signal zeitkontinuierlich, nach der Abtastung liegt eine zeitdiskrete Variante des Signals vor - \item kontinuierliches Signal (links), zeitdiskretes Signal (rechts) + \item kontinuierliches Signal (links) vs. zeitdiskretes Signal (rechts) \end{itemize*} - \begin{center} + \begin{center} \includegraphics[width=.3\linewidth]{Assets/Biosignalverarbeitung-kontinuierliches-signal.png} \includegraphics[width=.3\linewidth]{Assets/Biosignalverarbeitung-zeitdiskretes-signal.png} \end{center} @@ -970,79 +969,78 @@ \item EKG ist digitalisiert um es darstellen zu können \item stark instationäres Signal mit scharfer R-Zacke. Daraus ergibt sich ein periodisches Spektrum, da die R-Zacke einen nahezu impulsartigen Verlauf hat. \item Um Überlappung (Aliasing) der Spektren und Mehrdeutigkeiten zu vermeiden, muss gelten: $\frac{1}{T_A}\geq 2*f_{max}$ - \item Nyquist-Frequenz entspricht der halben Grundfrequenz der Abtastrate, sie begrenzt nach oben das bei Null beginnende sog. Basisband. + \item Nyquist-Frequenz entspricht der halben Grundfrequenz der Abtastrate, sie begrenzt nach oben das bei Null beginnende sog. Basisband. \item Basisband ist der Frequenzbereich, in dem man bei der Signalanalyse arbeitet \item damit gespiegelte Spektrum nicht bis ins Basisband reicht und dadurch das Vorhandensein real nichtexistenter Signalkomponenten vortäuscht, muss gewährleistet werden, dass die halbe AR höher liegt, als die höchste Frequenz des Signals, d.h. die AR muss mindestens doppelt so hoch sein, die die höchste vorhandene Frequenz \item periodische Wiederholung des Spektrums nach der Abtastung hat folgende praktische Bedeutung: Da sich das Spektrum mit jeder Harmonischen der AR wiederholt und gespiegelt wird, kann man ein bandbegrenztes Signal ins Basisband holen \end{itemize*} %\includegraphics[width=.5\linewidth]{Assets/Biosignalverarbeitung-rekonstruierter-sinus.png} - Oft wird nach der Abtastung und anschließenden Signalverarbeitung das Ergebnis im ursprünglichen Bereich des Signals benötigt. Dazu ist die sog. Rekonstruktion notwendig, also eine Übertragung aus dem Analyse-in den Originalbereich. Zur Rekonstruktion ist allgemein eine Interpolation zwischen den diskreten Punkten notwendig, oder aus Sicht der Filtertheorie die Anwendung eines Interpolationsfilters. Das einfachste Interpolationsfilter ist ein Tiefpass in Nähe der höchsten Signalfrequenz. Wurde bei der Abtastung das Abtasttheorem verletzt, so treten im rekonstruierten Signal Komponenten auf, die im Originalsignal nicht vorhanden waren, siehe oben. - - Dieser Effekt kann auch im Fernsehen beobachtet werden: Tragen die Sprecher z.B. ein Hemd mit einem sehr feinen Strichmuster, so reicht die Bildschirmauflösung -d.h. die räumliche Abtastrate -nicht aus, um das Muster richtig zu erfassen und am Fernsehmonitor entsteht sog. Moiree, d.h. großflächige Farbmuster, die es in der Realität nicht gibt. + \subsubsection{Rekonstruktion} + \begin{itemize*} + \item nach der Abtastung/Signalverarbeitung Ergebnis oft im ursprünglichen Bereich des Signals benötigt + \item Übertragung aus Analyse- in Originalbereich + \item Interpolation zwischen diskreten Punkten notwendig oder aus Sicht der Filtertheorie die Anwendung eines Interpolationsfilters + \item einfachste Interpolationsfilter ist Tiefpass in Nähe der höchsten Signalfrequenz. Wurde bei Abtastung das Abtasttheorem verletzt, treten im rekonstruierten Signal Komponenten auf, die im Originalsignal nicht vorhanden waren + \end{itemize*} + % Dieser Effekt kann auch im Fernsehen beobachtet werden: Tragen die Sprecher z.B. ein Hemd mit einem sehr feinen Strichmuster, so reicht die Bildschirmauflösung -d.h. die räumliche Abtastrate -nicht aus, um das Muster richtig zu erfassen und am Fernsehmonitor entsteht sog. Moiree, d.h. großflächige Farbmuster, die es in der Realität nicht gibt. %\includegraphics[width=.5\linewidth]{Assets/Biosignalverarbeitung-abtastung-rekonstruktion.png} \begin{itemize*} - \item $y(t)=x(t)*\sum^{\infty}_{n=-\infty} \sigma(t-nT_A) \Leftrightarrow Y(\omega)=X(\omega)*\sum^{\infty}_{n=-\infty}\sigma(\omega-n\frac{1}{T_A})$ + \item $y(t)=x(t)*\sum^{\infty}_{n=-\infty} \sigma(t-nT_A) \Leftrightarrow$ $Y(\omega)=X(\omega)*\sum^{\infty}_{n=-\infty}\sigma(\omega-n\frac{1}{T_A})$ \item Übergang aus dem kontinuierlichen Zeitbereich in eine Folge, d.h. Entkopplung von der Abtastperiode \item $y(n)=y(nT_A) \Leftrightarrow Y(k)=Y(k\omega_A/M)=\sum_{n=1}^M y(nT_A)^{-jkn/M}$ \item $Y(K)\Leftrightarrow FFT(y(n))\Rightarrow$ normierte Frequenz $\omega\in(0,2\pi)\vee f\in(0,1)$ \item Nyquist Frequenz $\omega_N=\pi$, $f_N=0,5$ - \item Nach der Abtastung und Digitalisierung hat das Signal die Form einer Zahlenfolge bzw. eines Vektors oder Matrix. Ist die Abtastrate unbekannt, so ist das Signal auch nicht mehr reproduzierbar. Da sich aber Analysen und digitale Filterung grundsätzlich auch ohne Kenntnis der Abtastrate durchführen lassen, wird die sog. normierte Frequenz eingeführt, die bei der Rekonstruktion durch eine reale Abtastrate ersetzt wird. + \item Nach Abtastung und Digitalisierung hat Signal Form einer Zahlenfolge/Vektors oder Matrix. Ist Abtastrate unbekannt, so ist Signal nicht mehr reproduzierbar + \item Da Analysen/digitale Filterung grundsätzlich ohne Kenntnis der Abtastrate durchführbar, wird sog. normierte Frequenz eingeführt, die bei Rekonstruktion durch reale Abtastrate ersetzt \end{itemize*} - \includegraphics[width=.5\linewidth]{Assets/Biosignalverarbeitung-amplitudenmodulation.png} + % \includegraphics[width=.5\linewidth]{Assets/Biosignalverarbeitung-amplitudenmodulation.png} + %\begin{itemize*} + % \item $EKG\_\{AM\}=EKG*sin(\omega\_c t)$ + % \item Beispiel zum Abtasttheorem: Das EKG wird für eine Kabelübertragung mit einem Träger bei 10kHz multipliziert, was nachrichtentechnisch einer Amplitudenmodulation entspricht. Das Spektrum spiegelt sich um den Träger herum, ähnlich wie bei der Abtastung. Hier gibt es allerdings nur eine Spiegelung, da der Träger eine Harmonische ist und somit im Spektrum nur eine Nadel darstellt. Es entstehen zwei Seitenbänder, das obere und das untere. Beide sind hinsichtlich des Informationsgehaltes völlig identisch. + % \item Die Frage ist nun zu beantworten, wie hoch die Abtastrate für ein solches Signal sein muss. + %\end{itemize*} + Abtasttheorem Kotelnikov, Channon ($T_A=1/2f_{max}$) \begin{itemize*} - \item $EKG\_\{AM\}=EKG*sin(\omega\_c t)$ - \item Beispiel zum Abtasttheorem: Das EKG wird für eine Kabelübertragung mit einem Träger bei 10kHz multipliziert, was nachrichtentechnisch einer Amplitudenmodulation entspricht. Das Spektrum spiegelt sich um den Träger herum, ähnlich wie bei der Abtastung. Hier gibt es allerdings nur eine Spiegelung, da der Träger eine Harmonische ist und somit im Spektrum nur eine Nadel darstellt. Es entstehen zwei Seitenbänder, das obere und das untere. Beide sind hinsichtlich des Informationsgehaltes völlig identisch. - \item Die Frage ist nun zu beantworten, wie hoch die Abtastrate für ein solches Signal sein muss. - \end{itemize*} - - Abtasttheorem Kotelnikov, Channon ($T\_A=1/2f\_\{max\}$) - - \begin{itemize*} - \item hinreichende/notwendige Bedingung \item hinreichend aber nicht notwendig $AR\geq 22ksps$ \item $Y(\omega)=X(\omega)*\sum^{\infty}_{n=-\infty} \sigma(\omega-n\frac{1}{T_A})$ \item notwendig und hinreichend: $AR\geq 2ksps$ %\item \includegraphics[width=.5\linewidth]{Assets/Biosignalverarbeitung-abtastung-kotelnikov.png} - \item Da die höchste Frequenz im Signal 11kHz beträgt, müsste die Abtastrate mind. 22ksps betragen. Dies ist eine hinreichende, aber keine notwendige Bedingung. Denn bezieht man sich auf die Wiederholung des Spektrums um jede Harmonische der Abtastrate, wird eine AR von 2ksps ausreichen. Damit passt eine Wiederholung des Spektrums in das Basisband, in wir ja analysieren. + \item höchste Frequenz im Signal 11kHz $\rightarrow$ Abtastrate mind. 22ksps. hinreichende, aber keine notwendige Bedingung + \item bezieht man sich auf Wiederholung des Spektrums um jede Harmonische der Abtastrate, wird eine AR von 2ksps ausreichen. Damit passt eine Wiederholung des Spektrums in das Basisband \end{itemize*} Pulsamplitudenmodulation (PAM) - \begin{itemize*} - \item nach Sample \& Hold - \begin{itemize*} - \item Zeit diskret - \item Pegel analog - \end{itemize*} + \item nach Sample \& Hold: Zeit diskret, Pegel analog %\item \includegraphics[width=.5\linewidth]{Assets/Biosignalverarbeitung-Pulsamplitudenmodulation.png} - \item Die Abtastung (sample \& hold) entspricht nachrichtentechnisch der PAM: die Werte treten in definierten Abständen entsprechend dem Pulsbreite (Abtastperiode) auf und haben einen kontinuierlichen Wertebereich. Allerdings spielt diese Modulationsart in der Nachrichtentechnik keine praktische Rolle, wichtig ist sie für die Theorie. + \item Abtastung entspricht nachrichtentechnisch der PAM: Werte treten in definierten Abständen entsprechend dem Pulsbreite (Abtastperiode) auf und haben kontinuierlichen Wertebereich + \item spielt in der Nachrichtentechnik keine praktische Rolle, wichtig für Theorie \end{itemize*} Mehrkanalsystem - Simultansampling - \begin{itemize*} - \item \includegraphics[width=.5\linewidth]{Assets/Biosignalverarbeitung-Mehrkanalsysteme.png} - \item Oft werden in der Technik und vor allem in der Medizin mehrkanalige Messsysteme benötigt. Für die Analyse ist von entscheidender Bedeutung, dass der zeitliche Zusammenhang der Kanäle identisch, oder zumindest bekannt ist. Beim echten Simultansampling werden alle Kanalsignale zum selben Zeitpunkt abgetastet und sequentiell digitalisiert. Damit reicht im Normalfall ein ADC für alle Kanäle. - %\item \includegraphics[width=.5\linewidth]{Assets/Biosignalverarbeitung-Mehrkanalsysteme2.png} - \begin{itemize*} - \item alle Kanäle im selben Augenblick abgetastet, echtes Simultansampling - \end{itemize*} - \item \includegraphics[width=.5\linewidth]{Assets/Biosignalverarbeitung-Mehrkanalsysteme3.png} - \begin{itemize*} - \item Um den HW-Aufwand zu minimieren, werden die Kanalsignale sequentiell abgetastet und digitalisiert, dabei kann die Einsparung an Elektronik beachtlich sein. Signalanalytisch kann sie jedoch problematisch sein: Aus der Signalsequenz wird das Simultansignal über die Laufzeitkorrektur in der FFT zurückgerechnet. Bei zeitkritischen Vorgängen ist diese Alternative zu verwerfen, da die durch die sequentielle Abtastung verlorengegangenen Signalteile durch Rückrechnung nicht mehr zu retten sind. - \end{itemize*} - %\item \includegraphics[width=.5\linewidth]{Assets/Biosignalverarbeitung-Mehrkanalsysteme4.png} - \begin{itemize*} - \item Versatz der Kanäle um $T_A/N$ - \item Rechnerische Korrektur der Abtastzeit (nicht-online-fähig) - \item $X^{}*(j\omega)=X(j\omega)^{j\omega T_A/N}$ - \end{itemize*} + \item Oft mehrkanalige Messsysteme benötigt + \item Für Analyse entscheidend, dass zeitlicher Zusammenhang der Kanäle identisch oder bekannt + \item Bei echtem Simultansampling werden alle Kanalsignale zum selben Zeitpunkt abgetastet und sequentiell digitalisiert. Im Normalfall reicht ein ADC für alle Kanäle + \item HW-Aufwand minimieren $\rightarrow$ Kanalsignale sequentiell abgetastet und digitalisiert + \item Signalanalytisch problematisch: Aus Signalsequenz wird Simultansignal über Laufzeitkorrektur in der FFT zurückgerechnet. Bei zeitkritischen Vorgängen ist dies zu verwerfen, da die durch die sequentielle Abtastung verlorengegangenen Signalteile durch Rückrechnung nicht mehr zu retten sind + \item parallele (links) vs sequentielle (rechts) Abtastung + \begin{center} + \includegraphics[width=.3\linewidth]{Assets/Biosignalverarbeitung-Mehrkanalsysteme.png} + \includegraphics[width=.3\linewidth]{Assets/Biosignalverarbeitung-Mehrkanalsysteme3.png} + \end{center} + \item simultan (links) vs versetzte (rechts) Kanäle\\ + \begin{center} + \includegraphics[width=.2\linewidth]{Assets/Biosignalverarbeitung-Mehrkanalsysteme2.png} + \includegraphics[width=.2\linewidth]{Assets/Biosignalverarbeitung-Mehrkanalsysteme4.png} + \end{center} + \item Versatz der Kanäle um $T_A/N$ + \item Rechnerische Korrektur der Abtastzeit (nicht-online-fähig) $X^*(j\omega)=X(j\omega)^{j\omega T_A/N}$ \end{itemize*} - \subsection{Digitalisierung}\label{digitalisierung} \subsubsection{Prinzipien der AD Wandlung}\label{prinzipien-der-ad-wandlung}