Reworked the header of each org-mode file

disturbance-ty-sr/index.org

@@ -1,4 +1,4 @@
-#+TITLE:Vibrations induced by both the translation stage and the slip-ring
+#+TITLE:Vibrations induced by simultaneous scans of the translation stage and of the slip-ring
 :DRAWER:
 #+STARTUP: overview
 
@@ -30,19 +30,111 @@
 :END:
 
 * Measurement description
-All the stages are OFF except Ty and Slip-Ring.
+** Setup                                                             :ignore:
+*Setup*:
+All the stages are OFF except the translation stage and the Slip-Ring.
 
-Each of the signal is amplified by voltage amplifier:
-- 40db!!!
-- AC
-- 1kHz
+Each of the signal is amplified by voltage amplifiers with the following settings:
+- Gain: 40dB
+- AC/DC option: AC
+- Low pass filter: 1kHz
 
-col1: marble
-col2: sample
+The slip-ring is rotating at 60rpm. At the same time, scans with the translation stage are done at 1Hz with an amplitude of 600000cnt (= 3mm).
 
-Rotating the slip-ring at 60rpm and scans with the Ty at 1Hz.
-Everything on the Z direction
+Two geophones are used to measure the motion in the vertical direction of the marble and of the sample.
 
-- meas 50: Slip-Ring at 1Hz, Ty OFF
-- meas 51: Slip-Ring at 1Hz, Ty ON (The current and cnt error of Ty is also registered)
-- meas 52: Slip-Ring at 1Hz, Ty 1Hz 600000cnt
+** Goal                                                              :ignore:
+*Goal*:
+- The goal is to estimate the vibrations induced by the simultaneous scans of the spindle (here the slip-ring is used as the spindle is not fully functional yet) and of the translation stage
+
+** Measurements                                                      :ignore:
+*Measurements*:
+
+Three measurements are done:
+| Measurement File   | Description                                                                  |
+|--------------------+------------------------------------------------------------------------------|
+| =mat/data_050.mat= | Slip-Ring at 1Hz, Ty OFF                                                     |
+| =mat/data_051.mat= | Slip-Ring at 1Hz, Ty ON (The current and cnt error of Ty is also registered) |
+| =mat/data_052.mat= | Slip-Ring at 1Hz, Ty 1Hz 600000cnt                                           |
+
+Each of the measurement =mat= file contains one =data= array with 3 columns:
+| Column number | Description                     |
+|---------------+---------------------------------|
+|             1 | Geophone on the marble          |
+|             2 | Geophone at the sample location |
+|             3 | Time                            |
+
+* Data Analysis
+  :PROPERTIES:
+  :header-args:matlab+: :tangle matlab/disturbance_ty_sr.m
+  :header-args:matlab+: :comments org :mkdirp yes
+  :END:
+  <<sec:disturbance_ty_sr>>
+
+** ZIP file containing the data and matlab files                     :ignore:
+#+begin_src bash :exports none :results none
+  if [ matlab/disturbance_ty_sr.m -nt data/disturbance_ty_sr.zip ]; then
+    cp matlab/disturbance_ty_sr.m disturbance_ty_sr.m;
+    zip data/disturbance_ty_sr \
+        mat/data_050.mat \
+        mat/data_051.mat \
+        mat/data_052.mat \
+        mat/Ty-when-Rz-1Hz-and-Ty-1Hz.csv \
+        mat/Ty-when-Rz-1Hz.csv \
+        disturbance_ty_sr.m
+    rm disturbance_ty_sr.m;
+  fi
+#+end_src
+
+#+begin_note
+  All the files (data and Matlab scripts) are accessible [[file:data/disturbance_ty_sr.zip][here]].
+#+end_note
+
+** Matlab Init                                              :noexport:ignore:
+#+begin_src matlab :tangle no :exports none :results silent :noweb yes :var current_dir=(file-name-directory buffer-file-name)
+  <<matlab-dir>>
+#+end_src
+
+#+begin_src matlab :exports none :results silent :noweb yes
+  <<matlab-init>>
+#+end_src
+
+** Load data
+#+begin_src matlab
+  ty_off = load('mat/data_050.mat', 'data'); ty_off = ty_off.data;
+  ty_on  = load('mat/data_051.mat', 'data'); ty_on  = ty_on.data;
+  ty_1hz = load('mat/data_052.mat', 'data'); ty_1hz = ty_1hz.data;
+#+end_src
+
+** Time domain plots
+#+begin_src matlab
+  figure;
+  hold on;
+  plot(ty_1hz(:, 3), ty_1hz(:, 1), 'DisplayName', 'Marble - Z');
+  plot(ty_on(:, 3),  ty_on(:, 1),  'DisplayName', 'Marble - Z');
+  plot(ty_off(:, 3), ty_off(:, 1), 'DisplayName', 'Marble - Z');
+  hold off;
+  xlabel('Time [s]'); ylabel('Voltage [V]');
+  xlim([0, 100]); ylim([-5, 5]);
+  legend('Location', 'northeast');
+#+end_src
+
+#+begin_src matlab
+  xlim([0, 1]);
+#+end_src
+
+#+begin_src matlab
+  figure;
+  hold on;
+  plot(ty_1hz(:, 3), ty_1hz(:, 2), 'DisplayName', 'Marble - Z');
+  plot(ty_on(:, 3),  ty_on(:, 2),  'DisplayName', 'Marble - Z');
+  plot(ty_off(:, 3), ty_off(:, 2), 'DisplayName', 'Marble - Z');
+  hold off;
+  xlabel('Time [s]'); ylabel('Voltage [V]');
+  xlim([0, 100]); ylim([-5, 5]);
+  legend('Location', 'northeast');
+#+end_src
+
+#+begin_src matlab
+  xlim([0, 1]);
+#+end_src