% Experimental data used in the case study (section 6) of 'On the 
% completeness of interface descriptions and the consistency of blocked 
% forces obtained in-situ'
close all
clear all
clc
%% Load data
rigidH = load('rigidH.mat'); % Acelerance FRF matrix measured on the rigidly
% coupled assembly (see figure 12a in the acompanying paper). Dimensions
% are 26 (response positions) x 29 (excitation positions) x 12801 (frequency bins) 
% Response positions 1:7 correspond to foot 1, 8:14 to foot 2, 15:21 to
% foot3, and 22:26 remote receiver DoFs.
% Excitation postions 1:21 correspond to the equivalent response positions,
% 22:29 are rtemote excitations on the source.
% The excitations and responses for each foot are ordered as so: z1, z2, z3, z4 (clockwise
% around each foot), x1, y1, y2.

resilientH = load('resilientH.mat'); % Acelerance FRF matrix measured on the rigidly
% coupled assembly (see figure 12b in the acompanying paper). Dimensions, excitation
% and response positions are the same as for rigidH.

rigidOppA = load('rigidOppA.mat'); % Operational time response measured on the rigidly
% coupled assembly. Dimensions are 26 (response positions) x 1671168 (time
% bins). Response positions corresponds to those of rigidH. Measured at a sample 
% rate of 25600 Hz

resilientOppA = load('resilientOppA.mat'); % Operational time response measured on the rigidly
% coupled assembly. Dimensions are 26 (response positions) x 1671168 (time
% bins). Response positions corresponds to those of resilientH. Measured at a sample 
% rate of 25600 Hz

freqV = load('freqVec.mat'); % Frequency vector for measured FRFs.