Team:Paris/Modeling/plot

From 2008.igem.org

Revision as of 15:45, 8 August 2008 by YannLC (Talk | contribs)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

To go back to the previous page : Clic here.


% Main function for the simulations

% Main function for the simulations

close all;

%% File to use to obtain :

% 1 - The ouput of the system

% 2 - A visualization of that output

% 3 - The value of a criteria describing whether there are oscillations

global Ncell;

global Next;

tmax=50;

t=[0:0.1:tmax];

t2=t.*20/log(2);

V0=zeros(1,Next);

for i=1:Ncell

V0(8+8*(i-1))=0;

end

%V0=[0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;1;0]; %Initial state

V=simul(t,V0);%Numerical simulation of the system

% Visualization of the output

Z1=[V(:,8)];

Z2=[V(:,2:4)];

% Z2=[V(:,9),V(:,10),V(:,11),V(:,12),V(:,13),V(:,14),V(:,15),V(:,16)];

% Z3=[V(:,17)];

%% tests oscill & FIFO

plot(t2,Z1);

legend('FlhDC');

figure;

plot(t2,Z2');

legend('FliL','FlgA','FlhB');

% title('Cell1');

%xlabel('time (min)');

% ylabel('concentrations');

%

% figure;

% plot(t,Z2);

% legend('FliA','FliL','FlgA','FlhB','AHL interieur','mTetR','TetR','flhDC');

% title('cell2');

% xlabel('time');

% ylabel('concentrations');

%

% figure;

% plot(t,Z3');

% legend('AHL ext');

% title('2 cells model');

% xlabel('time');

% ylabel('concentrations');

%% tests synchronisation

% figure;

% plot(t,[V(:,2),V(:,10)]');

% legend('fliL');

% title('Test Synchronisation');

% xlabel('time');

% ylabel('concentrations');

%

% figure;

% plot(t,[V(:,3),V(:,11)]);

% legend('flgA');

% title('Test Synchronisation');

% xlabel('time');

% ylabel('concentrations');

%

% figure;

% plot(t,[V(:,4),V(:,12)]);

% legend('flhB');

% title('Test Synchronisation');

% xlabel('time');

% ylabel('concentrations');