Add new test mat files
Before Width: | Height: | Size: 95 KiB After Width: | Height: | Size: 120 KiB |
Before Width: | Height: | Size: 22 KiB After Width: | Height: | Size: 23 KiB |
Before Width: | Height: | Size: 20 KiB After Width: | Height: | Size: 21 KiB |
Before Width: | Height: | Size: 32 KiB After Width: | Height: | Size: 33 KiB |
Before Width: | Height: | Size: 115 KiB After Width: | Height: | Size: 112 KiB |
Before Width: | Height: | Size: 28 KiB After Width: | Height: | Size: 29 KiB |
25
index.org
@ -46,6 +46,7 @@
|
|||||||
|
|
||||||
#+name: fig:test-bench-picture
|
#+name: fig:test-bench-picture
|
||||||
#+caption: Picture of the test bench. The Attocube and mirror are covered by a "bubble sheet"
|
#+caption: Picture of the test bench. The Attocube and mirror are covered by a "bubble sheet"
|
||||||
|
#+ATTR_ORG: :width 600px
|
||||||
[[file:figs/IMG-7865.JPG]]
|
[[file:figs/IMG-7865.JPG]]
|
||||||
|
|
||||||
** Matlab Init :noexport:ignore:
|
** Matlab Init :noexport:ignore:
|
||||||
@ -61,13 +62,14 @@
|
|||||||
The first measurement was made during ~17 hours with a sampling time of $T_s = 0.1\,s$.
|
The first measurement was made during ~17 hours with a sampling time of $T_s = 0.1\,s$.
|
||||||
|
|
||||||
#+begin_src matlab
|
#+begin_src matlab
|
||||||
load('./mat/long_test2.mat', 'x', 't')
|
load('./mat/long_test_plastic.mat', 'x', 't')
|
||||||
Ts = 0.1; % [s]
|
Ts = 0.1; % [s]
|
||||||
#+end_src
|
#+end_src
|
||||||
|
|
||||||
#+begin_src matlab :exports none
|
#+begin_src matlab :exports none
|
||||||
figure;
|
figure;
|
||||||
plot(t/60/60, 1e9*x)
|
plot(t/60/60, 1e9*x)
|
||||||
|
xlim([0, 17.5]);
|
||||||
xlabel('Time [h]'); ylabel('Displacement [nm]');
|
xlabel('Time [h]'); ylabel('Displacement [nm]');
|
||||||
#+end_src
|
#+end_src
|
||||||
|
|
||||||
@ -85,8 +87,8 @@ Let's fit the data with a step response to a first order low pass filter (Figure
|
|||||||
#+begin_src matlab
|
#+begin_src matlab
|
||||||
f = @(b,x) b(1)*(1 - exp(-x/b(2)));
|
f = @(b,x) b(1)*(1 - exp(-x/b(2)));
|
||||||
|
|
||||||
y_cur = x(t < 17*60*60);
|
y_cur = x(t < 17.5*60*60);
|
||||||
t_cur = t(t < 17*60*60);
|
t_cur = t(t < 17.5*60*60);
|
||||||
|
|
||||||
nrmrsd = @(b) norm(y_cur - f(b,t_cur)); % Residual Norm Cost Function
|
nrmrsd = @(b) norm(y_cur - f(b,t_cur)); % Residual Norm Cost Function
|
||||||
B0 = [400e-9, 2*60*60]; % Choose Appropriate Initial Estimates
|
B0 = [400e-9, 2*60*60]; % Choose Appropriate Initial Estimates
|
||||||
@ -99,7 +101,7 @@ The corresponding time constant is (in [h]):
|
|||||||
#+end_src
|
#+end_src
|
||||||
|
|
||||||
#+RESULTS:
|
#+RESULTS:
|
||||||
: 2.0576
|
: 2.0658
|
||||||
|
|
||||||
#+begin_src matlab :exports none
|
#+begin_src matlab :exports none
|
||||||
figure;
|
figure;
|
||||||
@ -107,6 +109,7 @@ The corresponding time constant is (in [h]):
|
|||||||
plot(t_cur/60/60, 1e9*y_cur);
|
plot(t_cur/60/60, 1e9*y_cur);
|
||||||
plot(t_cur/60/60, 1e9*f(B, t_cur));
|
plot(t_cur/60/60, 1e9*f(B, t_cur));
|
||||||
hold off;
|
hold off;
|
||||||
|
xlim([0, 17.5])
|
||||||
xlabel('Time [h]'); ylabel('Displacement [nm]');
|
xlabel('Time [h]'); ylabel('Displacement [nm]');
|
||||||
#+end_src
|
#+end_src
|
||||||
|
|
||||||
@ -125,7 +128,7 @@ We only select the data between =t1= and =t2=.
|
|||||||
The obtained displacement is shown in Figure [[fig:long_meas_time_domain_zoom]].
|
The obtained displacement is shown in Figure [[fig:long_meas_time_domain_zoom]].
|
||||||
|
|
||||||
#+begin_src matlab
|
#+begin_src matlab
|
||||||
t1 = 11; t2 = 17; % [h]
|
t1 = 10.5; t2 = 17.5; % [h]
|
||||||
|
|
||||||
x = x(t > t1*60*60 & t < t2*60*60);
|
x = x(t > t1*60*60 & t < t2*60*60);
|
||||||
x = x - mean(x);
|
x = x - mean(x);
|
||||||
@ -198,7 +201,7 @@ The time domain measurement is shown in Figure [[fig:short_meas_time_domain]].
|
|||||||
|
|
||||||
The Power Spectral Density of the measured displacement is computed
|
The Power Spectral Density of the measured displacement is computed
|
||||||
#+begin_src matlab
|
#+begin_src matlab
|
||||||
win = hann(ceil(length(x)/10));
|
win = hann(ceil(length(x)/20));
|
||||||
[p_2, f_2] = pwelch(x, win, [], [], 1/Ts);
|
[p_2, f_2] = pwelch(x, win, [], [], 1/Ts);
|
||||||
#+end_src
|
#+end_src
|
||||||
|
|
||||||
@ -230,6 +233,7 @@ The computed ASD of the two measurements are combined in Figure [[fig:psd_combin
|
|||||||
|
|
||||||
#+name: fig:picture-test-bench-aluminium-tube
|
#+name: fig:picture-test-bench-aluminium-tube
|
||||||
#+caption: Aluminium tube used to protect the beam path from disturbances
|
#+caption: Aluminium tube used to protect the beam path from disturbances
|
||||||
|
#+ATTR_ORG: :width 600px
|
||||||
[[file:figs/IMG-7864.JPG]]
|
[[file:figs/IMG-7864.JPG]]
|
||||||
|
|
||||||
** Matlab Init :noexport:ignore:
|
** Matlab Init :noexport:ignore:
|
||||||
@ -243,7 +247,7 @@ The computed ASD of the two measurements are combined in Figure [[fig:psd_combin
|
|||||||
|
|
||||||
** Aluminium Tube and Bubble Sheet
|
** Aluminium Tube and Bubble Sheet
|
||||||
#+begin_src matlab
|
#+begin_src matlab
|
||||||
load('./mat/long_test_plastic.mat');
|
load('./mat/short_test_plastic.mat');
|
||||||
Ts = 1e-4; % [s]
|
Ts = 1e-4; % [s]
|
||||||
#+end_src
|
#+end_src
|
||||||
|
|
||||||
@ -264,7 +268,7 @@ The computed ASD of the two measurements are combined in Figure [[fig:psd_combin
|
|||||||
|
|
||||||
** Only Aluminium Tube
|
** Only Aluminium Tube
|
||||||
#+begin_src matlab
|
#+begin_src matlab
|
||||||
load('./mat/long_test_alu_tube.mat');
|
load('./mat/short_test_alu_tube.mat');
|
||||||
Ts = 1e-4; % [s]
|
Ts = 1e-4; % [s]
|
||||||
#+end_src
|
#+end_src
|
||||||
|
|
||||||
@ -286,7 +290,7 @@ The time domain measurement is shown in Figure [[fig:short_meas_time_domain]].
|
|||||||
|
|
||||||
** Nothing
|
** Nothing
|
||||||
#+begin_src matlab
|
#+begin_src matlab
|
||||||
load('./mat/long_test_without_material.mat');
|
load('./mat/short_test_without_material.mat');
|
||||||
Ts = 1e-4; % [s]
|
Ts = 1e-4; % [s]
|
||||||
#+end_src
|
#+end_src
|
||||||
|
|
||||||
@ -319,11 +323,12 @@ The time domain measurement is shown in Figure [[fig:short_meas_time_domain]].
|
|||||||
hold off;
|
hold off;
|
||||||
set(gca, 'Xscale', 'log'); set(gca, 'Yscale', 'log');
|
set(gca, 'Xscale', 'log'); set(gca, 'Yscale', 'log');
|
||||||
ylabel('ASD [$m/\sqrt{Hz}$]'); xlabel('Frequency [Hz]');
|
ylabel('ASD [$m/\sqrt{Hz}$]'); xlabel('Frequency [Hz]');
|
||||||
|
xlim([1e-1, 5e3]);
|
||||||
legend('location', 'northeast');
|
legend('location', 'northeast');
|
||||||
#+end_src
|
#+end_src
|
||||||
|
|
||||||
#+begin_src matlab :tangle no :exports results :results file replace
|
#+begin_src matlab :tangle no :exports results :results file replace
|
||||||
exportFig('figs/asd_noise_comp_bubble_aluminium.pdf', 'width', 'wide', 'height', 'normal');
|
exportFig('figs/asd_noise_comp_bubble_aluminium.pdf', 'width', 'wide', 'height', 'tall');
|
||||||
#+end_src
|
#+end_src
|
||||||
|
|
||||||
#+name: fig:asd_noise_comp_bubble_aluminium
|
#+name: fig:asd_noise_comp_bubble_aluminium
|
||||||
|