Edinburgh/SBCode

From 2007.igem.org

(Difference between revisions)
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 STATES
+
********** MODEL MATLAB FUNCTIONS
-
d/dt(YFP) = YFPexp*ForwardPhase(time)-YFPdeg*YFP
+
function [result] = ForwardPhase(t)                                                                                                           
-
d/dt(GFP) = GFPexp*BackwardPhase(time)-GFPdeg*GFP
+
FlippingTime=10;                                                                                                                             
-
d/dt(YFP) =
+
DivisionTime=40;                                                                                                                             
-
 
+
temp=0;                                                                                                                                       
-
YFP(0) = 0
+
temp = mod(t,2*DivisionTime+2*FlippingTime);                                                                                                 
-
GFP(0) = 0
+
if temp < (DivisionTime)                                                                                                                     
-
 
+
-
********** 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;                                                                                                                                   
   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