#+TITLE:Vibrations induced by simultaneous scans of the translation stage and of the slip-ring
:DRAWER:
#+STARTUP: overview

#+LANGUAGE: en
#+EMAIL: dehaeze.thomas@gmail.com
#+AUTHOR: Dehaeze Thomas

#+HTML_LINK_HOME: ../index.html
#+HTML_LINK_UP: ../index.html

#+HTML_HEAD: <link rel="stylesheet" type="text/css" href="../css/htmlize.css"/>
#+HTML_HEAD: <link rel="stylesheet" type="text/css" href="../css/readtheorg.css"/>
#+HTML_HEAD: <link rel="stylesheet" type="text/css" href="../css/zenburn.css"/>
#+HTML_HEAD: <script type="text/javascript" src="../js/jquery.min.js"></script>
#+HTML_HEAD: <script type="text/javascript" src="../js/bootstrap.min.js"></script>
#+HTML_HEAD: <script type="text/javascript" src="../js/jquery.stickytableheaders.min.js"></script>
#+HTML_HEAD: <script type="text/javascript" src="../js/readtheorg.js"></script>

#+HTML_MATHJAX: align: center tagside: right font: TeX

#+PROPERTY: header-args:matlab  :session *MATLAB*
#+PROPERTY: header-args:matlab+ :comments org
#+PROPERTY: header-args:matlab+ :results none
#+PROPERTY: header-args:matlab+ :exports both
#+PROPERTY: header-args:matlab+ :eval no-export
#+PROPERTY: header-args:matlab+ :output-dir figs

#+PROPERTY: header-args:shell  :eval no-export
:END:

* Measurement description
** Setup                                                             :ignore:
*Setup*:
All the stages are OFF except the translation stage and the Slip-Ring.

Each of the signal is amplified by voltage amplifiers with the following settings:
- Gain: 40dB
- AC/DC option: AC
- Low pass filter: 1kHz

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).

Two geophones are used to measure the motion in the vertical direction of the marble and of the sample.

** 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