Skip to content
Snippets Groups Projects
Commit 470af0e0 authored by Claudio Mandrioli's avatar Claudio Mandrioli
Browse files

add run many

parent b6bb27ae
No related branches found
No related tags found
No related merge requests found
%"monteCarlo" example
%% Load data
disp('Loads data')
load('GNSSaidedINS_data.mat');
%% Load filter settings
disp('Loads settings')
settings=get_settings();
settings.P_treshold=0.1;
%run first simulation
out_data=GPSaidedINS(in_data,settings);
energy=out_data.energy;
error=out_data.error;
j=1;
for P_treshold=[3,6]%[2.3,2.8,3.3,3.8,4.3,4.8,5.3,5.8,6.3,6.8,7.3]
%clear energy error
settings.P_treshold=P_treshold;
% %run first simulation
% out_data=GPSaidedINS(in_data,settings);
% %energy=out_data.energy;
% %error=out_data.error;
% energy=[energy,out_data.energy];
% error=[error,out_data.error];
%run more simulations
for k=1:20
%% Run the GNSS-aided INS
disp('Running the GNSS-aided INS')
P_treshold
k
out_data=GPSaidedINS(in_data,settings);
energy=[energy,out_data.energy];
error=[error,out_data.error];
end
% %plot data
% figure(j)
% scatter(energy, error)
% xlabel('energy')
% ylabel('error')
% title(['P treshold = ', num2str(P_treshold)])
% grid
% j=j+1;
end
%plot data
figure(2)
scatter(energy, error)
xlabel('energy')
ylabel('error')
grid
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment