155TG2 Teoretická geodézie 2 / úloha 2: Porovnání verzí

Z GeoWikiCZ
mBez shrnutí editace
mBez shrnutí editace
Řádek 10: Řádek 10:
Podrobnější informace k převodu vektoru souřadnic z terestrického do nebeského systému zde budou ještě doplněny.
Podrobnější informace k převodu vektoru souřadnic z terestrického do nebeského systému zde budou ještě doplněny.
</span>
</span>
<!--
<!--
Na základě skutečných drah družic daných svými přesnými efemeridami v souborech formátu SP3, poskytovaných na serveru služby [https://cddis.nasa.gov/archive/gnss/products/ CDDIS] (pozor: nově nutná registrace) <!-- [ftp://cddis.gsfc.nasa.gov/gnss/products/ CDDIS] [http://igscb.jpl.nasa.gov/igscb/product IGS]-->, sledujte časový průběh oskulačních elementů. Použijte k tomu data pro den = den Vašeho narození + 15 let (resp. totožné datum, pouze rok se zvýší o 15, je-li tímto dnem neexistující 29.únor, použijte data z 1.března). Pokud by se v ročníku náhodou objevila dvojčata a neměla to štěstí mít rozdílné datum narození, zvolí si pro účel úlohy každé jiné datum nebo jinou družici. Výpočty provádějte pro družici PRN = 1. Pokud v daném souboru formátu SP3 není družice PRN = 1, použijte nejbližší vyšší PRN (2, 3, ... ).  
Na základě skutečných drah družic daných svými přesnými efemeridami v souborech formátu SP3, poskytovaných na serveru služby [https://cddis.nasa.gov/archive/gnss/products/ CDDIS] (pozor: nově nutná registrace) [ftp://cddis.gsfc.nasa.gov/gnss/products/ CDDIS] [http://igscb.jpl.nasa.gov/igscb/product IGS], sledujte časový průběh oskulačních elementů. Použijte k tomu data pro den = den Vašeho narození + 15 let (resp. totožné datum, pouze rok se zvýší o 15, je-li tímto dnem neexistující 29.únor, použijte data z 1.března). Pokud by se v ročníku náhodou objevila dvojčata a neměla to štěstí mít rozdílné datum narození, zvolí si pro účel úlohy každé jiné datum nebo jinou družici. Výpočty provádějte pro družici PRN = 1. Pokud v daném souboru formátu SP3 není družice PRN = 1, použijte nejbližší vyšší PRN (2, 3, ... ).  


;Stručný postup:
;Stručný postup:
Řádek 31: Řádek 32:
:Tz. musíte zadat polohu pólu (xp, yp), rozdíl časů UT1-UTC a rozdíl GPS-UTC, číslo družice (to je ta jednička v příkladu) a název souboru s přesnými efemeridami. Výstup potom přesměrovat do souboru. Výstup je tvořen řádky pro každou polohu družice, na řádku je postupně MJD, sekunda dne a souřadnice X, Y, Z v systému ICRF. Návod k použití se vypíše na obrazovku po spuštění prográmku bez parametrů.  
:Tz. musíte zadat polohu pólu (xp, yp), rozdíl časů UT1-UTC a rozdíl GPS-UTC, číslo družice (to je ta jednička v příkladu) a název souboru s přesnými efemeridami. Výstup potom přesměrovat do souboru. Výstup je tvořen řádky pro každou polohu družice, na řádku je postupně MJD, sekunda dne a souřadnice X, Y, Z v systému ICRF. Návod k použití se vypíše na obrazovku po spuštění prográmku bez parametrů.  
* ''Prográmek sp3crs byl vytvořen pro starší typ SP3 souboru, ve kterém nebyly střední chyby souřadnic a byly jinak číslované družice. U novějšího typu souboru (družice jsou číslované PG01-PG##, ne P  1 - P ## ; po zpracování tohoto novějšího typu souboru je výstupní soubor prázdný) musíte před jeho zpracováním všechny znaky "G0" nahradit dvěmi mezerami ("  "), znaky "G" následně nahradit jednou mezerou (" "), a umazat střední chyby (původně sloupce 6-9, po případném nahrazení znaků sloupce 7-10) - pro vysvětlení: v souboru se nacházejí sloupce PG#, x, y, z, oprava hodin, střed.ch.x, střed.ch.y, stř.ch.z, stř.ch.opravy. Pozor po vymazání dat z konců řádků tam nesmí zůstat ani žádné prázdné znaky (mezery).''
* ''Prográmek sp3crs byl vytvořen pro starší typ SP3 souboru, ve kterém nebyly střední chyby souřadnic a byly jinak číslované družice. U novějšího typu souboru (družice jsou číslované PG01-PG##, ne P  1 - P ## ; po zpracování tohoto novějšího typu souboru je výstupní soubor prázdný) musíte před jeho zpracováním všechny znaky "G0" nahradit dvěmi mezerami ("  "), znaky "G" následně nahradit jednou mezerou (" "), a umazat střední chyby (původně sloupce 6-9, po případném nahrazení znaků sloupce 7-10) - pro vysvětlení: v souboru se nacházejí sloupce PG#, x, y, z, oprava hodin, střed.ch.x, střed.ch.y, stř.ch.z, stř.ch.opravy. Pozor po vymazání dat z konců řádků tam nesmí zůstat ani žádné prázdné znaky (mezery).''
:''Pro vymazání sloupců doporučuji načtení do wordu, a označení sloupců se stisknutou klávesou ALT, ale určitě jsou i jiné způsoby (nebo si nástroj můžete naprogramovat sami). Za sloupcem s opravami hodin již nesmí být žádné další znaky (mezery) - při použití wordu možná budete muset sloupce vpravo od opravy hodin mazat nadvakrát.''
:''Pro vymazání sloupců doporučuji načtení do wordu, a označení sloupců se stisknutou klávesou ALT, ale určitě jsou i jiné způsoby (nebo si nástroj můžete naprogramovat sami). Za sloupcem s opravami hodin již nesmí být žádné další znaky (mezery) - při použití wordu možná budete muset sloupce vpravo od opravy hodin mazat nadvakrát.''
:''Pro odmítače MS Office doporučuji např. editor Emacs s možnosti výběru sloupce pomocí (C-x rt).''
:''Pro odmítače MS Office doporučuji např. editor Emacs s možnosti výběru sloupce pomocí (C-x rt).''
-->


 
Výpočet oskulačních elementů, derivace vektoru družice podle jednotlivých Keplerovských elementů:  
 
Výpočet oskulačních elementů <!--, derivace vektoru družice podle jednotlivých Keplerovských elementů-->:  
*[ftp://athena.fsv.cvut.cz/VG/VYG1/texty/poznamky_druzice.pdf poznamky_druzice.pdf]
*[ftp://athena.fsv.cvut.cz/VG/VYG1/texty/poznamky_druzice.pdf poznamky_druzice.pdf]
 
*[ftp://athena.fsv.cvut.cz/VG/VYG1/texty/derivace_dodatek.pdf derivace_dodatek.pdf]
 
-->
<!--
<!--
*[ftp://athena.fsv.cvut.cz/VG/VYG1/texty/derivace_dodatek.pdf derivace_dodatek.pdf]
Na základě skutečných drah družic, daných svými přesnými efemeridami sledujte časový průběh oskulačních elementů družice během daného dne. Přesné efemeridy pro zadaný den si dohledejte v souborech formátu SP3, poskytovaných na serveru služby [http://igscb.jpl.nasa.gov/igscb/product IGS]. Výpočty provádějte pro družici PRN = 1.  
Na základě skutečných drah družic, daných svými přesnými efemeridami sledujte časový průběh oskulačních elementů družice během daného dne. Přesné efemeridy pro zadaný den si dohledejte v souborech formátu SP3, poskytovaných na serveru služby [http://igscb.jpl.nasa.gov/igscb/product IGS]. Výpočty provádějte pro družici PRN = 1.  


Řádek 53: Řádek 47:


'''Derivace polohy družice podle Keplerových elementů jsou pro vás vypsány v dokumentu [ftp://athena.fsv.cvut.cz/VG/VYG1/texty/derivace_dodatek.pdf derivace_dodatek.pdf]'''
'''Derivace polohy družice podle Keplerových elementů jsou pro vás vypsány v dokumentu [ftp://athena.fsv.cvut.cz/VG/VYG1/texty/derivace_dodatek.pdf derivace_dodatek.pdf]'''
-->
-->


{{Teoretická geodézie}}
{{Teoretická geodézie}}

Verze z 15. 3. 2022, 12:36

Název úlohy

GPS - oskulační elementy dráhy družice

Zadání úlohy

Na základě skutečných drah družic daných svými přesnými efemeridami v souborech formátu SP3, poskytovaných na serveru služby CDDIS, sledujte časový průběh oskulačních elementů družice. Výpočet proveďte pro družici danou svým PRN a pro všechny epochy obsažené v zadaném jednodenním SP3-souboru.

Souřadnice družice v terestrickém systému ze soubru SP3 je nejprve nutno na základě známých parametrů orientace Země převést do nebeského systému. Dále postupně pro každé dvě následující polohy družice vypočtete oskulační Keplerovy parametry dráhy družice definované těmito dvěma body. Vypočtené časové řady oskulačních Keplerových parametrů vynesete vhodně do grafu (samostatně pro každý parametr).

Podrobnější informace k převodu vektoru souřadnic z terestrického do nebeského systému zde budou ještě doplněny.