Répondre à la discussion
Affichage des résultats 1 à 4 sur 4

resolution numérique equa diff d'ordre 2



  1. #1
    Heimdall

    resolution numérique equa diff d'ordre 2


    ------

    salut,

    J'essaie de simuler l'intéraction gravitationnelle entre N corps.

    Supposons un systeme S de N corps M_i de masse m_i en interaction gravitationnelle. Soit M_j un corps fixé dans ce système.




    Je dois résoudre pour trouver le mouvement de M_j, l'équation différentielle du second ordre suivante :




    je connais l'algorithme de Runge-Kutta d'ordre 4 pour résoudre les équations différentielles d'ordre 1, et je me demandais si je pouvais l'utiliser pour résoudre celle-ci d'ordre 2 ?

    J'ai pensé a scinder l'équation en un système de deux équations du premier ordre que je ferai passer l'une après l'autre dans runge kutta... est-ce une méthode qui marcherai ?

    le système serait :




    merci

    -----

  2. #2
    deep_turtle

    Re : resolution numérique equa diff d'ordre 2

    Oui, dédoubler le système d'équations en équations du premier ordre est ce qui se fait en général. Je te mets en garde contre un obstacle que tu vas rencontrer, le système que tu dois résoudre est fortement sensible aux erreurs numériques, à cause des interactions entre les particules qui s'approchent trop les unes des autres (ce qui donne des forces énormes et des dérivées énormes aussi).

    Dans le cas d'un système à deux corps, il y a une solution presque miracle qui consiste à faire un changement de variable astucieux dans les équations de départ (coordonnées de Levi-Civita), ça se généralise aux systèmes à N corps mais il faut bosser un peu...

    Je l'ai fait pour un système très simple : un corps attiré par deux autres fixes, et tu arrives comme ça à une précision relative de 10-10 sans trop de difficulté, là où le système de départ part dans le décor assez rapidement (erreurs relatives de l'ordre de l'unité dès qu'une interaction proche a lieu...).

  3. #3
    Heimdall

    Re : resolution numérique equa diff d'ordre 2

    merci deep_turtle pour ta réponse rapide, pour ton truc avec Levi-Civita j'ai tout de même l'impression que c'est un brin compliqué. Je vais essayer sans, et comme ça ça m'apprendra des trucs sur les erreurs numériques et tout, faut voir ça du bon côté.

  4. #4
    deep_turtle

    Re : resolution numérique equa diff d'ordre 2

    Yep tu as raison. Tu pourras toujours revenir si ça va pas comme tu veux !!

Sur le même thème :

Discussions similaires

  1. Equa diff linéaire d'ordre premier
    Par Uranie89 dans le forum Mathématiques du supérieur
    Réponses: 4
    Dernier message: 29/11/2007, 11h38
  2. Etude d'une équa diff d'ordre 1
    Par heperion dans le forum Mathématiques du supérieur
    Réponses: 0
    Dernier message: 04/11/2007, 17h05
  3. Equa diff compliquée (+résolution numérique sous Maple) ?
    Par RVmappeurCS dans le forum Mathématiques du supérieur
    Réponses: 5
    Dernier message: 22/03/2007, 21h27
  4. résolution equa. diff.
    Par Ram dans le forum Mathématiques du supérieur
    Réponses: 35
    Dernier message: 09/12/2005, 19h40
  5. Résolution- Equa Diff
    Par Darklingg dans le forum Mathématiques du supérieur
    Réponses: 9
    Dernier message: 14/04/2005, 12h07