26、化学反应工程中的MATLAB计算

化学反应工程中的MATLAB计算

1. MATLAB脚本 useadbpfr

1.1 脚本功能与代码

useadbpfr 脚本的主要功能是设置数据，并使用内置函数 ode45 来求解由 adbpfr 函数定义的微分方程。以下是具体的代码：

% useadbpfr.m
clear all;
Ca0 = 9.3; Fa0 = 146.7; T1 = 360; T2 = 333; k1 = 31.1; K2 = 3.03;
E = 65700; R = 8.314; dH = -6900; Vspan = [0 5]; X0 = 0;
[V X] = ode45(@adbpfr,Vspan,X0,[],Ca0,Fa0,T1,T2,k1,K2,E,R,dH);
T = 330 + 43.4265*X;
k = k1*exp(E*(1/T1 - 1./T)/R); Kc = K2*exp(dH*(1/T2 - 1./T)/R);
ra = -k*Ca0.*(1 - (1 + 1./Kc).*X); Xe = Kc./(1+Kc);
figure(1)
subplot(1, 2), plot(V,T), xlabel('V'), ylabel('T(K)')
subplot(1, 2), plot(V,-ra), xlabel('V'), ylabel('-r_A')
figure(2)
plot(V,X,V,Xe,'--'), xlabel('V'), ylabel('X,X_e'), legend('X'