Mathématicien du jour bonjour

Gomly · 9/06/2015 10:13:02

Nouveau problème , nouvelle publication.
Bonjour ,
j'essaie d'ajouter une fonction qui permette à un pilote novice de rejoindre un autre pilote en lui donnant un cap a suivre.

J'ai donc la longitude latitude de mes deux avions et j'applique la formule du cosinus à savoir (cos(t) = u.v / ||u||||v||)
avec u le vecteur qui joint les deux avions et v le vecteur unitaire des ordonnés.
Malheureusement après compilation ca ne me donne pas le résultat escompté. En effet si je suis le bon cap pour rejoindre l'avion le cap indicateur est totalement différent.

Je ne sais pas si c'est une erreur de math , une erreur de code ou une erreur de repère.

Je vous met le bout de code :

if(mPiloteName != "null"){
	  
	  for (unsigned int i = 0; i < 3; ++i) {
		stringstream myString;
		myString << "/ai/models/multiplayer[" << i ;
		string way(myString.str());
		std::string str;
		const char* waychar = way.c_str();
      
		const char* wayname = (way + "]/callsign").c_str();
		
		//const char* waynamechar = wayname.c_str();
		
		string name = fgGetString(wayname);
      
		if(mPiloteName == name){
			const char* latway = (way+"]/position/latitude-deg").c_str();
			double lat = fgGetDouble(latway);
			const char* longway = (way+"]/position/longitude-deg").c_str();
			double longi = fgGetDouble(longway);
			double mylat = fgGetDouble("/position/latitude-deg");
			double mylongi = fgGetDouble("/position/longitude-deg");
			
			double capasuivre = acos((lat - mylat) / (pow(pow(lat-mylat,2)+pow(longi-mylongi,2),0.5)))* 180.0 / M_PI;
			//double capasuivre = ((atan2(lat - mylat,longi - mylongi + M_PI/2)* 180.0 / M_PI)+360)%360;
			cout << lat << endl;
			cout << longi << endl;
			
			cout << "Mon cap   :   "<<fgGetDouble("/orientation/heading-deg")<<endl;
			cout << "Cap à suivre   :   "<<capasuivre << endl;
			cout << name << endl;
		}
      }  
  }

PS: Le code est pas hyper propre si vous avez des idées pour le simplifier hésitez pas !

Merci d'avance

Dernière modification par Gomly (9/06/2015 10:25:34)

ctesc356 · 9/06/2015 11:00:28

Bonjour,
en c++ et en math je ne vais pas aider beaucoup
mais peut-être une piste: il y a cette fonction dans ">Multiplayer >Pilot-List

regarde dans: "data/Nasal/multiplayer.nas" ligne 275 et+, func(id)
et: "data/geo.nas" ligne 184, "course_to"
Je pense que c'est là que ça se goupille

Gomly · 9/06/2015 14:12:03

J'ai as bien compris a quoi le brg correspond ?

ctesc356 · 9/06/2015 15:05:30

Gomly a écrit :

J'ai as bien compris a quoi le brg correspond ?

C'est la direction de cet avion par rapport au tien, donc le cap à prendre pour le rejoindre.

Gomly · 9/06/2015 15:15:00

Mmmmm ça peut être intéressant en effet ça pourra éviter mes calculs un peu foireux même si j'aimerai comprendre pourquoi ses calculs en question ne fonctionne pas ;p

J'attend la réponse des Matheux
Merci en tous cas !

ctesc356 · 9/06/2015 15:22:15

Gomly a écrit :

ça pourra éviter mes calculs un peu foireux même si j'aimerai comprendre pourquoi ses calculs en question ne fonctionne pas ;p

Je ne veux surtout pas te détourner des maths
Mais je pense qu'en décortiquant "geo.nas" et consorts tu trouveras comment c'est fait et peut-être l'erreur dans ton prog.

Gomly · 9/06/2015 15:48:45

J'ai fouillé un peu j'ai trouvé ca

call(func ret = math.mod(math.atan2(math.sin(dlon) * math.cos(dest._lat),
				math.cos(me._lat) * math.sin(dest._lat)
				- math.sin(me._lat) * math.cos(dest._lat)
				* math.cos(dlon)), 2 * math.pi) * R2D, nil, var err = []);

ce qui ressemble pas mal a un calcul d'angle j'ai essayé de reproduire ca dans mon code mais le math.mod me pose soucis je ne sais pas si c'est module , modulo ou autre chose ...
j'y connais pas grand chose en nasal

Clm76 · 10/06/2015 7:51:30

Bonjour,

le math.mod me pose soucis je ne sais pas si c'est module , modulo ou autre chose ...

C'est la fonction modulo en nasal. On la trouve aussi sous la forme math.fmod(). http://wiki.flightgear.org/Nasal_librar … fmod.28.29

Gomly · 10/06/2015 11:18:53

Oui j'ai trouvé la formule de FGFS fonctionne a 2pi pret donc j'ai modifié légèrement et ça me donne de bon résultat !!
Je me bas un peu avec les fichier XML et j'aurai réussi mon truc !

Merci pour vos aides

dany93 · 22/08/2015 12:54:17

Bonjour Gomly,

Tes calculs trigonométriques sont-ils valables avec cette géométrie non euclidienne ? On est ici sur un espace courbe (la terre) et les calculs d'angles à partir des coordonnées n'obéissent évidemment plus aux mêmes lois que la géométrie dans le plan. Sauf peut-être très localement (ce qui devrait te suffire) et pas trop loin de l'équateur.

Enfin, si j'ai bien compris, car je n'ai regardé tes formules que très rapidement... Et j'ai l’impression que c'est strictement de la géométrie plane.

Dany

Forum Flightgear France

Annonce

#1 9/06/2015 10:13:02

Mathématicien du jour bonjour

#2 9/06/2015 11:00:28

Re : Mathématicien du jour bonjour

#3 9/06/2015 14:12:03

Re : Mathématicien du jour bonjour

#4 9/06/2015 15:05:30

Re : Mathématicien du jour bonjour

#5 9/06/2015 15:15:00

Re : Mathématicien du jour bonjour

#6 9/06/2015 15:22:15

Re : Mathématicien du jour bonjour

#7 9/06/2015 15:48:45

Re : Mathématicien du jour bonjour

#8 10/06/2015 7:51:30

Re : Mathématicien du jour bonjour

#9 10/06/2015 11:18:53

Re : Mathématicien du jour bonjour

#10 22/08/2015 12:54:17

Re : Mathématicien du jour bonjour

Pied de page des forums