Team:Paris/Modeling/simul

From 2008.igem.org

(Difference between revisions)
 
(2 intermediate revisions not shown)
Line 5: Line 5:
To go back to the previous page : [[Team:Paris/Modeling#Corresponding codes|Clic here]].
To go back to the previous page : [[Team:Paris/Modeling#Corresponding codes|Clic here]].
 +
<html xmlns:v="urn:schemas-microsoft-com:vml"
<html xmlns:v="urn:schemas-microsoft-com:vml"
Line 16: Line 17:
<meta name=Generator content="Microsoft Word 11">
<meta name=Generator content="Microsoft Word 11">
<meta name=Originator content="Microsoft Word 11">
<meta name=Originator content="Microsoft Word 11">
-
<link rel=File-List href="simul_fichiers/filelist.xml">
+
<link rel=File-List href="simul2_fichiers/filelist.xml">
<title>% simul</title>
<title>% simul</title>
<!--[if gte mso 9]><xml>
<!--[if gte mso 9]><xml>
Line 23: Line 24:
   <o:LastAuthor> </o:LastAuthor>
   <o:LastAuthor> </o:LastAuthor>
   <o:Revision>1</o:Revision>
   <o:Revision>1</o:Revision>
-
   <o:TotalTime>0</o:TotalTime>
+
   <o:TotalTime>1</o:TotalTime>
-
   <o:Created>2008-08-06T17:08:00Z</o:Created>
+
   <o:Created>2008-08-08T15:37:00Z</o:Created>
-
   <o:LastSaved>2008-08-06T17:08:00Z</o:LastSaved>
+
   <o:LastSaved>2008-08-08T15:38:00Z</o:LastSaved>
   <o:Pages>1</o:Pages>
   <o:Pages>1</o:Pages>
-
   <o:Words>74</o:Words>
+
   <o:Words>76</o:Words>
-
   <o:Characters>409</o:Characters>
+
   <o:Characters>418</o:Characters>
   <o:Company> Supélec</o:Company>
   <o:Company> Supélec</o:Company>
   <o:Lines>3</o:Lines>
   <o:Lines>3</o:Lines>
   <o:Paragraphs>1</o:Paragraphs>
   <o:Paragraphs>1</o:Paragraphs>
-
   <o:CharactersWithSpaces>482</o:CharactersWithSpaces>
+
   <o:CharactersWithSpaces>493</o:CharactersWithSpaces>
   <o:Version>11.6408</o:Version>
   <o:Version>11.6408</o:Version>
  </o:DocumentProperties>
  </o:DocumentProperties>
</xml><![endif]--><!--[if gte mso 9]><xml>
</xml><![endif]--><!--[if gte mso 9]><xml>
  <w:WordDocument>
  <w:WordDocument>
-
  <w:GrammarState>Clean</w:GrammarState>
 
   <w:HyphenationZone>21</w:HyphenationZone>
   <w:HyphenationZone>21</w:HyphenationZone>
   <w:PunctuationKerning/>
   <w:PunctuationKerning/>
Line 68: Line 68:
font-family:"Times New Roman";
font-family:"Times New Roman";
mso-fareast-font-family:"Times New Roman";}
mso-fareast-font-family:"Times New Roman";}
-
span.GramE
 
-
{mso-style-name:"";
 
-
mso-gram-e:yes;}
 
@page Section1
@page Section1
{size:612.0pt 792.0pt;
{size:612.0pt 792.0pt;
Line 114: Line 111:
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
<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;
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
-
mso-ansi-language:EN-GB'>% <span class=GramE>simul.m :</span> simulation of the
+
mso-ansi-language:EN-GB'>% simul.m : simulation of the model</span><span
-
model</span><span lang=EN-GB style='font-family:"Courier New";mso-bidi-font-family:
+
lang=EN-GB style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
-
"Times New Roman";mso-ansi-language:EN-GB'><o:p></o:p></span></p>
+
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
Line 125: Line 122:
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
-
class=GramE><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
+
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
-
color:forestgreen;mso-ansi-language:EN-GB'>%%% to be launched.</span></span><span
+
mso-ansi-language:EN-GB'>%%% to be launched.</span><span lang=EN-GB
-
lang=EN-GB style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
+
style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
Line 148: Line 145:
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
-
class=GramE><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
+
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:blue;
-
color:blue;mso-ansi-language:EN-GB'>function</span></span><span lang=EN-GB
+
mso-ansi-language:EN-GB'>function</span><span lang=EN-GB style='font-size:10.0pt;
-
style='font-size:10.0pt;font-family:"Courier New";color:black;mso-ansi-language:
+
font-family:"Courier New";color:black;mso-ansi-language:EN-GB'>
-
EN-GB'> [y,t]=simul(t,xinit)</span><span lang=EN-GB style='font-family:"Courier New";
+
[y,t]=simul(t,xinit)</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>
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
<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
style='font-size:10.0pt;font-family:"Courier New";color:forestgreen'>%<span
-
style='mso-spacerun:yes'>  </span>[y<span class=GramE>,t</span>]=simul(t,xinit)</span><span
+
style='mso-spacerun:yes'>  </span>[y,t]=simul(t,xinit)</span><span
style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman"'><o:p></o:p></span></p>
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
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
-
class=GramE><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
+
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
-
color:forestgreen;mso-ansi-language:EN-GB'>%<span style='mso-spacerun:yes'>  
+
mso-ansi-language:EN-GB'>%<span style='mso-spacerun:yes'> </span>t : duration
-
</span>t</span></span><span lang=EN-GB style='font-size:10.0pt;font-family:
+
of the simulation</span><span lang=EN-GB style='font-family:"Courier New";
-
"Courier New";color:forestgreen;mso-ansi-language:EN-GB'> : duration of the
+
mso-bidi-font-family:"Times New Roman";mso-ansi-language:EN-GB'><o:p></o:p></span></p>
-
simulation</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
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
-
class=GramE><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
+
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
-
color:forestgreen;mso-ansi-language:EN-GB'>%<span style='mso-spacerun:yes'>  
+
mso-ansi-language:EN-GB'>%<span style='mso-spacerun:yes'> </span>xinit:
-
</span>xinit</span></span><span lang=EN-GB style='font-size:10.0pt;font-family:
+
initial state</span><span lang=EN-GB style='font-family:"Courier New";
-
"Courier New";color:forestgreen;mso-ansi-language:EN-GB'>: initial state</span><span
+
mso-bidi-font-family:"Times New Roman";mso-ansi-language:EN-GB'><o:p></o:p></span></p>
-
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
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
-
class=GramE><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
+
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
-
color:forestgreen;mso-ansi-language:EN-GB'>%<span style='mso-spacerun:yes'>  
+
mso-ansi-language:EN-GB'>%<span style='mso-spacerun:yes'> </span>y: output</span><span
-
</span>y</span></span><span lang=EN-GB style='font-size:10.0pt;font-family:
+
-
"Courier New";color:forestgreen;mso-ansi-language:EN-GB'>: output</span><span
+
lang=EN-GB style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
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>
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
Line 186: Line 177:
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
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
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:blue;
 +
mso-ansi-language:EN-GB'>global</span><span lang=EN-GB style='font-size:10.0pt;
 +
font-family:"Courier New";color:black;mso-ansi-language:EN-GB'> Ncell</span><span
lang=EN-GB style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
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>
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
Line 201: Line 199:
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
-
class=GramE><span lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";
+
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
-
color:black;mso-ansi-language:EN-GB'>options</span></span><span lang=EN-GB
+
mso-ansi-language:EN-GB'><span style='mso-spacerun:yes'> </span></span><span
-
style='font-size:10.0pt;font-family:"Courier New";color:black;mso-ansi-language:
+
lang=EN-GB style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
-
EN-GB'> = odeset(</span><span lang=EN-GB style='font-size:10.0pt;font-family:
+
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
-
"Courier New";color:#A020F0;mso-ansi-language:EN-GB'>'RelTol'</span><span
+
 
 +
<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;
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
-
mso-ansi-language:EN-GB'>,1e-4,</span><span lang=EN-GB style='font-size:10.0pt;
+
mso-ansi-language:EN-GB'>vectol=zeros(1,8*Ncell+1)+1e-4;</span><span
-
font-family:"Courier New";color:#A020F0;mso-ansi-language:EN-GB'>'AbsTol'</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;
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
-
mso-ansi-language:EN-GB'>,[1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4 1e-4]);</span><span
+
mso-ansi-language:EN-GB'>options = odeset(</span><span lang=EN-GB
 +
style='font-size:10.0pt;font-family:"Courier New";color:#A020F0;mso-ansi-language:
 +
EN-GB'>'RelTol'</span><span lang=EN-GB style='font-size:10.0pt;font-family:
 +
"Courier New";color:black;mso-ansi-language:EN-GB'>,1e-4,</span><span
 +
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:#A020F0;
 +
mso-ansi-language:EN-GB'>'AbsTol'</span><span lang=EN-GB style='font-size:10.0pt;
 +
font-family:"Courier New";color:black;mso-ansi-language:EN-GB'>,vectol);</span><span
lang=EN-GB style='font-family:"Courier New";mso-bidi-font-family:"Times New Roman";
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>
mso-ansi-language:EN-GB'><o:p></o:p></span></p>
Line 216: Line 224:
<p class=MsoNormal style='mso-layout-grid-align:none;text-autospace:none'><span
<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;
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:black;
-
mso-ansi-language:EN-GB'>[<span class=GramE>t,</span>y] = <span class=GramE>ode23(</span>@modele,t,xinit,options);</span><span
+
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'><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'>[t,y] = ode23(@modele,t,xinit,options);</span><span
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
lang=EN-GB style='font-size:10.0pt;font-family:"Courier New";color:forestgreen;
mso-ansi-language:EN-GB'>% This function computes the output of the ode
mso-ansi-language:EN-GB'>% This function computes the output of the ode
Line 229: Line 249:
</html>
</html>
 +
 +
|}<br style="clear:both" />
|}<br style="clear:both" />

Latest revision as of 15:46, 8 August 2008

To go back to the previous page : Clic here.


% simul

% simul.m : simulation of the model

%%% Commentary: This function is used in the premiermodele.m file that has

%%% to be launched.

%%% Code:

function [y,t]=simul(t,xinit)

% [y,t]=simul(t,xinit)

% t : duration of the simulation

% xinit: initial state

% y: output

global Ncell

% simulation

vectol=zeros(1,8*Ncell+1)+1e-4;

options = odeset('RelTol',1e-4,'AbsTol',vectol);

[t,y] = ode23(@modele,t,xinit,options);% This function computes the output of the ode previously defined.

 
Retrieved from "http://2008.igem.org/Team:Paris/Modeling/simul"