exp_240207_allanDeviation4.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
addpath(genpath('../../functions'));
datapath = getDataPath();
datapath = fullfile(datapath, "self_captured/240207_allan_deviation_4");
 
%% load data
data1 = jumpReadData(fullfile(datapath, "S-24-A-13026_000.BIN"), "time_correction",false);
data2 = jumpReadData(fullfile(datapath, "S-24-A-13036_000.BIN"), "time_correction",false);
data3 = jumpReadData(fullfile(datapath, "S-24-A-13169_000.BIN"), "time_correction",false);
 
% remove start because of temperature decrease and end because of me
% walking
data1.tt = data1.tt(4e5:end - 200 * 60,:);
data2.tt = data2.tt(4e5:end - 200 * 60,:);
data3.tt = data3.tt(4e5:end - 200 * 60,:);
 
 
%% parameter analysis
param14 = jumpNoiseAnalysis(data1, "plot", true);
param24 = jumpNoiseAnalysis(data2, "plot", true);
param34 = jumpNoiseAnalysis(data3, "plot", true);
 
%% average parameters
%accelerometer
acc.N = median([param14.acc.N, param24.acc.N, param34.acc.N]); % NoiseDensity [(m/s^2)/sqrt(Hz)]
acc.K = median([param14.acc.K, param24.acc.K, param34.acc.K]); % RandomWalk [m/s^2*sqrt(Hz)]
acc.B = median([param14.acc.B, param24.acc.B, param34.acc.B]); % BiasInstability [m/s^2]
 
%gyroscope
gyr.N = median([param14.gyr.N, param24.gyr.N, param34.gyr.N]); % NoiseDensity [°/s/sqrt(Hz)]
gyr.K = median([param14.gyr.K, param24.gyr.K, param34.gyr.K]); % RandomWalk [°/s*sqrt(Hz)]
gyr.B = median([param14.gyr.B, param24.gyr.B, param34.gyr.B]); % BiasInstability [°/s]
 
%magnetometer
mag.N = median([param14.mag.N, param24.mag.N, param34.mag.N]); % NoiseDensity [µT/sqrt(Hz)]
mag.K = median([param14.mag.K, param24.mag.K, param34.mag.K]); % RandomWalk [µT*sqrt(Hz)]
mag.B = median([param14.mag.B, param24.mag.B, param34.mag.B]); % BiasInstability [µT]
 
%high-g accelerometer
hig.N = median([param14.hig.N, param24.hig.N, param34.hig.N]); % NoiseDensity [(m/s^2)/sqrt(Hz)]
hig.K = median([param14.hig.K, param24.hig.K, param34.hig.K]); % RandomWalk [m/s^2*sqrt(Hz)]
hig.B = median([param14.hig.B, param24.hig.B, param34.hig.B]); % BiasInstability [m/s^2]
 
%temperature
tmp.N = median([param14.tmp.N, param24.tmp.N, param34.tmp.N]); % NoiseDensity [°C/sqrt(Hz)]
tmp.K = median([param14.tmp.K, param24.tmp.K, param34.tmp.K]); % RandomWalk [°C*sqrt(Hz)]
tmp.B = median([param14.tmp.B, param24.tmp.B, param34.tmp.B]); % BiasInstability [°C]
 
%pressure
prs.N = median([param14.prs.N, param24.prs.N, param34.prs.N]); % NoiseDensity [Pa/sqrt(Hz)]
prs.K = median([param14.prs.K, param24.prs.K, param34.prs.K]); % RandomWalk [Pa*sqrt(Hz)]
prs.B = median([param14.prs.B, param24.prs.B, param34.prs.B]); % BiasInstability [Pa]
 
%% save the average parameters to a file
fid = fopen('res_240207_allanResults.txt', 'w');
fprintf(fid, "*********** Results 2nd Allan Deviation Experiment ***********\n");
fprintf(fid, "Accelerometer:\n");
fprintf(fid, "NoiseDensity: %f [(m/s^2)/sqrt(Hz)]\n", acc.N);
fprintf(fid, "RandomWalk: %f [m/s^2 *sqrt(Hz)]\n", acc.K);
fprintf(fid, "BiasInstability: %f [m/s^2]\n\n", acc.B);
fprintf(fid, "Gyroscope:\n");
fprintf(fid, "NoiseDensity: %f [(°/s)/sqrt(Hz)]\n", gyr.N);
fprintf(fid, "RandomWalk: %f [°/s*sqrt(Hz)]\n", gyr.K);
fprintf(fid, "BiasInstability: %f [°/s]\n\n", gyr.B);
fprintf(fid, "Magnetometer:\n");
fprintf(fid, "NoiseDensity: %f [uT/sqrt(Hz)]\n", mag.N);
fprintf(fid, "RandomWalk: %f [uT*sqrt(Hz)]\n", mag.K);
fprintf(fid, "BiasInstability: %f [uT]\n\n", mag.B);
fprintf(fid, "High-g Accelerometer:\n");
fprintf(fid, "NoiseDensity: %f [(m/s^2)/sqrt(Hz)]\n", hig.N);
fprintf(fid, "RandomWalk: %f [m/s^2*sqrt(Hz)]\n", hig.K);
fprintf(fid, "BiasInstability: %f [m/s^2]\n\n", hig.B);
fprintf(fid, "Temperature:\n");
fprintf(fid, "NoiseDensity: %f [°C/sqrt(Hz)]\n", tmp.N);
fprintf(fid, "RandomWalk: %f [°C*sqrt(Hz)]\n", tmp.K);
fprintf(fid, "BiasInstability: %f [°C]\n\n", tmp.B);
fprintf(fid, "Pressure:\n");
fprintf(fid, "NoiseDensity: %f [mbar/sqrt(Hz)]\n", prs.N);
fprintf(fid, "RandomWalk: %f [mbar*sqrt(Hz)]\n", prs.K);
fprintf(fid, "BiasInstability: %f [mbar]\n", prs.B);
fprintf(fid, "*************************************************************\n");
fclose(fid);