<strong>DOTcvpSB</strong>: a <strong>Matlab</strong> <strong>Toolbox</strong> <strong>for</strong> <strong>Dynamic</strong> <strong>Optimization</strong> <strong>in</strong> <strong>Systems</strong> Biology11.5.6 MITSDOTcvp: dotcvp_mits_default.m1 function [data] = dotcvp_mits_default(data)23 data.options.mits.NLPtol = data.nlp.NLPtol; % NLP tolerance level4 data.options.mits.max_iter = data.nlp.MaxIter; % maximum number of iterations5 data.options.mits.MaxTime = data.nlp.MaxCPUTime; % maximum CPU time [sec]6 data.options.mits.fex = -<strong>in</strong>f; % stopp<strong>in</strong>g criteria: if f = fex --> stop7 data.options.mits.locstart = 0; % if true: MITS starts with a local solver run from the <strong>in</strong>itial guess89 end11.5.7 SRESDOTcvp: dotcvp_sres_default.m1 function [data] = dotcvp_sres_default(data)23 data.options.sres.mm = 'm<strong>in</strong>'; % ['max'|'m<strong>in</strong>'] (<strong>for</strong> maximization or m<strong>in</strong>imization)4 data.options.sres.lambda = 150; % population size (number of offspr<strong>in</strong>g) (100 to 200)5 data.options.sres.G = 100000000; % maximum number of generations6 data.options.sres.mu = round(data.options.sres.lambda/7); % parent number (mu/lambda usually 1/7), be<strong>for</strong>e: roundn(lambda/7,0)7 data.options.sres.pf = 0.45; % pressure on fitness <strong>in</strong> [0 0.5] try around 0.458 data.options.sres.varphi = 1; % expected rate of convergence (usually 1)9 data.options.sres.MaxTime = data.nlp.MaxCPUTime; % maximum CPU time [sec]10 data.options.sres.itermax = data.nlp.MaxIter; % maximum number of iterations11 data.options.sres.MaxFunEvals = <strong>in</strong>f; % maximum number of function evaluations12 data.options.sres.NLPtol = data.nlp.NLPtol; % NLP tolerance level13 data.options.sres.delta = 0.75; % allowed the time variation <strong>for</strong> time grid adaptation [0 < 0.75 < 1]1415 endPage – 53
Bibliography[1] E. Balsa-Canto, J. R. Banga, A. A. Alonso, and V. S. Vassiliadis. Efficient optimal control of bioprocessesus<strong>in</strong>g second-order <strong>in</strong><strong>for</strong>mation. Industrial and Eng<strong>in</strong>eer<strong>in</strong>g Chemistry Research, 39(11):4287–4295, 2000.[2] E. Balsa-Canto, J. R. Banga, A. A. Alonso, and V. S. Vassiliadis. Restricted second order <strong>in</strong><strong>for</strong>mation <strong>for</strong>the solution of optimal control problems us<strong>in</strong>g control vector parameterization. Journal of Process Control,12(2):243–255, 2 2002.[3] E. Balsa-Canto, V. S. Vassiliadis, and J. R. Banga. <strong>Dynamic</strong> optimization of s<strong>in</strong>gle- and multi-stage systemsus<strong>in</strong>g a hybrid stochastic-determ<strong>in</strong>istic method. Industrial and Eng<strong>in</strong>eer<strong>in</strong>g Chemistry Research, 44(5):1514–1523, 2005.[4] J. R. Banga, A. A. Alonso, and R. P. S<strong>in</strong>gh. Stochastic dynamic optimization of batch and semicont<strong>in</strong>uousbioprocesses. Biotechnology Progress, 13(3):326–335, 1997.[5] J. R. Banga, E. Balsa-Canto, E. G. Moles, and A. A. Alonso. <strong>Dynamic</strong> optimization of bioprocesses: Efficientand robust numerical strategies. Journal of Biotechnology, 117:407–419, 2005.[6] J. R. Banga, R. Irizarry-Rivera, and W. D. Seider. Stochastic optimization <strong>for</strong> optimal and model-predictivecontrol. Computers and Chemical Eng<strong>in</strong>eer<strong>in</strong>g, 22:603–612(10), 1998.[7] L. T. Biegler and I. E. Grossmann. Retrospective on optimization. Computers and Chemical Eng<strong>in</strong>eer<strong>in</strong>g,28(8):1169–1192, 7 2004.[8] A. E. Bryson and Y. Ch. Ho. Applied Optimal Control - <strong>Optimization</strong>, Estimation and Control. HemispherePublish<strong>in</strong>g corporation, Wash<strong>in</strong>gton, 1975.[9] C.-T. Chen and C. Hwang. Optimal control computation <strong>for</strong> differential-algebraic process systems withgeneral constra<strong>in</strong>ts. Chemical Eng<strong>in</strong>eer<strong>in</strong>g Communications, 97(1):9–26, 1 1990.[10] T. Coleman, M. A. Branch, and A. Grace. <strong>Optimization</strong> toolbox <strong>for</strong> use with matlab user’s guide version 2,1998.[11] S. Crescitelli and S. Nicoletti. Near optimal control of batch reactors. Chemical Eng<strong>in</strong>eer<strong>in</strong>g Science, 28:463–471, 1973.[12] S. A. Dadebo and K. B. McAuley. <strong>Dynamic</strong> optimization of constra<strong>in</strong>ed chemical eng<strong>in</strong>eer<strong>in</strong>g problemsus<strong>in</strong>g dynamic programm<strong>in</strong>g. Computers and Chemical Eng<strong>in</strong>eer<strong>in</strong>g, 19(5):513–525, 1995. 10: Diversen.[13] O. Exler, L. T. Antelo, J. A. Egea, A. A. Alonso, and J. R. Banga. A tabu search-based algorithm <strong>for</strong> mixed<strong>in</strong>tegernonl<strong>in</strong>ear problems and its application to <strong>in</strong>tegrated process and control system design. Computersand Chemical Eng<strong>in</strong>eer<strong>in</strong>g, 32(8):1877–1891, 9 2008.[14] O. Exler and K. Schittkowski. A trust region sqp algorithm <strong>for</strong> mixed-<strong>in</strong>teger nonl<strong>in</strong>ear programm<strong>in</strong>g. <strong>Optimization</strong>Letters, 1(3):269–280, 6 2007.[15] W. F. Feehery and P. I. Barton. <strong>Dynamic</strong> optimization with state variable path constra<strong>in</strong>ts. Computers andChemical Eng<strong>in</strong>eer<strong>in</strong>g, 22:1241–1256, 1998.[16] M. Fikar and M. A. Latifi. User’s guide <strong>for</strong> FORTRAN dynamic optimisation code DYNO. Technical Reportmf0201, LSGC CNRS, Nancy, France; SUT Bratislava, Slovak Republic, 2001.[17] A. C. H<strong>in</strong>dmarsh, P. N. Brown, K. E. Grant, S. L. Lee, R. Serban, D. E. Shumaker, and C. S. Woodward.Sundials: Suite of nonl<strong>in</strong>ear and differential/algebraic equation solvers,. ACM Transactions on MathematicalSoftware, 31(3):363–396, 2005.