Forum |  HardWare.fr | News | Articles | PC | S'identifier | S'inscrire | Shop Recherche
1986 connectés 

  FORUM HardWare.fr
  Programmation
  Python

  Intégration méthode d'Euler

 



 Mot :   Pseudo :  
 
Bas de page
Auteur Sujet :

Intégration méthode d'Euler

n°1938691
fusionfroi​de
Posté le 07-11-2009 à 23:18:00  profilanswer
 

Salut.
 
J'essaie de résoudre avec la méthode d'Euler (pour commencer) deux équations différentielles (imbriquées l'une dans l'autre) qui décrivent le mouvement d'un proton dans un cyclotron. Les équations différentielles que j'obtiens sont les suivantes : dVx=w²Vy et dVy=w²Vx  où est une constante. Vy : vitesse suivant y. Vx : vitesse suivant x. Donc dVx et dVy sont les dérivée de ces vitesses. condition initiale : Vy(0) = Vy0
 
Dans mon script, l'utilisation de la boucle for me pose le problème suivant :
 
normalement pour utiliser une boucle for il faut définir une séquence avant. Par exemple je définis une séquence de valeurs de Vx puis je l'utiliser dans une boucle for pour trouver les valeurs de Vy en faisant :for i in Vx ...ect. Le problème dans mon script c'est que je cherche à la fois les valeurs de Vy et celles de Vx.
 
Comment régler ce problème ? faut il s'en passer de la boucle for ? merci ;)
 

Code :
  1. from pylab import *
  2. import numpy as N
  3. import math as M
  4. q=1.602e-19
  5. B=1
  6. m=1.67e-27
  7. pas=input('saisir la valeur du pas=')
  8. def dVy(Vx):
  9. return (-(q*B)/m)*Vx
  10. def dVx(Vy):
  11. return ((q*B)/m)*Vy
  12. def Euler(pas):
  13.     Vyi=5
  14.     Vy=([Vyi])
  15.     for [quoi] in [quoi] :
  16.    
  17.         Vxi1=Vyi+dVx(Vyi)*pas
  18.         Vxi1=Vxi
  19.         Vyi1=Vxi+dVy(Vxi)*pas
  20.         Vyi1=Vyi
  21.        
  22.     Vy=N.append(Vy,Vyi1)
  23.     Vx=N.append(Vx,Vxi1)
  24.     Vy=N.append(Vy,10)
  25.     Vx=N.append(Vx,10)
  26.     print Vx
  27.     print Vy
  28.    
  29. Euler(pas)


 

mood
Publicité
Posté le 07-11-2009 à 23:18:00  profilanswer
 

n°1938704
el muchach​o
Comfortably Numb
Posté le 08-11-2009 à 09:56:54  profilanswer
 

Il ne faut pas utiliser for mais un do while ou équivalent, avec une condition sur la précision, par ex.


Aller à :
Ajouter une réponse
  FORUM HardWare.fr
  Programmation
  Python

  Intégration méthode d'Euler

 

Sujets relatifs
PHP5 / JavaBridge : IntégrationMéthode find et propriété after
La methode Range de l'objet Worksheets a échouébordure autour d'un lien (intégration e-mailing)
suppresion automatiquz d'une methodesuppresion automatique d'une methode de notre programme java
Entity Manager et Session beans methode d'updateintegration paypal
Integration de Birt dans le serveur d'application JBOSS GA 5.0.1Utilisation de la méthode Cells
Plus de sujets relatifs à : Intégration méthode d'Euler


Copyright © 1997-2022 Hardware.fr SARL (Signaler un contenu illicite / Données personnelles) / Groupe LDLC / Shop HFR