推進薬搭載率εと理想速度増分⊿Vの関係

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-lambda).*epsilon));
    plot(epsilon,deltaV);
end
legend(arrayfun(@(x) sprintf('lambda=%.3f',x), lambdas,UniformOutput=false));
title('単段式ロケットの推進薬搭載率と理想速度増分');
xlabel('推進薬搭載率ε');
ylabel('速度増分⊿V[m/s]');
big;