1
<?xml version="1.0" encoding="ISO-8859-1" standalone="no"?>
2
<!DOCTYPE MAN SYSTEM "../../manrev.dtd">
4
<LANGUAGE>eng</LANGUAGE>
6
<TYPE>Scilab Function</TYPE>
7
<DATE>April 1993</DATE>
8
<SHORT_DESCRIPTION name="systems"> a collection of dynamical system</SHORT_DESCRIPTION>
10
<CALLING_SEQUENCE_ITEM>[]=systems() </CALLING_SEQUENCE_ITEM>
14
A call to this function will load into Scilab a set of macros which
15
describes dynamical systems. Their parameters can be initiated by
16
calling the routine tdinit().
19
<SECTION label="Bioreact">
30
<SP>is the biomass concentration</SP>
33
<SP>is the sugar concentration</SP>
38
xdot(1)=mu_td(x(2))*x(1)- debit*x(1);
39
xdot(2)=-k*mu_td(x(2))*x(1)-debit*x(2)+debit*x2in;
43
where mu_td is given by
51
<SECTION label="Compet">
54
[xdot]=compet(t,x [,u])
58
a competition model. <VERB>x(1),x(2)</VERB> stands for two populations which grows on a same resource. <VERB>1/u</VERB> is the level of that resource ( 1 is the default value).
63
xdot(1) = ppr*x(1)*(1-x(1)/ppk) - u*ppa*x(1)*x(2) ,
64
xdot(2) = pps*x(2)*(1-x(2)/ppl) - u*ppb*x(1)*x(2) ,
68
"The macro <VERB>[xe]=equilcom(ue)</VERB>" computes an equilibrium point of the competition model and a fixed level of the resource ue ( default value is 1)
71
"The macro <VERB>[f,g,h,linsy]=lincomp([ue])</VERB>" gives the linearisation of the competition model ( with output y=x) around the equilibrium point xe=equilcom(ue). This macro returns [f,g,h] the three matrices of the linearised system. and linsy which is a Scilab macro [ydot]=linsy(t,x) which computes the dynamics of the linearised system
74
<SECTION label="Cycllim">
81
a model with a limit cycle
85
xdot=a*x+qeps(1-||x||**2)x
89
<SECTION label="Linear">
99
<SECTION label="Blinper">
106
a linear system with quadratic perturbations.
109
<SECTION label="Pop">
116
a fish population model
120
xdot= 10*x*(1-x/K)- peche(t)*x
124
<SECTION label="Proie">
126
a Predator prey model with external insecticide.
135
<SP>The prey ( that we want to kill )</SP>
138
<SP>the predator ( that we want to preserve )</SP>
141
<SP>mortality rate due to insecticide which destroys both prey and predator ( default value u=0)</SP>
146
xdot(1) = ppr*x(1)*(1-x(1)/ppk) - ppa*x(1)*x(2) - u*x(1);
147
xdot(2) = -ppm*x(2) + ppb*x(1)*x(2) - u*x(2);
151
The macro <VERB>[xe]=equilpp([ue])</VERB> computes the equilibrium point of the <VERB>p_p</VERB> system for the value <VERB>ue</VERB> of the command. The default value for <VERB>ue</VERB> is 0.</P>
154
xe(1) = (ppm+ue)/ppb;
155
xe(2) = (ppr*(1-xe(1)/ppk)-ue)/ppa;
159
<SECTION label="Lincom">
166
linear system with a feedback</P>