# Messtechnik HS2021 - Tutorial 4

## Aufgabe 1: Diskrete Fourier Transformation und Datenverarbeitung
---------
Ein Kohlenstoff NMR Signal ($^{13}{\rm C}$) wurde aufgenommen und mit einer Abtastrate von 15 kHz digitalisiert. Das komplexwertige FID wurde durch ein MATLAB Program digitalisiert, exportiert und ist in `fid.mat` zu finden.

---------

### 1a)
Generieren Sie aus dem komplexwertigen FID das dazugehörige Spektrum und plotten Sie dessen Realteil mit der korrekten Frequenzachse in kHz.

*Hinweis:* Um MATLAB Dateien in Python zu importieren, kann man die ``scipy.io.loadmat`` Funktion brauchen:

````
from scipy.io import loadmat 
data = loadmat('fid.mat')
fid = data['fid']
fid = np.squeeze(fid) 
````


```python

```

### 1b)
Durch einen Matched Filter lässt sich das Signal-zu-Rausch-Verhältnis (SNR) nachträglich verbessern. Das Signal kann beschrieben werden als
$$ y(t) = \exp(-i\omega_0 t) \exp(-\frac{t}{T_2}) $$
wobei $\omega_0$ die Larmorfrequenz ist und $T_2$ die transversale Relaxationszeit. Bei einem Matched Filter wird das Zeitsignal $ y(t) $ mit der Zerfallsfunktion $ y_e(t) = \exp(-\frac{t}{T_2}) $ multipliziert.
Zeigen Sie die Auswirkungen auf das Spektrum bei der Anwendung von Matched Filtern mit Zerfallszeiten $T_2 = [1,2,5,10,50,100]$.


```python

```

### 1c)
Finden Sie eine passende Zerfallszeit für den Matched Filter, sodass Sie eine gute Balance zwischen Auflösung und SNR erhalten.


```python

```

### 1d)
Nehmen Sie nun an die Aufnahmezeit des FID's ist beschränkt auf 20 ms. Erstellen Sie aus `fid.mat` das entsprechende abgeschnittene Signal, hängen Sie dann Nullen an (Zerofilling), sodass Sie am Ende ein Signal mit insgesamt 2048 Punkte erhalten. Was fällt Ihnen auf, wenn Sie die Spektren des abgeschnittenen und des originalen Signal mit dem vollständigen vergleichen? 


```python

```

### 1e)
Programmieren Sie ein Hamming Apodisierungsfenster 

$$ w_{\text{hamming}}(t) = 0.54 + 0.46 \cos(\frac{\pi t}{t_{\text{max}}})$$

für das bei 60 ms abgeschnittene FID. Wenden Sie das Hamming Window auf das abgeschnittene Signal an, zerofillen Sie auf 2048 Punkte und vergleichen Sie die Spektren mit und ohne Hamming Window.


```python

```

### 1f)
Zeigen Sie, dass ein Hamming Fenster die Auflösung im Frequenzbereich im Vergleich zu einem "hart" abgeschnittenen FID verschlechtert, jedoch die *ripples* besser unterdrückt. Verwenden Sie Zerofilling au 2048 Punkte. Vergleichen Sie dazu die entsprechende Faltungsfunktion im Frequenzbereich.


```python

```

## Aufgabe 2: Digitalisierung eines Zeitsignals
---------
Ein zeitlich abfallendes oszillierendes Signal hat die Form:
$\begin{equation}
	x(t) = 
	\begin{cases}
		 \exp \left( - \frac{t}{T_2} \right) \left\{ \cos(2\pi \nu_1 t) + 3\cos(2\pi \nu_2 t) \right\} & t\geq 0 \\
	 	0 & t < 0
	 \end{cases}
\end{equation}$
mit den beiden Frequenzkomponenten $ \nu_1 =10\,{\rm Hz}$ und $\nu_2 = 150\,{\rm Hz}$ sowie der Zerfallszeit $T_2  = 0.5\,{\rm s}$. Das Signal wird total durch 2048 Punkte aufgenommen und digitalisiert.

---------

### 2a)
Welche Abtastrate schlagen Sie für das Digitalisieren dieses Signals vor? Begründen Sie Ihre Antwort mit Hilfe des Abtasttheorems.
Berechnen Sie anhand der vorgeschlagenen Abtastrate das Zeitsignal und das dazugehörige Spektrum.



```python

```

### 2b)
Vergleichen Sie unterschiedlichen Abtastraten $f_s = 200,320,1000\,{\rm Hz}$ für das obengenannte Signal . Was für Erkenntnisse können Sie aus den entsprechenden Spektren gewinnen?


```python

```

a9fbf2c77ddaf418af642210a614aa5c32ce6cb1
ipynb
Jupyter Notebook
tutorial_6_signal_processing/tutorial_6_blank.ipynb
JeschkeLab/messtechnikpy
0ba4ecd8858edad77c2a1dae1daa0af25df75c58
MIT
tutorial_6_signal_processing/tutorial_6_blank.ipynb
JeschkeLab/messtechnikpy
0ba4ecd8858edad77c2a1dae1daa0af25df75c58
MIT
tutorial_6_signal_processing/tutorial_6_blank.ipynb
JeschkeLab/messtechnikpy
0ba4ecd8858edad77c2a1dae1daa0af25df75c58
MIT
Qwen/Qwen-72B
1. YES
2. YES

