Fermer X

Contents

1 L’attracteur de Lorenz

1.1 Script du fichier lorenz.sci

Le système dynamique est codé dans ce fichier, qui est appelé dans le fichier lorenz.sce.

  function ydot=lorenz(t,y,a)
    ydot=[sig*(y(2)-y(1));ro*y(1)-y(2)-y(1)*y(3);y(1)*y(2)-beta*y(3)];
  endfunction

1.2 Script du fichier lorenz.sce

Ce fichier permet de faire tourner le système dynamique que l’on vient de coder. Il est lancé par la commande exec lorenz.sce

  xbasc();
  pas=.01;
  t=0:pas:20;
  t0=0;
  y0=[-3,-6,12]';
  y=[];
  exec('lorenz.sce');
  a=12.5;
  y=ode(y0,t0,t,list(lorenz,a));
  param3d(y(1,:),y(2,:),y(3,:));
  read(%io(1),1,1,'(a)');
  xbasc();
  plot2d(y(2,:),y(3,:));
L'École des Ponts ParisTech est membre fondateur de

L'École des Ponts ParisTech est certifiée 

ParisTech ParisEst ParisTech

 

Copyright 2014
École des Ponts ParisTech
All rights reserved