Redo some figures for the huddle test

huddle-test-geophones/index.org

@@ -106,8 +106,8 @@ We load the data of the z axis of two geophones.
 ** Computation of the ASD of the measured voltage
 We first define the parameters for the frequency domain analysis.
 #+begin_src matlab :results none
-  win = hanning(ceil(length(x1)/100));
   Fs = 1/dt;
+  win = hanning(ceil(10*Fs));
 #+end_src
 
 #+begin_src matlab :results none
@@ -171,12 +171,12 @@ The ASD of the measured velocity is shown on figure [[fig:psd_velocity]].
   set(gca, 'xscale', 'log');
   set(gca, 'yscale', 'log');
   xlabel('Frequency [Hz]'); ylabel('PSD [m/s/sqrt(Hz)]')
-  xlim([2, 500]);
+  xlim([0.1, 500]);
 #+end_src
 
 #+NAME: fig:psd_velocity
 #+HEADER: :tangle no :exports results :results value raw replace :noweb yes
-#+begin_src matlab :var filepath="figs/psd_velocity.pdf" :var figsize="wide-tall" :post pdf2svg(file=*this*, ext="png")
+#+begin_src matlab :var filepath="figs/psd_velocity.pdf" :var figsize="full-tall" :post pdf2svg(file=*this*, ext="png")
   <<plt-matlab>>
 #+end_src
 
@@ -195,12 +195,12 @@ We also plot the ASD in displacement (figure [[fig:asd_displacement]]);
   hold off;
   set(gca, 'xscale', 'log'); set(gca, 'yscale', 'log');
   xlabel('Frequency [Hz]'); ylabel('PSD [m/s/sqrt(Hz)]')
-  xlim([2, 500]);
+  xlim([0.1, 500]);
 #+end_src
 
 #+NAME: fig:asd_displacement
 #+HEADER: :tangle no :exports results :results value raw replace :noweb yes
-#+begin_src matlab :var filepath="figs/asd_displacement.pdf" :var figsize="wide-tall" :post pdf2svg(file=*this*, ext="png")
+#+begin_src matlab :var filepath="figs/asd_displacement.pdf" :var figsize="full-tall" :post pdf2svg(file=*this*, ext="png")
   <<plt-matlab>>
 #+end_src
 
@@ -235,7 +235,7 @@ We also compute the coherence between the two signals (figure [[fig:coh_geophone
   xlabel('Frequency [Hz]'); ylabel('Phase');
 
   linkaxes([ax1,ax2],'x');
-  xlim([1, 500]);
+  xlim([0.1, 500]);
 #+end_src
 
 #+NAME: fig:tf_geophones
@@ -258,7 +258,7 @@ We also compute the coherence between the two signals (figure [[fig:coh_geophone
   plot(f, coh12);
   set(gca, 'xscale', 'log');
   xlabel('Frequency [Hz]'); ylabel('Coherence');
-  ylim([0,1]); xlim([1, 500]);
+  ylim([0,1]); xlim([0.1, 500]);
 #+end_src
 
 #+NAME: fig:coh_geophones
@@ -339,12 +339,12 @@ The instrumental noise is computed below. The result in V^2/Hz is shown on figur
   hold off;
   set(gca, 'xscale', 'log'); set(gca, 'yscale', 'log');
   xlabel('Frequency [Hz]'); ylabel('PSD [$V^2/Hz$]');
-  xlim([1, 500]);
+  xlim([0.1, 500]);
 #+end_src
 
 #+NAME: fig:intrumental_noise_V
 #+HEADER: :tangle no :exports results :results value raw replace :noweb yes
-#+begin_src matlab :var filepath="figs/intrumental_noise_V.pdf" :var figsize="wide-tall" :post pdf2svg(file=*this*, ext="png")
+#+begin_src matlab :var filepath="figs/intrumental_noise_V.pdf" :var figsize="full-tall" :post pdf2svg(file=*this*, ext="png")
   <<plt-matlab>>
 #+end_src
 
@@ -363,12 +363,12 @@ This is then further converted into velocity and compared with the ground veloci
   hold off;
   set(gca, 'xscale', 'log'); set(gca, 'yscale', 'log');
   xlabel('Frequency [Hz]'); ylabel('PSD [$m/s/\sqrt{Hz}$]');
-  xlim([1, 500]);
+  xlim([0.1, 500]);
 #+end_src
 
 #+NAME: fig:intrumental_noise_velocity
 #+HEADER: :tangle no :exports results :results value raw replace :noweb yes
-#+begin_src matlab :var filepath="figs/intrumental_noise_velocity.pdf" :var figsize="wide-tall" :post pdf2svg(file=*this*, ext="png")
+#+begin_src matlab :var filepath="figs/intrumental_noise_velocity.pdf" :var figsize="full-tall" :post pdf2svg(file=*this*, ext="png")
   <<plt-matlab>>
 #+end_src