"Model-based design of smart active packaging systems with antimicrobial activity" Matlab codes

This repository contains the Matlab codes used in the Results section of the paper  "Model-based design of smart active packaging systems with antimicrobial activity".

The different optimizations are performed using Differential Evolution (DE) algorithm. The Matlab codes can be freely downloaded from:

http://www1.icsi.berkeley.edu/~storn/code.html

In order to use DE for solving the problems described in the manuscript, users have two alternatives:
- Add the DE path to the Matlab path
- Copy the DeMat folder in the same directory as the folders containing the case studies codes.

There are three folders:

    case1_onlySafety --> It contains the codes for the first case study where only safety constraints are considered. The main file to perform the optimization is mainOptDesign.m. In this file the user may select the number of layers to be considered (line 20) and the configuration of the active packaging (lines 21-38). Other parameters (such as bounds on the decision variables or constraints) can be chosen/modified in the rest of the code. However, doing so would imply that the solution may differ from the one presented in the manuscript. Folder "results" contain the optimal design results presented in the manuscript.

    case2_safetyQuality --> It contains the codes for the second case study where both safety and quality constraints are considered. The main file to perform the optimization is mainOptDesign.m. In this file the user may select the number of layers to be considered (line 20) and the configuration of the active packaging (lines 21-38). Other parameters (such as bounds on the decision variables or constraints) can be chosen/modified in the rest of the code. However, doing so would imply that the solution may differ from the one presented in the manuscript. Folder "results" contain the optimal design results presented in the manuscript.

    case3_effectVariableTProfile --> It contains the codes for the third case study where the effect of perturbations of the storage temperature are studied. Two type of perturbations are considered: (i) a set of 50 randomly generated temperature-time profiles (mainRandomTprofiles.m). Function "createRandomTprofiles.m" can be used to generate random profiles in the storage temperature, however the ones used in the manuscript are already stored in "Tprofiles.mat". (ii) a set of 10 profiles with the same perturbation on the storage temperature introduced at different time instants (mainTprofilesDifferentTimes.m).