Last else if graph not plotting

Question

ANMOL BAKSHI 2020-11-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/634269-last-else-if-graph-not-plotting

评论： ANMOL BAKSHI 2020-11-3

采纳的回答： S. Walter

if i==3600

hold on;

plot(x,C);

elseif i==10800

plot(x,C);

elseif i==36000

plot(x,C);

elseif i==36300

plot(x,C);

hold off;

end

In the above program, if i varies from 0 to 36300 or more then last plot is not happening and if I add one more else if i ==36600 then that will not get plot.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

S. Walter 2020-11-2

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/634269-last-else-if-graph-not-plotting#answer_532034

编辑：S. Walter 2020-11-2

在 MATLAB Online 中打开

I tried the following:

% Clear the command window and close the figures
clc
close all
% Make a color vector for plots
col_vec = lines(4);
% Start a counter for the legend entry
cnt = 1;
% Loop over the entries
for i = 1 : 1 : 36300
    % Make random data
    x = 0:1:10;
    C = rand(size(x));
    % If/else statements - added some extra code for plotting clarity
    if i==3600
    	hold on;
    	plot(x,C,'-o','MarkerFaceColor',col_vec(cnt,:))
    	leg{cnt} = 'first';
    	cnt = cnt+1;
    elseif i==10800
    	plot(x,C,'-d','MarkerFaceColor',col_vec(cnt,:))
    	leg{cnt} = 'second';
    	cnt = cnt+1;
    elseif i==36000
    	plot(x,C,'-s','MarkerFaceColor',col_vec(cnt,:))
    	leg{cnt} = 'third';
    	cnt = cnt+1;
    elseif i==36300
    	plot(x,C,'-^','MarkerFaceColor',col_vec(cnt,:))
    	leg{cnt} = 'fourth';
    	hold off;
    end          
end
% Annotate the plot
legend(leg,'location','northeast')
axis([0 10 0 1.1])
grid on
xlabel('x')
ylabel('C')

When I do that, I get four lines with markers

However, as the other commenter said, placing hold on and hold off ahead of the if/else statements is more robust, because if the condition is not met, you might not get the right hold conditions.

3 个评论
显示 1更早的评论隐藏 1更早的评论

ANMOL BAKSHI 2020-11-2

clc; clear;

lengt = input('The calculation length domain (in m) = ');

points = input('No. of grid points = ');

deltaX = (lengt./(points-1));

deltat = input('Time step (in seconds) = ');

diff = input('Diffusivity = ');

t = input('Time (in s) upto which calculation is to be done = ');

x = 0:deltaX*100:lengt*100;

for i = 0:deltat:t

if i==0

for m = 1:points

cp0(m) = 0.15;

end

else

for n = 1:points

cp0(n) = C(n);

end

for j = 2:points-1

diffe(j) = diff;

diffw(j) = diff;

ae(j) = (diffe(j)./deltaX);

aw(j) = (diffw(j)./deltaX);

ap0(j) = (deltaX./deltat);

ap(j) = (ae(j) + aw(j) + ap0(j));

P(j) = -aw(j);

Q(j) = ap(j);

R(j) = -ae(j);

d(j) = (cp0(j).*ap0(j));

end

P(1) = 0;

Q(1) = 1;

R(1) = 0;

d(1) = 1;

P(points) = 0;

Q(points) = 1;

R(points) = 0;

d(points) = 0.2;

Rp(1) = R(1)/Q(1);

dp(1) = d(1)/Q(1);

for k=2:points-1

Rp(k) = R(k)./(Q(k)-(P(k).*Rp(k-1)));

dp(k) = (d(k)-(P(k).*dp(k-1)))./(Q(k)-(P(k).*Rp(k-1)));

end

dp(points) = (d(points)-(P(points).*dp(points-1)))./(Q(points)-(P(points).*Rp(points-1)));

C(points) = dp(points);

for k=points-1:-1:1

C(k) = dp(k)-(Rp(k).*C(k+1));

end

if i==3600

hold on

plot(x,C);

C3600 = C;

elseif i==10800

plot(x,C);

C10800 = C;

elseif i==36000

plot(x,C);

C36000 = C;

elseif i==36300

plot(x,C);

C36300 = C;

hold off

end

Not working here.

S. Walter 2020-11-2

在 MATLAB Online 中打开

I tried to run your code, setting the following variables:

% The calculation length domain (in m)
lengt = 1;
% No. of grid points 
points = 100;
deltaX = (lengt./(points-1));
% Time step (in seconds)
deltat = 1;
% Diffusivity
diff = 1e-5;
% Maximum time (in seconds)
t = 370000;

And I get the following result

It's a bit difficult to see but the result of i = 36000 and i = 36300 are basically on top of each other. Below is a picture of the result zoomed in to axis([49 51 0.57 0.59]):

Just as a note, you may want to evaluate how you are doing your time step. You currently have your

and iteration in time as the same thing but your iterations must be integers. This line:

i = 0:deltat:t

ANMOL BAKSHI 2020-11-3

Yes, got it.

Really thanks for the suggestions and help.

请先登录，再进行评论。

Answer 2

KSSV 2020-11-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/634269-last-else-if-graph-not-plotting#answer_532014

在 MATLAB Online 中打开

Try like this:

hold on 
if i==3600
               
                plot(x,C);
                
              elseif i==10800
                plot(x,C);
                
              elseif i==36000
                plot(x,C);
               
              elseif i==36300
                plot(x,C);
              
                     
end          
hold off

Though I don't understand..what is the purpose to plot like this.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

ANMOL BAKSHI 2020-11-2

hold on

if i==3600

plot(x,C);

C3600 = C;

elseif i==10800

plot(x,C);

C10800 = C;

elseif i==36000

plot(x,C);

C36000 = C;

elseif i==36300

plot(x,C);

C36300 = C;

end

hold off

Tried that. Now, here C36300 is getting stored but plot for last else if is not happening.

请先登录，再进行评论。

Last else if graph not plotting

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

3 个评论
显示 1更早的评论隐藏 1更早的评论

更多回答（1 个）

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

另请参阅

类别

标签

产品

Community Treasure Hunt

Last else if graph not plotting

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

3 个评论 显示 1更早的评论隐藏 1更早的评论

更多回答（1 个）

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

另请参阅

类别

标签

产品

Community Treasure Hunt

WeChat

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

3 个评论
显示 1更早的评论隐藏 1更早的评论

1 个评论
显示 -1更早的评论隐藏 -1更早的评论