% 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');
"
