From 5e3a78ac2f85361b7515ea73880048d34294384e Mon Sep 17 00:00:00 2001 From: RicCada <riccardo.cadamuro@skywarder.eu> Date: Sun, 8 Oct 2023 17:46:03 +0200 Subject: [PATCH] Removed errored simulation --- .../postprocessing/AFD/mainPostprocess.m | 20 ++++--------------- 1 file changed, 4 insertions(+), 16 deletions(-) diff --git a/RoccarasoFlight/postprocessing/AFD/mainPostprocess.m b/RoccarasoFlight/postprocessing/AFD/mainPostprocess.m index 2788f77..5584079 100644 --- a/RoccarasoFlight/postprocessing/AFD/mainPostprocess.m +++ b/RoccarasoFlight/postprocessing/AFD/mainPostprocess.m @@ -61,15 +61,7 @@ if opt.flagCA settings.motor.expThrust = Thrust; settings.tb = settings.motor.expTime(end); settings.tControl = settings.tb; - - settings.mDotTotalTime = [diff(settings.motor.expM)./diff(settings.motor.expTime) 0]; - - [TfMDOT, YfMDOT, TaMDOT, YaMDOT, bound_valueMDOT] = stdRun(settings); - load ascent_plot.mat; - data_ascentMDOT = data_ascent; - clear data_ascent; - settings.mDotTotalTime = zeros(size(settings.motor.expThrust)); [Tf, Yf, Ta, Ya, bound_value] = stdRun(settings); load ascent_plot.mat; @@ -81,36 +73,33 @@ if opt.flagCA hold on plot(res.time, res.altitude.alt_NasSim); plot(Ta, -Ya(:, 3)) - plot(TaMDOT, -YaMDOT(:, 3)) xlabel('Time [s]'); ylabel('Altitude [m]'); title('Altitude'); grid on - legend('NAS - flight', 'NAS - simulated', 'Simulation','Simulation mDot', 'Location','southoutside', 'Orientation', 'horizontal'); + legend('NAS - flight', 'NAS - simulated', 'Simulation', 'Location','southoutside', 'Orientation', 'horizontal'); figure plot(res.time, res.velocity.vel_NasFlight); hold on plot(res.time, res.velocity.vel_NasSim); plot(Ta, vecnorm(Ya(:, 4:6), 2, 2)); - plot(TaMDOT, vecnorm(YaMDOT(:, 4:6), 2, 2)); xlabel('Time [s]'); ylabel('Velocity [m/s]'); title('Total Velocity'); grid on - legend('NAS - flight', 'NAS - simulated', 'Simulation', 'Simulation mDot','Location','southoutside', 'Orientation', 'horizontal'); + legend('NAS - flight', 'NAS - simulated', 'Simulation', 'Location','southoutside', 'Orientation', 'horizontal'); figure plot(res.time, res.velocity.velVert_NasFlight); hold on plot(res.time, res.velocity.velVert_NasSim); plot(Ta, -data_ascent.velocities(3, :)); - plot(TaMDOT, -data_ascentMDOT.velocities(3, :)); xlabel('Time [s]'); ylabel('Velocity [m/s]'); title('Vertical Velocity'); grid on - legend('NAS - flight', 'NAS - simulated', 'Simulation', 'Simulation mDot', 'Location','southoutside', 'Orientation', 'horizontal'); + legend('NAS - flight', 'NAS - simulated', 'Simulation','Location','southoutside', 'Orientation', 'horizontal'); figure @@ -119,12 +108,11 @@ if opt.flagCA plot(res.time(indexCO:end),res.CA.CA_NasSim(indexCO:end)); plot(res.time(indexCO:end), [res.CA.CA0(indexCO:indexOpeningARB) res.CA.CA100(indexOpeningARB+1:end)]); plot(data_ascent.integration.t, data_ascent.coeff.CA); - plot(data_ascentMDOT.integration.t, data_ascentMDOT.coeff.CA); xlabel('Time [s]'); ylabel('CA [-]'); title('CA from telemetry'); grid on - legend('NAS - flight', 'NAS - simulated', 'MD', 'Simulation', 'Simulation mDot', 'Location','southoutside', 'Orientation', 'horizontal'); + legend('NAS - flight', 'NAS - simulated', 'MD', 'Simulation', 'Location','southoutside', 'Orientation', 'horizontal'); end -- GitLab