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;