exp_240508_allanDeviationDifferentAxes.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;
 
%% add functions path and get data path
datapath = getDataPath();
datapath = fullfile(datapath, "self_captured/240508_allan_deviation_different_axes");
 
%% load data
 
% check if already a .mat file exists
if ~exist(fullfile(datapath, "data.mat"), 'file')
 
    disp("Loading data from BIN files and saving as .mat file...");
 
    % load data
    data(1) = jumpReadData(fullfile(datapath, "S-24-A-13026_000.BIN"));
    data(2) = jumpReadData(fullfile(datapath, "S-24-A-13036_000.BIN"));
    data(3) = jumpReadData(fullfile(datapath, "S-24-A-13169_000.BIN"));
    
    % save data
    save(fullfile(datapath, "data.mat"), 'data');
 
    disp("Data loaded and saved as .mat file.");
else
    % load data
 
    disp("Loading data from .mat file...");
 
    load(fullfile(datapath, "data.mat"));
 
    disp("Data loaded.");
end
 
%% plot acceleration data
figure;
for i = 1:3
    subplot(3,1,i);
    plot(data(i).tt.t, data(i).tt.acc);
    title("Accel. data " + num2str(i));
    xlabel("Time [s]");
    ylabel("Acceleration [m/s^2]");
end
 
%% parameter analysis
param(1) = jumpNoiseAnalysis(data(1), "plot", true, "toAnalyze","acc");
param(2) = jumpNoiseAnalysis(data(2), "plot", true, "toAnalyze","acc");
param(3) = jumpNoiseAnalysis(data(3), "plot", true, "toAnalyze","acc");