152VYG1 Vyšší geodézie 1 - úloha 4: Porovnání verzí
→Numerické zadání: pridano numericke zadani |
|||
(Není zobrazeno 38 mezilehlých verzí od 6 dalších uživatelů.) | |||
Řádek 1: | Řádek 1: | ||
==Termín odevzdání== | ==Termín odevzdání== | ||
13.12.2012 | |||
==Název úlohy== | ==Název úlohy== | ||
GPS - oskulační elementy dráhy družice | GPS - oskulační elementy dráhy družice | ||
==Zadání úlohy== | ==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 [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í + 20 let (resp. totožné datum, pouze rok se zvýsí o 20, je-li tímto dnem neexistující 29.únor, použijte data z 1.dubna). 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: | |||
# Souřadnice dané družice převeďte ze systému ITRF do nebeského rámce ICRF pomocí programu [http://athena.fsv.cvut.cz/TEG1.2008/zadani/sp3crs.exe sp3crs.exe]. Vstupem kromě souboru ve formátu SP3 jsou také souřadnice pólu a hodnoty UT1 - UTC, GPS - UTC (v souboru ERP). Pozor na typ souboru SP3 - viz poznámky k výpočtu. | |||
# Postupně pro každé dvě následující polohy družice vypočtěte keplerovy parametry (oskulační) dráhy definované těmito dvěma body. | |||
# Vypočtené časové řady parametrů vyneste vhodně do grafu, samostatně pro každý parametr. | |||
;Poznámky k výpočtu | |||
* Nejprve převeďte datum (datum Vašeho narození + 20 let) na číslo GPS týdne (GPSW - GPS week) a číslo dne v týdnu (DOW - day of week). Číslování dnů v týdnu je následující: 0. den je neděle, 6. den je sobota. K převodu časů můžete použít [ftp://athena.fsv.cvut.cz/VG/VYG1/SW/prevod_casu.zip jednoduché funkce pro Matlab], ve kterých jsou vyřešeny převody mezi různými časovými škálami. Funkce jsou vytvořeny pro Octave, ale měly by fungovat i pod MatLabem. Jinou a možná rychlejší alternativou je použít [http://www.geod.nrcan.gc.ca/tools-outils/tools_e.php GPS kalendář]. | |||
* Pro dané datum vyhledejte na serveru IGS soubory s přesnými drahami družic (nápověda: dráhy nejdete někde pod odkazem DATA & PRODUCTS). Soubor, který hledáte, musí mít následující název: | |||
:*IGS{GPSW}{DOW}.SP3.* | |||
:tz. jméno je tvořeno zkratkou IGS, číslem GPS týdne, číslem dne, příponou SP3 a případně další příponou označující typ komprese. Kromě tohoto souboru budete potřebovat i soubor s parametry rotace Země, který má podobný název IGS{GPSW}7.ERP.* a vyskytuje se ve stejném adresáři. V něm najdete rozdíl časů UT1-UTC a polohu pólu. Je-li přípona souboru .Z jde o kompresi pomoci přikazu compress. WinZip tento formát rozbalit neumí, ve Windows použijte programy [http://www.7-zip.org/ 7-zip], [ftp://athena.fsv.cvut.cz/VG/VYG1/SW/compress.exe compress.exe] ([ftp://athena.fsv.cvut.cz/VG/VYG1/SW/compress.txt návod], [ftp://athena.fsv.cvut.cz/VG/VYG1/SW/decompr.bat decompr.bat]), nebo [ftp://athena.fsv.cvut.cz/VG/VYG1/SW/gzip.exe gzip.exe]. | |||
* Nyní budete potřebovat najít rozdíl mezi časem TAI a UTC, tzv. number of leap seconds, pro Váš konkrétní datum. Tento údaj je k dispozici na serveru služby [http://www.iers.org IERS] nebo přímo v tomto [http://hpiers.obspm.fr/eop-pc/earthor/utc/TAI-UTC_tab.html přehledu]. Vztahy pro převod časů jsou následující: | |||
:*TAI = GPS + 19 sec (konstantní rozdíl) | |||
:*TAI = UTC + leap_sec (variabilní rozdíl) | |||
* Máte-li potřebné informace, můžete souřadnice družice v souborech SP3 v systému ITRF transformovat do systému ICRF. Použijte k tomu prográmek [http://athena.fsv.cvut.cz/TEG1.2008/zadani/sp3crs.exe sp3crs.exe]. Použití prográmku je velmi jednoduché, stačí na příkazovou řádku napsat: | |||
:* sp3crs -xp 80630 -yp 175490 -ut 741016 -gp 13 1 IGS07600.SP3 > vystup.txt | |||
:* pozn. v linuxových distribucích snad půjde soubor sp3crs.exe také spustit (zkoušeno na Fedoře 11) | |||
: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. Formát lze převést ručně, pomocí vlastního prográmku nebo použít již nějaký vytvořený.'' | |||
:* ''Manuálně to lze udělat takto: U novějšího typu souboru jsou družice číslované PG01-PG##, ne P 1 - P ## a musíte tedy 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. Editaci sloupců podporuje ve Windows např. [http://www.pspad.com/ pspad], v Linuxových distribucích určite [http://www.gnu.org/software/emacs/ emacs] (výběr sloupce pomocí (C-x rt)), [http://www.geany.org/ geany] a jistě mnoho dalších.'' | |||
:* ''Ke konverzi můžete využít např. tento perlovský prográmek [ftp://athena.fsv.cvut.cz/VG/VYG1/SW/sp3c_2_sp3a.pl sp3c_2_sp3a.pl] pomocí příkazu:'' perl sp3c_2_sp3a.pl < igs15991.sp3 > output.sp3, ''kde perl získáte např. tady [http://www.perl.org/get.html download perl]'' | |||
*'''Pokud by vám tento postup opakovaně nefungoval, použijte o 5 let starší sp3 soubor.''' | |||
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/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. | ||
Dané civilní datum si převeďte na datum GPS pomocí programu [ftp://athena.fsv.cvut.cz/vg/VYG1/zadani/uloha4/gpstime.exe gpstime.exe]. Souřadnice dané družice převeďte ze systému ITRF do nebeského rámce ICRF pomocí programu [ftp://athena.fsv.cvut.cz/vg/VYG1/zadani/uloha4/sp3crs.exe sp3crs.exe]. Vstupem kromě souboru formátu SP3 jsou také souřadnice pólu a hodnoty UT1 - UTC, GPS – UTC, které zjistíte ze souboru ERP. | Dané civilní datum si převeďte na datum GPS pomocí programu [ftp://athena.fsv.cvut.cz/vg/VYG1/zadani/uloha4/gpstime.exe gpstime.exe]. Souřadnice dané družice převeďte ze systému ITRF do nebeského rámce ICRF pomocí programu [ftp://athena.fsv.cvut.cz/vg/VYG1/zadani/uloha4/sp3crs.exe sp3crs.exe]. Vstupem kromě souboru formátu SP3 jsou také souřadnice pólu a hodnoty UT1 - UTC, GPS – UTC, které zjistíte ze souboru ERP. | ||
Dále pak postupně z každých dvou po sobě následujících epoch vypočítejte Keplerovy elementy (oskulační) dráhy definované těmito dvěma body. Vypočtené časové řady elementů vyneste vhodně do grafu, samostatně pro každý element. | Dále pak postupně z každých dvou po sobě následujících epoch vypočítejte Keplerovy elementy (oskulační) dráhy definované těmito dvěma body. Vypočtené časové řady elementů vyneste vhodně do grafu, samostatně pro každý element. | ||
'''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]''' | |||
==Numerické zadání== | ==Numerické zadání== | ||
{| class = "border" | {| class = "border" | ||
Řádek 134: | Řádek 172: | ||
| 58|| Vyhnalová Klára || 23.11.1999 | | 58|| Vyhnalová Klára || 23.11.1999 | ||
|} | |} | ||
--> | |||
<!--'''NOVÉ - nefungující sp3crs.exe pro mladší ročníky:''' | |||
Omlouváme se všem mladším ročníkům. Program sp3crs.exe byl vytvořen pro starší formát sp3 a pro novější nefunguje. Máte 2 možnosti: | |||
* stáhnout si o 5 let mladší sp3 soubor | |||
* upravit nefungující soubor s daty podle návodu níže na této stránce--> | |||
[ | [http://geo.fsv.cvut.cz/wiki/index.php/Vy%C5%A1%C5%A1%C3%AD_geod%C3%A9zie_1 Vyšší geodézie 1] | ||
{{Teoretická geodézie}} |
Aktuální verze z 26. 9. 2012, 08:21
Termín odevzdání
13.12.2012
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 IGS, sledujte časový průběh oskulačních elementů. Použijte k tomu data pro den = den Vašeho narození + 20 let (resp. totožné datum, pouze rok se zvýsí o 20, je-li tímto dnem neexistující 29.únor, použijte data z 1.dubna). 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
- Souřadnice dané družice převeďte ze systému ITRF do nebeského rámce ICRF pomocí programu sp3crs.exe. Vstupem kromě souboru ve formátu SP3 jsou také souřadnice pólu a hodnoty UT1 - UTC, GPS - UTC (v souboru ERP). Pozor na typ souboru SP3 - viz poznámky k výpočtu.
- Postupně pro každé dvě následující polohy družice vypočtěte keplerovy parametry (oskulační) dráhy definované těmito dvěma body.
- Vypočtené časové řady parametrů vyneste vhodně do grafu, samostatně pro každý parametr.
- Poznámky k výpočtu
- Nejprve převeďte datum (datum Vašeho narození + 20 let) na číslo GPS týdne (GPSW - GPS week) a číslo dne v týdnu (DOW - day of week). Číslování dnů v týdnu je následující: 0. den je neděle, 6. den je sobota. K převodu časů můžete použít jednoduché funkce pro Matlab, ve kterých jsou vyřešeny převody mezi různými časovými škálami. Funkce jsou vytvořeny pro Octave, ale měly by fungovat i pod MatLabem. Jinou a možná rychlejší alternativou je použít GPS kalendář.
- Pro dané datum vyhledejte na serveru IGS soubory s přesnými drahami družic (nápověda: dráhy nejdete někde pod odkazem DATA & PRODUCTS). Soubor, který hledáte, musí mít následující název:
- IGS{GPSW}{DOW}.SP3.*
- tz. jméno je tvořeno zkratkou IGS, číslem GPS týdne, číslem dne, příponou SP3 a případně další příponou označující typ komprese. Kromě tohoto souboru budete potřebovat i soubor s parametry rotace Země, který má podobný název IGS{GPSW}7.ERP.* a vyskytuje se ve stejném adresáři. V něm najdete rozdíl časů UT1-UTC a polohu pólu. Je-li přípona souboru .Z jde o kompresi pomoci přikazu compress. WinZip tento formát rozbalit neumí, ve Windows použijte programy 7-zip, compress.exe (návod, decompr.bat), nebo gzip.exe.
- Nyní budete potřebovat najít rozdíl mezi časem TAI a UTC, tzv. number of leap seconds, pro Váš konkrétní datum. Tento údaj je k dispozici na serveru služby IERS nebo přímo v tomto přehledu. Vztahy pro převod časů jsou následující:
- TAI = GPS + 19 sec (konstantní rozdíl)
- TAI = UTC + leap_sec (variabilní rozdíl)
- Máte-li potřebné informace, můžete souřadnice družice v souborech SP3 v systému ITRF transformovat do systému ICRF. Použijte k tomu prográmek sp3crs.exe. Použití prográmku je velmi jednoduché, stačí na příkazovou řádku napsat:
- sp3crs -xp 80630 -yp 175490 -ut 741016 -gp 13 1 IGS07600.SP3 > vystup.txt
- pozn. v linuxových distribucích snad půjde soubor sp3crs.exe také spustit (zkoušeno na Fedoře 11)
- 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. Formát lze převést ručně, pomocí vlastního prográmku nebo použít již nějaký vytvořený.
- Manuálně to lze udělat takto: U novějšího typu souboru jsou družice číslované PG01-PG##, ne P 1 - P ## a musíte tedy 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. Editaci sloupců podporuje ve Windows např. pspad, v Linuxových distribucích určite emacs (výběr sloupce pomocí (C-x rt)), geany a jistě mnoho dalších.
- Ke konverzi můžete využít např. tento perlovský prográmek sp3c_2_sp3a.pl pomocí příkazu: perl sp3c_2_sp3a.pl < igs15991.sp3 > output.sp3, kde perl získáte např. tady download perl
- Pokud by vám tento postup opakovaně nefungoval, použijte o 5 let starší sp3 soubor.
Výpočet oskulačních elementů, derivace vektoru družice podle jednotlivých Keplerovských elementů: