Discussion :

[kbenjamin]

Bonjour,

Je suis entrain d'écrire un code Matlab pour simuler les orbites de la Terre et de Mars autour du soleil.

Pour cela j'ai travaillé en polaire en 2D:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
%Pour la Terre
e=.0167;
a=149.6*10^6;
Re=(a.*(1-e.^2))./(1-e*cos(theta));
polar(theta,Re,'b . -'); 
%Même structure pour Mars

J'obtiens deux jolies ellipses coplanaires qui correspondent à la trajectoire de la terre et de Mars autour du Soleil. Jusqu'à pas de problème.

Cependant maintenant j'aimerai travailler en 3D en rajoutant l'angle de déphasage i entre l'ellipse idéale coplanaire et l'ellipse réelle (voire figure). Et la je suis bloqué je sais que je dois utiliser la donction plot3(x,y,z) pour obtenir l'ellipse en 3D je peux aussi utiliser la fonction pol2cart pour obtenir les coordonnées de x et y à partir de mes calculs. Mais je n'arrive pas à obtenir ce déphasage.

Est ce que quelqu'un a la solution?

Merci pour votre aide.

Benjamin

[duf42]

Bonjour,

Je comprends pas bien, ton problème est pour calculer la nouvelle elipse ou pour la tracer?

Bonne journée,
Duf

[kbenjamin]

Bonjour duf42,

Je ne sais pas ni la calculer ni la tracer. Ce qui me bloque c'est que l'angle i est l'angle entre l'ellipse coplanaire et l'ellipse réelle et non pas l'angle entre l'ellipse réelle et le centre de l'ellipse.

A cause de cela je ne peux pas utiliser la fonction sph2cart...

Une petite idée?