From 7368a30d66e8ac155bb4b8f8de205dd809e0d921 Mon Sep 17 00:00:00 2001
From: jaseg <git-bigdata-wsl-arch@jaseg.de>
Date: Wed, 13 May 2020 13:22:17 +0200
Subject: modulation, freq estimation: add some blurb, citations

---
 ma/safety_reset.tex | 71 ++++++++++++++++++++++++++++++++++++-----------------
 1 file changed, 48 insertions(+), 23 deletions(-)

(limited to 'ma/safety_reset.tex')

diff --git a/ma/safety_reset.tex b/ma/safety_reset.tex
index 3462740..303b519 100644
--- a/ma/safety_reset.tex
+++ b/ma/safety_reset.tex
@@ -1011,30 +1011,45 @@ controllable load:
     \item[Modulation amplitude] proportionally related to modulation power. In a practical setup we might realize a
         modulation power up to a few hundred \si{\mega\watt} which would yield maybe a few tens of \si{\milli\hertz} of
         frequency amplitude.
-    \item[Modulation pre-emphasis and slew-rate control]. Pre-emphasis might be necessary to ensure an adequate SNR at
-        the receiver. Slew-rate control and other shaping measures might be necessary to reduce the impact of these
-        sudden load changes on the transmitter's primary function (say, aluminium smelting) and to prevent disturbances
-        to grid components.
+    \item[Modulation pre-emphasis and slew-rate control]. Pre-emphasis might be necessary to ensure an adequate
+        Signal-to-Noise ratio (SNR) at the receiver. Slew-rate control and other shaping measures might be necessary to
+        reduce the impact of these sudden load changes on the transmitter's primary function (say, aluminium smelting)
+        and to prevent disturbances to grid components.
     \item[Modulation frequency]. For a practical implementation a careful study would be necessary to determine an
         optimal frequency band for operation. On one hand we need to prevent disturbances to the grid such as through
-        excitation of some local or inter-area modes. On the other hand we need to optimize SNR and data rate to achieve
-        optimal latency between transmission start and successful reception and to reduce the overall burden on
-        transmitter and grid.
+        excitation of some local or inter-area modes. On the other hand we need to optimize Signal-to-Noise ratio (SNR)
+        and data rate to achieve optimal latency between transmission start and successful reception and to reduce the
+        overall burden on transmitter and grid.
     \item[Further modulation parameters]. The modulation itself has numerous parameters that are discussed in sec.\
         \ref{mod_params} below.
 \end{description}
 
-\section{From grid frequency to a reliable communications channel}
-% FIXME
+\section{From grid frequency to a reliable communication channel}
 
 \subsection{Channel properties}
-% FIXME
+In this section we will explore how we can construct a reliable communication channel from the analog primitive we
+outline in the previous section. Our load control approach to grid frequency modulation leads to a channel with the
+following properties.
 
+\begin{description}
+    \item[Slow-changing.] Accurate grid frequency measurements need several periods of the mains sine wave. Faster
+        sampling rates can be achieved with more complex specialized synchrophasor estimation algorithms but this will
+        result in a tradeoff between sampling rate and accuracy\cite{belega01}.
+    \item[Analog.] Grid frequency is an analog signal.
+    \item[Noisy.] While stable over long periods of time thanks to Load-Frequency Control\cite{entsoe04} it shows
+        considerable random short-term variations. In addition our modulation amplitude is limited by technical and
+        economic constraints so we have to find a system that will work at poor SNRs.
+    \item[Polarized.] Grid frequency measurements have an inherent sense of \emph{up} (higher frequencies). We can use
+        this in a polarized modulation scheme to encode information without first transmitting some reference signal to
+        establish this polarization.
+\end{description}
 
 \subsection{Modulation and its parameters}
-
 \label{mod_params}
 
+In this section we will consider how to select a good set of parameters for a modulation scheme fitting grid frequency
+modulation.
+
 The sensitivity of the grid to oscillation at particular frequencies described above means we should avoid any
 modulation technique that would concentrate a lot of energy in a small bandwidth. Taking this principle to its extreme
 provides us with a useful pointer towards techniques that might work well: Spread-spectrum techniques. By employing
@@ -1285,12 +1300,12 @@ required precision for manageable averaging times--we would need either a ADC sa
 for a reconstruction through interpolated readings an impractically high ADC resolution.
 
 Detail on the inner workings of commercial phasor measurement units is scarce but given their essential role to SCADA
-systems there is a large amount of academic research on such algorithms\cite{narduzzi01,derviskadic01}. A popular
-approach to these systems is to perform a Short-Time Fourier Transform (STFT) on ADC data sampled at high sampling rate
-(e.g. \SI{10}{\kilo\hertz}) and then perform some analysis on the frequency-domain data to precisely locate the strong peak
-around \SI{50}{\hertz}. A key observation here is that FFT bin size is going to be much larger than required frequency
-resolution. This fundamental limitiation follows from the nyquist criterion %FIXME maybe cite?  and if we had to process
-an \emph{arbitrary} signal this would highly limit our practical measurement accuracy
+systems there is a large amount of academic research on such algorithms\cite{narduzzi01,derviskadic01,belega01}. A
+popular approach to these systems is to perform a Short-Time Fourier Transform (STFT) on ADC data sampled at high
+sampling rate (e.g. \SI{10}{\kilo\hertz}) and then perform some analysis on the frequency-domain data to precisely
+locate the strong peak around \SI{50}{\hertz}. A key observation here is that FFT bin size is going to be much larger
+than required frequency resolution. This fundamental limitiation follows from the nyquist criterion %FIXME maybe cite?
+and if we had to process an \emph{arbitrary} signal this would highly limit our practical measurement accuracy
 \footnote{
     Some software packages providing FFT or STFT primitives such as scipy\cite{virtanen01} allow the user to
     super-sample FFT output by specifying an FFT width larger than input data length, padding the input data with zeros
@@ -1351,12 +1366,22 @@ resolution and despite numerous distortions.
 Published grid frequency estimation algorithms such as \textcite{narduzzi01} or \textcite{derviskadic01} are rather
 sophisticated and use a combination of techniques to reduce numerical errors in FFT calculation and peak fitting. Given
 that we do not need reference standard-grade accuracy for our application we chose to start with a very basic algorithm
-instead. We chose to use a general approach developed by experimental physicists at CERN that is described by
-\textcite{gasior01}. This approach assumes a general sinusoidal signal superimposed with harmonics and broadband noise.
-Applicable to a wide spectrum of practical signal analysis tasks it is a reasonable first-degree approximation of the
-much more sophisticated estimation algorithms developed specifically for power systems.  Some algorithms have components
-such as kalman filters\cite{narduzzi01} that require a phyiscal model.  As a general algorithm from \textcite{gasior01}
-does not require this kind of application-specific tuning, eliminating one source of error.
+instead. We chose to use a general approach to estimate the precise fundamental frequency of an arbitrary signal that
+was developed by experimental physicists at CERN and that is described by \textcite{gasior01}. This approach assumes a
+general sinusoidal signal superimposed with harmonics and broadband noise.  Applicable to a wide spectrum of practical
+signal analysis tasks it is a reasonable first-degree approximation of the much more sophisticated estimation algorithms
+developed specifically for power systems.  Some algorithms have components such as kalman filters\cite{narduzzi01} that
+require a phyiscal model.  As a general algorithm from \textcite{gasior01} does not require this kind of
+application-specific tuning, eliminating one source of error.
+
+The \textcite{gasior01} algorithm passes the windowed input signal through a DFT, then interpolates the signal's
+fundamental frequency by fitting a wavelet such as a gaussian to the largest peak in the DFT results. The bias parameter
+of this curve fit is an accurate estimation of the signal's fundamental frequency. This algorithm is similar to the
+simpler interpolated DFT algorithm used as a reference in much of the synchrophasor estimation
+literature\cite{borkowski01}. The three-term variant of the maximum sidelobe decay window often used there is a blackman
+window with parameter $\alpha = \frac{1}{4}$. Analysis has shown\cite{belega01} that the interpolated DFT algorithm is
+worse than algorithms involving more complex models under some conditions but that there is \emph{no free lunch} meaning
+that more complex perform worse when the input signal deviates from their models.
 
 \subsection{Frequency sensor hardware design}
 \label{sec-fsensor}
-- 
cgit 

