Why do the lines not plot to completion

2 次查看（过去 30 天）
Mary Jean Savitsky
Mary Jean Savitsky 2020-12-2
评论： Mary Jean Savitsky 2020-12-3
Any ideas on why the red and yellow line are stopping? I want them to go all the way down but I cant figure out what the issue is. It is a graph of Effectiveness over observable modulus and I copied my code below. I tried removing the xlim but that didnt work. Would love to hear any ideas! Thank you!
%% first order line
observ_mod=logspace(-2,1);
n=tanh(observ_mod)./observ_mod;
figure(2)
loglog(observ_mod,n)
xlim([-2 10])
ylim([-1 1])
hold on
%% Beta=2 line
F=0.5
b_2=2
syms theta(x)
equation_41b_2=odeToVectorField(diff(theta,2) == (F*theta)/(1+b_2*theta)); % Solves the equation 10.6.41b from the textbook
vector2=matlabFunction(equation_41b_2, 'vars', {'x','Y'}); % Determines the vector
solution2=ode45(vector2, [0 1], [.39 0]);
factor2=((1+b_2)/(F))*(solution2.y(2,:)); % Effective Factor First Order 1
observ_mod2=(factor2*F)/(1+b_2); % Observable Modulus First Order 1
%% Beta=5 line
b_5=5
syms theta(x)
equation_41b_5=odeToVectorField(diff(theta,2) == (F*theta)/(1+b_5*theta)); % Solves the equation 10.6.41b from the textbook
vector5=matlabFunction(equation_41b_5, 'vars', {'x','Y'}); % Determines the vector
solution5=ode45(vector5, [0 1], [.39 0]);
factor5=((1+b_5)/(F))*(solution5.y(2,:)); % Effective Factor First Order 1
observ_mod5=(factor5*F)/(1+b_5); % Observable Modulus First Order 1
% observ_mod=logspace(-2,1);
% n=tanh(observ_mod)./observ_mod;
%% Plotting
% loglog(factor1,1-observ_mod1)
loglog(factor2,1-observ_mod2)
loglog(factor5,1-observ_mod5)
% xlim([10^-2 10])
ylim([0.1 1])
title('Effectiveness Factor vs. Observable Modulus')
xlabel('Observable Modulus')
ylabel('Effectiveness Factor, n')

请先登录，再回答此问题。

采纳的回答

Ive J
Ive J 2020-12-2
Increase integration interval in ode45:
solution2=ode45(vector2, [0 25], [.39 0]);
solution5=ode45(vector5, [0 25], [.39 0]);

更多回答（1 个）

Walter Roberson
Walter Roberson 2020-12-2
solution2=ode45(vector2, [0 1], [.39 0]);
factor2=((1+b_2)/(F))*(solution2.y(2,:)); % Effective Factor First Order 1
so factor2 will be set according to the output of the ode -- not according to the input of the ode.
loglog(factor2,1-observ_mod2)
and that factor2 will be used as the independent variable for the second plot, even though it is obviously a dependent variable. There is no solid reason ahead of time to expect that the output of the ode will have any particular range.
observ_mod2=(factor2*F)/(1+b_2); % Observable Modulus First Order 1
In practice your factor2 appear to be in strictly increasing order. Your observe_mod2 is calculated from that with a simple transformation. You use your factor2 as the x axis for your loglog() plot. That means that as long as factor2 is finely-enough sampled that you do not get visual artifacts, that you might as well instead just substitute factor2 values that you make up with linspace() or logspace() without ever having done the ode45(). In context, the only purpose of the ode45 is to establish the range of values to plot over, and you do not even want to use that range, so you might as well not even call ode45() .
  1 个评论
Mary Jean Savitsky
Mary Jean Savitsky 2020-12-3
I've tried both reccomendations, Im trying to get it to look like the plot below but for some reason my Beta=2 line keeps ending below x=10^1. It's not that big of an issue but just curious to see what more I could do to get it closer.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Numerical Integration and Differential Equations 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by