Edinburgh/SBCode
From 2007.igem.org
< Edinburgh(Difference between revisions)
Luca.Gerosa (Talk | contribs) |
Luca.Gerosa (Talk | contribs) |
||
(6 intermediate revisions not shown) | |||
Line 1: | Line 1: | ||
+ | This is the [http://www.sbtoolbox.org/ SBToolBox] code we developed for simulation of the ODEs model. | ||
+ | == The SBToolBox Code == | ||
+ | |||
+ | ********** MODEL NAME | ||
+ | Model of the Division PoPer Construct of Edinburgh iGEM 2007 Team | ||
+ | ********** MODEL NOTES | ||
+ | The model represents the Division PoPper construct with the addition of a reporter protein placed downstream | ||
+ | ********** MODEL STATES | ||
+ | d/dt(YFP) = YFPexp*ForwardPhase(time)-YFPdeg*YFP | ||
+ | d/dt(GFP) = GFPexp*BackwardPhase(time)-GFPdeg*GFP | ||
+ | YFP(0) = 0 | ||
+ | GFP(0) = 0 | ||
+ | ********** MODEL PARAMETERS | ||
+ | YFPexp = 0.0001 | ||
+ | GFPexp = 0.0001 | ||
+ | GFPdeg = 0.05 | ||
+ | YFPdeg = 0.05 | ||
+ | ********** MODEL VARIABLES | ||
+ | ********** MODEL REACTIONS | ||
+ | ********** MODEL FUNCTIONS | ||
+ | ********** MODEL EVENTS | ||
+ | |||
+ | ********** MODEL MATLAB FUNCTIONS | ||
+ | function [result] = ForwardPhase(t) | ||
+ | FlippingTime=15; | ||
+ | DivisionTime=80; | ||
+ | temp=0; | ||
+ | temp = mod(t,2*DivisionTime+2*FlippingTime); | ||
+ | if temp < (DivisionTime) | ||
+ | result=1; | ||
+ | else result=0; | ||
+ | end | ||
+ | return | ||
+ | |||
+ | function [result] = BackwardPhase(t) | ||
+ | FlippingTime=15; | ||
+ | DivisionTime=80; | ||
+ | temp=0; | ||
+ | temp = mod(t,2*DivisionTime+2*FlippingTime); | ||
+ | if temp > (DivisionTime+FlippingTime) | ||
+ | if temp < (2*DivisionTime+FlippingTime) | ||
+ | result=1; | ||
+ | else result=0; | ||
+ | end | ||
+ | else result=0; | ||
+ | end | ||
+ | return |
Latest revision as of 15:34, 24 October 2007
This is the [http://www.sbtoolbox.org/ SBToolBox] code we developed for simulation of the ODEs model.
The SBToolBox Code
********** MODEL NAME Model of the Division PoPer Construct of Edinburgh iGEM 2007 Team ********** MODEL NOTES The model represents the Division PoPper construct with the addition of a reporter protein placed downstream ********** MODEL STATES d/dt(YFP) = YFPexp*ForwardPhase(time)-YFPdeg*YFP d/dt(GFP) = GFPexp*BackwardPhase(time)-GFPdeg*GFP YFP(0) = 0 GFP(0) = 0 ********** MODEL PARAMETERS YFPexp = 0.0001 GFPexp = 0.0001 GFPdeg = 0.05 YFPdeg = 0.05 ********** MODEL VARIABLES ********** MODEL REACTIONS ********** MODEL FUNCTIONS ********** MODEL EVENTS
********** MODEL MATLAB FUNCTIONS function [result] = ForwardPhase(t) FlippingTime=15; DivisionTime=80; temp=0; temp = mod(t,2*DivisionTime+2*FlippingTime); if temp < (DivisionTime) result=1; else result=0; end return function [result] = BackwardPhase(t) FlippingTime=15; DivisionTime=80; temp=0; temp = mod(t,2*DivisionTime+2*FlippingTime); if temp > (DivisionTime+FlippingTime) if temp < (2*DivisionTime+FlippingTime) result=1; else result=0; end else result=0; end return