From 470af0e059338c385ac60e63a3f9108724dc85c9 Mon Sep 17 00:00:00 2001
From: Claudio Mandrioli <claudio.mandrioli@control.lth.se>
Date: Mon, 8 Oct 2018 09:23:03 +0200
Subject: [PATCH] add run many
---
matlab/run_many.m | 53 +++++++++++++++++++++++++++++++++++++++++++++++
1 file changed, 53 insertions(+)
create mode 100644 matlab/run_many.m
diff --git a/matlab/run_many.m b/matlab/run_many.m
new file mode 100644
index 0000000..b20d4c3
--- /dev/null
+++ b/matlab/run_many.m
@@ -0,0 +1,53 @@
+%"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
--
GitLab