Team:DTU-Denmark/Matlab

From 2011.igem.org

(Difference between revisions)
(Simulation)
(Simulation)
Line 7: Line 7:
[[File:DTU_Model8.txt]]
[[File:DTU_Model8.txt]]
 +
 +
<pre>
 +
********** MODEL NAME
 +
Dimensionless form. Catalytical.
 +
 +
********** MODEL NOTES
 +
Kinetic model of trap-RNA system.
 +
Parameters are estimated from literature.
 +
 +
********** MODEL STATES
 +
d/dt(m) = 1 - m - k_s*alpha_m*m*s/(beta_m*beta_s)
 +
d/dt(s) = (beta_s/beta_m)*(alpha_s/alpha_m - s - k_r*alpha_m * s * r /(beta_s*beta_r))
 +
d/dt(r) = (beta_r/beta_m)*(alpha_r/alpha_m - r)
 +
 +
m(0) = 1
 +
s(0) = 0
 +
r(0) = 0                                   
 +
 +
********** MODEL PARAMETERS
 +
alpha_m = 10
 +
alpha_s = 0
 +
alpha_r = 0
 +
beta_m = 0.0257
 +
beta_s = 0.0257
 +
beta_r = 0.0257
 +
k_s = 0.00082
 +
k_r = 0.0082
 +
 +
 +
********** MODEL VARIABLES
 +
 +
********** MODEL REACTIONS
 +
 +
 +
********** MODEL FUNCTIONS
 +
 +
 +
********** MODEL EVENTS
 +
event = gt(time,1), alpha_s, 40
 +
event = gt(time,3), alpha_r, 200
 +
event = gt(time,6), alpha_r, 0
 +
********** MODEL MATLAB FUNCTIONS
 +
</pre>
The script running simulation and generating figures.
The script running simulation and generating figures.
{{:Team:DTU-Denmark/Templates/Standard_page_end}}
{{:Team:DTU-Denmark/Templates/Standard_page_end}}

Revision as of 02:25, 22 September 2011

Matlab code

Simulation

Temporal simulation is performed using the Systems Biology Toolbox 2 http://www.sbtoolbox.org/ environment with numerical integration using ode45. The used models.

File:DTU Model7.txt

File:DTU Model8.txt

********** MODEL NAME
Dimensionless form. Catalytical.

********** MODEL NOTES
Kinetic model of trap-RNA system.
Parameters are estimated from literature.

********** MODEL STATES
d/dt(m) = 1 - m - k_s*alpha_m*m*s/(beta_m*beta_s)
d/dt(s) = (beta_s/beta_m)*(alpha_s/alpha_m - s - k_r*alpha_m * s * r /(beta_s*beta_r))
d/dt(r) = (beta_r/beta_m)*(alpha_r/alpha_m - r)

m(0) = 1
s(0) = 0
r(0) = 0                                    

********** MODEL PARAMETERS
alpha_m = 10
alpha_s = 0
alpha_r = 0
beta_m = 0.0257
beta_s = 0.0257
beta_r = 0.0257
k_s = 0.00082
k_r = 0.0082


********** MODEL VARIABLES

********** MODEL REACTIONS
	

********** MODEL FUNCTIONS


********** MODEL EVENTS
event = gt(time,1), alpha_s, 40
event = gt(time,3), alpha_r, 200
event = gt(time,6), alpha_r, 0
********** MODEL MATLAB FUNCTIONS

The script running simulation and generating figures.