Intégration d'équation différentielle sous python pour un TIPE

[Tim6714]

Bonjour, je suis actuellement étudiant en prépa MP et je rencontre des problèmes sur mon TIPE. J'ai pour but de modéliser la trajectoire d'astéroïdes sous python. J'ai déjà réussi pour une modélisation assez simple en utilisant la méthode de runge kutta 4. exemple mon code.PNGexemple model.PNG Voici mon premier code pour vous donner une idée de mon niveau si cela peut vous aider( Je sais que ce n'est pas très compact mais je ne m'intéresse pas à ca pour l'instant).

J'ai voulu alors modéliser mon problème d'une manière plus compliquée en prenant en compte l'action de certaines autres planètes sur mon astéroïde. En réutilisant la méthode de runge kutta comme avant je fais face à un problème. J'ai pris en compte l'action de la Terre, Jupiter et Saturne ensemble. Je veux d'abord pouvoir représenter comme sur la deuxième photo leurs trajectoires en 3D. Mon problème est que dans ce modèle mes équations sont couplées et je n'arrive pas à écrire un programme sur la méthode de runge kutta qui me permet d'obtenir leurs positions comme sur la première photo. Auriez vous une piste pour m'aider ?
-----

[polo974]

Bonujour,
Il n'y a (presque) rien de pire que du code sous forme d'image.
Pour envoyer du code en gardant l'indentation, il faut "aller en mode avancé" et utiliser la balise code (bouton #)

exemple
    Indenté

Ne pas hésiter à utiliser "prévisualisation du message...

Ah, aussi, mettre les import, même si on a compris que np c'est numpy...

[umfred]

comme leur mouvement est relatif entre elles, il serait peut-être plus simple de prendre le soleil comme point de référence pour dessiner leurs trajectoires ?

[Brinicle]

Bonjour,

Je suppose que pour l'instant il n'y a que le soleil et l'astéroïde, non ?

Déjà, il y a trois fonctions, alors qu'une seule suffirait, cela simplifie déjà pas mal le problème (et évite de calculer plusieurs fois la même chose).

Ensuite, je vois qu'il y a de grosses difficultés avec numpy... il n'y a rien de pire que les listes d'arrays, il faudrait utiliser un seul array (de taille (N+1,6)), ce qui raccourcit aussi énormément le code (ta boucle fait peur). A mon avis, tu devrais d'abord commencer par cela pour pouvoir faire plus compliqué par la suite.

Et je penses que tu te méprends sur le fonctionnement de insert().

[A voir en vidéo sur Futura]