aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authorNao Pross <np@0hm.ch>2022-08-31 20:13:46 +0200
committerNao Pross <np@0hm.ch>2022-08-31 20:13:46 +0200
commit1a8e52ded09496a745756882253c1d3fd1e76996 (patch)
tree0f1c7f772bb57e318382d4bd621b5f7181ef0e3f
parentMerge remote-tracking branch 'upstream/master' (diff)
downloadSeminarSpezielleFunktionen-1a8e52ded09496a745756882253c1d3fd1e76996.tar.gz
SeminarSpezielleFunktionen-1a8e52ded09496a745756882253c1d3fd1e76996.zip
kugel: Feedback and minor changes, add reference
-rw-r--r--buch/papers/kugel/applications.tex40
-rw-r--r--buch/papers/kugel/references.bib12
-rw-r--r--buch/papers/kugel/spherical-harmonics.tex6
3 files changed, 35 insertions, 23 deletions
diff --git a/buch/papers/kugel/applications.tex b/buch/papers/kugel/applications.tex
index b527ebd..f8f3edd 100644
--- a/buch/papers/kugel/applications.tex
+++ b/buch/papers/kugel/applications.tex
@@ -18,7 +18,8 @@ the most interesting applications we came across during our research.
\subfigure[Gauss' Law \label{kugel:fig:eeg-flux}]%
{\includegraphics[width=.4\linewidth]{papers/kugel/figures/flux}}
\caption{
- Electroencephalography.
+ Courtesy of C. Hope \cite{sheerman-chase_volunteer_2012} for picture (a),
+ and Wikimedia \cite{maschen_english_2013} for (b).
\label{kugel:fig:eeg}
}
\end{figure}
@@ -34,6 +35,7 @@ relate to the spherical harmonics, we will first quickly recap a bit of physics,
electrodynamics to be precise.
\subsubsection{Electrodynamics}
+\nocite{griffiths_introduction_2015}
In section \ref{kugel:sec:construction:eigenvalue} we have shown that the
spherical harmonics arise from the surface spherical Laplacian operator, whose
@@ -46,23 +48,23 @@ electric potential $\phi(x, y, z)$:
\nabla^2 \phi
= \nabla \cdot \nabla \phi
= \nabla \cdot \mathbf{E}
- = \rho / \varepsilon,
+ = \frac{1}{\varepsilon} \rho,
\quad \text{or} \quad
\iiint_\Omega \nabla \cdot \mathbf{E} \, dv
= \iint_{\partial \Omega} \mathbf{E} \cdot d\mathbf{s}
- = \Phi / \varepsilon.
+ = \frac{1}{\varepsilon} \Phi.
\end{equation*}
Put into words: on the left we have the differential form, where we recall that
the Laplacian (which is a second derivative) is the divergence of the gradient.
Unpacking the notation we first see that we have the gradient of the potential,
which is just the electric field $\mathbf{E}$, and then the divergence of said
electric field is proportional to the charge density $\rho$. So, the Laplacian
-of the electric potential is the charge density! For those that are more
-familiar with the integral form of Maxwell's equation, we have also included an
-additional step using the divergence theorem, which brings us to the electric
-Flux, which by Gauss' law (shown in the iconic\footnote{Every electrical
-engineer has seen this picture so many times that is probably burnt in their
-eyes.} figure \ref{kugel:fig:eeg-flux}) equals the net electric charge.
+of the electric potential is proportional to the charge density! For those that
+are more familiar with the integral form of Maxwell's equation, we have also
+included an additional step using the divergence theorem, which brings us to the
+electric Flux $\Phi$, which by Gauss' law (shown in the iconic\footnote{Every
+electrical engineer has seen this picture so many times that is probably burnt
+in their eyes.} figure \ref{kugel:fig:eeg-flux}) equals the net electric charge.
Now, an important observation is that if we switch to spherical coordinates, the
physics does not change. So, the spherical Laplacian $\sphlaplacian$ of the
@@ -96,7 +98,7 @@ finite linear combination of spherical harmonics:
= \sum_{n=1}^N \sum_{m=-n}^n a_{m,n} Y^m_n(\vartheta, \varphi),
\end{equation*}
where the values $a_{m,n}$ are the unknowns of our interpolation problem. Now to
-the measurements: we let $\phi_1, \phi_2, \ldots, p_M$ be the measured voltages
+the measurements: we let $\phi_1, \phi_2, \ldots, \phi_M$ be the measured voltages
at points in space $p_1, p_2, \ldots, p_M$ (position of the electrodes). To
simplify, we will assume that the electrodes are reasonably evenly distributed,
which means that we have no points that are on top of each other or at wildly
@@ -132,10 +134,10 @@ that (hint: eigenvalues)
\end{equation*}
So that when substituted into \eqref{kugel:eqn:eeg-min} results in
\begin{align*}
- \int_{\partial S} \left|
+ \int_{\partial S} \biggl|
\sum_{n=1}^N \sum_{m=-n}^n n(n+1) a_{m,n}
Y^m_n(\vartheta, \varphi)
- \right|^2 ds
+ \biggr|^2 ds
= \sum_{m, m'} \sum_{n, n'} a_{m',n'} \overline{a_{m,n}}
n'(n'+1) n(n+1)
\underbrace{\int_{\partial S} Y^{m'}_{n'} \overline{Y^m_n} \, ds}_{
@@ -165,15 +167,15 @@ discuss a few interesting implications and problems.
The most interesting perhaps unforeseen fact is that with this method we are
getting a free (!) spectral analysis, since the coefficients $a_{m,n}$ are the
-spectrum of the interpolated electric field $V(\vartheta, \varphi)$. However,
-like in the non spherical Fourier transformation, we only get a \emph{finite}
-resolution since our measurement are spatially discrete. In fact, if we know the
-mean angular inter-electrode distance $\gamma$ we can actually formulate a
-Nyquist frequency just like in the usual Fourier theory:
+spectrum of the interpolated electric potential $V(\vartheta, \varphi)$.
+However, like in the non spherical Fourier transformation, we only get a
+\emph{finite} resolution since our measurement are spatially discrete. In fact,
+if we know the mean angular inter-electrode distance $\gamma$ we can actually
+formulate a Nyquist frequency just like in the usual Fourier theory:
\begin{equation}
f_N = \frac{\pi}{2T}
\iff
- n_N = \left\lfloor \frac{\pi}{2\gamma} \right\rfloor.
+ n_N = \biggl\lfloor \frac{\pi}{2\gamma} \biggr\rfloor.
\end{equation}
Before concluding this overview of EEG, we should point out that in practice
@@ -195,7 +197,7 @@ we will not discuss it here, since this is getting too long already. Another
important issue is that in the real world, we cannot ``evenly distribute'' the
electrodes on our head. As shown in the image, most of the electrodes are on a
cap, and then there are just a few on the face, and almost none near the jawline
-and chin. This not something that can be ignored, and in fact, makes the
+and chin. This is not something that can be ignored, and in fact, makes the
analysis much more difficult. Finally, the most obvious problem is that human
heads are not perfect spheres. Here too, it is possible to account for this fact
and model the head with a more complex shape at the cost of making the math
diff --git a/buch/papers/kugel/references.bib b/buch/papers/kugel/references.bib
index 07e4d5d..1405d7c 100644
--- a/buch/papers/kugel/references.bib
+++ b/buch/papers/kugel/references.bib
@@ -342,4 +342,14 @@ This photo was taken during studies that resulted in the publication: Hope, C, S
note = {Publisher: Taylor \& Francis
\_eprint: https://doi.org/10.3109/00207458808986175},
keywords = {cross-validated current source density (csd) estimation, cross-validated functional brain mapping, cross-validated laplacian estimation, {ROLANDO} {BISCAY}-{LIRIO}, source derivations, spatial analysis, spherical harmonic fourier expansion},
-} \ No newline at end of file
+}
+
+@artwork{maschen_english_2013,
+ title = {English: Divergence theorem in {EM}},
+ url = {https://commons.wikimedia.org/wiki/File:Divergence_theorem_in_EM.svg},
+ shorttitle = {English},
+ author = {{Maschen}},
+ urldate = {2022-08-31},
+ date = {2013-05-12},
+ file = {Wikimedia Snapshot:/Users/npross/Zotero/storage/MCML5F28/FileDivergence_theorem_in_EM.html:text/html},
+}
diff --git a/buch/papers/kugel/spherical-harmonics.tex b/buch/papers/kugel/spherical-harmonics.tex
index b540531..fb5a144 100644
--- a/buch/papers/kugel/spherical-harmonics.tex
+++ b/buch/papers/kugel/spherical-harmonics.tex
@@ -871,8 +871,8 @@ computational cost lower by a factor of six \cite{davari_new_2013}.
The goal of this subsection's part is to apply the recurrence relations of the
$P^m_n(z)$ functions to the Spherical Harmonics. With some little adjustments
we will be able to have recursion equations for them too. As previously written
-the most of the work is already done. Now it is only a matter of minor
-mathematical operations/rearrangements. We can start by listing all of them:
+most of the work is already done. Now it is only a matter of minor mathematical
+operations/rearrangements. We can start by listing all of them:
\begin{subequations}
\begin{align}
Y^m_n(\vartheta, \varphi) &= \dfrac{1}{(2n+1)\cos \vartheta} \left[
@@ -899,7 +899,7 @@ mathematical operations/rearrangements. We can start by listing all of them:
\begin{proof}[Proof of \eqref{kugel:eqn:rec-sph-harm-1}]
We can multiply both sides of equality in \eqref{kugel:eqn:rec-leg-1} by $e^{im
\varphi}$ and perform the substitution $z=\cos \vartheta$. After a few simple
- algebraic steps, we will obtain the relation we are looking for
+ algebraic steps, we will obtain the relation we are looking for.
\end{proof}
\begin{proof}[Proof of \eqref{kugel:eqn:rec-sph-harm-2}]