Ho aggiunto anche gli ltri 2 risultati...è probabile che alcuni risultino doppi oppure non sia possibile calcolarli (adesso faccio le prove con DErive, se lo trovo, perchè con Matlab fa piuttosto schifo la cosa...)
Tranquillo![]()
Twitter @veke87
Adesso ci sono...(lo studio logora...)
Per lat2 ottengo 4 valori:
>> syms k distanza lat1 x long1 long2
>> eq= k * acos( sin(lat1) * sin(x) + cos(lat1) * cos(x) * cos(long2- long1))-distanza
>> solve(eq)
ans =
atan((-cos(lat1)*(-cos(lat1)*cos(-long2+long1)*cos(distanza/k)+(cos(lat1)^2*cos(-long2+long1)^2*cos(distanza/k)^2+sin(distanza/k)^2-cos(lat1)^2-cos(lat1)^2*sin(-long2+long1)^2*sin(distanza/k)^2+cos(lat1)^4*sin(-long2+long1)^2)^(1/2))/(-1+cos(lat1)^2*sin(-long2+long1)^2)*cos(-long2+long1)+cos(distanza/k))/sin(lat1)
(-cos(lat1)*cos(-long2+long1)*cos(distanza/k)+(cos(lat1)^2*cos(-long2+long1)^2*cos(distanza/k)^2+sin(distanza/k)^2-cos(lat1)^2-cos(lat1)^2*sin(-long2+long1)^2*sin(distanza/k)^2+cos(lat1)^4*sin(-long2+long1)^2)^(1/2))/(-1+cos(lat1)^2*sin(-long2+long1)^2))
atan((cos(lat1)*(cos(lat1)*cos(-long2+long1)*cos(distanza/k)+(cos(lat1)^2*cos(-long2+long1)^2*cos(distanza/k)^2+sin(distanza/k)^2-cos(lat1)^2-cos(lat1)^2*sin(-long2+long1)^2*sin(distanza/k)^2+cos(lat1)^4*sin(-long2+long1)^2)^(1/2))/(-1+cos(lat1)^2*sin(-long2+long1)^2)*cos(-long2+long1)+cos(distanza/k))/sin(lat1)
-(cos(lat1)*cos(-long2+long1)*cos(distanza/k)+(cos(lat1)^2*cos(-long2+long1)^2*cos(distanza/k)^2+sin(distanza/k)^2-cos(lat1)^2-cos(lat1)^2*sin(-long2+long1)^2*sin(distanza/k)^2+cos(lat1)^4*sin(-long2+long1)^2)^(1/2))/(-1+cos(lat1)^2*sin(-long2+long1)^2))
E questo per long2:
>> syms k distanza lat1 lat2 long1 x
>> eq= k * acos( sin(lat1) * sin(lat2) + cos(lat1) * cos(lat2) * cos(x - long1))-distanza
eq =
k*acos(sin(lat1)*sin(lat2)+cos(lat1)*cos(lat2)*cos (x-long1))-distanza
>> solve(eq)
ans =
long1+acos((-sin(lat1)*sin(lat2)+cos(distanza/k))/cos(lat1)/cos(lat2))
Perdona la rincoglionitaggine ma k ?
Twitter @veke87
Hmm long2 non sembra funzionare ...
Twitter @veke87
Download Derive
Questo potrebbe esserti moooooolto d'aiuto; il "problema" è che ci sono tutte le condizioni da mettere, che poi andranno a escludere certi risultati.
Per usare questo programma, all'inizio imposta subito:
InputMode := Word
(così prende le "parole" come variabili, non singole letetre)
scrivi l'equazione nella barra in basso
(es:k·ACOS(SIN(lat1)·SIN(lat2) + COS(lat1)·COS(lat2)·COS(long2 - long1)) - distanza = 0)
La selezioni, vai nel menu "Solve" -> Expression -> selezioni la variabile che ti interessa e clicki su "solve", così ti dà i risultati condizionati
In effetti se ci pensi bene l'arcocoseno assume valori da 0 a Pigreco quindi da 0 a 3.14 circa...
La longitudine va da -180 a 180 quindi non va bene quella funzione.. Correggetemi se sbaglio.
Twitter @veke87