Newsgroup: it.scienza.astronomia.amatoriale
Titolo: Calcolo del (JD)algoritmi e differenza SOLEX10.1L

Qualche giorno fa' mentre controllavo l'algoritmo per il calcolo del JD (Julian date) 
o giorno giuliano, con questo algoritmo.
Formula excel .
=(TRONCA(365,25*(Q50+4716))+TRONCA(30,6*(S50+1))+B6+T50-1524,5+((B9/24)+(B10/1440)+(B11/(24*3600))))
variabili di riporto :
Q50=SE(B7>2;B8;B8-1)
S50=SE(B7>2;B7;B7+12)
T50 = SE((B8+B7/100+B6/10000)<1582,1015;0;2-TRONCA(Q50/100)+TRONCA(TRONCA(Q50/100)/4))
Variabili
B6=gg  |  B7=mese  |  B8=anno  |  B9= ora  |  B10=minuti  | B11=secondi
==========================================================================
Controllavo per curiosita' andando indietro per raffrontare la posizione del Sole e 
della luna, con quelli calcolati da me' e ottenuti con JMeeus e Chapront, dove sono 
riuscito ad ottenere un ottimo risultato di precisione, nel range (-4800ac. +6500bc. 
(foto : 
http://yfrog.com/0oscostamentij
 ) high accuracy.
Ora per esempio controllando con questo algoritmo la data del  01/01/-10.000ac ho 
JD= -1931441,5 e un calcolo del Tempo effemeridi di (TE)= -118,999630390144
TE=((JD-2415020)/36525)
Inizialmente imposto la data sul Solex10.l -10000.0101 , poi entro nel programma e 
visualizzavo questi
TT:-10000/01/01: 00:00:00
JD: -1931442,5  Saturday
step= 1,000   d
Ora con l'agoritmo che ho usato in excel ho :
Algor.Sopra: JD= -1931441,5 domenica
Progr.Solex: JD= -1931442,5 sabato
Scarto 1 g !
Potrei sapere gentilmente quale e' l'algoritmo usato sul Solex??? [per intero e non 
frammenti... :)  ]
Dove e' l'errore ???
So' che il calcolo del giorno giuliano parte dalla data 1 Gennaio 4713 ac. e quindi 
per le date indietro non sarebbe corretto, ma per le effemeridi usare il TE 
(ephemerides time) e' fondamentale .
Ora gentilmente a cui dirotto anche la stessa domanda al Prof.Vitagliano nel caso 
legga e saluto cordialmente, sarebbe possibile sapere come mai tale errore ??
Qualcuno lo sa' ??
Grazie e Saluti.

 Grazie e Saluti.
Dal Manuale in dotazione:
About Time reference, the formulas (*New10.1*) used for computing
DeltaT (DT-UT) are taken From Espenak & Meeus (2005, Canon of Solar
Eclipses) , see 
http://eclipse.gsfc.nasa.gov/SEcat5/deltatpoly.html
A Wikipedia e da tanti mesi è stata inserita nelle "cose da fare"
quella dedicata alla stesura della voce: Solex chissà che ne pensa il
padre :-) . Nessuno per ora ha scritto nulla, ahimè :-/

"Eliseo" ha scritto nel messaggio

Dal Manuale in dotazione:

A Wikipedia e da tanti mesi � stata inserita nelle "cose da fare"
quella dedicata alla stesura della voce: Solex chiss� che ne pensa il
padre :-) . Nessuno per ora ha scritto nulla, ahim� :-/
========================================================
Grazie ma conosco gia' questo link e' per il calcolo del deltaT, io parlo da come vedi 
del Julian date (JD) ovvero calcolo del giorno giuliano corretto e per ora ho solo 
visto un errore, non so' se e del mio algoritmo, che per'altro ho gia' controllato in 
giro e sembra corretto, mentre vedo discordanza sul programma Solex10.1L .
Se venissero messi anche in linea e in una pagina anche gli algoritmi matematici, le 
persone imparano almeno qualcosa e vedono lo sviluppo di quello che fa' poi un 
programma compilato !
Cosi' si vede solo il valore calcolato e la gente si affida a quello senza capire come 
voglio anche io i passi e capire perche' deve venire anche un certo valore !
E' mai possibile che si deve solo prendere un listato, compilarlo e eseguirlo per 
sentirsi contenti ???
Come del resto non sono il solo che vuole vedere come avviene un calcolo e non mi 
sembra tanto anche chiedere di IMPARARE come avviene.
Pertanto mi sembra anche lecito chiedere della discordanza.
Del resto i fondamentali calcoli non devono essere rimandati a manuali poi da pagare, 
perche' sarebbe il colmo anche quello !
Cribbio e' mai possibile che l'astronomia e la meccanica celeste la si deve imparare 
da un eseguibile e non la si deve imparare dalle fondamenta ???
Come i relativi passi per calcolare AR e declinazione del Sole ecc...
Mi piacerebbe inserire qui i passi affinche molti sappiano almeno come si calcolano le 
posizioni ... per il resto attendo risposta alla mia domanda iniziale.
Ad esempio per spiegarti il grafico, vedi a questo link i calcoli che ho fatto e i 
termini addittivi per la luna..
http://img193.yfrog.com/img193/6233/calcolosoleluna.jpg
  (con IE8 ingrandire a 125% 
per visualizzare meglio )
Questo e' quello che intendo per calcolo .
Cordiali Saluti Eliseo.
Marco 

Possibile che nessuno sia in grado di dare una risposta attinente al problema in cui 
ho notato ???
Differenza jd sul solex e rispetto al valore che ho io, (algoritmo postato)
Cordiali Saluti e cieli serenissimi. 


Scusami ma forse sono io a non capire...
La risposta te la potrebbe dare solo chi ha libero accesso al codice
sorgente del Solex, il prof Aldo Vitagliano.
Ma c'è un'altra situazione sulla quale occorre "fare mente locale",
così riassumibile:
"Dove non esistono osservazioni, non vi possono essere valori certi
del DeltaT." Esistono formule che si basano sull'accelerazione
secolare della Luna, ma sono stime di essa più o meno accurate ma
sempre stime. Con le stime si fanno le estrapolazioni che sono tanto
più "rischiose" a mano a mano che ci allontaniamo dal periodo in cui
abbiamo osservazioni pratiche. il valore del DeltaT, piaccia o no, lo
si ricava solo da osservazioni astronomiche. Potrei scansionarti un
testo di circa 40 anni fa, che prevedeva per questo anno 2009 un
DeltaT nell'ordine di +80 secondi. e sono solo 40 anni fa. Immàginati
l'incertezza quando si parla di migliaia o svariate decine di migliaia
di anni dall'epoca attuale. Gli scarti potrebbero assommare a ore,
forse a giorni!
Eliseo

 Eliseo
Da quello che mi pare di aver letto, la sua domanda non riguarda
minimamente il delta-T ma il calcolo del Giorno Giuliano partendo dalla
data del calendario e relativa "inversa".
Inoltre, se non ricordo male, sul sito di Vitagliano c'� la possibilit� di
scaricare il "sorgente" basic della prima versione del programma (il Solex
1) e quindi, se l'autore non ha modificato quella parte del programma,
esaminare gli algoritmi utilizzati per le conversioni Giorno Giuliano <-->
data.
D'altra parte non sappiamo neppure quali sono i limiti di utilizzo degli
algoritmi utilizzati "Marco", quindi � possibile che gli uni e/o gli altri
non siano utilizzati correttamente.
La prima cosa che viene in mente � che uno dei suddetti algoritmi "mal
digerisca" Giorni Giuliani (G.G.) negativi; comunque non sono sicuro che
la giusta sequenza dei G.G. sia
-3, -2, -1, 0, 1, 2, 3, ecc. ecc.
e non invece
-3, -2, -1, 1, 2, 3, ecc. ecc.
Una volta appurata quale sia quella corretta, basta fare un po' di
conversioni di data <--> G.G. attorno allo "zero" per verificare chi
sbaglia e chi no!
Zeta Cancri
P.S. Vogliamo poi conoscere i risultati delle verifiche!
-- 
questo articolo e` stato inviato via web dal servizio gratuito 
http://www.newsland.it/news
 segnala gli abusi ad 
ab...@newsland.it

Dimenticavo!
Se dal punto di vista "matematico" pu� anche essere corretto, da un punto
di vista "logico" non � che abbia molto senso calcolare l'esatta data
corrispondente a valori negativi del Giorno Giuliano: al quell'epoca il
calendario (quello Giuliano, visto che quello Gregoriano � entrato in
vigore solo pochi secoli fa!) non esisteva, tanto meno le settimane e i
giorni della settimana...
Zeta Cancri


> On 11 Dic, 16:00, "Marco"

Gentile Zeta .C
Si' e esatto mi riferisco alla discordanza del calcolo del giorno giuliano, che come 
specificato ovviamente digerisce male i numeri prima del 4713ac, come specificato.
Il codice dell'algoritmo che ho usato e' nel primo messaggio postato, a disposizione 
di tutti.
Se ti interessa quello che hai visto ti posto appena terminato alcune modifiche, il 
foglio calcolo completo e disponibile, scaricato tempo fa' dal link postato da un 
certo gioacchino, un ottimo lavoro ha fatto, e che porto avanti anche io come 
progetto, ho trovato una buona precisione che ho migliorato. In cui corrispondenti a 
molti esempi del JMeeus Astronomical Algorithms ecc...
Sempre ringraziando potresti fornirmi il link del listato ??? cosi' lo scarico e trovo 
il problema ???
Ti ringrazio cordialmente e ti porgo i miei cordiali Saluti.
Marco.


Si' ma ha senso se mi serve calcolare il Tempo effemeridi, altrimenti acquisti un 
programma che lo fa', allora mi toglierebbe il gusto di fare alcuni calcoli in excel.
Sicuramente saprai che per risolvere l'algoritmo, devi avere il T  tempo effemeridi 
per tali formule ...
Longitudine Media geometrica Sole
=(RESTO(279,6966778+(36000,76892*T)+(0,0003025*T^2);360))
Tempo effemeridi =
((JD0 + 2415020,0)/36525)
JD0 = JD + Correzione deltaT (dal delta T passi al discorso dell'esimio Eliseo.)
Tanto per gradirti !
Cordiali Saluti.
Marco.

 Marco.
Prego...
Zeta Cancri
