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

Verze z 21. 2. 2022, 02:54

Název úlohy

GPS - výpočet drah družic

Zadání úlohy

Příklad 1.

Dráhy družic jsou zadány svými Keplerovými elementy $a$ , $e$ , $\omega$ , $i$ , $\Omega$ v inerciálním nebeském geocentrickém systému. Dále uvažujte, že se každá družice v čase $t_{0}=0^{s}$ nachází v perigeu své eliptické dráhy (pohyb družice tedy bude sledován v relativních časových jednotkách od průchodu perigeem) a družice vykonává nerušený pohyb.

Družice	$a$ [km]	$e$	$\omega$ [°]	$i$ [°]	$\Omega$ [°]
GPS (a)	26560	0.01	0	55	0
GPS (b)	26560	0.01	0	55	60
GOCE	6629	0.004	144.2	96.6	257.7
MOLNIJA	26554	0.7	270	63	245

Pro každou ze zadaných družic určete při zanedbání vlivu precese, nutace, pohybu pólu, variace v délce dne a za předpokladu, že v čase $t_{0}=0^{s}$ mají nebeský systém a terestrický systém souhlasně orientované osy v prostoru:

vypočtěte oběžnou dobu družice a udejte ji v [h min s];
vykreslete dráhy družic v nebeském systému (3D) v intervalu 1 dne;
vykreslete dráhy družic v terestrickém systému (3D) v intervalu 1 dne;
vykreslete stopu dráhy družice na zemském povrchu (Ground Track Plot) v intervalu 3 dnů.

Pro výpočet uvažujte geocentrickou gravitační konstantu $GM=398600.5km^{3}s^{-2}$ a rotaci terestrického systému úhlovou rychlostí $\omega _{E}=7292115.1467*10^{-11}s^{-1}$ . Výpočty drah proveďte v dostatečně jemném kroku výpočetních epoch (např. po $30^{s}$ ). Dosažené výsledky komentujte, různé typy drah vzájemně srovnejte.

Příklad 2.

K dispozici máte soubor vysílaných (broadcast) efemerid ve formátu RINEX (tzv. navigační zprávu, označení koncovkou 'n') pro stanici GOPE sítě EUREF EPN vytvořený dne x.x.2021.

Pro GPS družici danou PRN kódem spočítejte pomocí vysílaných efemerid její polohu v čase $t$ dne x.x.2021. Souřadnice družice GPS určete v systému WGS84, ve kterém jsou zadány i dráhové elementy. Pro výpočet Keplerových elementů v zadaném čase $t$ použijte vysílané efemeridy dané družice v zadaném souboru RINEX v epoše $t_{0}$ , která je nejbližší nižší zadanému času $t$ .

Pro zadanou GPS družici a stanovisko pozorovatele dané zeměpisnými souřadnicemi nakreslete následující grafy pro celý den x.x.2021:

graf znázorňující polohu družice vůči pozorovateli v místním topocentrickém systému (Sky Plot);
graf znázorňující závislost elevačního úhlu na čase (Elevation Plot).

Polohy GPS družice během zadaného dne vypočítejte na základě Keplerových elementů družice v čase $t_{0}$ .

Dokumenty ke stažení

Poznámky k nerušenému a rušenému pohybu družic: poznamky_druzice.pdf.

Parametry elipsoidu WGS84 můžete nalézt zde.

Obrázek k transformaci z geocentrické na topocentrickou souřadnicovou soustavu: topoc.jpg.

Datový soubor pro Matlab obsahující pobřežní linie kontinentů: coast.mat (proměnná lat obsahuje vektor zeměpisných šířek a long vektor zeměpisných délek pobřežní linie).

Funkce pro Matlab na tvorbu Sky Plotu: skyplot.m.

Pro kontrolu správnosti vašich výpočtů můžete využít grafické výstupy pro řešení 1.příkladu (nerušený pohyb) vyhotovené pro první 3 hodiny pohybu družic: nebeský systém, terestrický systém, Ground Track Plots: GPS (a), GPS (b), GOCE, MOLNIJA.

@@ Řádek 30: / Řádek 30: @@
 '''Příklad 2.'''
-(bude doplněn ve 2.týdnu výuky)
+K dispozici máte soubor vysílaných (broadcast) efemerid ve formátu RINEX (tzv. navigační zprávu, označení koncovkou 'n') pro stanici GOPE sítě EUREF EPN vytvořený dne x.x.2021.
-<!--
-K dispozici máte soubor vysílaných efemerid ve formátu RINEX (tzv. navigační zprávu, označeni koncovkou 'n') pro stanici GOPE sítě EUREF EPN vytvořený dne 1.9.2016.
-Pro družici danou PRN kódem spočítejte pomocí vysílaných efemerid její polohu v čase <math>t</math> dne 1.9.2016. Souřadnice družice určete v systému WGS84, ve kterém jsou zadány i dráhové elementy.
+Pro GPS družici danou PRN kódem spočítejte pomocí vysílaných efemerid její polohu v čase <math>t</math> dne x.x.2021. Souřadnice družice GPS určete v systému WGS84, ve kterém jsou zadány i dráhové elementy. Pro výpočet Keplerových elementů v zadaném čase <math>t</math> použijte vysílané efemeridy dané družice v zadaném souboru RINEX v epoše <math>t_0</math>, která je nejbližší nižší zadanému času <math>t</math>.
-Pro výpočet Keplerových elementů v čase <math>t</math> použijte vysílané efemeridy v epoše nejbližší nižší zadanému času.
-Pro tuto družici nakreslete následující grafy:
-* graf znázorňující závislost elevačního úhlu na čase,
-* a graf znázorňující polohu družice vůči pozorovateli v místním topocentrickém systému (skyplot).
-Oba grafy sestavte pro celý celý den 1.9.2016. Polohy družice během dne vypočítejte na základě Keplerových elementů družice v čase <math>t_0</math>. Jako zeměpisné souřadnice pozorovatele použijte souřadnice stanice uvedené na stránkách sítě CZEPOS (ETRS89).
+Pro zadanou GPS družici a stanovisko pozorovatele dané zeměpisnými souřadnicemi nakreslete následující grafy pro celý den x.x.2021:
+* graf znázorňující polohu družice vůči pozorovateli v místním topocentrickém systému (Sky Plot);
+* graf znázorňující závislost elevačního úhlu na čase (Elevation Plot).
+Polohy GPS družice během zadaného dne vypočítejte na základě Keplerových elementů družice v čase <math>t_0</math>.
+<!--
 Soubory RINEX (navigační efemeridy) je možné si stáhnout na stránkach [ftp://igs.ensg.ign.fr/pub/igs/data/ IGS].
 ==Numerické zadání==
-Numerické zadání s elipsoidickými souřadnicemi identických bodů v obou soustavách S1 a S2 naleznete v adresáři http://athena.fsv.cvut.cz:8000/TG2/helmert/zadani v souboru '''tg2_2022_u1_xx.m''', kde '''xx''' je číslo zadání. Číslo zadání studenta odpovídá číslování uvedenému na stránkách cvičení TGD2.
+Numerické zadání pro 2. příklad naleznete v adresáři http://athena.fsv.cvut.cz:8000/TG2/orbits/zadani v souboru '''tg2_2022_u1_xx.m''', kde '''xx''' je číslo zadání. Číslo zadání studenta odpovídá číslování uvedenému na stránkách cvičení TGD2.
 -->