Edinburgh/SBCode
From 2007.igem.org
(Difference between revisions)
Luca.Gerosa (Talk | contribs) |
Luca.Gerosa (Talk | contribs) |
||
Line 1: | Line 1: | ||
********** MODEL NAME | ********** MODEL NAME | ||
Model of the Division PoPer Construct of Edinburgh iGEM 2007 Team | Model of the Division PoPer Construct of Edinburgh iGEM 2007 Team | ||
- | |||
********** MODEL NOTES | ********** MODEL NOTES | ||
- | |||
The model represents the Division PoPper construct with the addition of a reporter protein placed downstream | 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=10; | |
- | + | DivisionTime=40; | |
- | + | temp=0; | |
- | + | temp = mod(t,2*DivisionTime+2*FlippingTime); | |
- | + | if temp < (DivisionTime) | |
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | ********** MODEL MATLAB FUNCTIONS | + | |
- | function [result] = ForwardPhase(t) | + | |
- | FlippingTime=10; | + | |
- | DivisionTime=40; | + | |
- | temp=0; | + | |
- | temp = mod(t,2*DivisionTime+2*FlippingTime); | + | |
- | if temp < (DivisionTime) | + | |
result=1; | result=1; | ||
- | else result=0; | + | else result=0; |
- | end | + | end |
- | return | + | return |
- | function [result] = BackwardPhase(t) | + | function [result] = BackwardPhase(t) |
- | FlippingTime=10; | + | FlippingTime=10; |
- | DivisionTime=40; | + | DivisionTime=40; |
- | temp=0; | + | temp=0; |
- | temp = mod(t,2*DivisionTime+2*FlippingTime); | + | temp = mod(t,2*DivisionTime+2*FlippingTime); |
- | if temp > (DivisionTime+FlippingTime) | + | if temp > (DivisionTime+FlippingTime) |
if temp < (2*DivisionTime+FlippingTime) | if temp < (2*DivisionTime+FlippingTime) | ||
result=1; | result=1; | ||
Line 55: | Line 41: | ||
end | end | ||
else result=0; | else result=0; | ||
- | end | + | end |
- | return | + | return |
Revision as of 14:18, 24 October 2007
********** 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=10; DivisionTime=40; temp=0; temp = mod(t,2*DivisionTime+2*FlippingTime); if temp < (DivisionTime) result=1; else result=0; end return function [result] = BackwardPhase(t) FlippingTime=10; DivisionTime=40; 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