clear; close all; clc; figure(); hold on; a0_g = linspace(1,3.5,50); Lambdas = [2 2.5 3 3.5 4]; %質量比の逆数 y_pks = []; a0_g_peaks = []; for idx=1:length(Lambdas) Lambda = Lambdas(idx); y_g_Isp2 = 1./a0_g.*(1+... -1/Lambda*(1+log(Lambda)…

clear; close all; clc; C = 4413; %m/s lambdas = [0.025, 0.05]; %ペイロード比 2.5%, 5% epsilon = linspace(0.85,0.95,10); %推進薬搭載率 85%~95% figure(); hold on; for idx=1:length(lambdas) lambda = lambdas(idx); deltaV = C*log(1./(1-(1-lambd…