Universite de Nice L3MASS annee Departement de Mathematiques NOM Date PRENOM Salle et heure
4 pages
Français

Découvre YouScribe en t'inscrivant gratuitement

Je m'inscris

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

Découvre YouScribe en t'inscrivant gratuitement

Je m'inscris
Obtenez un accès à la bibliothèque pour le consulter en ligne
En savoir plus
4 pages
Français
Obtenez un accès à la bibliothèque pour le consulter en ligne
En savoir plus

Description

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


Sujets

Informations

Publié par
Nombre de lectures 18
Langue Français

Exrait

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)
  • Accueil Accueil
  • Univers Univers
  • Ebooks Ebooks
  • Livres audio Livres audio
  • Presse Presse
  • BD BD
  • Documents Documents