Skip to content
Snippets Groups Projects
Commit 632ebd77 authored by Pier Francesco A. Bachini's avatar Pier Francesco A. Bachini
Browse files

[ADA] Small changes after rebase

parent 6d85a84c
No related branches found
No related tags found
1 merge request!50Ada majority voting development
......@@ -11,8 +11,8 @@ settings.fault_sim.N_faulty_sensors = -1; % if set to -1 it will go to manual fa
settings.fault_sim.fault_type = ["no fault", "no fault", "no fault"];
if settings.fault_sim.N_faulty_sensors == -1
settings.fault_sim.selected_sensors = [1];
settings.fault_sim.fault_type = ["freezing", "no fault", "no fault"];
settings.fault_sim.selected_sensors = [];
settings.fault_sim.fault_type = ["no fault", "no fault", "no fault"];
end
% fault parameters
......
......@@ -89,7 +89,7 @@ end
if ~exist("../commonFunctions/graphics/general-utilities/","dir")
warning('To export file you need to download the repository, read the README file in the folder')
end
std_plots(simOutput,settings,contSettings,mission,rocket,environment)
std_plots(simOutput,settings,contSettings,mission,environment)
sensor_plots(simOutput, environment, rocket, settings);
% report_plots(simOutput,settings,contSettings)
......
function std_plots(simOutput, settings,contSettings,mission,rocket,environment)
function std_plots(simOutput, settings,contSettings,mission,environment)
if ~exist("report_images\"+mission.name,"dir")
mkdir("report_images\"+mission.name)
......@@ -70,7 +70,7 @@ end
drawnow
%% ada
figures.ada = figure('Position',[100,100,600,400]);
figures.ada = figure('Name', 'ADA vs Trajectory');
hold on; grid on;
for ii = 1:contSettings.ADA_N_instances
plot(simOutput.sensors.ada.time, simOutput.sensors.ada.data{ii}.xv(:,1),'DisplayName',strcat('$ADA\_', num2str(ii), '_{z}$'));
......@@ -95,7 +95,7 @@ drawnow
if contSettings.run_old_ada
figure('Position', [100,100,600,400]);
figures.ADAComp = figure('Name', 'ADA Comparisons');
hold on; grid on;
plot(simOutput.sensors.ada.time, simOutput.sensors.old_ada.xv(:,1), 'LineWidth', 1.5, 'DisplayName', "$run\_ADA_{z}$");
......@@ -106,23 +106,25 @@ if contSettings.run_old_ada
end
legend("Interpreter", "latex");
title("Comparison between run\_ADA and majority voting ADA");
drawnow
figure('Position', [100,100,600,400]);
figures.ADAErr = figure('Name', 'ADA Absolute error wrt run_ADA');
subplot(2,1,1); hold on; grid on;
subplot(2,1,2); hold on; grid on;
for ii = 1:contSettings.ADA_N_instances
subplot(2,1,1);
plot(simOutput.sensors.ada.time, simOutput.sensors.old_ada.xv(:,1) - simOutput.sensors.ada.data{ii}.xv(:,1), 'DisplayName', strcat('$ADA\_', num2str(ii), '_{z}$'));
disp("ADA " + num2str(ii) + " mean z error: " + num2str(mean(simOutput.sensors.old_ada.xv(:,1) - simOutput.sensors.ada.data{ii}.xv(:,1))) + " m");
disp("ADA " + num2str(ii) + " std z error: " + num2str(std(simOutput.sensors.old_ada.xv(:,1) - simOutput.sensors.ada.data{ii}.xv(:,1))) + " m");
% disp("ADA " + num2str(ii) + " mean z error: " + num2str(mean(simOutput.sensors.old_ada.xv(:,1) - simOutput.sensors.ada.data{ii}.xv(:,1))) + " m");
% disp("ADA " + num2str(ii) + " std z error: " + num2str(std(simOutput.sensors.old_ada.xv(:,1) - simOutput.sensors.ada.data{ii}.xv(:,1))) + " m");
subplot(2,1,2);
plot(simOutput.sensors.ada.time, simOutput.sensors.old_ada.xv(:,2) - simOutput.sensors.ada.data{ii}.xv(:,2), 'DisplayName', strcat('$ADA\_', num2str(ii), '_{vz}$'));
disp("ADA " + num2str(ii) + " mean vz error: " + num2str(mean(simOutput.sensors.old_ada.xv(:,2) - simOutput.sensors.ada.data{ii}.xv(:,2))) + " m/s");
disp("ADA " + num2str(ii) + " std vz error: " + num2str(std(simOutput.sensors.old_ada.xv(:,2) - simOutput.sensors.ada.data{ii}.xv(:,2))) + " m/s");
% disp("ADA " + num2str(ii) + " mean vz error: " + num2str(mean(simOutput.sensors.old_ada.xv(:,2) - simOutput.sensors.ada.data{ii}.xv(:,2))) + " m/s");
% disp("ADA " + num2str(ii) + " std vz error: " + num2str(std(simOutput.sensors.old_ada.xv(:,2) - simOutput.sensors.ada.data{ii}.xv(:,2))) + " m/s");
end
subplot(2,1,1); legend("Interpreter", "latex");
subplot(2,1,2); legend("Interpreter", "latex");
sgtitle("Absolute error between run\_ADA and majority voting ADA instances");
drawnow
end
%% reference
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment