Code covered by the BSD License  

Highlights from
MPCC Solution of Genetic Regulatory Circuit Design Problems

image thumbnail

MPCC Solution of Genetic Regulatory Circuit Design Problems

by

James Allison

 

Demonstration of a novel MPCC-based technique for solving gene circuit design problems.

constraint(xd,p)
function C = constraint(xd,p)
[tout,Xout] = ode15s(@(t,X)sys(t,X,xd,p),[p.t0 p.tf],p.xinit);
O1_index = length(find(tout<p.t1));
t_peak = p.t1 + 5;

Opeak_index = length(find(tout<t_peak));
Opeak = Xout(Opeak_index,3);
t_peak_end = p.t1 + 5 *2;
Opeak_end_index = length(find(tout<t_peak_end));

Opeak_end = Xout(Opeak_end_index,3);

O2 = Xout(end,3);
O1 = Xout(O1_index,3);
sensitivity_num = (Opeak-O1)/O1;

sensitivity_den = (p.I2-p.I1)/p.I1;

G(1) = 1.2*O1 - Opeak ;
G(2) = O1 - Opeak + 0.01;
G(3) = Opeak_end - Opeak + 0.01;
G(4) = -O1 + 0.01;
G(5) = -Opeak + 0.01;
C = G;
end

Contact us