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

 

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