Team:Paris/Modeling/plot

From 2008.igem.org

(Difference between revisions)
(New page: {{Paris/Menu}} {|cellspacing="5" cellpadding="10" style="background:#649CD7; width: 965px;" |-valign="top" |style="background:#ffffff"| |}<br style="clear:both" />)
Line 3: Line 3:
|-valign="top"
|-valign="top"
|style="background:#ffffff"|
|style="background:#ffffff"|
 +
 +
<html xmlns:v="urn:schemas-microsoft-com:vml"
 +
xmlns:o="urn:schemas-microsoft-com:office:office"
 +
xmlns:w="urn:schemas-microsoft-com:office:word"
 +
xmlns="http://www.w3.org/TR/REC-html40">
 +
 +
<head>
 +
<meta http-equiv=Content-Type content="text/html; charset=windows-1252">
 +
<meta name=ProgId content=Word.Document>
 +
<meta name=Generator content="Microsoft Word 11">
 +
<meta name=Originator content="Microsoft Word 11">
 +
<link rel=File-List href="premiermodele_fichiers/filelist.xml">
 +
<title>% Main function for the simulations</title>
 +
<!--[if gte mso 9]><xml>
 +
<o:DocumentProperties>
 +
  <o:Author> </o:Author>
 +
  <o:LastAuthor> </o:LastAuthor>
 +
  <o:Revision>1</o:Revision>
 +
  <o:TotalTime>1</o:TotalTime>
 +
  <o:Created>2008-08-06T17:08:00Z</o:Created>
 +
  <o:LastSaved>2008-08-06T17:09:00Z</o:LastSaved>
 +
  <o:Pages>1</o:Pages>
 +
  <o:Words>94</o:Words>
 +
  <o:Characters>522</o:Characters>
 +
  <o:Company> Supélec</o:Company>
 +
  <o:Lines>4</o:Lines>
 +
  <o:Paragraphs>1</o:Paragraphs>
 +
  <o:CharactersWithSpaces>615</o:CharactersWithSpaces>
 +
  <o:Version>11.6408</o:Version>
 +
</o:DocumentProperties>
 +
</xml><![endif]--><!--[if gte mso 9]><xml>
 +
<w:WordDocument>
 +
  <w:GrammarState>Clean</w:GrammarState>
 +
  <w:HyphenationZone>21</w:HyphenationZone>
 +
  <w:PunctuationKerning/>
 +
  <w:ValidateAgainstSchemas/>
 +
  <w:SaveIfXMLInvalid>false</w:SaveIfXMLInvalid>
 +
  <w:IgnoreMixedContent>false</w:IgnoreMixedContent>
 +
  <w:AlwaysShowPlaceholderText>false</w:AlwaysShowPlaceholderText>
 +
  <w:Compatibility>
 +
  <w:BreakWrappedTables/>
 +
  <w:SnapToGridInCell/>
 +
  <w:WrapTextWithPunct/>
 +
  <w:UseAsianBreakRules/>
 +
  <w:DontGrowAutofit/>
 +
  </w:Compatibility>
 +
  <w:BrowserLevel>MicrosoftInternetExplorer4</w:BrowserLevel>
 +
</w:WordDocument>
 +
</xml><![endif]--><!--[if gte mso 9]><xml>
 +
<w:LatentStyles DefLockedState="false" LatentStyleCount="156">
 +
</w:LatentStyles>
 +
</xml><![endif]-->
 +
<style>
 +
<!--
 +
/* Style Definitions */
 +
p.MsoNormal, li.MsoNormal, div.MsoNormal
 +
{mso-style-parent:"";
 +
margin:0cm;
 +
margin-bottom:.0001pt;
 +
mso-pagination:widow-orphan;
 +
font-size:12.0pt;
 +
font-family:"Times New Roman";
 +
mso-fareast-font-family:"Times New Roman";}
 +
span.GramE
 +
{mso-style-name:"";
 +
mso-gram-e:yes;}
 +
@page Section1
 +
{size:612.0pt 792.0pt;
 +
margin:70.85pt 70.85pt 70.85pt 70.85pt;
 +
mso-header-margin:36.0pt;
 +
mso-footer-margin:36.0pt;
 +
mso-paper-source:0;}
 +
div.Section1
 +
{page:Section1;}
 +
-->
 +
</style>
 +
<!--[if gte mso 10]>
 +
<style>
 +
/* Style Definitions */
 +
table.MsoNormalTable
 +
{mso-style-name:"Tableau Normal";
 +
mso-tstyle-rowband-size:0;
 +
mso-tstyle-colband-size:0;
 +
mso-style-noshow:yes;
 +
mso-style-parent:"";
 +
mso-padding-alt:0cm 5.4pt 0cm 5.4pt;
 +
mso-para-margin:0cm;
 +
mso-para-margin-bottom:.0001pt;
 +
mso-pagination:widow-orphan;
 +
font-size:10.0pt;
 +
font-family:"Times New Roman";
 +
mso-ansi-language:#0400;
 +
mso-fareast-language:#0400;
 +
mso-bidi-language:#0400;}
 +
</style>
 +
<![endif]--><!--[if gte mso 9]><xml>
 +
<o:shapedefaults v:ext="edit" spidmax="2050"/>
 +
</xml><![endif]--><!--[if gte mso 9]><xml>
 +
<o:shapelayout v:ext="edit">
 +
  <o:idmap v:ext="edit" data="1"/>
 +
</o:shapelayout></xml><![endif]-->
 +
</head>
 +
 +
<body lang=FR style='tab-interval:35.4pt'>
 +
 +
<div class=Section1>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
 +
mso-ansi-language:EN-GB'>% Main function for the simulations</span><span
 +
lang=EN-GB style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
 +
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'>clear, close </span><span lang=EN-GB style='font-size:
 +
10.0pt;font-family:"Courier New";color:#A020F0;mso-ansi-language:EN-GB'>all</span><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'>;</span><span lang=EN-GB style='font-family:"Courier New";
 +
mso-bidi-font-family:"Times New Roman";mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'><span style='mso-spacerun:yes'> </span></span><span
 +
lang=EN-GB style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
 +
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
 +
mso-ansi-language:EN-GB'>%% File to use to obtain :</span><span lang=EN-GB
 +
style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
 +
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
 +
mso-ansi-language:EN-GB'>%<span style='mso-spacerun:yes'>  </span>1 - The
 +
ouput of the system</span><span lang=EN-GB style='font-family:"Courier New";
 +
mso-bidi-font-family:"Times New Roman";mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
 +
mso-ansi-language:EN-GB'>%<span style='mso-spacerun:yes'>  </span>2 - A
 +
visualization of that output</span><span lang=EN-GB style='font-family:"Courier New";
 +
mso-bidi-font-family:"Times New Roman";mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
 +
mso-ansi-language:EN-GB'>%<span style='mso-spacerun:yes'>  </span>3 - The
 +
value of a criteria describing whether there are oscillations</span><span
 +
lang=EN-GB style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
 +
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
 +
mso-ansi-language:EN-GB'><span style='mso-spacerun:yes'> </span></span><span
 +
lang=EN-GB style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
 +
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'>tmax=200;</span><span lang=EN-GB style='font-family:
 +
"Courier New";mso-bidi-font-family:"Times New Roman";mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'>t=[0:0.1:tmax];</span><span lang=EN-GB
 +
style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
 +
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
 +
mso-ansi-language:EN-GB'>%Prop=0.3;</span><span lang=EN-GB style='font-family:
 +
"Courier New";mso-bidi-font-family:"Times New Roman";mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
 +
mso-ansi-language:EN-GB'><span style='mso-spacerun:yes'> </span></span><span
 +
lang=EN-GB style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
 +
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'>V0=[1;0;0;0;0;0;0;0;0];</span><span lang=EN-GB
 +
style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;mso-ansi-language:
 +
EN-GB'>%Initial state</span><span lang=EN-GB style='font-family:"Courier New";
 +
mso-bidi-font-family:"Times New Roman";mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'>V=simul(t,V0);</span><span lang=EN-GB
 +
style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;mso-ansi-language:
 +
EN-GB'>%Numerical simulation of the system</span><span lang=EN-GB
 +
style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
 +
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
 +
mso-ansi-language:EN-GB'><span style='mso-spacerun:yes'> </span></span><span
 +
lang=EN-GB style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
 +
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
 +
mso-ansi-language:EN-GB'>% Visualization of the output</span><span lang=EN-GB
 +
style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
 +
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'>Z=[V(:,1),V(:,5),V(:,6),V(:,9)];</span><span
 +
lang=EN-GB style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
 +
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'>figure;</span><span lang=EN-GB style='font-family:
 +
"Courier New";mso-bidi-font-family:"Times New Roman";mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'>plot(t,V');</span><span lang=EN-GB style='font-family:
 +
"Courier New";mso-bidi-font-family:"Times New Roman";mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'>legend(</span><span lang=EN-GB style='font-size:10.0pt;
 +
font-family:"Courier New";color:#A020F0;mso-ansi-language:EN-GB'>'FlhDC'</span><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'>,</span><span lang=EN-GB style='font-size:10.0pt;
 +
font-family:"Courier New";color:#A020F0;mso-ansi-language:EN-GB'>'FliA'</span><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'>,</span><span lang=EN-GB style='font-size:10.0pt;
 +
font-family:"Courier New";color:#A020F0;mso-ansi-language:EN-GB'>'FliL'</span><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'>,</span><span lang=EN-GB style='font-size:10.0pt;
 +
font-family:"Courier New";color:#A020F0;mso-ansi-language:EN-GB'>'FlgA'</span><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'>,</span><span lang=EN-GB style='font-size:10.0pt;
 +
font-family:"Courier New";color:#A020F0;mso-ansi-language:EN-GB'>'FlhB'</span><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'>,</span><span lang=EN-GB style='font-size:10.0pt;
 +
font-family:"Courier New";color:#A020F0;mso-ansi-language:EN-GB'>'AHL
 +
interieur'</span><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
 +
color:black;mso-ansi-language:EN-GB'>,</span><span lang=EN-GB style='font-size:
 +
10.0pt;font-family:"Courier New";color:#A020F0;mso-ansi-language:EN-GB'>'AHL
 +
ext'</span><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
 +
color:black;mso-ansi-language:EN-GB'>,</span><span lang=EN-GB style='font-size:
 +
10.0pt;font-family:"Courier New";color:#A020F0;mso-ansi-language:EN-GB'>'mTetR'</span><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'>,</span><span lang=EN-GB style='font-size:10.0pt;
 +
font-family:"Courier New";color:#A020F0;mso-ansi-language:EN-GB'>'TetR'</span><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
 +
mso-ansi-language:EN-GB'>);</span><span lang=EN-GB style='font-family:"Courier New";
 +
mso-bidi-font-family:"Times New Roman";mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
 +
mso-ansi-language:EN-GB'>%legend('FlhDC','FlhB','AHL int','TetR');</span><span
 +
lang=EN-GB style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
 +
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
style='font-size:10.0pt;font-family:"Courier New";color:forestgreen'>% <span
 +
class=GramE>xlabel(</span>'Time');</span><span style='font-family:"Courier New";
 +
mso-bidi-font-family:"Times New Roman"'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
 +
style='font-size:10.0pt;font-family:"Courier New";color:forestgreen'>% <span
 +
class=GramE>ylabel(</span>'Concentrations')</span><span style='font-family:
 +
"Courier New";mso-bidi-font-family:"Times New Roman"'><o:p></o:p></span></p>
 +
 +
<p class=MsoNormal><o:p>&nbsp;</o:p></p>
 +
 +
</div>
 +
 +
</body>
 +
 +
</html>
|}<br style="clear:both" />
|}<br style="clear:both" />

Revision as of 17:27, 6 August 2008

% Main function for the simulations

% Main function for the simulations

clear, 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

tmax=200;

t=[0:0.1:tmax];

%Prop=0.3;

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

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

% Visualization of the output

Z=[V(:,1),V(:,5),V(:,6),V(:,9)];

figure;

plot(t,V');

legend('FlhDC','FliA','FliL','FlgA','FlhB','AHL interieur','AHL ext','mTetR','TetR');

%legend('FlhDC','FlhB','AHL int','TetR');

% xlabel('Time');

% ylabel('Concentrations')