exp_240213_pressureMembrane.m

Go to the documentation of this file.

%
% For LICENSE / TERMS OF USE, see the file LICENSE in the root directory of the repository.
 
clear; close all; clc;
 
%% get data path and files
datapath = getDataPath();
datapath = fullfile(datapath, "self_captured/240213_pressure_membrane");
 
% get all the .BIN files in the directory
files = dir(fullfile(datapath, "*.BIN"));
 
% import the data
for i   = 1:length(files)
data(i) = jumpReadData(fullfile(datapath, files(i).name));
end
 
% downsample the data to 5 Hz (actual pressure data is 5 Hz)
for  i     = 1:length(data)
data(i).tt = retime(data(i).tt, "regular", "nearest", "SampleRate", 5);
end
 
%% find risetimes
% sensor 1 (indices found by hand)
start_idx_1 = [  1, 163, 326, 518, 693,  892, 1061, 1267];
stop_idx_1  = [146, 289, 453, 670, 868, 1044, 1227, 1347];
 
% sensor 2 (indices found by hand)
start_idx_2 = [  3, 157, 296, 490, 685,  883, 1057, 1262];
stop_idx_2  = [147, 281, 448, 676, 871, 1043, 1217, 1345];
 
% sensor 3 (indices found by hand)
start_idx_3 = [  6, 158, 298, 506, 691,  885, 1056, 1262];
stop_idx_3  = [148, 281, 448, 673, 871, 1046, 1224, 1350];
 
start_idx = [start_idx_1; start_idx_2; start_idx_3];
stop_idx  = [stop_idx_1; stop_idx_2; stop_idx_3];
 
% find risetimes
 
debug = false;
 
for i = 1:length(data)
    for j=1:length(start_idx)
        risetimes(i,j) = risetime(data(i).tt.prs(start_idx(i,j):stop_idx(i,j)),5);
        if debug
            risetime(data(i).tt.prs(start_idx(i,j):stop_idx(i,j)),5);
            title(sprintf("Risetime %0.2f s", round(risetimes(j),2)));
            ylabel("Level (mbar)")
        end
    end
end
 
risetimes = reshape(risetimes, [1, numel(risetimes)]);
 
stem(risetimes);
 
% result: it seems that the risetimes are not very consistent, but the
% average risetime is around 0.52 s. for the model I'm assuming a mean risetime
% of 0.6 s with a standard deviation of 0.2 s.