add simulink file for IFF

This commit is contained in:
Thomas Dehaeze 2021-06-10 11:38:37 +02:00
parent ac7e7b9ef4
commit fd8a62c819
4 changed files with 18 additions and 6 deletions

Binary file not shown.

View File

@ -82,20 +82,20 @@ V_sin = generateSinIncreasingAmpl('Ts', 1/Fs, ...
'smooth_ends', true);
%% Zero Excitation
% Trec_start = 10; % Start time for Recording [s]
% Trec_dur = 10; % Recording Duration [s]
%
% Tsim = 2*Trec_start + Trec_dur; % Simulation Time [s]
Trec_start = 10; % Start time for Recording [s]
Trec_dur = 10; % Recording Duration [s]
Tsim = 2*Trec_start + Trec_dur; % Simulation Time [s]
V_zero = generateShapedNoise('Ts', 1/Fs, ...
'V_mean', 3.25, ...
'V_mean', 1.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_hf(2,:), V_noise_hf(1,:));
V_exc = timeseries(V_zero(2,:), V_zero(1,:));
%% Plot
figure;

BIN
matlab/iff_measure.slx Normal file

Binary file not shown.

12
matlab/iff_setup.m Normal file
View File

@ -0,0 +1,12 @@
run('identif_setup.m');
%%
Kiff_g1 = (1/(s + 2*pi*40))*... % Low pass filter (provides integral action above 40Hz)
(s/(s + 2*pi*30))*... % High pass filter to limit low frequency gain
(1/(1 + s/2/pi/500))*... % Low pass filter to be more robust to high frequency resonances
eye(6); % Diagonal 6x6 controller
Kiff = c2d(Kiff_g1, Ts, 'tustin');
%%
save('./frf_data.mat', 'Fs', 'Ts', 'Tsim', 'Trec_start', 'Trec_dur', 'V_exc', 'Kiff');