Merge remote-tracking branch 'origin/master'

37 files changed, 316 insertions, 58 deletions

diff --git a/doc/thesis/Fading.tex b/doc/thesis/Fading.tex
index f78c8db..3a9eaca 100644
--- a/doc/thesis/Fading.tex
+++ b/doc/thesis/Fading.tex
@@ -53,7 +53,7 @@
 	arrows,
 	patterns,
 }
-\tikzexternalize[mode=graphics if exists]{Fading}
+\tikzexternalize[prefix=figures/externalized/]{Fading}
 \pgfdeclarelayer{background}
 \pgfdeclarelayer{foreground}
 \pgfsetlayers{background,main,foreground}
diff --git a/doc/thesis/Makefile b/doc/thesis/Makefile
index 0afb1f6..8e0d1da 100644
--- a/doc/thesis/Makefile
+++ b/doc/thesis/Makefile
@@ -1,5 +1,5 @@
 TEX := xelatex
-TEXARGS := -halt-on-error -interaction=nonstopmode -shell-escape=1
+TEXARGS := --halt-on-error --interaction=nonstopmode --shell-escape=1 --enable-write18
 
 SOURCES := \
 	Fading.tex \
@@ -28,13 +28,15 @@ SOURCES := \
 	figures/tikz/packet-frame.tex \
 	figures/tikz/qpsk-simulations-static.tex \
 	figures/tikz/qpsk-simulations-dynamic.tex \
-	figures/tikz/qpsk-hardware.tex
+	figures/tikz/qam-simulations-dynamic.tex \
+	figures/tikz/hardware.tex
 
 DATA := \
 	figures/data/multipath_frequency_response.dat \
 	figures/data/rice_pdf.dat \
 	\
 	figures/data/qpsk_hw_lab.dat \
+	figures/data/qam_hw_lab.dat \
 	\
 	figures/data/qpsk_sim_nochannel.dat \
 	\
@@ -46,7 +48,12 @@ DATA := \
 	figures/data/qpsk_sim_dynamic_grblock_los_etsi_doppler5_k5.dat \
 	figures/data/qpsk_sim_dynamic_grblock_los_etsi_doppler70_k5.dat \
 	figures/data/qpsk_sim_dynamic_grblock_nlos_etsi_doppler5.dat \
-	figures/data/qpsk_sim_dynamic_grblock_nlos_etsi_doppler70.dat
+	figures/data/qpsk_sim_dynamic_grblock_nlos_etsi_doppler70.dat \
+	\
+	figures/data/qam_sim_dynamic_grblock_los_etsi_doppler5_k5.dat \
+	figures/data/qam_sim_dynamic_grblock_los_etsi_doppler70_k5.dat \
+	figures/data/qam_sim_dynamic_grblock_nlos_etsi_doppler5.dat \
+	figures/data/qam_sim_dynamic_grblock_nlos_etsi_doppler70.dat
 
 
 # Get the main file from the file
@@ -65,6 +72,7 @@ all: $(PDF)
 	/usr/bin/env python3 $<
 
 %.pdf: %.tex $(SOURCES) $(DATA)
+	mkdir -p figures/externalized
 	$(TEX) $(TEXARGS) $<
 	# bibliography
 	biber $(basename $(MAIN))
@@ -78,6 +86,7 @@ clean:
 	@rm -vf tex/*.aux
 	@rm -vf chapters/*.aux
 	@rm -vf figures/tikz/*.aux
+	@rm -vfrd figures/externalized
 
 cleanall: clean
 	@rm -vf $(PDF)
diff --git a/doc/thesis/chapters/implementation.tex b/doc/thesis/chapters/implementation.tex
index 3150b66..c22f6d3 100644
--- a/doc/thesis/chapters/implementation.tex
+++ b/doc/thesis/chapters/implementation.tex
@@ -293,24 +293,6 @@ Thus, they will be distributed among the other whole numbers. A window function
 	
 \end{lstlisting}
 
-%
-%\begin{figure}
-%	\centering
-%	\input{figures/tikz/qpsk-sim-constellations-static}
-%	\caption{
-%		Constellation diagrams for a simulated link using QPSK with AWGN and Rayleighan fading.
-%	}
-%\end{figure}
-%
-% \begin{figure}
-% 	\centering
-% 	\input{figures/tikz/qpsk-sim-constellations-dynamic}
-% 	\caption{
-% 		Constellation diagrams for a simulated link using QPSK with AWGN and Rayleighan fading. The paramters are: frequency offset of 0.2 \%, \SI{100}{\milli\volt} noise, dopper shift for \(v = \SI{2}{\meter\per\second}\), and a NLOS urban PDP.
-% 	}
-% 	\label{fig:dynamic-exp}
-% \end{figure}
-
 \subsection{Fading with statistical model}
 
 In order to represent the effect of the multipaht fading not only statically, a second model was created using the Frequency Selective Fading Model from Gnu Radio, according to \ref{sec:statistical-model},which was implemented after the algorithm from the paper \cite{Alimohammad2009}, with the help of the sum-of sinusoid principal (SOS). The algorithm in this block is implemented with the aim that only a small number of sinusoids are needed. 
@@ -370,16 +352,6 @@ The numbers of tags used in this case are similar to the number of given values.
 	\caption{Extended Typical Urban model (ETU) ETSI Standard PDP values for multipath fading propagation conditions \cite{ETSI}. \label{tab:etsi-tap-values}}
 \end{table}
 
-% \begin{figure}
-% 	\centering
-% 	\input{figures/tikz/qpsk-sim-constellations-dynamic-exp-NLOS-5}
-% 	\caption{
-% 		Constellation diagrams for a simulated link using QPSK and Rayleighan fading. With the ETU model and a Doppler frequency of \(\SI{5}{\hertz}\).
-% 	}
-% 	\label{fig:dynamic-exp-real}
-% \end{figure}
-
-
 \subsection{Measurements/Demonstration}
 
 \skelpar[5]{
@@ -447,15 +419,6 @@ As in \ref{sec:GUI} described the GUI was implemented, but unfortunately the par
 
 The second part which is missing is to be able to change the timing plot for the different scattering plots.
 
-%\begin{figure}
-%	\centering
-%	\label{fig:GUI}
-%	\input{figures/screenshots/gui_screenshot.png}
-%	\caption{
-%		Screenshot from the GUI
-%	}
-%\end{figure}
-
 % TODO : Piczure of the setup
 %TODO: Plots from the Hardware 
 \subsection{Incomplete parts}
@@ -484,18 +447,26 @@ Without those only the amplitudes could be seen in the Plots, with all the noise
 \newpage
 \begin{figure}
 	\centering
-	\label{fig:qpsk-simulations-dynamic}
 	\input{figures/tikz/qpsk-simulations-dynamic}
 	\caption{
 		Simulations with a dynamic fading channel model using PDP values of the Extended Typical Urban model (ETU) of the ETSI standard normative Annex B.2 in \cite{ETSI}. The color gradient represents progression in time.
+		\label{fig:qpsk-simulations-dynamic}
 	}
 \end{figure}
 \newpage
 \begin{figure}
 	\centering
-	\input{figures/tikz/qpsk-hardware}
+	\input{figures/tikz/qam-simulations-dynamic}
+	\caption{
+		TODO QAM simulation
+	}
+\end{figure}
+\begin{figure}
+	\centering
+	\input{figures/tikz/hardware}
 	\caption{
 		TODO QPSK hardware
 	}
 \end{figure}
+\newpage
 \restoregeometry
diff --git a/doc/thesis/chapters/theory.tex b/doc/thesis/chapters/theory.tex
index c73f94d..dfea46d 100644
--- a/doc/thesis/chapters/theory.tex
+++ b/doc/thesis/chapters/theory.tex
@@ -69,7 +69,7 @@ Before explaining how the two carrier signals are generated, we first need to di
 			= 0, \text{ and } \label{eqn:orthogonal-condition} \\
 		\langle \phi_k, \phi_k \rangle
 			&= \int_T \phi_k \phi_k^*  \,dt = 1,
-			\text{ where } k \text{ is either } i \text{ or } q. \label{eqn:orthonormal-condition}
+			\text{ where } k \text{ is either } i \text{ or } q \,\text{\cite{Gallager}}. \label{eqn:orthonormal-condition}
 	\end{align}
 \end{subequations}
 Provided these rather abstract conditions, let's define a new signal 
diff --git a/doc/thesis/figures/data/flowgraphs/channel_qam_hw_lab.dat b/doc/thesis/figures/data/flowgraphs/channel_qam_hw_lab.dat
new file mode 100644
index 0000000..4a04026
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/channel_qam_hw_lab.dat
diff --git a/doc/thesis/figures/data/flowgraphs/channel_qam_sim_dynamic_grblock_los_etsi_doppler5_k5.dat b/doc/thesis/figures/data/flowgraphs/channel_qam_sim_dynamic_grblock_los_etsi_doppler5_k5.dat
new file mode 100644
index 0000000..49070d5
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/channel_qam_sim_dynamic_grblock_los_etsi_doppler5_k5.dat
diff --git a/doc/thesis/figures/data/flowgraphs/channel_qam_sim_dynamic_grblock_los_etsi_doppler70_k5.dat b/doc/thesis/figures/data/flowgraphs/channel_qam_sim_dynamic_grblock_los_etsi_doppler70_k5.dat
new file mode 100644
index 0000000..13a6201
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/channel_qam_sim_dynamic_grblock_los_etsi_doppler70_k5.dat
diff --git a/doc/thesis/figures/data/flowgraphs/channel_qam_sim_dynamic_grblock_nlos_etsi_doppler5.dat b/doc/thesis/figures/data/flowgraphs/channel_qam_sim_dynamic_grblock_nlos_etsi_doppler5.dat
new file mode 100644
index 0000000..207213c
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/channel_qam_sim_dynamic_grblock_nlos_etsi_doppler5.dat
diff --git a/doc/thesis/figures/data/flowgraphs/channel_qam_sim_dynamic_grblock_nlos_etsi_doppler70.dat b/doc/thesis/figures/data/flowgraphs/channel_qam_sim_dynamic_grblock_nlos_etsi_doppler70.dat
new file mode 100644
index 0000000..3a01a87
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/channel_qam_sim_dynamic_grblock_nlos_etsi_doppler70.dat
diff --git a/doc/thesis/figures/data/flowgraphs/equalized_qam_hw_lab.dat b/doc/thesis/figures/data/flowgraphs/equalized_qam_hw_lab.dat
new file mode 100644
index 0000000..fa4d8e5
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/equalized_qam_hw_lab.dat
diff --git a/doc/thesis/figures/data/flowgraphs/equalized_qam_sim_dynamic_grblock_los_etsi_doppler5_k5.dat b/doc/thesis/figures/data/flowgraphs/equalized_qam_sim_dynamic_grblock_los_etsi_doppler5_k5.dat
new file mode 100644
index 0000000..c235697
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/equalized_qam_sim_dynamic_grblock_los_etsi_doppler5_k5.dat
diff --git a/doc/thesis/figures/data/flowgraphs/equalized_qam_sim_dynamic_grblock_los_etsi_doppler70_k5.dat b/doc/thesis/figures/data/flowgraphs/equalized_qam_sim_dynamic_grblock_los_etsi_doppler70_k5.dat
new file mode 100644
index 0000000..59d13f2
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/equalized_qam_sim_dynamic_grblock_los_etsi_doppler70_k5.dat
diff --git a/doc/thesis/figures/data/flowgraphs/equalized_qam_sim_dynamic_grblock_nlos_etsi_doppler5.dat b/doc/thesis/figures/data/flowgraphs/equalized_qam_sim_dynamic_grblock_nlos_etsi_doppler5.dat
new file mode 100644
index 0000000..db59982
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/equalized_qam_sim_dynamic_grblock_nlos_etsi_doppler5.dat
diff --git a/doc/thesis/figures/data/flowgraphs/equalized_qam_sim_dynamic_grblock_nlos_etsi_doppler70.dat b/doc/thesis/figures/data/flowgraphs/equalized_qam_sim_dynamic_grblock_nlos_etsi_doppler70.dat
new file mode 100644
index 0000000..2f0b72c
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/equalized_qam_sim_dynamic_grblock_nlos_etsi_doppler70.dat
diff --git a/doc/thesis/figures/data/flowgraphs/locked_qam_hw_lab.dat b/doc/thesis/figures/data/flowgraphs/locked_qam_hw_lab.dat
new file mode 100644
index 0000000..f3c9b92
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/locked_qam_hw_lab.dat
diff --git a/doc/thesis/figures/data/flowgraphs/locked_qam_sim_dynamic_grblock_los_etsi_doppler5_k5.dat b/doc/thesis/figures/data/flowgraphs/locked_qam_sim_dynamic_grblock_los_etsi_doppler5_k5.dat
new file mode 100644
index 0000000..ecd61c9
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/locked_qam_sim_dynamic_grblock_los_etsi_doppler5_k5.dat
diff --git a/doc/thesis/figures/data/flowgraphs/locked_qam_sim_dynamic_grblock_los_etsi_doppler70_k5.dat b/doc/thesis/figures/data/flowgraphs/locked_qam_sim_dynamic_grblock_los_etsi_doppler70_k5.dat
new file mode 100644
index 0000000..7dc562e
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/locked_qam_sim_dynamic_grblock_los_etsi_doppler70_k5.dat
diff --git a/doc/thesis/figures/data/flowgraphs/locked_qam_sim_dynamic_grblock_nlos_etsi_doppler5.dat b/doc/thesis/figures/data/flowgraphs/locked_qam_sim_dynamic_grblock_nlos_etsi_doppler5.dat
new file mode 100644
index 0000000..3a48c8a
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/locked_qam_sim_dynamic_grblock_nlos_etsi_doppler5.dat
diff --git a/doc/thesis/figures/data/flowgraphs/locked_qam_sim_dynamic_grblock_nlos_etsi_doppler70.dat b/doc/thesis/figures/data/flowgraphs/locked_qam_sim_dynamic_grblock_nlos_etsi_doppler70.dat
new file mode 100644
index 0000000..e6ead95
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/locked_qam_sim_dynamic_grblock_nlos_etsi_doppler70.dat
diff --git a/doc/thesis/figures/data/flowgraphs/synchronized_qam_hw_lab.dat b/doc/thesis/figures/data/flowgraphs/synchronized_qam_hw_lab.dat
new file mode 100644
index 0000000..cdbbe25
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/synchronized_qam_hw_lab.dat
diff --git a/doc/thesis/figures/data/flowgraphs/synchronized_qam_sim_dynamic_grblock_los_etsi_doppler5_k5.dat b/doc/thesis/figures/data/flowgraphs/synchronized_qam_sim_dynamic_grblock_los_etsi_doppler5_k5.dat
new file mode 100644
index 0000000..32ed9d0
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/synchronized_qam_sim_dynamic_grblock_los_etsi_doppler5_k5.dat
diff --git a/doc/thesis/figures/data/flowgraphs/synchronized_qam_sim_dynamic_grblock_los_etsi_doppler70_k5.dat b/doc/thesis/figures/data/flowgraphs/synchronized_qam_sim_dynamic_grblock_los_etsi_doppler70_k5.dat
new file mode 100644
index 0000000..f642f4f
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/synchronized_qam_sim_dynamic_grblock_los_etsi_doppler70_k5.dat
diff --git a/doc/thesis/figures/data/flowgraphs/synchronized_qam_sim_dynamic_grblock_nlos_etsi_doppler5.dat b/doc/thesis/figures/data/flowgraphs/synchronized_qam_sim_dynamic_grblock_nlos_etsi_doppler5.dat
new file mode 100644
index 0000000..a1418dd
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/synchronized_qam_sim_dynamic_grblock_nlos_etsi_doppler5.dat
diff --git a/doc/thesis/figures/data/flowgraphs/synchronized_qam_sim_dynamic_grblock_nlos_etsi_doppler70.dat b/doc/thesis/figures/data/flowgraphs/synchronized_qam_sim_dynamic_grblock_nlos_etsi_doppler70.dat
new file mode 100644
index 0000000..842c36a
--- /dev/null
+++ b/doc/thesis/figures/data/flowgraphs/synchronized_qam_sim_dynamic_grblock_nlos_etsi_doppler70.dat
diff --git a/doc/thesis/figures/data/qam_hw_lab.py b/doc/thesis/figures/data/qam_hw_lab.py
new file mode 100644
index 0000000..fde109a
--- /dev/null
+++ b/doc/thesis/figures/data/qam_hw_lab.py
@@ -0,0 +1,35 @@
+#!/usr/bin/env python3
+
+import utils
+import numpy as np
+
+# get array of samples [channel, synchronized, equalized, locked]
+samples = utils.load_samples(__file__)
+
+# range of samples we want to show
+start = 25e3
+end = start +800
+
+# select every second samples
+select_samples = lambda arr: arr[int(start):int(end):8]
+values = map(select_samples, samples)
+
+# split into imaginary and real parts
+get_parts = lambda v: (np.real(v), np.imag(v))
+parts = [p for v in map(get_parts, values) for p in v]
+
+# add 'samplenr' metadata
+parts += [np.arange(0, len(parts[0]))]
+
+# zip data and add header
+data = np.array(list(zip(*parts)))
+headers = [
+    "channel_re",      "channel_im",
+    "synchronized_re", "synchronized_im",
+    "equalized_re",    "equalized_im",
+    "locked_re",       "locked_im",
+    "samplenr",
+]
+
+# save to file
+utils.save_to_file(__file__, data, headers)
diff --git a/doc/thesis/figures/data/qam_sim_dynamic_grblock_los_etsi_doppler5_k5.py b/doc/thesis/figures/data/qam_sim_dynamic_grblock_los_etsi_doppler5_k5.py
new file mode 100644
index 0000000..a644402
--- /dev/null
+++ b/doc/thesis/figures/data/qam_sim_dynamic_grblock_los_etsi_doppler5_k5.py
@@ -0,0 +1,35 @@
+#!/usr/bin/env python3
+
+import utils
+import numpy as np
+
+# get array of samples [channel, synchronized, equalized, locked]
+samples = utils.load_samples(__file__)
+
+# range of samples we want to show
+start = 10e3
+end = start +40000
+
+# select every second samples
+select_samples = lambda arr: arr[int(start):int(end):200]
+values = map(select_samples, samples)
+
+# split into imaginary and real parts
+get_parts = lambda v: (np.real(v), np.imag(v))
+parts = [p for v in map(get_parts, values) for p in v]
+
+# add 'samplenr' metadata
+parts += [np.arange(0, len(parts[0]))]
+
+# zip data and add header
+data = np.concatenate([np.array(list(zip(*parts)))])
+headers = [
+    "channel_re",      "channel_im",
+    "synchronized_re", "synchronized_im",
+    "equalized_re",    "equalized_im",
+    "locked_re",       "locked_im", 
+    "samplenr",
+]
+
+# save to file
+utils.save_to_file(__file__, data, headers)
diff --git a/doc/thesis/figures/data/qam_sim_dynamic_grblock_los_etsi_doppler70_k5.py b/doc/thesis/figures/data/qam_sim_dynamic_grblock_los_etsi_doppler70_k5.py
new file mode 100644
index 0000000..8f35d50
--- /dev/null
+++ b/doc/thesis/figures/data/qam_sim_dynamic_grblock_los_etsi_doppler70_k5.py
@@ -0,0 +1,35 @@
+#!/usr/bin/env python3
+
+import utils
+import numpy as np
+
+# get array of samples [channel, synchronized, equalized, locked]
+samples = utils.load_samples(__file__)
+
+# range of samples we want to show
+start = 25e3
+end = start +1e3
+
+# select every second samples
+select_samples = lambda arr: arr[int(start):int(end):8]
+values = map(select_samples, samples)
+
+# split into imaginary and real parts
+get_parts = lambda v: (np.real(v), np.imag(v))
+parts = [p for v in map(get_parts, values) for p in v]
+
+# add 'samplenr' metadata
+parts += [np.arange(0, len(parts[0]))]
+
+# zip data and add header
+data = np.array(list(zip(*parts)))
+headers = [
+    "channel_re",      "channel_im",
+    "synchronized_re", "synchronized_im",
+    "equalized_re",    "equalized_im",
+    "locked_re",       "locked_im",
+    "samplenr"
+]
+
+# save to file
+utils.save_to_file(__file__, data, headers)
diff --git a/doc/thesis/figures/data/qam_sim_dynamic_grblock_nlos_etsi_doppler5.py b/doc/thesis/figures/data/qam_sim_dynamic_grblock_nlos_etsi_doppler5.py
new file mode 100644
index 0000000..e434809
--- /dev/null
+++ b/doc/thesis/figures/data/qam_sim_dynamic_grblock_nlos_etsi_doppler5.py
@@ -0,0 +1,35 @@
+#!/usr/bin/env python3
+
+import utils
+import numpy as np
+
+# get array of samples [channel, synchronized, equalized, locked]
+samples = utils.load_samples(__file__)
+
+# range of samples we want to show
+start = 20e3
+end = start +10000
+
+# select every second samples
+select_samples = lambda arr: arr[int(start):int(end):100]
+values = map(select_samples, samples)
+
+# split into imaginary and real parts
+get_parts = lambda v: (np.real(v), np.imag(v))
+parts = [p for v in map(get_parts, values) for p in v]
+
+# add 'samplenr' metadata
+parts += [np.arange(0, len(parts[0]))]
+
+# zip data and add header
+data = np.array(list(zip(*parts)))
+headers = [
+    "channel_re",      "channel_im",
+    "synchronized_re", "synchronized_im",
+    "equalized_re",    "equalized_im",
+    "locked_re",       "locked_im",
+    "samplenr"
+]
+
+# save to file
+utils.save_to_file(__file__, data, headers)
diff --git a/doc/thesis/figures/data/qam_sim_dynamic_grblock_nlos_etsi_doppler70.py b/doc/thesis/figures/data/qam_sim_dynamic_grblock_nlos_etsi_doppler70.py
new file mode 100644
index 0000000..5b93750
--- /dev/null
+++ b/doc/thesis/figures/data/qam_sim_dynamic_grblock_nlos_etsi_doppler70.py
@@ -0,0 +1,31 @@
+#!/usr/bin/env python3
+
+import utils
+import numpy as np
+
+# get array of samples [channel, synchronized, equalized, locked]
+samples = utils.load_samples(__file__)
+
+# range of samples we want to show
+start = 1e3
+end = start +800
+
+# select every second samples
+select_samples = lambda arr: arr[int(start):int(end):8]
+values = map(select_samples, samples)
+
+# split into imaginary and real parts
+get_parts = lambda v: (np.real(v), np.imag(v))
+parts = [p for v in map(get_parts, values) for p in v]
+
+# zip data and add header
+data = np.array(list(zip(*parts)))
+headers = [
+    "channel_re",      "channel_im",
+    "synchronized_re", "synchronized_im",
+    "equalized_re",    "equalized_im",
+    "locked_re",       "locked_im"
+]
+
+# save to file
+utils.save_to_file(__file__, data, headers)
diff --git a/doc/thesis/figures/data/qpsk_hw_lab.py b/doc/thesis/figures/data/qpsk_hw_lab.py
index 741f32f..f6f816b 100644
--- a/doc/thesis/figures/data/qpsk_hw_lab.py
+++ b/doc/thesis/figures/data/qpsk_hw_lab.py
@@ -18,13 +18,17 @@ values = map(select_samples, samples)
 get_parts = lambda v: (np.real(v), np.imag(v))
 parts = [p for v in map(get_parts, values) for p in v]
 
+# add 'samplenr' metadata
+parts += [np.arange(0, len(parts[0]))]
+
 # zip data and add header
 data = np.array(list(zip(*parts)))
 headers = [
     "channel_re",      "channel_im",
     "synchronized_re", "synchronized_im",
     "equalized_re",    "equalized_im",
-    "locked_re",       "locked_im"
+    "locked_re",       "locked_im",
+    "samplenr",
 ]
 
 # save to file
diff --git a/doc/thesis/figures/data/qpsk_sim_dynamic_grblock_los_etsi_doppler5_k5.py b/doc/thesis/figures/data/qpsk_sim_dynamic_grblock_los_etsi_doppler5_k5.py
index f621fd0..819800f 100644
--- a/doc/thesis/figures/data/qpsk_sim_dynamic_grblock_los_etsi_doppler5_k5.py
+++ b/doc/thesis/figures/data/qpsk_sim_dynamic_grblock_los_etsi_doppler5_k5.py
@@ -18,13 +18,17 @@ values = map(select_samples, samples)
 get_parts = lambda v: (np.real(v), np.imag(v))
 parts = [p for v in map(get_parts, values) for p in v]
 
+# add 'samplenr' metadata
+parts += [np.arange(0, len(parts[0]))]
+
 # zip data and add header
 data = np.array(list(zip(*parts)))
 headers = [
     "channel_re",      "channel_im",
     "synchronized_re", "synchronized_im",
     "equalized_re",    "equalized_im",
-    "locked_re",       "locked_im"
+    "locked_re",       "locked_im",
+    "samplenr"
 ]
 
 # save to file
diff --git a/doc/thesis/figures/data/qpsk_sim_dynamic_grblock_los_etsi_doppler70_k5.py b/doc/thesis/figures/data/qpsk_sim_dynamic_grblock_los_etsi_doppler70_k5.py
index c82fbd1..d2147df 100644
--- a/doc/thesis/figures/data/qpsk_sim_dynamic_grblock_los_etsi_doppler70_k5.py
+++ b/doc/thesis/figures/data/qpsk_sim_dynamic_grblock_los_etsi_doppler70_k5.py
@@ -18,13 +18,17 @@ values = map(select_samples, samples)
 get_parts = lambda v: (np.real(v), np.imag(v))
 parts = [p for v in map(get_parts, values) for p in v]
 
+# add 'samplenr' metadata
+parts += [np.arange(0, len(parts[0]))]
+
 # zip data and add header
 data = np.array(list(zip(*parts)))
 headers = [
     "channel_re",      "channel_im",
     "synchronized_re", "synchronized_im",
     "equalized_re",    "equalized_im",
-    "locked_re",       "locked_im"
+    "locked_re",       "locked_im",
+    "samplenr"
 ]
 
 # save to file
diff --git a/doc/thesis/figures/data/qpsk_sim_dynamic_grblock_nlos_etsi_doppler5.py b/doc/thesis/figures/data/qpsk_sim_dynamic_grblock_nlos_etsi_doppler5.py
index c82fbd1..15e8fd0 100644
--- a/doc/thesis/figures/data/qpsk_sim_dynamic_grblock_nlos_etsi_doppler5.py
+++ b/doc/thesis/figures/data/qpsk_sim_dynamic_grblock_nlos_etsi_doppler5.py
@@ -18,13 +18,17 @@ values = map(select_samples, samples)
 get_parts = lambda v: (np.real(v), np.imag(v))
 parts = [p for v in map(get_parts, values) for p in v]
 
+# add 'samplenr' metadata
+parts += [np.arange(0, len(parts[0]))]
+
 # zip data and add header
 data = np.array(list(zip(*parts)))
 headers = [
     "channel_re",      "channel_im",
     "synchronized_re", "synchronized_im",
     "equalized_re",    "equalized_im",
-    "locked_re",       "locked_im"
+    "locked_re",       "locked_im",
+    "samplenr",
 ]
 
 # save to file
diff --git a/doc/thesis/figures/data/qpsk_sim_dynamic_grblock_nlos_etsi_doppler70.py b/doc/thesis/figures/data/qpsk_sim_dynamic_grblock_nlos_etsi_doppler70.py
index eae2afc..5194407 100644
--- a/doc/thesis/figures/data/qpsk_sim_dynamic_grblock_nlos_etsi_doppler70.py
+++ b/doc/thesis/figures/data/qpsk_sim_dynamic_grblock_nlos_etsi_doppler70.py
@@ -18,13 +18,17 @@ values = map(select_samples, samples)
 get_parts = lambda v: (np.real(v), np.imag(v))
 parts = [p for v in map(get_parts, values) for p in v]
 
+# add 'samplenr' metadata
+parts += [np.arange(0, len(parts[0]))]
+
 # zip data and add header
 data = np.array(list(zip(*parts)))
 headers = [
     "channel_re",      "channel_im",
     "synchronized_re", "synchronized_im",
     "equalized_re",    "equalized_im",
-    "locked_re",       "locked_im"
+    "locked_re",       "locked_im",
+    "samplenr"
 ]
 
 # save to file
diff --git a/doc/thesis/figures/tikz/qpsk-hardware.tex b/doc/thesis/figures/tikz/hardware.tex
index 68a31a6..d13926b 100644
--- a/doc/thesis/figures/tikz/qpsk-hardware.tex
+++ b/doc/thesis/figures/tikz/hardware.tex
@@ -13,7 +13,8 @@
 				#5
 			]
 
-			\addplot[only marks, scatter] table [x index = #2, y index = #3] {#1};
+			\addplot[only marks, scatter, point meta = \thisrow{samplenr}]
+				table [x index = #2, y index = #3] {#1};
 		\end{axis}
 	\end{tikzpicture}
 	\hfill
@@ -28,15 +29,15 @@
 		text width = 3em,
 		align = right,
 	},
-	title = {LOS, Indoors (Lab)},
+	title = {QPSK LOS, Indoors (Lab)},
 }
 %
 \makeplot{figures/data/qpsk_hw_lab.dat}{0}{1}{.25}{
-	title = {LOS, Indoors (Lab)},
+	title = {TODO},
 }
 %
-\makeplot{figures/data/qpsk_hw_lab.dat}{0}{1}{.25}{
-	title = {LOS, Indoors (Lab)},
+\makeplot{figures/data/qam_hw_lab.dat}{0}{1}{2}{
+	title = {QAM LOS, Indoors (Lab)},
 }
 \newline
 
@@ -51,7 +52,7 @@
 	},
 }
 \makeplot{figures/data/qpsk_hw_lab.dat}{2}{3}{.25}{}
-\makeplot{figures/data/qpsk_hw_lab.dat}{2}{3}{.25}{}
+\makeplot{figures/data/qam_hw_lab.dat}{2}{3}{2}{}
 \newline
 
 % -------------------------------------------------------------------
@@ -65,7 +66,7 @@
 	},
 }
 \makeplot{figures/data/qpsk_hw_lab.dat}{4}{5}{2}{}
-\makeplot{figures/data/qpsk_hw_lab.dat}{4}{5}{2}{}
+\makeplot{figures/data/qam_hw_lab.dat}{4}{5}{2}{}
 \newline
 
 % -------------------------------------------------------------------
@@ -79,5 +80,5 @@
 	},
 }
 \makeplot{figures/data/qpsk_hw_lab.dat}{6}{7}{2}{}
-\makeplot{figures/data/qpsk_hw_lab.dat}{6}{7}{2}{}
+\makeplot{figures/data/qam_hw_lab.dat}{6}{7}{2}{}
 \newline
diff --git a/doc/thesis/figures/tikz/qam-simulations-dynamic.tex b/doc/thesis/figures/tikz/qam-simulations-dynamic.tex
new file mode 100644
index 0000000..88bea4a
--- /dev/null
+++ b/doc/thesis/figures/tikz/qam-simulations-dynamic.tex
@@ -0,0 +1,85 @@
+% vim: set ts=2 sw=2 noet:
+
+\newcommand{\makeplot}[5]{%
+	\hfill
+	\tikzset{external/force remake}
+	\begin{tikzpicture}
+		\begin{axis}[
+				width  = {\linewidth / 3.3},
+				height = {\linewidth / 3.3},
+				grid = major,
+				xmin = {-#4}, xmax = {#4},
+				ymin = {-#4}, ymax = {#4},
+				colormap name=viridis,
+				#5
+			]
+
+			\addplot[only marks, scatter, point meta = \thisrow{samplenr}]
+				table [x index = #2, y index = #3] {#1};
+		\end{axis}
+	\end{tikzpicture}
+	\hfill
+}
+
+% -------------------------------------------------------------------
+
+\noindent
+\makeplot{figures/data/qam_sim_dynamic_grblock_nlos_etsi_doppler5.dat}{0}{1}{4}{
+	ylabel = {Channel with ISI},
+	yticklabel style = {
+		text width = 3em,
+		align = right,
+	},
+	title = {NLOS, \(\Delta f = \SI{5}{\hertz}\)},
+}
+%
+\makeplot{figures/data/qam_sim_dynamic_grblock_los_etsi_doppler5_k5.dat}{0}{1}{4}{
+	title = {LOS, \(\Delta f = \SI{5}{\hertz}, K = 5\)},
+}
+%
+\makeplot{figures/data/qam_sim_dynamic_grblock_los_etsi_doppler70_k5.dat}{0}{1}{4}{
+	title = {LOS, \(\Delta f = \SI{70}{\hertz}, K = 5\)},
+}
+\newline
+
+% -------------------------------------------------------------------
+
+\noindent
+\makeplot{figures/data/qam_sim_dynamic_grblock_nlos_etsi_doppler5.dat}{2}{3}{4}{%
+	ylabel = {Synchronized},
+	yticklabel style = {
+		text width = 3em,
+		align = right,
+	},
+}
+\makeplot{figures/data/qam_sim_dynamic_grblock_los_etsi_doppler5_k5.dat}{2}{3}{4}{}
+\makeplot{figures/data/qam_sim_dynamic_grblock_los_etsi_doppler70_k5.dat}{2}{3}{4}{}
+\newline
+
+% -------------------------------------------------------------------
+
+\noindent
+\makeplot{figures/data/qam_sim_dynamic_grblock_nlos_etsi_doppler5.dat}{4}{5}{2}{%
+	ylabel = {Equalized},
+	yticklabel style = {
+		text width = 3em,
+		align = right,
+	},
+}
+\makeplot{figures/data/qam_sim_dynamic_grblock_los_etsi_doppler5_k5.dat}{4}{5}{2}{}
+\makeplot{figures/data/qam_sim_dynamic_grblock_los_etsi_doppler70_k5.dat}{4}{5}{2}{}
+\newline
+
+% -------------------------------------------------------------------
+
+\noindent
+\makeplot{figures/data/qam_sim_dynamic_grblock_nlos_etsi_doppler5.dat}{6}{7}{2}{%
+	ylabel = {Locked},
+	yticklabel style = {
+		text width = 3em,
+		align = right,
+	},
+}
+\makeplot{figures/data/qam_sim_dynamic_grblock_los_etsi_doppler5_k5.dat}{6}{7}{2}{}
+\makeplot{figures/data/qam_sim_dynamic_grblock_los_etsi_doppler70_k5.dat}{6}{7}{2}{}
+\newline
diff --git a/doc/thesis/figures/tikz/qpsk-simulations-dynamic.tex b/doc/thesis/figures/tikz/qpsk-simulations-dynamic.tex
index 4d01c92..0764c2a 100644
--- a/doc/thesis/figures/tikz/qpsk-simulations-dynamic.tex
+++ b/doc/thesis/figures/tikz/qpsk-simulations-dynamic.tex
@@ -13,7 +13,8 @@
 				#5
 			]
 
-			\addplot[only marks, scatter] table [x index = #2, y index = #3] {#1};
+			\addplot[only marks, scatter, point meta = \thisrow{samplenr}]
+				table [x index = #2, y index = #3] {#1};
 		\end{axis}
 	\end{tikzpicture}
 	\hfill