Figure 1: Picture of the Setup
+
-
Figure 2: Zoom on the APA
-
1 Setup
+
+
-
1 Setup
-
1.1 Parameters
+
+
-
1.1 Parameters
Ts = 1e-4; @@ -85,8 +85,8 @@
-
1.2 Filter White Noise
+
+
-
1.2 Filter White Noise
Glpf = 1/(1 + s/2/pi/500); @@ -98,13 +98,13 @@ Gz = c2d(Glpf, Ts, 'tustin');
-
2 Run Experiment and Save Data
+
+
-
2 Run Experiment and Save Data
-
2.1 Load Data
+
+
-
2.1 Load Data
data = SimulinkRealTime.utils.getFileScopeData('data/apa95ml.dat').data;
@@ -113,8 +113,8 @@ Gz = c2d(Glpf, Ts, 'tustin');
-
2.2 Save Data
+
+
-
2.2 Save Data
u = data(:, 1); % Input Voltage [V] @@ -131,16 +131,16 @@ t = data(:, 3); % Time [s]
-
3 Huddle Test
+
+
-
3 Huddle Test
-
3.1 Time Domain Data
+
+
3.1 Time Domain Data
-
+
-
Figure 3: Measurement of the Mass displacement during Huddle Test
@@ -148,8 +148,8 @@ t = data(:, 3); % Time [s]
-
3.2 PSD of Measurement Noise
+
+
-
3.2 PSD of Measurement Noise
Ts = t(end)/(length(t)-1); @@ -165,7 +165,7 @@ win = hanning(ceil(1*Fs));
+
-
Figure 4: Amplitude Spectral Density of the Displacement during Huddle Test
@@ -174,16 +174,16 @@ win = hanning(ceil(1*Fs));
-
4 Transfer Function Estimation with m=5kg
+
+
-
4 Transfer Function Estimation with m=5kg
-
4.1 Time Domain Data
+
+
4.1 Time Domain Data
-
+
-
Figure 5: Time domain signals during the test
@@ -191,8 +191,8 @@ win = hanning(ceil(1*Fs));
-
4.2 Comparison of the PSD with Huddle Test
+
+
-
4.2 Comparison of the PSD with Huddle Test
Ts = t(end)/(length(t)-1); @@ -209,7 +209,7 @@ win = hanning(ceil(1*Fs));
+
-
Figure 6: Comparison of the ASD for the identification test and the huddle test
@@ -217,8 +217,8 @@ win = hanning(ceil(1*Fs));
-
4.3 Compute TF estimate and Coherence
+
+
-
4.3 Compute TF estimate and Coherence
win = hann(ceil(1/Ts)); @@ -229,30 +229,25 @@ win = hanning(ceil(1*Fs));
+
-
Figure 7: Coherence
+
-
-
Figure 8: Estimation of the transfer function from input voltage to displacement
-
Glpf = 1/(1 + s/2/pi/500); --
-
diff --git a/index.org b/index.org
index 828562a..aa5688a 100644
--- a/index.org
+++ b/index.org
@@ -290,29 +290,3 @@
#+caption: Estimation of the transfer function from input voltage to displacement
#+RESULTS:
[[file:figs/apa95ml_5kg_10V_tf.png]]
-
-#+begin_src matlab
- Glpf = 1/(1 + s/2/pi/500);
-#+end_src
-
-#+begin_src matlab :exports none
- figure;
- ax1 = subplot(2, 1, 1);
- hold on;
- plot(f, abs(tf_est./squeeze(freqresp(Glpf, f, 'Hz'))), 'k-')
- set(gca, 'Xscale', 'log'); set(gca, 'Yscale', 'log');
- ylabel('Amplitude'); xlabel('Frequency [Hz]');
- hold off;
-
- ax2 = subplot(2, 1, 2);
- hold on;
- plot(f, 180/pi*(angle(tf_est./squeeze(freqresp(Glpf, f, 'Hz')))), 'k-')
- set(gca, 'Xscale', 'log'); set(gca, 'Yscale', 'lin');
- ylabel('Phase'); xlabel('Frequency [Hz]');
- hold off;
- ylim([-180, 180]);
- yticks([-180, -90, 0, 90, 180]);
-
- linkaxes([ax1,ax2], 'x');
- xlim([10, 5e3]);
-#+end_src
Created: 2020-07-20 lun. 13:17
+Created: 2020-07-20 lun. 13:22