Commit 42b7b5c2 authored by Claudio Mandrioli's avatar Claudio Mandrioli
Browse files

polished, commented and introduced print statements in run_many_cycling file

parent f3feffa7
%"monteCarlo" example
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Script for example of energy-accuracy trade-off evaluatio with the
% cycling trace
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Load data
disp('Loads data')
......@@ -12,29 +15,26 @@ settings=get_settings_cycling();
color=[[0.7,0.0,0.0];[0.0,0.7,0.0];[0.7,0.0,0.7];[0.0,0.0,0.7];[0.0,0.7,0.7];[0.7,0.7,0.0]];
%vector of refrences for P
P_treshold=[4.0,5.0,6.0,8.0,10.0,12.0];
%nubmer of runs per each P
n_run=1;
%nubmer of simulations per each P
n_run=3;
%allocate memory for output
energy=zeros(n_run,length(P_treshold));
error=zeros(n_run,length(P_treshold));
for j=1:length(P_treshold)
%set P reference
settings.P_treshold=P_treshold(j);
%run simulations
for k=1:n_run
%% Run the GNSS-aided INS
disp('Running the GNSS-aided INS')
P_treshold(j)
k
out_data=GPSaidedINS_cycling(in_data,settings);
energy(k,j)=out_data.energy;
error(k,j)=out_data.error;
%% Iterate over the threshold for tr(P)
for j=1:length(P_treshold)
disp(['Starting simulations with threshold for P equal to:', num2str(P_treshold(j))]) %display reached treshold
settings.P_treshold=P_treshold(j); %set P reference
for k=1:n_run %iterate over the number of simulations per threshold
disp(['Simulation number:',num2str(k)]) %display simulation number
out_data=GPSaidedINS_cycling(in_data,settings); % Run the GNSS-aided INS
energy(k,j)=out_data.energy; % store energy
error(k,j)=out_data.error; % store positioning error
end
end
%plot data
%% Plot data
figure
hold on
xlabel('energy')
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment