example.m

模拟化工生产精馏开工过程,运行软件为MATLAB

%用函数ode45()来解微分方程，并画图计算结果(注：ode45()采用龙格库塔算法求解)
%example.m
%初始化x1，x2,...,x22。x1即xD,x22即xB
[t,x]=ode45(@ivpodefun,[0 500],[0.55;0.545;0.54;0.535;0.53;0.525;0.52;0.515;0.51;0.505;0.5;0.495;0.49;0.485;0.48;0.475;0.47;0.465;0.46;0.455;0.45;0.445]);
%画图
plot(t,x(:,1),'-',t,x(:,22),'--')
xlabel('t/min');
ylabel('x');
legend('xD','xB');