How to iterate xa from 0 to 1 with steps of 0.05 in the function r=F*(k*C*(​1-xa))/(1+​K*C*(1-xa)​). All variables are known except for r.

5 次查看(过去 30 天)
I have to iterate the conversion of a reaction (x) in order to get the different values of catalyst mass (m).
The conversion 'x' would go from 0 to 1 with spaces of 0.05.
The function I need to iterate is: m=F*(k_kin*C*(1-x))/(1+k_ads*C*(1-x))
I know the value of all the variables except for m, which is the mass that I'd like to calculate.
F=0.11574; k_kin=1.3*10^-6; C=1.086; k_ads=0.9986
Please, help me. I haven't used matlab in years and I've forgotten how to make a 'for' loop work.
Thank you in advance.

采纳的回答

Star Strider
Star Strider 2019-5-23
Everything except ‘x’ are scalars, so you can do this using a vectorization approach without the loop:
F=0.11574; k_kin=1.3*10^-6; C=1.086; k_ads=0.9986;
x = 0:0.05:1;
m=F*(k_kin*C*(1-x))./(1+k_ads*C*(1-x));
figure
plot(x, m)
grid
xlabel('Reaction')
ylabel('Catalyst Mass')
... including the plot! (Note the (./) denoting element-wise division.)

更多回答(2 个)

Geoff Hayes
Geoff Hayes 2019-5-23
编辑:Geoff Hayes 2019-5-23
Marina - try
m = zeros(20,1);
k = 1;
for x = 0:0.05:1
m(k) = F*(k_kin*C*(1-x))/(1+k_ads*C*(1-x))
k = k + 1;
end
edit to above: repalce == with = in assignment to m
  8 个评论
Marina Gonzalez
Marina Gonzalez 2019-5-23
Thank you! That has really helped. I think I'm doing something wrong because for a conversion of 0, the mass of the catalyst is way bigger than for the total conversion.
Captura.PNG
I'm doing something wrong but I don't know what. I'll think about it.
Thank you once again!

请先登录,再进行评论。


Marina Gonzalez
Marina Gonzalez 2019-5-23
It's me again- sorry!
What if the function instead of being the one I said before was m(k) = F*(1+k_ads*C*(1-x))/(k_kin*C*(1-x))
Would anything change? Because when I plot this new equation nothing shows up on the graph
  3 个评论

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Data Preprocessing 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by