作业帮matlab带循环的反应器程序报错初学matlab,我前几天看到一个带循环操作的平推流反应器计算最优循环比下反应器最佳体积的例子,想要自己试着做一下,但是总是报错.想请教各位怎样修改这个程
来源：学生作业帮助网 编辑：作业帮 时间：2024/05/08 20:59:27
matlab带循环的反应器程序报错初学matlab,我前几天看到一个带循环操作的平推流反应器计算最优循环比下反应器最佳体积的例子,想要自己试着做一下,但是总是报错.想请教各位怎样修改这个程
matlab带循环的反应器程序报错
初学matlab,我前几天看到一个带循环操作的平推流反应器计算最优循环比下反应器最佳体积的例子,想要自己试着做一下,但是总是报错.想请教各位怎样修改这个程序.
自催化反应A+R——2R,其反应速率方程式为-rA=KcAcR,在70℃下等温的进行反应,在此温度下k=1.512m3/kmol.h,其他数据如下：cA0=0.99kmol/m3,cR0=0.01kmol/m3,vo=10m3/h,哟求反应的转化率xA=0.99,若采用带有循环操作的平推流反应器,求最佳循环比下的反应器体积.
程序如下：
function V=objfunc(beita)
global CA0 CR0 k XAf
CA0=0.99;CR0=0.01;k=1.512;XAf=0.99;v0=10;
CAf=CA0*(1-XAf);CA1=(CA0+beita*CAf)/(1+beita);
sol=quad(@func,CAf,CA1,[],[],beita);
V=(1+beita)*v0*sol;
%-----------------------------------------------------------------------
function y=func(CA,beita)
global CA0 CR0 k XAf
CRf=CR0+CA0*XAf;CR1=(CR0+beita*CRf)/(1+beita);
CA1=(CA0+beita*CAf)/(1+beita);
CR=CR1+CA1-CA;
rate=k*CA.*CR;
y=1./rate;
beita=fminsearch(@objfunc,1);
disp(beita)
但是一运行就出现：
Input argument "beita" is undefined.
Error in ==> objfunc at 4
CAf=CA0*(1-XAf);CA1=(CA0+beita*CAf)/(1+beita);
流程示意图
beita=fminsearch(@objfunc,1);
disp(beita)
最后这两句有问题
matlab带循环的反应器程序报错初学matlab,我前几天看到一个带循环操作的平推流反应器计算最优循环比下反应器最佳体积的例子,想要自己试着做一下,但是总是报错.想请教各位怎样修改这个程
beita = beita;
CA0=0.99;CR0=0.01;k=1.512;XAf=0.99;v0=10;