Exponential Plot not acting correct

조회 수: 1 (최근 30일)
Christina Reid
Christina Reid . 2021년 5월 16일
편집: the cyclist . 2021년 5월 16일
I have a script that is running, but the exponential plot is not acting correctly. As the x axis increases, the y axis should as well and it seems like it is doing the oppoiste of that. I have attached the equation that is needed, in the 'ideal_rocket_equation' attachment. Does anyone know how I can fix the plot
%% Coursework Number 1:
% Basing on the "Ideal Rocket Flight Performance" theory compute and show
% through suitable curves the optimal (maximum) total payload ratio
% lamda_tot of a multisage rocket with identical stages (c_i = c and k_si =
% k_s) as a function of the number of stages N.
% delta_V/c equal to 0.5, 1.0, 2.0, 3.0, and 4.0
%k_s equal to 0.1
% For N = 1
ks_1 = 0.1;
%x = linspace(0,5.0);
x = [ 0.5, 1.0, 2.0, 3.0, 4.0];
delta_V = x;
line_color = ['b' 'g' 'y' 'c' 'm'];
stages = cell(1, length(line_color));
%n = 1:5;
for i = 1: length(line_color);
stages{i} = sprintf('Stage %d' ,i);
hold on
%get new values
lamda_totmax = ((exp(-delta_V./i) - ks_1)/(1-ks_1)).^i
%lamda_totmax = ((exp(-delta_V./i)- ks_1)/(1 - ks_1)).^i;
semilogy ( x, lamda_totmax, '-', 'Color', line_color(i), 'LineWidth', 2)
title (' Total payload ratio as a function of demensionless ideal velocity and Number of stages', 'FontSize', 16)
ylabel ( 'Payload Ratio (lamda_tot)')
xlabel ( 'Dimensionless Ideal Velocity (delta_V/Isp*g_0')
legend( stages, 'Location', 'southwest')

답변 (1개)

the cyclist
the cyclist 2021년 5월 16일
편집: the cyclist 님. 2021년 5월 16일
The plot behavior looks correct to me, according to the equation.
If delta_V == 0, then lambda == 1. (You don't calculate/plot this point, but that is the result you get.)
In the limit that delta_V gets large, lambda gets smaller, and goes to
ks_1 = 0.1;
-ks_1 / (1 - ks_1) % when N == 1
ans = -0.1111
This seems to be what your plot does. So, I don't see a problem. But maybe I am not understanding something.


Help CenterFile Exchange에서 Programming에 대해 자세히 알아보기

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by