Compute sensor/actuator constants

matlab/frf_setup.m

@@ -27,7 +27,7 @@ V_noise = generateShapedNoise('Ts', 1/Fs, ...
 %% Sweep Sine
 gc = 0.1;
 xi = 0.5;
-wn = 2*pi*92.7;
+wn = 2*pi*94.3;
 
 % Notch filter at the resonance of the APA
 G_sweep = 0.2*(s^2 + 2*gc*xi*wn*s + wn^2)/(s^2 + 2*xi*wn*s + wn^2);
@@ -41,15 +41,6 @@ V_sweep = generateSweepExc('Ts',      Ts, ...
                            'sweep_type',   'log', ...
                            'V_exc', G_sweep*1/(1 + s/2/pi/500));
 
-V_sweep_lf = generateSweepExc('Ts',      Ts, ...
-                           'f_start', 0.1, ...
-                           'f_end',   10, ...
-                           'V_mean',  3.25, ...
-                           't_start', Trec_start, ...
-                           'exc_duration', Trec_dur, ...
-                           'sweep_type',   'log', ...
-                           'V_exc', 0.2);
-
 %% High Frequency Shaped Noise
 [b,a] = cheby1(10, 2, 2*pi*[300 2e3], 'bandpass', 's');
 wL = 0.005*tf(b, a);
@@ -70,19 +61,6 @@ V_sin = generateSinIncreasingAmpl('Ts', 1/Fs, ...
                                   't_start', Trec_start, ...
                                   'smooth_ends', true);
 
-%% Zero Excitation
-% Trec_start = 10;  % Start time for Recording [s]
-% Trec_dur   = 40; % Recording Duration [s]
-% 
-% Tsim = 2*Trec_start + Trec_dur;  % Simulation Time [s]
-
-V_zero = generateShapedNoise('Ts', 1/Fs, ...
-                             'V_mean', 3.25, ...
-                             't_start', Trec_start, ...
-                             'exc_duration', Trec_dur, ...
-                             'smooth_ends', true, ...
-                             'V_exc', tf(0));
-
 %% Select the excitation signal
 V_exc = timeseries(V_noise(2,:), V_noise(1,:));
 
@@ -91,12 +69,12 @@ figure;
 tiledlayout(1, 2, 'TileSpacing', 'Normal', 'Padding', 'None');
 
 ax1 = nexttile;
-plot(V_exc.Time, squeeze(V_exc.Data));
+plot(V_exc(1,:), V_exc(2,:));
 xlabel('Time [s]'); ylabel('Amplitude [V]');
 
 ax2 = nexttile;
-win = hanning(floor(length(squeeze(V_exc.Data))/8));
-[pxx, f] = pwelch(squeeze(V_exc.Data), win, 0, [], Fs);
+win = hanning(floor(length(V_exc)/8));
+[pxx, f] = pwelch(V_exc(2,:), win, 0, [], Fs);
 plot(f, pxx)
 xlabel('Frequency [Hz]'); ylabel('Power Spectral Density [$V^2/Hz$]');
 set(gca, 'xscale', 'log'); set(gca, 'yscale', 'log');