How to plot more than one graph on a single tile using tiledlayout?

Question

Ron_S 2022-4-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1688659-how-to-plot-more-than-one-graph-on-a-single-tile-using-tiledlayout

评论： Star Strider 2022-4-7

I'm solving 17 species through ode45 (using a for loop with three different initial values). I want to plot all three solutions of a single species on one tile using tiledlayout. With the code I currently have it only plots the third iteration of each species. Could you please provide your expert advise? Thank you!

% Numerically finding the all-positive equilibrium point

% Remove the contents of the workspace

clear

% Set the initial time domain for integration

tstart=0;tend=150;

tspan=[tstart,tend];

y0 = [1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21];

for k = 1:length(y0)

[t,y] = ode45(@Chol_model3,tspan,[y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k)]);

end

tiledlayout(2,2)

nexttile

plot(t,y(:,3),'g','LineWidth', 1,'LineStyle','-');

title('B')

xlabel('Time')

ylabel('Species concentration')

nexttile

plot(t,y(:,4),'k','LineWidth', 1,'LineStyle','-');

nexttile

plot(t,y(:,5),'r','LineWidth', 1,'LineStyle','-');

nexttile

plot(t,y(:,6),'b','LineWidth', 1,'LineStyle','-');

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2022-4-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1688659-how-to-plot-more-than-one-graph-on-a-single-tile-using-tiledlayout#answer_934754

在 MATLAB Online 中打开

I cannot run the posted code, however it is definitely possible to plot more than one series on a specific tile using the hold function —

t = 0:10;

tiledlayout(1,2)

nexttile

plot(t, rand(1,11))

hold on

plot(t, randn(1,11))

hold off

nexttile

plot(t, rand(1,11))

hold on

plot(t, randn(1,11))

hold off

.

13 个评论
显示 11更早的评论隐藏 11更早的评论

Ron_S 2022-4-5

My apologies, the ode45 solver uses the Chol_model3.m file with code posted below. (I have not been able to combine these two files into one script file yet).

function dydt = Chol_model3(t,y)

%Read in the parameters from a file

%Chol_parameters

theta=10;

eta=0.01;

mu=5;

gamma=1;

p1=1;

p2=0.8;

k1=0.08;

k2=0.05;

k3=0.1;

k4=0.5;

k5=0.08;

k6=0.01;

k7=0.01;

k8=0.05;

k9=0.1;

k10=10;

k11=0.1;

k12=1;

k13=1;

k14=0.5;

k15=1;

k16=1;

k17=0.01;

k18=1;

R1=0.1;

R2=0.1;

C=y(1);Bc=y(2);B=y(3);Ce=y(4);S=y(5);L=y(6);Ldl=y(7);H=y(8);I=y(9);Cp=y(10);

Cm=y(11);E=y(12);O=y(13);Io=y(14);X=y(15);Xo=y(16);A=y(17);D=y(18);

eqn1=theta*Ce-p1*B*C;

eqn2=p1*B*C-eta*Bc*I;

eqn3=mu-p1*B*C-p2*B;

eqn4=k1*Cp+k3*E-k2*Ce-k4*Ce-theta*Ce+k5*H;

eqn5=p2*B-k9*S;

eqn6=k9*S-k6*L-k15*L*D;

eqn7=k6*L-k7*Ldl;

eqn8=k9*S-k5*H-k16*Io*H;

eqn9=k9*S-k11*I*O-k17*I-eta*Bc*I;

eqn10=k2*Ce+k7*Ldl-k1*Cp-k8*Cp-R1;

eqn11=k8*Cp-k10*Cm - R2;

eqn12=k4*Ce-k3*E-k18*A*E;

eqn13=k10*Cm-k11*I*O-k12*X*O;

eqn14=k11*I*O-k16*Io*H;

eqn15=gamma-k12*X*O;

eqn16=k12*X*O-k13*Xo-k14*Xo;

eqn17=k13*Xo-k18*A*E;

eqn18=k14*Xo-k15*L*D;

dydt=[eqn1; eqn2; eqn3; eqn4; eqn5; eqn6; eqn7; eqn8; eqn9; eqn10; eqn11; eqn12;

eqn13; eqn14; eqn15; eqn16; eqn17; eqn18];

Ron_S 2022-4-5

编辑：Ron_S 2022-4-5

Without the tiledlayout code, I can get the three graphs to show:

y0 = [1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21];

for k = 1:length(y0)

[t,y] = ode45(@Chol_model3,tspan,[y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k)]);

plot(t,y(:,3),'g','LineWidth', 1,'LineStyle','-');

title('B')

xlabel('Time')

ylabel('Species concentration')

end

If I add the tiledlayout code it only plots one graph. I'm not sure why?

y0 = [1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21];

for k = 1:length(y0)

[t,y] = ode45(@Chol_model3,tspan,[y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k)]);

tiledlayout (2,2)

nexttile

plot(t,y(:,3),'g','LineWidth', 1,'LineStyle','-');

hold on

title('B')

xlabel('Time')

ylabel('Species concentration')

end

hold off

Star Strider 2022-4-5

在 MATLAB Online 中打开

I am still getting used to tiledlayout since I do not use it often.

The first of these draws a plot in each tile, however when I run them together, the second tiledlayout call clears the first set, so I leave the first set to you to experiment with. (This is unlike the subplot function that can create separate subplots in a single figure window. There might be a way to do that with tiledlayout, however I cannot find any documentation for it.) Run them separately to see the results of each.

The problem is that your original code only plots one integration in each tile, similar to my first set here. So the hold call has no effect.

The second creates a single tile and plots all the differential equation solutions in it using hold.

I am not certain what you are doing, however it is definitely possible to plot multiple series in a single tile using hold. It would certainly be possible to do that here, except the code you posted onlly plots one column of the differential equation integration in each iteration of the for loop.

% With 'tiledlayout' & Three Axes —

tspan = [0 150];

% y0 = [1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21];

tiledlayout(3,1)

y0 = 1:10:21;

for k = 1:length(y0)

y0k = ones(1,18)*y0(k);

% [t,y] = ode45(@Chol_model3,tspan,[y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k)]);

[t,y] = ode45(@Chol_model3,tspan,y0k);

nexttile

plot(t,y(:,3),'g','LineWidth', 1,'LineStyle','-');

grid

title('B')

xlabel('Time')

ylabel('Species concentration')

end

% With 'tiledlayout' & One Axes —

tspan = [0 150];

% y0 = [1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21 1:10:21];

tiledlayout(1,1)

nexttile

hold on

y0 = 1:10:21;

for k = 1:length(y0)

y0k = ones(1,18)*y0(k);

% [t,y] = ode45(@Chol_model3,tspan,[y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k) y0(k)]);

[t,y] = ode45(@Chol_model3,tspan,y0k);

plot(t,y(:,3),'g','LineWidth', 1,'LineStyle','-');

grid

title('B')

xlabel('Time')

ylabel('Species concentration')

end

hold off

function dydt = Chol_model3(t,y)

%Read in the parameters from a file

%Chol_parameters

theta=10;

eta=0.01;

mu=5;

gamma=1;

p1=1;

p2=0.8;

k1=0.08;

k2=0.05;

k3=0.1;

k4=0.5;

k5=0.08;

k6=0.01;

k7=0.01;

k8=0.05;

k9=0.1;

k10=10;

k11=0.1;

k12=1;

k13=1;

k14=0.5;

k15=1;

k16=1;

k17=0.01;

k18=1;

R1=0.1;

R2=0.1;

C=y(1);Bc=y(2);B=y(3);Ce=y(4);S=y(5);L=y(6);Ldl=y(7);H=y(8);I=y(9);Cp=y(10);

Cm=y(11);E=y(12);O=y(13);Io=y(14);X=y(15);Xo=y(16);A=y(17);D=y(18);

eqn1=theta*Ce-p1*B*C;

eqn2=p1*B*C-eta*Bc*I;

eqn3=mu-p1*B*C-p2*B;

eqn4=k1*Cp+k3*E-k2*Ce-k4*Ce-theta*Ce+k5*H;

eqn5=p2*B-k9*S;

eqn6=k9*S-k6*L-k15*L*D;

eqn7=k6*L-k7*Ldl;

eqn8=k9*S-k5*H-k16*Io*H;

eqn9=k9*S-k11*I*O-k17*I-eta*Bc*I;

eqn10=k2*Ce+k7*Ldl-k1*Cp-k8*Cp-R1;

eqn11=k8*Cp-k10*Cm - R2;

eqn12=k4*Ce-k3*E-k18*A*E;

eqn13=k10*Cm-k11*I*O-k12*X*O;

eqn14=k11*I*O-k16*Io*H;

eqn15=gamma-k12*X*O;

eqn16=k12*X*O-k13*Xo-k14*Xo;

eqn17=k13*Xo-k18*A*E;

eqn18=k14*Xo-k15*L*D;

dydt=[eqn1; eqn2; eqn3; eqn4; eqn5; eqn6; eqn7; eqn8; eqn9; eqn10; eqn11; eqn12;

eqn13; eqn14; eqn15; eqn16; eqn17; eqn18];

end

.

Ron_S 2022-4-7

Hopefully it will be accepted! I'll keep you posted.

Star Strider 2022-4-7

Thank you!

If you also send it to me as an email through my profile page, please provide a context so I’ll know what it refers to, in addition to including a link to this thread. If it’s just a link, I’ll probably not recognise it, and I’m averse to clicking on uinknown links.

请先登录，再进行评论。