MATLAB在化工领域的应用 举例

完整的MATLAB程序如下所列：

function decayinCSTR



global ct0 yA0 tao k kd

ct0=1.0；

yA0＝0.8；

tao＝0.02；

k＝45；

kd＝9；

c0＝［0.8 1］；

t0＝0；t1＝0.5；

[t,c]＝ode45（@decay，［t0 t1］,c0）；

nt＝length(c(:,1));

for i=1:nt

x(i)=1-(1+yA0)/(1+c(i,1)/ct0)*c(i,1)/ct0/yA0;

end

plot(t,c(:,1),’b-‘,t,c(:2),’b--‘,t,x,’b*’)

xlabel(’time(h)’)

ylabel(‘cA,a,x’)

legend(‘cA’,’a’.’x’)



function dcdt=decay(t,c)



global ct0 yA0 tao k kd

dcdt(1)=ct0*yA0/tao-((1+yA0)/(1+c(1)/ct0/yA0)+c(2)*tao*k)/tao*c(1);

dcdt(2)=-kd*c(1)*c(2);

dcdt=[dcdt(1);dcdt(2)];







运行上述程序得反应器出口原料油浓度、转化率和反应器中催化剂活性随

时间的变化如下图所示。

.......想学学minitab.......没看懂.......

现在已经看不懂了，以前还学过的