why is it not graphing my plot? im getting the data but no data points on my graph

Question

Ivan Morales 2018-7-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/411484-why-is-it-not-graphing-my-plot-im-getting-the-data-but-no-data-points-on-my-graph

评论： Star Strider 2018-7-22

function prob8=project8(~,V)
Fa=V(1);
Fb=V(2);
Fc=V(3);
Fd=V(4);
Fe=V(5);
Fg=V(6);
Fw=V(7);
k1=0.014;
k2=0.007;
k3=0.014;
k4=0.45;
% molar flow rates in mol/s
Fao=10;
Fbo=5;
%volumetric flow rate dm^3/s
vo=100;
%total concentration in mol/dm^3
Cto=0.147;
%volume in dm^3
V=1000;
Ft=Fao+Fbo;
%Total molar flow%
Ft=Fa+Fb+Fc+Fd+Fe+Fg+Fw;
%defining each concentration%
Ca=Cto*(Fa/Ft);
Cb=Cto*(Fb/Ft);
Cc=Cto*(Fc/Ft);
Cd=Cto*(Fd/Ft);
Ce=Cto*(Fe/Ft);
Cg=Cto*(Fg/Ft);
Cw=Cto*(Fw/Ft);
ra=(-k1*Ca*(Cb^(1/2)))-(k2*Ca^2);
rb=-0.5*k1*Ca*(Cb^(1/2));
rc=(k1*Ca*(Cb^(1/2)))-(k3*Cc)-(k4*Cd);
rd=(0.5*k2*((Ca^2)))-(k4*Cd);
re=k3*Cc;
rg=k4*Cd;
rw=k3*Cc;
prob8(1)=ra;
prob8(2)=rb;
prob8(3)=rc;
prob8(4)=rd;
prob8(5)=re;
prob8(6)=rg;
prob8(7)=rw;
prob8=prob8';
%calling the functiin project8 to solve for each value%
prob8o=[10 5 0 0 1*10^(-16) 0 1*10^(-16)];
Vspan=[0 1000];
[V,y]=ode45('project8',Vspan, prob8o)
plot(V,prob8o(:,1));
hold on 
plot(V,prob8o(:,2));
hold on
plot(V,prob8o(:,3));
hold on
plot(V,prob8o(:,4));
hold on
plot(V,prob8o(:,5));
hold on 
plot(V,prob8o(:,6));
hold on
plot(V,prob8o(:,7));
legend('Fa','Fb','Fc','Fd','Fe','Fw','Fg')
title ('Flowrate Vs Volume of Reactor')
xlabel('Volume(dm3)')
ylabel('Flowrate (mol/s)')

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2018-7-22

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/411484-why-is-it-not-graphing-my-plot-im-getting-the-data-but-no-data-points-on-my-graph#answer_329775

在 MATLAB Online 中打开

First, you need to call the function correctly, so use the ‘@’ operator to create a function handle in your ode45 call.

Second, you need to plot the correct variables. Those are the ‘V’ vector and the ‘y’ matrix your ode45 call returns.

Your ‘project8’ function works correctly, so you only need to change the call to it, and the variables you plot.

Try this:

% calling the function project8 to solve for each value %
prob8o=[10 5 0 0 1*10^(-16) 0 1*10^(-16)];
Vspan=[0 1000];
[V,y]=ode45(@project8,Vspan, prob8o);
plot(V,y)
legend('Fa','Fb','Fc','Fd','Fe','Fw','Fg')
title ('Flowrate Vs Volume of Reactor')
xlabel('Volume(dm3)')
ylabel('Flowrate (mol/s)')