[C] programme sur rotation d'un satellite autour de la Terre

programme sur rotation d'un satellite autour de la Terre [C] - C - Programmation

Marsh Posté le 20-05-2006 à 22:40:34    

Salut à tous. J'aurais simplement eu besoin d'un petit coup de main pour mon programme qui compile bien mais qui ne fonctionne pas comme je le souhaiterais.
 
En fait, je veux juste tracer la trajectoire d'un satellite à partir d'une vitesse et d'une position de départ. J'utilise les coordonnées cartésiennes.
 
Je me demande si l'erreur ne provient pas des fonctions vitessex-vitessey: est-ce que si dans une fonction par exemple je mets vx=vx+1, à l'appel de cette fonction dans une boucle "for" le vx de droite (donc le vx à l'itération précédente) va  pouvoir être utilisé? Je ne sais pas si vous voyez ce que je veux dire, allez voir les fonctions vitessex-vitessey pour mieux comprendre. Merci :)
 
Voici le code:  
 

Code :
  1. [EDIT] erreure de ma part, voir plus bas le bon code que j'ai reposté

Message cité 1 fois
Message édité par PoLuxR le 20-05-2006 à 23:08:08

---------------
Bandit 400
Reply

Marsh Posté le 20-05-2006 à 22:40:34   

Reply

Marsh Posté le 20-05-2006 à 22:47:50    

PoLuxR a écrit :

Salut à tous. J'aurais simplement eu besoin d'un petit coup de main pour mon programme qui compile bien mais qui ne fonctionne pas comme je le souhaiterais.


Tu ne respectes pas les regles de codage du langage C. Le comportement est indéfni.


Compiling: main.c
main.c: In function `main_':
main.c:27: error: implicit declaration of function `vitessex'
main.c:27: warning: nested extern declaration of `vitessex'
main.c:27: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c:27: error: implicit declaration of function `accelerationx'
main.c:27: warning: nested extern declaration of `accelerationx'
main.c:29: error: implicit declaration of function `vitessey'
main.c:29: warning: nested extern declaration of `vitessey'
main.c:29: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c:29: error: implicit declaration of function `accelerationy'
main.c:29: warning: nested extern declaration of `accelerationy'
main.c: At top level:
main.c:37: warning: declaration of 'y0' shadows a global declaration
C:/Program Files/CodeBlocks/include/math.h:251: warning: shadowed declaration is here
main.c:37: warning: no previous prototype for 'vitessex'
main.c: In function `vitessex':
main.c:39: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c:39: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c: At top level:
main.c:44: warning: declaration of 'y0' shadows a global declaration
C:/Program Files/CodeBlocks/include/math.h:251: warning: shadowed declaration is here
main.c:44: warning: no previous prototype for 'vitessey'
main.c: In function `vitessey':
main.c:46: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c:46: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c: At top level:
main.c:52: warning: declaration of 'y0' shadows a global declaration
C:/Program Files/CodeBlocks/include/math.h:251: warning: shadowed declaration is here
main.c:52: warning: no previous prototype for 'accelerationx'
main.c: In function `accelerationx':
main.c:55: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c:55: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c: At top level:
main.c:61: warning: declaration of 'y0' shadows a global declaration
C:/Program Files/CodeBlocks/include/math.h:251: warning: shadowed declaration is here
main.c:61: warning: no previous prototype for 'accelerationy'
main.c: In function `accelerationy':
main.c:64: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c:64: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
Process terminated with status 1 (0 minutes, 0 seconds)
4 errors, 26 warnings


Le plus simple, pour commencer, est de suivre la regle 'définir avant d'utiliser... (petit coup de réindentation au passage...)

Code :
  1. #include<stdio.h>
  2. #include<stdlib.h>
  3. #include<math.h>
  5. #define G 6.67e-11
  6. #define M 5.9736e24
  8. double dt = 1; //déclaration du pas temporel
  10. int vitessex(double x0, double y0, double vxp){        //fonction vitesse selon x (vxp=vx précédent)
  11.    double vx0;
  12.    vx0 = vxp + dt * ( -(G * M * x0) / (pow((pow(x0, 2) + pow(y0, 2)), 1.5)));    //vx0=vx0+dt*ax0
  14.    return (vx0);
  15. }
  17. int vitessey(double x0, double y0, double vyp){        //fonction vitesse selon y (vyp = vy précédent)
  18.    double vy0;
  19.    vy0 = vyp + dt * ( -(G * M * y0) / (pow((pow(x0, 2) + pow(y0, 2)), 1.5)));    //vy0=vy0+dt*ay0
  21.    return (vy0);
  22. }
  25. int accelerationx(double x0, double y0){                    //fonction acceleration selon x
  26.    double ax0;
  28.    ax0 = -(G * M * x0) / (pow((pow(x0, 2) + pow(y0, 2)), 1.5));
  30.    return (ax0);
  31. }
  34. int accelerationy(double x0, double y0){                    //fonction acceleration selon y
  35.    double ay0;
  37.    ay0 = -(G * M * y0) / (pow((pow(x0, 2) + pow(y0, 2)), 1.5));
  39.    return (ay0);
  40. }
  42. int main(){
  44.    double x, y, vx, vy, ax, ay; //déclaration position-vitesse-acceleration
  45.    double t; //déclaration variable temporelle
  47.    x = 35784000; //position initiale du satellite (metres)
  48.    y = 0; //
  49.    vx = 0; //vitesse initiale du satellite (m/s)
  50.    vy = 2.648865e8; // perimetre/période de révolution
  51.    ax = -3.357201e-2; //acceleration initiale du satellite=-v²/r (m/s²)
  52.    ay = 0; //
  54.    FILE *PositionS = fopen("PositionS.out", "w" );
  56.    for (t = 0;t <= 86160;t = t + dt)
  57.    { //developpement limité abscisses & ordonnées
  58.       fprintf(PositionS, "%f %f\n", x, y);       //création du fichier affichant la position du satellite
  60.       x = x + dt * vitessex(x, y, vitessex(x, y)) + ((pow(dt, 2)) / 2) * accelerationx(x, y);    //calcul abscisse position
  62.       y = y + dt * vitessey(x, y) + ((pow(dt, 2)) / 2) * accelerationy(x, y); //calcul ordonnée position
  64.    }
  67.    return 0;
  68. }


Léger mieux, mai c'est pas encore ça...


Compiling: main.c
main.c:10: warning: declaration of 'y0' shadows a global declaration
C:/Program Files/CodeBlocks/include/math.h:251: warning: shadowed declaration is here
main.c:10: warning: no previous prototype for 'vitessex'
main.c: In function `vitessex':
main.c:12: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c:12: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c: At top level:
main.c:17: warning: declaration of 'y0' shadows a global declaration
C:/Program Files/CodeBlocks/include/math.h:251: warning: shadowed declaration is here
main.c:17: warning: no previous prototype for 'vitessey'
main.c: In function `vitessey':
main.c:19: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c:19: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c: At top level:
main.c:25: warning: declaration of 'y0' shadows a global declaration
C:/Program Files/CodeBlocks/include/math.h:251: warning: shadowed declaration is here
main.c:25: warning: no previous prototype for 'accelerationx'
main.c: In function `accelerationx':
main.c:28: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c:28: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c: At top level:
main.c:34: warning: declaration of 'y0' shadows a global declaration
C:/Program Files/CodeBlocks/include/math.h:251: warning: shadowed declaration is here
main.c:34: warning: no previous prototype for 'accelerationy'
main.c: In function `accelerationy':
main.c:37: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c:37: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c: In function `main_':
main.c:60: error: too few arguments to function `vitessex'
main.c:60: warning: passing arg 3 of `vitessex' as floating rather than integer due to prototype
main.c:60: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
main.c:62: error: too few arguments to function `vitessey'
main.c:62: warning: passing arg 2 of `pow' as floating rather than integer due to prototype
Process terminated with status 1 (0 minutes, 0 seconds)
2 errors, 23 warnings


Le plus gros problème est que tu ne passes pas le bon nombre de paramètres aux fonctions vitessex() et vitessey()... Je ne sais pas comment corriger, car je ne connais pas la partie 'métier' de ton code. (Au feeling, je mets vx, vy...)
 
J'ai passé les types retournés en double (je doute fort que tu ais voulu des int !)
 
Ceci compile et fonctionne (un fichier est crée ...)

Code :
  1. #include<stdio.h>
  2. #include<stdlib.h>
  3. #include<math.h>
  5. #define G 6.67e-11
  6. #define M 5.9736e24
  8. //fonction vitesse selon x (vxp=vx précédent)
  9. static double vitessex(double x, double y, double vxp, double dt)
  10. {
  11.    double vx0;
  12.    vx0 = vxp + dt * ( -(G * M * x) / (pow((pow(x, 2.0) + pow(y, 2.0)), 1.5)));    //vx0=vx0+dt*ax0
  14.    return vx0;
  15. }
  17. //fonction vitesse selon y (vyp = vy précédent)
  18. static double vitessey(double x, double y, double vyp, double dt)
  19. {
  20.    double vy0;
  21.    vy0 = vyp + dt * ( -(G * M * y) / (pow((pow(x, 2.0) + pow(y, 2.0)), 1.5)));    //vy0=vy0+dt*ay0
  23.    return vy0;
  24. }
  26. //fonction acceleration selon x
  27. static double accelerationx(double x, double y)
  28. {
  29.    double ax0;
  31.    ax0 = -(G * M * x) / (pow((pow(x, 2.0) + pow(y, 2.0)), 1.5));
  33.    return ax0;
  34. }
  36. //fonction acceleration selon y
  37. static int accelerationy(double x, double y)
  38. {
  39.    double ay0;
  41.    ay0 = -(G * M * y) / (pow((pow(x, 2.0) + pow(y, 2.0)), 1.5));
  43.    return ay0;
  44. }
  46. int main()
  47. {
  48.    /* invariants */
  49.    double const dt = 1; //déclaration du pas temporel
  50.    double const vx = 0; //vitesse initiale du satellite (m/s)
  51.    double const vy = 2.648865e8; // perimetre/période de révolution
  52.    double const ax = -3.357201e-2; //acceleration initiale du satellite=-v²/r (m/s²)
  53.    double const ay = 0; //
  54.    /* variables */
  55.    double x = 35784000; //position initiale du satellite (metres)
  56.    double y = 0; //
  57.    double t; //déclaration variable temporelle
  59.    FILE *PositionS = fopen("PositionS.out", "w" );
  61.    if (PositionS != NULL)
  62.    {
  63.       for (t = 0;t <= 86160;t = t + dt)
  64.       {
  65.          //developpement limité abscisses & ordonnées
  66.          //création du fichier affichant la position du satellite
  67.          fprintf(PositionS, "%f %f\n", x, y);
  68.          //calcul abscisse position
  69.          x = x + dt * vitessex(x, y, vitessex(x, y, vx, dt), dt)
  70.              + ((pow(dt, 2.0)) / 2.0) * accelerationx(x, y);
  71.          //calcul ordonnée position
  72.          y = y + dt * vitessey(x, y, vy, dt) + ((pow(dt, 2.0)) / 2.0) * accelerationy(x, y);
  73.       }
  74.       fclose (PositionS) , PositionS = NULL;
  75.    }
  76.    return 0;
  77. }


Message édité par Emmanuel Delahaye le 20-05-2006 à 23:17:31

---------------
Des infos sur la programmation et le langage C: http://www.bien-programmer.fr Pas de Wi-Fi à la maison : http://www.cpl-france.org/
Reply

Marsh Posté le 20-05-2006 à 23:05:58    

Ah mince! Je vous ai donné le mauvais programme, un truc faux que j'avais essayé, du coup c'était dur pour vous de bien comprendre la question!! Voici le "bon"code:
 

Code :
  1. #include<stdio.h>
  2. #include<stdlib.h>
  3. #include<math.h>
  5. #define G 6.67e-11
  6. #define M 5.9736e24
  8. double dt=1; //déclaration du pas temporel
  10. int main(){
  12. double x,y,vx,vy,ax,ay; //déclaration position-vitesse-acceleration
  13. double t; //déclaration variable temporelle
  15. x=35784000;//position initiale du satellite (metres)
  16. y=0; //
  17. vx=0; //vitesse initiale du satellite (m/s)
  18. vy=2.648865e8; // perimetre/période de révolution
  19. ax=-3.357201e-2; //acceleration initiale du satellite=-v²/r (m/s²)
  20. ay=0; //
  22. FILE *PositionS = fopen("PositionS.out","w" );
  24. for(t=0;t<=86160;t=t+dt){ //developpement limité abscisses & ordonnées
  25.     fprintf(PositionS,"%f %f\n",x,y);       //création du fichier affichant la position du satellite
  26.    
  27.     x=x+dt*vitessex(x,y)+((pow(dt,2))/2)*accelerationx(x,y);    //calcul abscisse position
  28.    
  29.     y=y+dt*vitessey(x,y)+((pow(dt,2))/2)*accelerationy(x,y);//calcul ordonnée position
  30.    
  31. }
  34. return 0;
  35. }
  37. int vitessex(double x0,double y0){        //fonction vitesse selon x
  38.     double vx0;
  39.     vx0=vx0+dt*(-(G*M*x0)/(pow((pow(x0,2)+pow(y0,2)),1.5)));    //vx0=vx0+dt*ax0
  40.    
  41.     return(vx0);
  42. }
  44. int vitessey(double x0,double y0){        //fonction vitesse selon y  
  45.     double vy0;
  46.     vy0=vy0+dt*(-(G*M*y0)/(pow((pow(x0,2)+pow(y0,2)),1.5)));    //vy0=vy0+dt*ay0
  47.    
  48.     return(vy0);
  49. }
  52. int accelerationx(double x0,double y0){                    //fonction acceleration selon x
  53.     double ax0;
  54.    
  55.     ax0=-(G*M*x0)/(pow((pow(x0,2)+pow(y0,2)),1.5));
  56.    
  57.     return(ax0);
  58. }
  61. int accelerationy(double x0,double y0){                    //fonction acceleration selon y
  62.     double ay0;
  63.    
  64.     ay0=-(G*M*y0)/(pow((pow(x0,2)+pow(y0,2)),1.5));
  65.    
  66.     return(ay0);
  67. }


 

Message cité 1 fois

---------------
Bandit 400
Reply

Marsh Posté le 20-05-2006 à 23:18:42    

PoLuxR a écrit :

Ah mince! Je vous ai donné le mauvais programme, un truc faux que j'avais essayé, du coup c'était dur pour vous de bien comprendre la question!! Voici le "bon"code:


Merci, je viens de passe 30 minutes à corriger ton ancien code... Je te laisse faire la transposition...
 


---------------
Des infos sur la programmation et le langage C: http://www.bien-programmer.fr Pas de Wi-Fi à la maison : http://www.cpl-france.org/
Reply

Marsh Posté le 21-05-2006 à 15:36:52    

Merci beaucoup Emmanuel. J'ai fais la transposition et j'obtiens un fichier de sortie:
35784000.000000 0.000000
-1.#QNAN0         -1.#QNAN0
-1.#QNAN0         -1.#QNAN0
-1.#QNAN0         -1.#QNAN0
etc...
 
Mais comme je l'ai dit en début de topic, mon problème n'est pas un problème de compilation, mais un problème en sortie de programme car le fichier m'affiche des valeurs fausses (surtout les ordonnées qui stagnent)! Donc si quelqu'un a une idée....  
 
Merci :)


---------------
Bandit 400
Reply

Marsh Posté le 21-05-2006 à 17:18:09    

Code :
  1. #include<stdio.h>
  2. #include<stdlib.h>
  3. #include<math.h>
  5. #define G 6.67e-11
  6. #define M 5.9736e24
  8. double dt=1; //déclaration du pas temporel
  10. int main(){
  12. double x,y,vx,vy,ax,ay; //déclaration position-vitesse-acceleration
  13. double t; //déclaration variable temporelle
  15. x=35784000;//position initiale du satellite (metres)
  16. y=0; //
  17. vx=0; //vitesse initiale du satellite (m/s)
  18. vy=2.648865e8; // perimetre/période de révolution
  19. ax=-3.357201e-2; //acceleration initiale du satellite=-v²/r (m/s²)
  20. ay=0; //
  22. FILE *PositionS = fopen("PositionS.out","w" );
  24. for(t=0;t<=86160;t=t+dt){ //developpement limité abscisses & ordonnées
  25.     fprintf(PositionS,"%f %f\n",x,y);       //création du fichier affichant la position du satellite
  27.     vx=vitessex(x,y);
  28.     ax=accelerationx(x,y);
  29.     x=x+dt*vx+((pow(dt,2))/2)*ax;    //calcul abscisse position
  30.    
  31.     vy=vitessey(x,y);
  32.     ay=accelerationy(x,y);
  33.     y=y+dt*vy+((pow(dt,2))/2)*ay;//calcul ordonnée position
  34.    
  35. }
  38. return 0;
  39. }
  41. int vitessex(double x0,double y0){        //fonction vitesse selon x
  42.     double vx0;
  43.     vx0=vx0+dt*(-(G*M*x0)/(pow((pow(x0,2)+pow(y0,2)),1.5)));    //vx0=vx0+dt*ax0
  44.    
  45.     return(vx0);
  46. }
  48. int vitessey(double x0,double y0){        //fonction vitesse selon y  
  49.     double vy0;
  50.     vy0=vy0+dt*(-(G*M*y0)/(pow((pow(x0,2)+pow(y0,2)),1.5)));    //vy0=vy0+dt*ay0
  51.    
  52.     return(vy0);
  53. }
  56. int accelerationx(double x0,double y0){                    //fonction acceleration selon x
  57.     double ax0;
  58.    
  59.     ax0=-(G*M*x0)/(pow((pow(x0,2)+pow(y0,2)),1.5));
  60.    
  61.     return(ax0);
  62. }
  65. int accelerationy(double x0,double y0){                    //fonction acceleration selon y
  66.     double ay0;
  67.    
  68.     ay0=-(G*M*y0)/(pow((pow(x0,2)+pow(y0,2)),1.5));
  69.    
  70.     return(ay0);
  71. }

Message cité 2 fois

---------------
Bandit 400
Reply

Marsh Posté le 21-05-2006 à 18:15:28    

PoLuxR a écrit :

Code :
  1. #include<stdio.h>
  2. #include<stdlib.h>
  3. #include<math.h>
  5. #define G 6.67e-11
  6. #define M 5.9736e24
  8. double dt=1; //déclaration du pas temporel
  10. int main(){
  12. <...>



Tu t'obstines à ne pas respecter les regles du langage. Je t'ai indiqué comment faire. Fait le ou je ne peux plus rien pour toi.


---------------
Des infos sur la programmation et le langage C: http://www.bien-programmer.fr Pas de Wi-Fi à la maison : http://www.cpl-france.org/
Reply

Marsh Posté le 21-05-2006 à 22:29:51    

PoLuxR a écrit :

Code :
  1. #include<stdio.h>
  2. #include<stdlib.h>
  3. #include<math.h>
  5. #define G 6.67e-11
  6. #define M 5.9736e24
  8. double dt=1; //déclaration du pas temporel
  10. int main(){
  12. double x,y,vx,vy,ax,ay; //déclaration position-vitesse-acceleration
  13. double t; //déclaration variable temporelle
  15. x=35784000;//position initiale du satellite (metres)
  16. y=0; //
  17. vx=0; //vitesse initiale du satellite (m/s)
  18. vy=2.648865e8; // perimetre/période de révolution
  19. ax=-3.357201e-2; //acceleration initiale du satellite=-v²/r (m/s²)
  20. ay=0; //
  22. FILE *PositionS = fopen("PositionS.out","w" );
  24. for(t=0;t<=86160;t=t+dt){ //developpement limité abscisses & ordonnées
  25.     fprintf(PositionS,"%f %f\n",x,y);       //création du fichier affichant la position du satellite
  27.     vx=vitessex(x,y);
  28.     ax=accelerationx(x,y);
  29.     x=x+dt*vx+((pow(dt,2))/2)*ax;    //calcul abscisse position
  30.    
  31.     vy=vitessey(x,y);
  32.     ay=accelerationy(x,y);
  33.     y=y+dt*vy+((pow(dt,2))/2)*ay;//calcul ordonnée position
  34.    
  35. }
  38. return 0;
  39. }
  41. int vitessex(double x0,double y0){        //fonction vitesse selon x
  42.     double vx0;
  43.     vx0=vx0+dt*(-(G*M*x0)/(pow((pow(x0,2)+pow(y0,2)),1.5)));    //vx0=vx0+dt*ax0
  44.    
  45.     return(vx0);
  46. }
  48. int vitessey(double x0,double y0){        //fonction vitesse selon y  
  49.     double vy0;
  50.     vy0=vy0+dt*(-(G*M*y0)/(pow((pow(x0,2)+pow(y0,2)),1.5)));    //vy0=vy0+dt*ay0
  51.    
  52.     return(vy0);
  53. }
  56. int accelerationx(double x0,double y0){                    //fonction acceleration selon x
  57.     double ax0;
  58.    
  59.     ax0=-(G*M*x0)/(pow((pow(x0,2)+pow(y0,2)),1.5));
  60.    
  61.     return(ax0);
  62. }
  65. int accelerationy(double x0,double y0){                    //fonction acceleration selon y
  66.     double ay0;
  67.    
  68.     ay0=-(G*M*y0)/(pow((pow(x0,2)+pow(y0,2)),1.5));
  69.    
  70.     return(ay0);
  71. }



Les fonctions "vitesseX" et "vitesseY" sont très similaires. Même remarque pour "accelerationX" et "accelerationY"
Peut-être serait-il pratique de les regrouper, quitte à y inclure un paramètre supplémentaire pour faire le choix "X" ou "Y".
Petit pb: tes fonctions de type "int" retournent des doubles...
 
 


---------------
Vous ne pouvez pas apporter la prospérité au pauvre en la retirant au riche.
Reply

Marsh Posté le 09-06-2006 à 01:33:19    

Merci à tous pour votre aide. Finalement j'ai réussi à terminer ce programme, avec simplement une erreure à un scalaire multiplicatif près. Voici les lignes, si ça peut aider quelqu'un qui fait une recherche à l'avenir:
 

Code :
  1. #include<stdio.h>
  2. #include<stdlib.h>
  3. #include<math.h>
  5. double GM=3.98439e8; //G en N.km².kg-², à la puissance -17 car nous sommes en km pour éviter les calculs trop gros  
  7. double dt=1; //déclaration du pas temporel
  10. double vitessex(double x0,double y0,double vx){        //fonction vitesse selon x
  11.     double vx0;
  12.     vx0=vx+dt*(-(GM*x0)/(pow((pow(x0,2)+pow(y0,2)),1.5)));    //vx0=vx0+dt*ax0
  13.    
  14.     return(vx0);
  15. }
  17. double vitessey(double x0,double y0, double vy){        //fonction vitesse selon y  
  18.     double vy0;
  19.     vy0=vy+dt*(-(GM*y0)/(pow((pow(x0,2)+pow(y0,2)),1.5)));    //vy0=vy0+dt*ay0
  20.    
  21.     return(vy0);
  22. }
  25. double accelerationx(double x0,double y0){                    //fonction acceleration selon x
  26.     double ax0;
  27.    
  28.     ax0=-(GM*x0)/(pow((pow(x0,2)+pow(y0,2)),1.5)); //ax0=(-GM*x0)/((x²+y²)^1.5)
  29.    
  30.     return(ax0);
  31. }
  34. double accelerationy(double x0,double y0){                    //fonction acceleration selon y
  35.     double ay0;
  36.    
  37.     ay0=-(GM*y0)/(pow((pow(x0,2)+pow(y0,2)),1.5)); //ay0=(-GM*y0)/((x²+y²)^1.5)
  38.    
  39.     return(ay0);
  40. }
  45. int main(){
  47. double x,y,vx,vy,ax,ay; //déclaration position-vitesse-acceleration
  48. double t; //déclaration variable temporelle
  50. printf("\t Bonjour, choisissez une valeur d'abscisse (x) pour la position initiale du satellite (en km)\n" );//choix de x
  51. scanf("%lf",&x);
  53. printf("\t Choisissez une valeur d'ordonnée (y) pour la position initiale du satellite (en km)\n" );//choix de y
  54. scanf("%lf",&y);
  56. printf("\t Choisissez une valeur pour la composante de la vitesse du satellite selon les abscisses (vx en km/s)\n" );//vx
  57. scanf("%lf",&vx);
  59. printf("\t Choisissez une valeur pour la composante de la vitesse du satellite selon les ordonnées (vy en km/s)\n" );//vy
  60. scanf("%lf",&vy);
  63. ax=0; //acceleration initiale du satellite
  64. ay=0; //
  66. FILE *PositionS = fopen("orbite.out","w" );
  68. for(t=0;t<=86160;t=t+dt){ //developpement limité abscisses & ordonnées
  69.     fprintf(PositionS,"%f %f\n",x,y);       //création du fichier affichant la position du satellite
  71.     vx=vitessex(x,y,vx);  // assimilation des appels de fonctions à un nom (vx,vy,ax,ay) pour pouvoir utiliser  
  72.     vy=vitessey(x,y,vy);  // la condition initiale
  73.     ax=accelerationx(x,y);//
  74.     ay=accelerationy(x,y);//
  77.     x=x+dt*vx+((pow(dt,2))/2.)*ax;    //calcul abscisse position
  78.     y=y+dt*vy+((pow(dt,2))/2.)*ay; //calcul ordonnée position
  80.     if (sqrt(x*x+y*y)<6374) break;    //si racine de (x²+y²) inférieur au rayon de la Terre, on stoppe l'itération
  82. }
  84. fclose(PositionS);
  86. FILE *graphe = fopen("graphe.out","w" );
  87. fprintf(graphe,"set terminal gif\n" );
  88. fprintf(graphe,"set output 'Orbite.gif'\n" );
  89. fprintf(graphe,"plot 'orbite.out' u 1:2\n" );
  90. fprintf(graphe,"q\n" );
  91. fclose(graphe);
  93. system("gnuplot graphe.out" );
  95. system("konqueror Orbite.gif" );
  96. return 0;
  97. }

Message cité 1 fois

---------------
Bandit 400
Reply

Marsh Posté le 09-06-2006 à 13:55:33    

PoLuxR a écrit :

Code :
  1. double accelerationx(double x0,double y0){                    //fonction acceleration selon x
  2.     double ax0;
  3.    
  4.     ax0=-(GM*x0)/(pow((pow(x0,2)+pow(y0,2)),1.5)); //ax0=(-GM*x0)/((x²+y²)^1.5)
  5.    
  6.     return(ax0);
  7. }
  10. double accelerationy(double x0,double y0){                    //fonction acceleration selon y
  11.     double ay0;
  12.    
  13.     ay0=-(GM*y0)/(pow((pow(x0,2)+pow(y0,2)),1.5)); //ay0=(-GM*y0)/((x²+y²)^1.5)
  14.    
  15.     return(ay0);
  16. }



 

Code :
  1. double accelerationxy(double x0, double y0, double axe){                    //fonction acceleration selon x ou y au choix
  2.      return (-GM * axe)/(pow((pow(x0,2)+pow(y0,2)),1.5));
  3. }


Et tu appelles "accelerationxy()" en lui passant en 3° paramètre "x0" si tu veux une accélération selon X, ou "y0" si tu veux une accélération selon Y.
On peut même remplacer cette fonction par un "#define"...
 
Même méthode pour "vitessex" et "vitessey"
 :sol:


Message édité par Sve@r le 09-06-2006 à 22:39:28

---------------
Vous ne pouvez pas apporter la prospérité au pauvre en la retirant au riche.
Reply

Sujets relatifs:

Leave a Replay

Make sure you enter the(*)required information where indicate.HTML code is not allowed