Cet ouvrage fait partie de la bibliothèque YouScribe
Obtenez un accès à la bibliothèque pour le lire en ligne
En savoir plus

Universite de Nice L3MASS annee Departement de Mathematiques NOM Date PRENOM Salle et heure

4 pages
Universite de Nice L3MASS, annee 2011-2012 Departement de Mathematiques NOM : Date : . PRENOM : Salle et heure : . Feuille-question du TP 5 Systeme de Lotka-Voltera, son linearise, sa perturbee 1 Le systeme de Lotka-Voltera Considerons le systeme d'equations differentielles { x? = x(1? 2y) y? = ?y(1? 2x) (1) Le code Scilab suivant permet d'en esquisser le champ et en tracer des trajectoires. clear; //////////////// Lotka-Voltera xset(window,0); function xprim=f(x,y); xprim=x*(1-2*y); endfunction; function yprim=g(x,y); yprim=-y*(1-2*x); endfunction; function vprim=www(t,v); vprim=[f(v(1),v(2)),g(v(1),v(2))]'; endfunction; xMin=0;xMax=+1;yMin=0;yMax=+1; fchamp(www,0,xMin :0.1 :xMax,yMin :0.1 :yMax); a=gca(); a.data bounds=[xMin,yMin;xMax,yMax]; Tmax=10;M0=[1.1,1.1]; N=100;petitpas=Tmax/N; t=0 :petitpas :Tmax; for numerotraj

  • function yprim

  • points stationnaires

  • scilab

  • systeme de lotka-voltera

  • solution au systeme

  • scilab pour le linearise

  • linearise du systeme au voisinage du point


Voir plus Voir moins
Universit´edeNice De´partementdeMathe´matiques NOM : PRENOM :
Date : Salle et heure :
L3MASS,anne´e2011-2012 . .
Feuille-question du TP 5 Syste`medeLotka-Voltera,sonlin´earis´e,saperturbe´e
1Lesyst`emedeLotka-Voltera Conside´ronslesyste`med´equationsdie´rentielles 0 x=x(12y) 0 y=y(12x) Le codeScilabsuivant permet d’en esquisser le champ et en tracer des trajectoires. clear ; //////////////// Lotka-Voltera xset("window",0) ; function xprim=f(x,y); xprim=x*(1-2*y); endfunction; function yprim=g(x,y); yprim=-y*(1-2*x); endfunction; function vprim=www(t,v); vprim=[f(v(1),v(2)),g(v(1),v(2))]’ ; endfunction ; xMin=0 ;xMax=+1 ;yMin=0 ;yMax=+1 ; fchamp(www,0,xMin :0.1 :xMax,yMin :0.1 :yMax); a=gca() ;a.data bounds=[xMin,yMin;xMax,yMax] ; Tmax=10 ;M0=[1.1,1.1] ; N=100 ;petitpas=Tmax/N ; t=0 :petitpas :Tmax; for numerotraj=1 :10 M0=[rand(),rand()]’ ; M=ode(M0,0,t,www) ; x=M(1, :);y=M(2, :); plot(x,y) ; end ; 1. Esquissezdes solutions. Quel est le comportement des solutions?
2.Ou`etcommentonte´te´choisieslesconditionsinitiales?
1
(1)
Un pour Un
Permettre à tous d'accéder à la lecture
Pour chaque accès à la bibliothèque, YouScribe donne un accès à une personne dans le besoin