Module 8 - Biochemical Reactors

bio_ss.m Solves for the steady-state of the biochemical reactor by using fsolve.

Note: The fsolve function must be available. This is standard with MATLAB 4.x, but requires the optimization toolbox for MATLAB 5.x. To run this file, simply enter

  x0 = [1;2];   % initial guess for biomass and substrate concentrations
  x = fsolve('bio_ss',x0);

at the MATLAB prompt. The solution obtained for these sets of parameters is

  x = 
      0.9951
      1.5122

bio.m Solves for the dynamics of the biochemical reactor by using ode45.

To reproduce Figure M8.3, use:

   x0 = [0.75;2];
   [t,x] = ode45('bio',[0 30],x0); % MATLAB 5.x
   [t,x] = ode45('bio',0,30,x0);   % MATLAB 4.x

bio_phas.m Generates the phase-plane plot shown in Figure M8.4.

   bio_phas

Please note that you need to change the call to ode45 within bio_phas.m if you are using MATLAB 5.x.