Programma per Calcolo Effemeridi e altri dati.

[lucianob]

Ciao @CristiaPi
-Concordo punto 1
-Concordo punto 2
.. concordossssssss

Grazie per il ragionamento su l grafico e le DE430 .
Io non ci capisco tanto di statistica e non ho fatto l'universita' o il liceo, la mia matematica e la conoscenza me la sono fatta leggendo, sappi questo, purtroppo non ho molto tempo, ma quello che dici mi interessa perché va a riempire le mie lacune statistiche, avanzate di universita', pensa che ho problemi con i vari matriciali integrali e via dicendo .
Quindi ............................

Visto che tu hai utilizzato le DE 430, io ho utilizzato VSOP87 il calcolo per la determinazione avviene con la formula del quale indicata in "1988.03_Planetary_VSOP87.Solutions.pdf" pag.310 Rif.form.(3)

A*100000000*cos[(B)+(C * JDE) che poi sarebbe tradotta...
\alphaT a¨1 *A*Cos(B+CT)-T \alpha A*C sin(B+CT)
(Le formule non mi riescono diamine)
VSOP87B.foglio
--------------------------------------------------------------------

Premetto che
1) I file sono per tutti chi ne avesse voglia.
2) non sono un laureato e quindi, Statistica, Integrali, Sommatorie matriciali o terminali non me ne intendo. Indi non ditemi "vai a vedere", debbo andare an studiare interi libri e non me la sento piu'. Serve una mano per la precisione. !
3) Non conosco VB o C++ anche se da giovane mi ci ero messo ma forse pensavo a qualche donzella di troppo.
4) Riesci a tirare giu' i termini per un calcolo preciso visto che parli di statistica e termini che a me' sono ostici,
ad esempio formule come quella che ho risolto per la VSOP87B

Ma invece prendendo quella che mi hai parlato e per il quale sei riuscito a fare determinazioni ?
Vedo che e' dal MJD50000 ovvero circa ~10/10/1995 al ~23/12/2014 MJD 57013 come mai hai scelto questo periodo, quando magari si prende es. 01/01/1900 al 01/01/1901 , 365g. di dati ?

Visto che tu ci sei riuscito a decifrare la Teoria, senza documentazione, sono sicuro che magari hai trovato forse qualcosa, ma condividere sarebbe utile.
Mi risolvi il problema ? delle line di risoluzione per Longitudine, Latitudine, Raggio vettore ? ,. perché sono quelli che a me' servono ! Vorrei utilizzarli sul foglio di calcolo ?

Non capisco una cosa perdonami, sulla tua proporzione ricordo ancora qualcosa
il prodotto dei medi è uguale al prodotto degli estremi.
a memoria mi risulta che . Es.. 10:x=2:4 in questo caso x =0,543 ! qui sbagli qualcosa
per la soluzione ricordo che (10*4)/2 = il tuo x !
Lo inserisci tu ! Non e' calcolato.

Se tu vuoi sapere a cosa corrisponde 0,543, basta fare la proporzione:
0,543 : 1 = x : 27,1783 ottenendo
x= 27,1783 * 0,543 / 1 = 14,7578 giorni, cioè il semiasse maggiore lunare è soggetto ad una perturbazione (periodica) pari a 3400 km ogni 2 settimane circa (cioè 14,7578 giorni). A cosa può essere dovuta questa perturbazione così grande? Sicuramente al fatto che la Luna ogni mezza orbita si trova in opposizione e in congiunzione col Sole, che contribuisce per il 50% all'accelerazione a cui la Luna è sottoposta (6 mm/s2), mentre l'altro 50% è dovuto alla Terra.

Il 27,1783 che mi ha riferito, il Periodo Orbitale Medio, ma quanti solo 225 gg..?.
Ti spiego meglio . considerando che La luna gira mantenendo la sua stessa faccia in maniera Geocentrica, nel modello reale e Tolemaico.

Il mese sidereo è l'intervallo tra due successivi allineamenti della Luna e di una stella su un meridiano e dura 27 giorni, 7 ore, 43 minuti e 11 secondi.
Questo valore indica il tempo reale di rivoluzione che corrisponde anche al periodo vero di rotazione. Il mese sinodico o lunare è l'intervallo tra due allineamenti uguali Luna-Sole (opposizioni o congiunzioni) e dura 29 giorni, 12 ore, 44 minuti e 3 secondi. Le differenze tra mese sidereo e sinodico derivano dal fatto che mentre la Luna ruota intorno alla Terra, questa si è spostata per il suo moto di rivoluzione intorno al Sole. Per ritrovare l'allineamento con Terra e Sole, è necessario che la Luna compia un ulteriore arco intorno alla Terra, corrispondente all'arco di circa 27° che la Terra ha compito in 27 giorni. Per percorrere questo arco la Luna impiega circa 2 giorni, per questo il mese sinodico è più lungo di 2 giorni.

Mi passi l'algoritmo di come e' calcolato ?????? esempio, quantità dei termini, periodo, Quanto era il semiasse maggiore orbitale, la distanza media ? Mi servono gli algoritmi ... E capire da dove sono stati trovati !
mi c'e un errore di fondo.

Poi i 6 mm/s².. di accelerazione ?????
Io ho bisogno di dati e calcolo !!!!! Come lo calcoli il valore del semiasse cosi' a spanne ?
Le integrazioni si' ma condizionate a dati reali non "jx" o -1y esempio
Cioe' cosi' non va bene.
3400 km ogni 2 settimane ?.. Christian esiste anche la librazione lunare ?. Il Movimento lunare e' uno dei più ostici da calcolare, proprio per i suoi correttivi che sono reali e sono dei dati ! Non polinomiali, questo sono per luna.
Poi non parliamo di Plutone, che e' perturbato, ma molto variabile in orbita, per'altro non ci sono ancora in giro i dati polinomiali, ma probabilmente debbono ancora essere elaborati, ma aspettiamo.

Una curiosità il libro di Jean Meeus 1 e 2 li hai aperti ? Hai visto le tabelle ? .
Quelli sulle sommatorie e sulle correzioni lunari in Longitudine, in Raggio Vettore, in Parallasse ????, addirittura ci sono delle tabelle polinomiali per il calcolo reale della velocità, che io invece uso a rigore di riferimento quello un po' differente (es. nel foglio calcolo)

Ti ripeto ho bisogno di riferimenti e formule reali che possa capire, perché siamo su un terreno minato ! Ma tanto, fatto di libri dove ho preso dati.

Se ti interessa e sono sicuro visto che mi serve lo sviluppo di DE430, (tra l'altro mi sorprende perche' deriva dalla INPOP10a, del quale non ne ero entusiasta.... )come vedi a pag.4 del Pdf Top2013, trovi, termini come [ a, lambda, k, h, q. p. L. B. R ]

Ti ho caricato sul server dove ho il sito questi file che ti consiglio di dare un'occhiata .
Che mi sono serviti per i calcoli della posizione dei pianeti !, Che vedi come ho inserito ben 35.080 termini della teoria ! Usando i temini " VSOP87B" non la vers. A o C. per un motivo che esula dallo spiegare qui al momento. Abbastanza lungo e articolato.
1988.03_Planetary_VSOP87.Solutions.pdf
2002.11_Pos.Vel_VSOP87.pdf
-------------------------------------------------------------------------------------

I file nuovi che cercavi e che so dove trovare.
http://astronomiaexcel.superweb.ws/dwnl/Top2013.pdf
http://astronomiaexcel.superweb.ws/d...P2013_file.pdf

Vorrei inserire i 70.000 mi pare della nuova, tu usi i dati della DE430, quindi hai trovato un modo.
Facciamo un passo alla volta mi debbo sincronizzare anche io con i tuoi dati, la testa tra le nuvole per problemi del foglio, poi i periodi orbitali lunari.
Fai attenzione a non confondere i vari periodi o perturbazioni, o le varie cose, e' un vespaio e aprendolo alla larga ci si punge e ci si fa male.
Comunque posso dirti complimenti ma gli algoritmi e i dati a me' servono come il burro di arachidi su alcuni piatti .

Rimango in attesa e starò assente alcuni giorni .
Cordiali.

[CristiaPi]

Scusami @lucianob, ma mi sto un po' perdendo tra tutte quelle domande. Sicuramente non ci capiamo bene perché a te serve qualcosa di diverso da quello che ti ho proposto, ma non c'è nessun problema; ho solo provato a proporti una cosa che però stona con il resto del tuo file perché in realtà non serve a calcolare la posizione dei corpi celesti; è un'aggiunta molto marginale. Come non detto!
Se posso aiutarti tra teorie planetarie, DE/LE, ecc lo faccio volentieri.

Buona pedalata!

[lucianob]

Tranquillo, non preoccuparti ci si puo' aiutare..
Perché purtroppo per mettere dei dati debbo andare a leggere come questi dati sono presi e che veridicità hanno.
Apprezzo e sono contento che esistano " " Soluzioni " " . anzi lo so'.
tante Domande ho che aspettano risposte...

Le variazioni e i grafici sono veramente interessato, intanto prima sistemo il codice VBA che automatizza e poi passo ai grafici e alle problematiche, contando di mettermi a studiare qualcosa su di essi !
Rimane indietro la VSOP2013/TOP2013, i cui una parte di termini sono gia' inseriti in calcolo in ampliamento nel foglio !
Sono in grado di calcolare
e=V(h²+k²)
h=e*sen (w + W)
k=e *cos (w + W)
λ
p=tan (i/2)*sen(W) -- i* sen (W)
Q= tan (i/2)*cos(W) -- i* cos (W)

Ma a prendere i dati, farne l'algoritmo per poter avere un risultato finale, come minimo mi verrà un mal di testa, un infarto e ahimé, tanto nervoso ...

Per grafici seri intendo che non mi fa' dormire tranquillo quello esposto al post #51 ultimo link !
Ma sto lavorando anche per il Sorgere e tramonto dei pianeti che ho visto fattibile, da un'altro foglio che mi e' stato passato da un americano, a cui ho aperto il codice del mio. Ma solo in questi casi, condividi bene che io do' il massimo .

Saluti che mi bevo una bella bevanda redBull che mi aiuta...

PS. Sono in ferie fino al 27/8 come diciamo ... a Speruma.-.

[CristiaPi]

Tranquillo, non preoccuparti ci si puo' aiutare..
Perché purtroppo per mettere dei dati debbo andare a leggere come questi dati sono presi e che veridicità hanno.

Sacrosanto! E' giustissimo accertarsi di ciò che si sta usando.

Ma a prendere i dati, farne l'algoritmo per poter avere un risultato finale, come minimo mi verrà un mal di testa, un infarto e ahimé, tanto nervoso ...

Aspetta, aspetta! Se devi rischiare un infarto, ti prometto di scriverti il codice VBA per le VSOP2013. Sarà un lavorone, perché prima dovrò scriverlo in C++, poi lo dovrò tradurre in VBA e verificare...

Saluti che mi bevo una bella bevanda redBull che mi aiuta...

Salute!

PS. Sono in ferie fino al 27/8 come diciamo ... a Speruma.-.

Continuerai a frequentare il forum in questo periodo?

Ciao
Cristiano

[CristiaPi]

Ciao @lucianob,

Ho finito di scrivere e testare il codice C++; riporto i risultati per Plutone e per EMB, così puoi verificarli.
Ho riportato solo la prima riga del file usato per la verifica del risultato:
ftp://ftp.imcce.fr/pub/ephem/planets...n/VSOP2013.ctl
perché è quella che ricavi dalle VSOP2013; gli altri valori li ottieni dalla prima riga.

Ho anche indicato con t_calc il tempo di calcolo in secondi; con Excel dovrai moltiplicarlo non di poco. Un altro problema che potrebbe venir fuori con Excel è relativo all'allocazione di memoria, sono necessari circa 30 MB per ogni corpo caricato.

Risultati:

EMB

t_calc= 0.543, Julian Date JD= 2411545.000000 (TDB)
1.00000963581671 4.81887776423271 -0.00362423603660404 0.0163481299890599 0.000124873024391335 -1.10805144206119e-005

t_calc= 0.537, Julian Date JD= 2415545.000000 (TDB)
1.00000514351976 4.51233656382664 -0.00363537003815904 0.016361939190527 0.000112741225225195 -1.00803820155472e-005

t_calc= 0.535, Julian Date JD= 2419545.000000 (TDB)
1.00000737602512 4.20581952066578 -0.00366953954432214 0.016349557707845 9.87435814643221e-005 -8.3659242103311e-006

t_calc= 0.537, Julian Date JD= 2423545.000000 (TDB)
1.00001520076846 3.89926266200131 -0.0037015196550015 0.0163336036318818 8.66253683305102e-005 -7.76820051593722e-006

t_calc= 0.537, Julian Date JD= 2427545.000000 (TDB)
1.00001983068475 3.59270693959904 -0.00369726367347217 0.0162977176757242 7.44490107197197e-005 -5.60240767508277e-006

t_calc= 0.538, Julian Date JD= 2431545.000000 (TDB)
1.00001148285321 3.2861200051585 -0.00368798180900027 0.0163256237200869 6.20847477877878e-005 -4.66827974499357e-006

t_calc= 0.537, Julian Date JD= 2435545.000000 (TDB)
1.0000163384048 2.97957797278067 -0.00372887243956091 0.0163382014578303 5.00080734490766e-005 -4.68056689175179e-006

t_calc= 0.534, Julian Date JD= 2439545.000000 (TDB)
1.00000339818452 2.67302769646763 -0.00372124730759786 0.0162970035684459 3.79132023264989e-005 -3.21353633729501e-006

t_calc= 0.539, Julian Date JD= 2443545.000000 (TDB)
0.999991568903723 2.36648455224777 -0.00368988198199365 0.0162713070394488 2.4477419590655e-005 -2.50452148461629e-006

t_calc= 0.535, Julian Date JD= 2447545.000000 (TDB)
0.999991821705856 2.05993749983886 -0.0036985132391304 0.0162829234456087 1.20454896041833e-005 -1.88442538322343e-006

t_calc= 0.531, Julian Date JD= 2451545.000000 (TDB)
0.999996427318172 1.75341273413234 -0.00373390341672845 0.0162796458462242 -6.05480558893857e-007 6.26788324827015e-007


---------------------------------------------------------------------------------------

PLUTONE

t_calc= 0.209, Julian Date JD= 2411545.000000 (TDB)
39.4227219158885 1.39102609612849 -0.177750842251796 -0.175166170874624 -0.0515906602811809 0.140092456459022

t_calc= 0.204, Julian Date JD= 2415545.000000 (TDB)
39.3129146523696 1.67019970408671 -0.178361242924535 -0.177992462587011 -0.051622819362119 0.140264629713917

t_calc= 0.203, Julian Date JD= 2419545.000000 (TDB)
39.3965378946634 1.95130318491159 -0.179005270247882 -0.175863087279178 -0.0517317679335813 0.140037026721951

t_calc= 0.204, Julian Date JD= 2423545.000000 (TDB)
39.4657239282022 2.23341919437779 -0.180710328746584 -0.175325853257892 -0.0517880153175245 0.139910582730333

t_calc= 0.204, Julian Date JD= 2427545.000000 (TDB)
39.5593694964467 2.50968470762029 -0.181221132092434 -0.173462738880123 -0.0517677830883876 0.139798620905566

t_calc= 0.203, Julian Date JD= 2431545.000000 (TDB)
39.7847670022264 2.78480463419741 -0.184594183459701 -0.169997266272932 -0.0515878838831906 0.139565740495477

t_calc= 0.203, Julian Date JD= 2435545.000000 (TDB)
39.8105744753874 3.06040192388532 -0.18528774873777 -0.170873669134016 -0.0515453634713328 0.139620115170303

t_calc= 0.204, Julian Date JD= 2439545.000000 (TDB)
39.6397923114637 3.33406831955637 -0.182272089159026 -0.171173464871989 -0.051446597195315 0.139712587060577

t_calc= 0.205, Julian Date JD= 2443545.000000 (TDB)
39.5355974283669 3.61058648361179 -0.180834532811945 -0.171194321899969 -0.0514051941266498 0.139836233117331

t_calc= 0.203, Julian Date JD= 2447545.000000 (TDB)
39.3577302458111 3.89185496356853 -0.177361562383721 -0.171212877558505 -0.0516130546695339 0.139859403110761

t_calc= 0.201, Julian Date JD= 2451545.000000 (TDB)
39.264854264792 4.17260457759312 -0.175864116726266 -0.170123414288863 -0.0517015913689731 0.139865451440449

[lucianob]

Ciao @CristiaPi
Grazie ma aspetta .. che la riga di codice a questo indirizzo
VSOP2013 PLUTO Code.
Ti spiego quello che conosco -

Allora EMB per ora non serve per tutta l'elaborazione, perché non mi interessa inoltre non trovo i polinomiali di calcolo e neppure sono in grado di calcolarli senza dati.
Poi andrei a complicarmi la vita tra coordinate Topocentriche, Geocentriche e heliocentriche e sinceramente faccio fatica con le traslazioni, anche perché quelle Planetarie sono diverse dai libri di scuola e universitari .

Poi , apri i file per cortesia. VSOP2013.pdf e Top2013.pdf, nel link che ti avevo indicato.

"The terms of series are given under the form => T(alfa) * (S*sen(phi) + C* cos(phi) )
In pratica di darebbe già la soluzione, ma ahimè .-.... non ci sono le serie come nella VSOP87B (S, C e altre 2) cioè i Coefficienti ..
Il T e' il tempo effemeridi, diverso dai secondi che esprime il TDT, che e' quello del foglio di calcolo TE/TDT/TT
.Usato per il calcolo dei termini di serie. Il valore che hai messo tu non serve a nulla mi dispiace.
Ti ho già scritto che si calcola prima con, TE/TDT/ =[( JD -2415020 ) / 36525] qua' senza le minime correzioni dovute a UT1, UT2 e amenità piccole variazioni dovuti a motivi planetari.

TI DICE ANCHE
In a VSOP2013 file, series are put in order of the variables <a,,k,h,q,p>
and, for each variable, in order of the time power ().

Poi i dati nel file LINEA
========================================
Julian Date JD 2411545.0 (TDB) => sai che significa ?
Giorno Giuliano Tempo Baricentro terrestre in data > in data 26/06/1890 alle ore 12TU (tempo universale)
39.4227219159 1.3910260961 -0.1777508423 -0.1751661709 -0.0515906603 0.1400924565
17.6212944088 43.6567983387 -9.7715435955 -0.0020295942 0.0006425559 0.0005199070
17.6213054610 43.9412232438 8.4004533071 -0.0020295941 0.0003827270 0.0007325993

Ci sono ben 18 serie di dati, li conosci ??????????
Perché io riconosco in ordine; a, l, k, h, q, p (sopra)
Ma sotto a malapena capisco 17,6 che e' l'inclinazione variabile, gli altri per me' sono arabo !
Il file dice anche .:
Each VSOP2013 file corresponds to a planet and contains trigonometric
series, functions of time (Periodic series and Poisson series), that
represent the 6 elliptic elements of the planet:

contengono elementi trigonometrici della serie, nella funzione di tempo nella serie di Poisson, questi rappresentano 6 elementi ellitici dei pianeti !
Ma il resto non e' specificato come vorrei capire.
Se qualcuno non lo scrive e riesco a prendere i termini per una sommatoria dei termini ok, altrimenti così a grandi spanne e inserendo un " t_calc= 0.209, " con tutta la tua buona volontà che apprezzo, ma e' sbagliato, cosi' dico non serve a nulla............. perdiamo tempo a scrivere post, a spiegare a cercare di capire io e faticare tu,
RINGRAZIANDOTI , ma devi conoscere come prendere i dati e come riconoscerli, anche se sei in gamba per le statistiche e con il C++ .

Questo il motivo mancando alcuni elementi chiari di calcolo e non potendo prendere il toro per le corna, devo abdicare dall'utilizzare una migliore precisione della VSOP2013 !


Tu hai messo .....................

*PLUTONE*
t_calc= 0.209, Julian Date JD= 2411545.000000 (TDB)
39.4227219158885 1.39102609612849 -0.177750842251796 -0.175166170874624 -0.0515906602811809 0.140092456459022

Che si riferiscono ad: a, l, k, h, q, p ! Se riuscito a estrapolarli, ma sono uguali a quelli che ci sono già !

I termini che si utilizzano sono i termini di Delaunay, e gli elementi besseliani, su questo argomento non c'e' poi tanto di Universitario spiegato dal JPL, cos' come le traslazioni RKS per le coordinate.
Neppure sul libro Explanatory_Supplementary_Almanac, o Astronomical Almanac 2017 and Companion .
A cui ti invito caldamente di leggerli ! Trovi qualcosa su Google Books, che grazie a dio qualcosa leggi, ma molte pagine sono in copyright
Per capire gli elementi !

Mi spiace che leggendo il file VSOP2013, o un altro inerente, io in fortran non ci capisco una maz*, basic manco.
Poi per calcolare quel k, h, q, p, devono avere il valore di eccentricità che nel tuo estrapolare non esiste.
Quindi prima si dovrebbe prendere la formula e andare a ritroso a essa.

Io ho imparato che per calcolare la posizione di un pianeta, debbo avere dei valori che sono certi
(a)Semiasse maggiore (UA), (e) Eccentricità, poi i polinomiali dell'argomento del perielio, il calcolo della longitudine del nodo ascendente, il valore di Anomalia media, con i termini per poter partire correttamente.
Poi i POLINOMIALI che ti avevo gia' scritto !!!!!!!!!!!! Cose che normalmente vengono messi a diffusione per gli astronomi, ma vedo che negli ultimi anni una chiusura mentale, non solo cerebralmente lesa di quella parte di sistema che invece i dati te li mette a pagamento, o sei di una Universita' e allora puoi accederne a una parte, oppure i dati veri che vorresti avere non riesci a reperirli ! PUNTO !

E avevo anche scritto che sulla VSOP87B, ero riuscito a calcolare vedendo una formala di Sommatoria e con un po' di fatica sono riuscito a capire, grazie al libro di Jean Meeus Astronomical Algorithms, 1 e 2.
A riuscire a calcolare il tutto.
QUA' SIAMO IN UN MONDO DIVERSO !!!!!!

Qua si ragione partendo sempre da T calcolato come detto, le integrazioni non le conosco e Vedendo i file VSOP2013, che sinceramente e con tutto il cuore non mi ci sono ancora messo a vedere bene, dovrei cercare in rete se trovo qualcosa, ma guardando a grandi spanne, se vedi quel file e se vedi invece quelli che ti ho messo io .........
VSOP87 VERSION B2 MERCURY VARIABLE 1 (LBR) *T**0 1583 TERMS HELIOCENTRIC DYNAMICAL ECLIPTIC AND EQUINOX J2000
2110 2 1 0 0 0 0 0 0 0 0 0 0 0 0.09025255625 -0.39983457843 0.40989414977 1.48302034195 26087.90314157420
Qua almeno avevi dei coefficienti e come utilizzarli.

Dall'altra parte invece mi trovo Radianti, mi trovo valori in cui per utilizzarli devo avere Spice che gia' mi crea problemi, poi cercando spice mi trovo ad avere problemi quando lo trovo vedo solo delle utility .!
Insomma problemi su problemi e senza riuscire a risolvere come voglio io .

E gia' abbastanza complesso il foglio che oggi tanto per raccontartene una,
Inserisco data, inserisco selezionando un test su Latitudine oltre i 55°, mi trovo una serie di celle con errore " #num " errore, cerco di capire e mi viene in mente che se salgo ... in LAT
Inteso Inizio Crepuscolo Astronomico 18° Sun is 18° Below Horizon, leggo errore, ................. mentre devo dovrebbe essere visualizzato > " --:--:-- "
Poi salendo in lat aumentando vedo Crepuscolo nautico e poi il crepuscolo civile, Sun is 12° Below Horizon,
e mi chiedo ma come mai ........? dove' l'errore ?

A questo punto sto immettendo un "IF" "se" condizionale se salgo sul valore calcolato, deve visualizzarsi un assente . Quindi per farlo cerco un valore di latitudine decimale, per inserire un condizionale e verificare di proseguire se Latitudine scende.
Perché la durata del giorno e della notte durano quasi 6 mesi e variano in base a latitudine.
Cioe' non me ne ero mai reso conto, dell'errore. e mentre lo sto rifacendo ex novo, voglio trovare soluzioni e risolvere una volta per tutti i problemi !

Mipiange il cuore non inserire la VSOP2013 , ma ora dico ORA non ho tempo e per farlo o iniziare dobbiamo prima leggere dalla VSOP87, riaprire qualche libro, leggere anche tu i libri che ti ho menzionato, devi conoscere come si calcola TE, dal giorno giuliano, ATTENTO CRISTIAN, che in rete ci sono JD che non vanno oltre il 1900 indietro .
Devi prendere quello mio . Che arriva a oltre -12000 bc in JD ! Anche se a te' parrebbe inutile, per me' no'.
(Per'altro sono stato l'unico essere pirla, a dire al Prof. Vitagliano, che in una versione vecchia del Solex c'era un errore del tempo giuliano, (io non avevo detto che la routine l'avevo trovata su un sorgente in una pagina web da un calcolo del NOAA), questo chiesi e feci presente in un gruppo di astronomia sui newsgroup, circa 6 anni or, sono, mi rimando' a un file basic il primo iniziale e poi neppure un grazie. Alla fine la versione successiva era corretto .) !

Quindi non fare integrazioni se prima non sai i dati che debbono utilizzarsi, il tempo in astronomia e' uno solo parte dal TDB Tempo dinamico baricentro terrestre, in arXiv trovi una marea di documenti !

Ti mando un grosso Saluto
Ringraziandoti, al momento lascia tutto cosi', che per fare questo post ci ho perso quasi 90' per verificare i dati, leggere il file celermente, andare a controllare le linee sulla VSOP87 e andare a cercare qualche post precedente.

E per favore cerca quei libri e vedi come sono poste le cose.
Apri per favore e leggi.
Jean Meeus Astronomical Algorithms_1Ed e 2 Ed. (io l'ho in pdf)
Jean Meeus Formulae for Calculators_4ed (l'ho in pdf)
Explanatory Supplement Astronomical Almanac (lo trovi su Google books)
NUMERICAL INTEGRATION.pdf (del Prof. Vitagliano) (regalo sotto)

Numerical_integration.pdf
Aggiungo una info, l'integrazione usata da Vitagliano e' riportata su un articolo che lessi di una rivista astronomica per pochi, cui parla di Feymann il matematico a cui si ispiro' ! .. Vedi un po' tu, magari ce la facciamo ...


CIAO

[CristiaPi]

Ciao @lucianob,
Siccome ho perso più di mezza giornata per scrivere il codice C++ che utilizza le VSOP2013 (che a me non servono assolutamente a niente), voglio provare un ultimo tentativo, poi eventualmente chiudo qui la faccenda.

Ti chiedo: sei interessato ad usare le VSOP2013?

Se vuoi rispondermi, fallo, per cortesia, soltanto mettendo una croce sul sì oppure sul no, senza aggiungere altro:

[ ] SI [ ] NO.

[lucianob]

Ecco la foto del progetto spero possa uscire in Ottobre 2017, se un paio di cose vanno come dico, ma come sempre ci sono sempre ostentamenti e problemi da risolvere.
Altrimenti lo faccio slittare a Dicembre 2017 --------------

Diciamo che il 50% e' in codice VBA . Trasformato dal sottoscritto e routine VALIDE .
Non commerciali e non free, che sono riuscito ad elaborare. Che saranno protette.

Aiuti sono giunti inizialmente da @mbrio76 a cui ringrazio e un altro utente esterno, altri al momento non ricordo perche' sono gli aiuti come il vento, passa pulisce e non lo risenti piu...
Idee e spunti, ma tutto il resto sono andato avanti da solo imparando parte di quello che mi serve.
Il Codice VBA al momento consta di ben 49 Costanti, 186 funzioni diverse di calcolo !
Diciamo che ho fatto moltissimo, "come sempre da solo " !!!

Posto la foto che spero piaccia l'interfaccia piu' segnata e decisa, fatta gia' per una risoluzione oramai quasi di tutti in uso di 1920 x 1080 (io uso una GTX-970 OC/ed. nVidia ASUS come scheda Video) -
Ma posso assicurare che rispetto a prima in VBA ho al momento migliorato la velocità di calcolo, visto che vengono fatte una marea di operazioni nei vari fogli !
E nell'ultima con le nuove aggiunte si piantava quando cambiavo di netto la longitudine ...
Ora vedo che non lo fà piu'.
C'e' e ammetto una moltitudine di celle che superano le 40.000 !!!!

Spero e mi auguro ci siano persone che aiutino, altrimenti il progetto tra Febbraio / Maggio 2018 , arriva al culmine e chiude. Inoltre scadrà perche' e' temporizzato per via di aggiornamenti che debbono eseguirsi per alcune routine.
Ho voluto e spero di trovare la soluzione per l'esatta posizione del Polo nord e del suo spostamento, calcolo complesso tra i vati UTC/UT e variezzzz... parte inserito ma il calcolo complesso al momento no' .
Ci sarebbero variazioni di aggiunta e calcolo.

Al momento i diagrammi saranno validi e funzionali solamente dal 1900 / 3000 dc. Per via delle eccezioni limitazioni di excel che, non consente di andare entro il 1900 !!!!!
Perche' ho preso parte di alcuni calcoli aggiunti in questi 3 mesi, da un foglio di calcolo, pubblicato dal NOAA come dimostrazione che si puo' avere precisione !
Grazie a loro che hanno condiviso, ho potuto portare avanti qualcosa...

Comunque questa e' solo la foto di quello che e' ora, pulito, piu chiare le sezioni, ma sempre un portento di tantissimi dati .
Qualche errore corretto grazie a qualcuno che ogni tanto spunta al capoli

Godetevi le foto ! che penso siano piu' eloquenti di quanto dico,.
PS Sono stati aggiunti il Rise dei pianeti e l'angolo che spero qualcuno avrà modo di controllare, trovato il calcolo su un foglio che mi hanno mandato, come sempre condiviso da un americano !

Foto.Effemeridi2017.10-1 ------ Foto.Effemeridi2017.10-2

Questo e' tutto .
Spero sia pronto per il 1 Ottobre 2017 -
Un Saluto a voi tutti .

ATTENDO CONSIGLI E PROPOSTE GRAZIE !
30/08/2017

Parte iniziale codice VBA
--------------------------------------------------------------------

Attribute VB_Name = "astroDreaker"
'Astro-functions Dreaker Rev.2017.10
'Editied and added to by LucianoB (dreaker project 1997, VBA 2017, RA and time UT, and for all angles
'measured in degrees, all instants defined by the decimal fraction of days after J2000.0.
...................
.................
.... ... ' [01] Funzione trigonometrica - RADIANTI SENO di Angolo (RADIANTI)
' Return: Ritorna il coseno di un angolo. SIN(RAD) - =DegSIN(xx) SIN(radianti)
' ************************************************** ******** OK *****/
...

[mbrio76]

Gran lavoro @lucianob!

Sempre felice di poterti dare qualche dritta, compatibilmente con gli impegni prioritari

[lucianob]

Un Cordiale Saluto a tutta la Comunità di Astronomia.com
Un Cordiale Saluto alla zona Admin & Mod @Stefano Simoni, @etruscastro, @garmau, @Red Hanuman
Ai molti utenti che hanno commentato, messo “Grazie” ecc..
@alemin, @Alfiere, @Andrea86, @Anjiloh, @Beppe, @Borja, @Cagnaccio, @CristiaPi, @Enrico Corsaro, @fabpc, @Gaetano M., @givi, @marcolunardo, @mbrio76, @Morimondo, @ncrndr, @nicola66, @nullasfugge, @Pierluigi Panunzi, @primaluce, @Starfighter, @SVelo, @vega41

Questi solo una parte, mentre tutti quelli che non sono menzionati a cui va il resto dei saluti (ne ho ancora tanti), per gli utenti e coloro che seguono il progetto che, in questi 2 anni mi hanno scritto, sia dal forum che altri privatamente a cui hanno trovato nel foglio proprieta' pdf il mio email.
Arriviamo al punto :

PREVIEW Foglio di Calcolo Effemeridi2017.10, progetto che aspettava da circa 2 anni boni .
Ci sono molte novità sulla finale e aspetto molti di voi, di scaricare il file e dare un giudizio, una critica costruttiva, un parere (non le tavanate bestiali), un consiglio, appunti e quant'altro nel futuro del progetto Common copyright sotto Freeware e Open (in parte), di licenza, cui versa questo “Foglio di Calcolo Effemeridi in Excel” ............... Con il solo menzionare l'attribuzione stessa, ma per usi NON-COMMERCIALI, come previsto in licenza.
Come potete notare i miei allegati utilizzo Informativi che accompagnano il foglio sono non solo protetti, visto che servono solo alla lettura, ma hanno all’interno proprio I marchi di Licenza gratuita Open, Freeware espressi ! - Ma anche nelle proprietà del file trovate la mia email per scrivermi.
Per doverosa correttezza, che chi lo scarica non si trovi da pagare NULLA A NESSUNO !
Un Grazie vi sarei grato a tutti !

NOVITA’
Parecchie esempio come accennato, il progetto andrà verso il VBA . Al momento uso Office 2010, mah.. mi sa che, passerò al 2016, anche se molti hanno ancora il MS.off-excel 2010. Ma vedo alcune limitazioni nel codice VBA dal excel 2010 al 2016, nel senso che su uno funziona e sull’altro no’, perché si è evoluto, ho girato la rete e sui Newsgroup dove seguo excel, sono iscritto al Forum excel (invitato e ringrazio), sito pregiato.
Forumexcel.it, ma anche a Archi.forumup . per seguire i post e cercare di imparare il VBA, postando le mie esperienze e quanto concerne. Ma anche pesantissimo con tante cose belle che si possono fare. Senza dimenticare OpenOffice a cui invece ti apre il foglio di calcolo !
E a cui io invece penso che se il file e’ chiuso nel progetto dei miei dati che servono alla struttura portante, lasciando il resto aperto, ovvio che scelgo sia eseguibile il file e limitato ! Con le dovute protezioni.

Le migliorie le elenco all’uscita del file spero fine mese, primi Ottobre.!
Alla Apertura del File eseguibile, vi troverete un MessageBox – che vi avverte di
1) abilitare le Macro OBBLIGO PER AVERE iL FOGLIO FUNZIONANTE E, NON AVERE PROBLEMI .
2) Menu di excel ridotto all’osso e minimale, con note da leggere-
3) Troverete 4 Fogli: Abbreviazioni (Utili per sapere un po’ di cose a riguardo le nozioni), FotoMeccanica (foto di elementi astronomici importanti, utili e soprattutto fanno capire determinati concetti), NOTE-LEGGERE (che si attiva appena aprite l’eseguibile), infine “Foglio Effemeridi”
4) Dovete solo inserire, Giorno , Mese, Anno, Ora, Minuti, Secondi, Nuovo inserimento : TEMPERATURA Luogo, e Pressione atmosferica, tutto nelle celle bianche .
E poi Selezionate la città vicino alla vostra, click sopra si apre un menu a discesa e selezionate scorrendo la barra laterale, il foglio calcolerà la variazione !
Se la VS città vicina non c’e’ Richiedetemela, Scrivetemi, che l’aggiungo !!!!!!
Datemi riscontro perfavore !!!!!!!!!!!!

5) Inseriti i parametri il foglio calcola il tutto e vi fornisce i dati e dico tantissimi dati . Spero che alcuni di astronomia, che alla stessa non basta avere un telescopio perché l’ha il vicino, occorre anche avere elementi e sapere i dati che riguardano la Terra, la Luna, Il Sole, il Sistema Solare, altrimenti non avete nulla di concreto .
6) La semplicita’ e’ data per il fatto che, molti non sanno bene excel, qui invece basta inserire pochi dati
PS…………… OBBLIGO DI INFORMARMI a parte gli scherzi per eventuali errori, anche io sbaglio !
Nel contempo ricordo
I LINK

DOWNLOAD
Dentro i file zip trovate questi file sotto.
=> 01_Algoritmi_generali_2017.10.pdf
=> 02_Spiegazioni_Foglio_2017.10.pdf
=> 03_Riferimenti-Libri_v.2017.10.pdf
=> Astronomical_Constants_2017.pdf
=> eccentric.Terra.2mill.y86-99-2004.gif
=> obliq.Terra_laskar86-99.gif
=> tabella_costanti.pdf
=> Effemeridi2017.10.exe <= NOVITA’ in eseguibile ! testatelo per favore !

In arrivo appena completato e sistemato come dico io che sarà aggiunto
“ 04_Librerie.VBA_v.2017.docx “

LINK DIRETTO, nel caso abbiate ancora problemi
PRW-Effemeridi2017.10.rar
PRW-Effemeridi2017.10.zip
-------------
Altro sito
AS_PRW-Effemeridi2017.10.rar
AS_PRW-Effemeridi2017.10.zip
--------------
Altro sito
AE_PRW-Effemeridi2017.10.rar
AE_PRW-Effemeridi2017.10.zip

PRW sta per PREVIEW ! non per la finale.

GRAZIE A TUTTI E SPERIAMO TRA 2 SETTIMANE A FARE USCIRE IL RESTO. E Chiedo Riscontri.

PS. Tutti i file della vecchia versione sono cancellati, non piu disponibili !

PS DIMENTICAVO AGGIUNGERE, USARE EXCEL VERSIONE 14(2010) IN POI, altrimenti si possono generare errori . Stavamo verificando io e assieme a Borja per un errore.