diff --git a/RoccarasoFlight/postprocessing/AFD/mainPostprocess.m b/RoccarasoFlight/postprocessing/AFD/mainPostprocess.m index 2788f77efc30fbd6e88a1b6edc5a6b64e8226041..55840790e67aa00f0863951ade240612623aef62 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