Poster un nouveau sujet Poster une réponse
Géodésie
Auteur Message
ld
Avatar

Autorisation : Membre
Nb de messages : 373
Inscrit le : Jeu 23 Sep 2010, 16:05
Posté le : Dim 15 Fév 2015, 19:55   Citer 

EDIT
Bonsoir les filles !
Voilà de quoi calculer la distance entre deux points de latitude/longitude connues sur notre planète.

Code
:Degré
:Disp "LATITUDE°
:Prompt A,B
:Disp "LONGITUDE°
:Prompt C,D
:Disp 6378arccos(sin(A)sin(B)+cos(A)cos(B)cos(D-C))pi/180,"KM"


Mes sources :
http://web.archive.org/web/20250310103330/https://geodesie.ign.fr/contenu/fichiers/Distance_longitude_latitude.pdf
https://integraledesmaths.pages-perso.free.fr/idm/PagePrincipale.htm#https://integraledesmaths.pages-perso.free.fr/idm/GeoAEDistSph.htm
https://www.dcode.fr/distance-geographique
Et Linkakro évidemment ;D

----------------------
Note Éphéméride Chrono Bombrman BlueFishEncryptor Labygénérateur Nord BandeauDeTexte LunarPhase Syracuse Papy-Jump Convertisseur ! SolveurDeSudoku VérificateurDeCodesBanquaire PrgmADiviser CryptageDeVigenere MiniSolv Reed-Solomon Code128Géodésie
 Adresse email Haut de page Bas de page 
 
linkakro



Autorisation : Membre
Nb de messages : 3738
Inscrit le : Lun 19 Oct 2009, 21:25
Posté le : Dim 15 Fév 2015, 23:02   Citer 


Citer
6378000arccos(sin(phiA)sin(phiB)+cos(phiA)cos(phiB)cos(lamdaB-lamdaA))


edit: le message précédent mentionnait des problèmes de calculs, mais cela a été résolu plus tard et le message modifié

Le produit rayon*angle nécessite les radians. Les calculs intermédiaires peuvent utiliser autant degré que radian pourvu qu'on soit cohérent.
Tu peux forcer les unités indépendamment du mode avec les symboles du menu angle en suffixes.

J'ai saisi dans une ti84plus le calcul du cas particulier du demi équateur et le résultat est cohérent (à moins d'1% près) : 20_037_077 mètres.
Ta saisie doit être erronée.
Je songe que arccos() en français est nommée cos^-1() en anglais dans les TI. Phi est la latitude, lamda la longitude.

----------------------
ti82statfr: 2008, inscrit: 2009, ti84pocketfr: noël2011, ti30xbmultiview: iut 2012-2014
Perfectionniste, manque tact. Pas le temps de tout publier depuis 2011. Répond toujours aux questions. (rédigé juin 2014)

Pour tout le monde et surtout les débutants, quelques-uns des articles courants :
*Traductions Algorithmie/Ti-Basic.
*Caractères spéciaux sur Tout82
Les défauts du TI-Basic : Goto_versus_algo et DelVar/End/Lbl/guillemet/store
 Adresse email Haut de page Bas de page 
 
Google_AdSense



Publicité    



 
ld
Avatar

Autorisation : Membre
Nb de messages : 373
Inscrit le : Jeu 23 Sep 2010, 16:05
Posté le : Dim 15 Fév 2015, 23:30   Citer 

Ok c'est bon je viens de comprendre tout seul.
Il fallait mettre PI au lieu de 180.
En fait faut d'abord convertir la latitude et la longitude en radian.
Dsl, j'suis nul on m'a donner mon bac S avec 4/20 en maths, c'est pas une blague ! cry.gif
Pitié enterez ce post moodforlove.gif

----------------------
Note Éphéméride Chrono Bombrman BlueFishEncryptor Labygénérateur Nord BandeauDeTexte LunarPhase Syracuse Papy-Jump Convertisseur ! SolveurDeSudoku VérificateurDeCodesBanquaire PrgmADiviser CryptageDeVigenere MiniSolv Reed-Solomon Code128Géodésie
 Adresse email Haut de page Bas de page 
 
linkakro



Autorisation : Membre
Nb de messages : 3738
Inscrit le : Lun 19 Oct 2009, 21:25
Posté le : Dim 15 Fév 2015, 23:41   Citer 

menu angle 2nde+matrice.
^r et ° représentent radian et degrés. Les coefficients 180/pi et pi/180 sont appliqués selon le mode.
^r applique 180/pi pour (convertir radian en degré) si le mode est degré, applique 1 en mode radian.
° applique pi/180 (pour convertir degré en radian) si le mode est radian, applique 1 en mode degré.

On peut donc écrire une formule comme celle-ci pour l'opération directe.

Code
 // variables A,B,C,D en degrés !
 // expression qui s'adapte au mode de la calculatrice
sin(A°)sin(B°)+ cos(A°)cos(B°)cos((C-D)°)


-- Problème :
On ne peut pas simplement convertir à volonté, mais cela tient bien compte du mode de la calculatrice.
On peut facilement adapter à l'antécédent des fonction ainsi, mais c'est plus compliqué pour l'image des réciproques.
Exemple dans notre cas : avec des entrées en degrés et le mode degré, arccos produira des degrés, alors que le produit rayon*angle nécessite des radians. Puisque les opérateurs ne permettent pas dans ce mode degré de convertir en radian, il faut le faire nous-même, et on perd la possibilité d'utiliser le mode radian avec la même formule. Ainsi j'impose le mode degré.


Code
Disp "°LAT
Prompt A,B
Disp "°LONG
Prompt C,D
Degree // mode degré impératif pour cette expression !
 // c'est impératif à cause de l'image du arccos
Disp 6378000*arccos( sin(A)sin(B)+ cos(A)cos(B)cos(D-C) ) *pi/180

edit : corrigé C-D en D-C (signalé message suivant)

----------------------
ti82statfr: 2008, inscrit: 2009, ti84pocketfr: noël2011, ti30xbmultiview: iut 2012-2014
Perfectionniste, manque tact. Pas le temps de tout publier depuis 2011. Répond toujours aux questions. (rédigé juin 2014)

Pour tout le monde et surtout les débutants, quelques-uns des articles courants :
*Traductions Algorithmie/Ti-Basic.
*Caractères spéciaux sur Tout82
Les défauts du TI-Basic : Goto_versus_algo et DelVar/End/Lbl/guillemet/store
 Adresse email Haut de page Bas de page 
 
ld
Avatar

Autorisation : Membre
Nb de messages : 373
Inscrit le : Jeu 23 Sep 2010, 16:05
Posté le : Lun 16 Fév 2015, 0:41   Citer 

J'aime l'idée, je prends !
Par contre c'est (D-C) et pas (C-D).

----------------------
Note Éphéméride Chrono Bombrman BlueFishEncryptor Labygénérateur Nord BandeauDeTexte LunarPhase Syracuse Papy-Jump Convertisseur ! SolveurDeSudoku VérificateurDeCodesBanquaire PrgmADiviser CryptageDeVigenere MiniSolv Reed-Solomon Code128Géodésie
 Adresse email Haut de page Bas de page 
 
Poster un nouveau sujet Poster une réponse





  Page générée en 6 requêtes