Dráha planetky-nalezení-efemerida-amatérsky

[grajda]

Ahoj- Vesta by se každopádně hodila! Co se týče vícero měření a "proložení křivky", tak na upřesňování dráhy jsou odborníci Jindřich Žižka, pravděpodobně i KpS a vaclav1. Pokud máš na mysli závislost afa, delta na čase, asi by pomohl nhějaký polynom, nalezený buď přímo metodou nejmenších čtverců, nebo snad i šíleným rozvojem nějakých těch rovnic kolem- vždycky by u extrapolace ale byla omezená přesnost, zhoršující se s přibývajícím časem.

PS: Mám od Jindřicha Žižky bakalářku kolem ccd astrometrie, fláknul jsem to na naše stránky- zajímavá, ale nemám čas. Kromě toho poslal vaclav1 makra pro excel k výpočtu drah, jsou to různé metody, víc zatím nevím, páč jsem se snažil do začátku druhýho semestru dohnat svůj skluz v práci a osobním životě a teď začíná zase příprava do školy. Oběma děkuju maximálně za ochotu a vstřícný přístup. Tákže, až si najdu chvilku, zařadím soubory do wiki. Nebo to někdo udělejte za mě, přeposlat soubory je mnohem lehčí- napište na gracias(zav)ipm.cz, ale až v pondělí!!!

Bohužel, času moc nemám a ani nebudu mít, protože se ke škole tentokrát hodlám postavit zodpovědněji- máme známky A až F (F<=50%) a průměrná známka C za první semestr mě nevyhovuje. Kombinované studium však žádá příliš mnoho času, čím víc, tím lepší známka...

[KpS]

A napadlo mě, jestli by se těmi mými nepřesně změřenými polohami Vesty nedala proložit nějaká křivka, která by ty odchylky eliminovala.Z Vašich pěti měření Vesty jsou pro určení dráhy vhodná jen tři. Měření z 18.10. a 3.11. by dráhu zhoršila. Elementy určené ze zbývajících tří měření nejsou špatné, ale není čím je vylepšit. Předpověď z nich plynoucí pro 1.1.2009 se od reality liší o 13'. Jakousi představu o vztahu mezi přesností měření, přesností z nich získaných elementů a odtud plynoucích předpovědí může pro zde diskutovaný případ dát tato tabulka.
Prokládat měření by mohlo mít smysl v případě, že máme více časově blízkých měření, která bychom pak nahradili jedním fiktivním. Tomu podobný postup byl užitečný v době, kdy byla důležitá každá ušetřená aritmetická operace. Dnes je jednodušší najít přibližné elementy a k nim pak hledat opravy, které nejlépe vystihnou všechna měření, která máme k dispozici.

[grajda]

Na naše stránky jsem přidal program a makra od Vaclav1. Je to zatím neošetřeno, jak to uložit na počítač a rozjet je tam napsáno. Díky za to. Podle zběžnýho pohledu uno perfekto. Jen je tam ještě jeden prográmek, u kterého Vaclav1 neověřoval správnost výpočtu okamžiku průchodu periheliem, resp. střední anomálie- bacha na to.

Buď to někdo vytřiďte za mě, nebo to nechejte na mě, ale pak to může trvat i pár týdnů... [smiley=vrolijk_26.gif]

A přidávám další z úžasných nových smajlíků:

[smiley=vrolijk_1.gif] [smiley=tekst-toppie.gif] [smiley=2vrolijk_08.gif] [smiley=2vrolijk_08.gif] [smiley=2vrolijk_08.gif]

[vaclav1]

To grajda:
Už mám opraveno M a T v aplikaci Gauss pod Visual Studio 2005 v C# a je to poměrně přesné na to že počítám dráhu asteroidu Pallas ze 3. pozic a MPC pro asteroid Pallas má dráhu vypočtenu pro 5482 pozic a od roku 1839 do roku 1993.Jsou rozdíly v elementech +- 0.001 proti přesným hodnotám.
Orbital Elements
The following orbital elements are taken from MPC 24084:
    (2) Pallas
Epoch 2007 Oct. 27.0 TT = JDT 2454400.5  (M-c)    Bowell
M 242.48990  (2000.0) P   Q
n   0.21357240     Peri.  310.31268     -0.56753236     -0.82051565
a   2.7718686 Node   173.13654     +0.81429696     -0.54709512
e   0.2307662 Incl.   34.84076     -0.12176814     +0.16565326
P   4.61     H    4.13    G   0.11
From 5482 observations at 63 oppositions, 1839-1993, mean residual 0".55.

Výpočet pomocí vlastního software Gauss
Planetka n  q  a  e  M  T Incl.     Node     Peri.
2 Pallas 0,2135752 2,1240752 2,7718445 0,2336961 328,19135 2444833,91247 34,79457 173,34361 309,92517
Epoch TT = 2446370,56742 JD vztaženo k ekvinokciu 2000.0

Planetka P[x] P[y] P[z] Q[x] Q[y] Q[z] kontrolaP kontrolaQ
2 Pallas -0,5644612 0,8162306 -0,1230903 -0,8228052 -0,5444065 0,1631355 1 1

[grajda]

Díky za tu opravu, teď máme perfektní nástroj pro kontrolu vlastních výpočtů a úžasná věc na tom je, že je to amatérská práce. Dostal jsem email s opraveným prográmkem a k večeru ho budu moct uploadnout taktéž... Je hezký, že v okolí máte o tohle vlákno, nemající popravdě úroveň, kterou by si zasloužilo, zájem. Přes víkend na to snad taky mrknu. Ale budu rád, když stihnu opravit, co jsem napáchal a snad se i něco málo naučit...

[grajda]

Dneska přišel od vaclav1 email s dalším programem:

Dobrý den,

V příloze zasílám program Gauss pro výpočet předběžných elementů dráhy asteroidu.

Pro instalaci je potřeba rozbalit adresář a nakopírovat tak jak je do C:.

Je nutné mít na PC nainstalován doplněk Microsoft Net Framework verze 2.0 a vyšší.

V případě eliptické dráhy s excentricitou v intervalu 0.95 až 1.05 je použit rozvoj
pro elitické dráhy blízké parabole ze Subbotina - Úvod do teoretické astronomie.

V případě hyperboly s excentricitou > 1.05 je řešena Keplerova rovnice(použita
funkce hyperbolický sinus) .

V případě eliptické dráhy s excentrcitou < 0.95 řešena klasická Keplerova rovnice.

Můžete nyní vesele testovat.

Díky, až bude nálada na inteligentnější činnosti, než jen přispívání do jiných vláken, všechno uploadnu, popíšu a zařadím do našich stránek.

[Jindrich Zizka]

Dobry den,
dodelal jsem program v Matlabu na vypocty Gaussovou metodou a shodl se ve vysledcich s Vaclavem1. Ted pracuji na programu pro vypocet poruch, tedy evoluce drahovych elementu planetek v case, cili sekularni teorie. Pristupuji k nemu naprosto obecne a resim numericky patricnou soustavu pohybovych diferenialnich rovnic (lze postupovat i tzv. Cowellovou metodou). Pokud nekdo ma o skript v Matlabu zajem, necht me kontaktuje na adrese zizka@hvezdarna.cz. Je k tomu udelane i graficke rozhrani. Zdravi jindrich Zizka

[Jindrich Zizka]

Tak mam hotovy skript v Matlabu na urcovani elementu Gaussovou metodou vcetne grafickeho rozhrani. Poslal jsem vam to do mailu. jedna se o soubory gaussgrafika.m, gaussgrafika.fig a qfce.m. Otevrenim gaussgrafika.m a stisknutim F5 dojde k prelozeni a spusteni skriptu. Kdyby byly nejake potize obratte se na me mailem. Budu rad za otestovani dat, ale mela by byt v poradku, testoval jsem je z vice zdroju. zdravi Zizka

[grajda]

Dnes jsem měl až cca do 22h plno práce, proto skript rád pomůžu otestovat, ale nejdřív až zítra v práci.

Předpokládám použití jednotek rad, AU.

[KpS]

K výpočtu kruhové dráhy jsem nyní přidal popis výpočtu dráhy ze tří pozorování. Text je o něco delší, je to však způsobeno především nutností použít některé ze vztahů platných pro pohyb po kuželosečce. Tato komplikace by však neměla zastínit svou podstatou jednoduché myšlenky, na kterých je celý postup založen. Text je doplněn přikladem výpočtu elementů planetky Vesta z měření, která provedl Tulen a tabulkou v Excelu. Tyto tři soubory jsou součástí e3.zip.

[grajda]

Pánové... Díky, díky, díky. Vy mě zabijete... Chodí mě smršť emailů od vaclav1 a jindrich zizka, kolem nebmech, k tomu od mnoha jiných lidí kolem jiných věcí, nestíhám na to ani pořádně odpovídat. Kolega phobos je kdoví kde a já jsem slíbil kolegovi jindrich zizka otestovat jeho skript s gaussovou metodou výp. el. el. dráhy a ještě ani to nestihnul. Tohhle není o tom, že si chvilku sednu k AF a píšu ptákoviny. Z práce jsem se vrátil teprve před chvilkou (a často se vracím i později) a od začátku semestru jsem se ještě ani nepodíval na svoji vlastní školu, natož na nebmech. A potřebuji se taky občas z toho shonu vylízat nicneděláním. PROSÍM, zařaďte to někdo, máte- li chuť, na naše stránky za mě! Nejsu mašina! Díky.

[grajda]

Osobně jsem si stáhnul soubory na počítač, s uploadem to však nebude hnedle... Mám žalostně pomalý mobilní internet a málo času. Opět je to krátký text, nečetl jsem ho, ale protože je od KpS, tak ho můžu dopředu vřele doporučit.

[grajda]

Nemohl jsem to vydržet a text jsem si prolítnul až po výpočet elementů- to už je jen nějaká geometrie kolem. Vždycky jsem si myslel, že určit elementy eliptické dráhy je nesmírně složitý problém, ale KpS to dokázal tak srozumitelně a stručně podat, že mám chuť si započítat- "vždyť je to v základu tak prosté". Ale to bohužel nepůjde. KpS: Gigantická pochvala!

Tak proto se v ZNM píše o těch plochách trojúhelníků... : )

[grajda]

Neodolal jsem a uploadnul to, ale teď už se fakt musím věnovat práci.

http://astro.ceskyserver.cz/wiki/index. ... femeridy_p...

Viz. užitečný software, literatura, uskut. výpočty.

[vaclav1]

To grajda:
Velká omluva, trošku jsem se zdržel s programováním pohybových rovnic pro metodu Cowell(perturbace elementů přes souřadnice x, y, z a složky rychlostí x´, y´, z´).Diferenciální rovnice řeším metodou Runge - Kutta 5 řádu přesnosti.To jest, že pokud je krok intervalu 1 den tak přesnost určení souřadnic je řádu max. 10 na -5.Což asi bude plně vyhovovat pro období cca 1-2 roky pro přesnější výpočty deklinací a rektascenzí s přesností na zlomky vteřin a sekund.Jinak poruchy od velkých planet mohou hýbat se souřadnicemi asteroidu až v řádu 10 na -3! a chyby při výpočtu rektascenze mohou dosahovat i 5 sec a v deklinaci 30".