Rozdział 2
Poło˙zenia ciał na sferze
Streszczenie
Kierunki do ciał niebieskich mo˙zna okre´sla´c równnie˙z przy zało˙zeniu, ˙ze ciała te znajduj ˛
a si˛e na
sferze o promieniu równym jedno´sci. Geometria sfery jest geometri ˛
a na zakrzywionej powierz-
chni dwuwymiarowej i w wielu przypadkach ró˙zni si˛e od euklidesowej geometrii na płaszczy´znie.
Elementami sfery wykorzystywanymi w astronomii s ˛
a koła wielkie i skonstruowane z ich po-
moc ˛
a dwuk ˛
aty i trójk ˛
aty sferyczne. W geometrii na sferze koła wielkie pełni ˛
a rol˛e analogiczn ˛
a do
prostych w planimetrii. Suma k ˛
atów wewn˛etrznych w trójk ˛
acie sferycznym jest zawsza wi˛eksza
od
. Na sferze mog ˛
a istnie´c trójk ˛
aty, w których wszystkie k ˛
aty s ˛
a k ˛
atami prostymi. Elementy
trójk ˛
ata sferycznego spełniaj ˛
a kilka grup równa´n pozwalaj ˛
acych na rozwi ˛
azywanie wielu zagad-
nie´n z zakresu astronomii sferycznej. Do najcz˛e´sciej wykorzystywanych wzorów nale˙z ˛
a wzory
sinusów i cosinusów. Problemy rozwi ˛
azywane w ramach trygonometrii sferycznej daj ˛
a si˛e tak˙ze
uj ˛
a´c w formali´zmie wektorowym.
Poło˙zenie ciała na sferze ustalone jest z pomoc ˛
a dwóch k ˛
atów, azymutalnego
i polarnego
.
Ka˙zdy układ współrz˛ednych sferycznych wymaga okre´slenia bieguna układu, wzgl˛edem którego
mierzony jest k ˛
at
, oraz koła wielkiego pełni ˛
acego rol˛e płaszczyzny odniesienia, słu˙z ˛
acej jako
pocz ˛
atek rachuby dwu´sciennego k ˛
ata
.
W astronomii wykorzystuje si˛e wiele ró˙znych układów współrz˛ednych, st ˛
ad konieczna jest
umiej˛etno´s´c przeliczania współrz˛ednych pomi˛edzy dwoma układami.
Transformacje współ-
rz˛ednych dokonuje si˛e drog ˛
a rozwi ˛
azania odpowiedniego trójk ˛
ata sferycznego albo za pomoc ˛
a
macierzy obrotu, w szczególno´sci z wykorzystaniem k ˛
atów Eulera.
Słowa kluczowe: sfera niebieska, koło małe, koło wielkie, bieguny koła wielkiego, k ˛
at sferyczny,
dwuk ˛
at sferyczny, trójk ˛
at paralaktyczny, nadmiar sferyczny, współrz˛edne sferyczne, współrz˛edne
prostok ˛
atne, triada ortogonalna, transformacje współrz˛ednych, k ˛
aty Eulera
10
Poło˙zenia ciał na sferze
2.1
Definicja poło˙ze ´n punktów w przestrzeni
Umawiamy si˛e, ˙ze otaczaj ˛
aca nas przestrze´n jest przestrzeni ˛
a euklidesow ˛
a, trójwymiarow ˛
a, ˙ze
w tej przestrzeni dany jest układ prostok ˛
atnych osi współrz˛ednych
(x;
y
;
z
)
rozpi˛ety na trójce
wersorów
i;
j;
k
.
1
Pomi˛edzy wersorami obowi ˛
azuj ˛
a nast˛epuj ˛
ace zale˙zno´sci
i
=
j
k
j
=
k
i
k
=
i
j
i
T
i
=
j
T
j
=
k
T
k
=
1
(2.1)
Trójk˛e
i;
j;
k
mo˙zemy uj ˛
a´c w formie macierzy
3
3
zwanej triad ˛
a
R
R
=
[i;
j;
k℄
(2.2)
Jej elementy s ˛
a cosinusami kierunkowymi kierunków
i;
j;
k
. Transpozycja triady
R
ma posta´c
R
T
=
2
4
i
T
j
T
k
T
3
5
(2.3)
Iloczyn
R
T
R
=
2
4
i
T
i
i
T
j
i
T
k
j
T
i
j
T
j
j
T
k
k
T
i
k
T
j
k
T
k
3
5
=
2
4
1
0
0
0
1
0
0
0
1
3
5
(2.4)
Co poci ˛
aga
R
1
=
R
T
, czyli triada
R
jest macierz ˛
a ortogonaln ˛
a.
Stwierdzenie, ˙ze wektor
r
opisuje poło˙zenie ciała niebieskiego oznacza, ˙ze do naszej dys-
pozycji s ˛
a trzy składowe
(x;
y
;
z
)
, czyli trzy liczby wyznaczone wzgl˛edem triady
R
. Na pi´smie
wyra˙zamy to stwierdzenie za pomoc ˛
a zapisu
r
=
R
2
4
x
y
z
3
5
(2.5)
2.2
Elementy geometryczne na sferze
Sfera.
Sfera jest to powierzchnia, której punkty s ˛
a równo odległe od punktu zwanego ´srodkiem sfery.
Je´sli w ´srodku sfery umie´scimy pocz ˛
atek układu współrz˛ednych, to dla punktów poło˙zonych na
sferze jednostkowej b˛edzie
r
T
r
=
1
Sfera jest powierzchni ˛
a dwuwymiarowa, sko ´nczon ˛
a ale nieograniczon ˛
a. Geometria sferyczna jest
geometri ˛
a na powierzchni zakrzywionej, jest to geometria dwuwymiarowa znacznie ró˙zni ˛
aca si˛e
1
A w jaki sposób ustalono orientacj˛e trójki wersorów
i;
j;
k
? Pytanie to wykracza poza ramy tego wykładu, bowiem
po´srednio dotyczy spraw ostatecznych podobnie jak pytanie — czy najpierw było jajo, czy kura? Poniewa˙z astronomia
nie zajmuje si˛e takimi problemami, dlatego poprzestajemy na tym, ˙ze elementy tej trójki s ˛
a znamymi cosinusami kierunk-
owymi prostych, wzdłu˙z których le˙z ˛
a wersory
i;
j;
k
.
2.2 Elementy geometryczne na sferze
11
X
A
B
P
Q
O
Rysunek 2.1: Ilustracja elementów sfery:
O
jest ´srodkiem sfery, okr ˛
ag przechodz ˛
acy przez punkty
A;
X
;
B
tradycyjnie nazwany jest kołem wielkim. Punkty
P
i
Q
s ˛
a biegunami koła wielkiego
AX
B
.
od euklidesowej geometrii na płaszczy´znie. W szczególno´sci, na sferze nie istniej ˛
a linie proste,
ich rol˛e graj ˛
a okr˛egi tradycyjnie zwane kołami wielkimi.
Koło wielkie.
Ka˙zde przeci˛ecie sfery płaszczyzn ˛
a jest okr˛egiem. Przeci˛ecie sfery płaszczyzn ˛
a przechodz ˛
ac ˛
a
przez ´srodek sfery jest kołem wielkim, np. koło
AX
z rysunku 2.1. Promie´n koła wielkiego
jest równy promieniowi sfery, a wektory
r
opisuj ˛
ace poło˙zenia punktów koła wielkiego spełniaj ˛
a
równanie
n
T
r
=
0
(2.6)
gdzie
n
jest wektorem definiuj ˛
acym jeden z biegunów koła wielkiego. Ko ´nce ´srednicy prostopadłej
do koła wielkiego sfery nazywamy biegunami tego koła. Na rysunku 2.1 biegunami koła wielkiego
AX
s ˛
a punkty
P
i
Q
, punkty poło˙zone diametralnie. Zauwa˙zmy, ˙ze dowolne koło wielkie prze-
chodz ˛
ace przez jeden z biegunów
P
, musi tak˙ze przechodzi´c przez drugi biegun
Q
.
Dwa punkty sfery, które nie s ˛
a punktami diametralnymi jak np.
A
i
X
, wyznaczaj ˛
a koło
wielkie jednoznacznie, bowiem ł ˛
acznie ze ´srodkiem sfery (punktem
O
) jednoznacznie okre´slaj ˛
a
płaszczyzn˛e, której przeci˛ecie ze sfer ˛
a jest kołem wielkim. Na punktach
A
,
X
rozpi˛ete s ˛
a dwa łuki,
mniejszy z nich nazywany lini ˛
a geodezyjn ˛
a, jest najkrótsz ˛
a krzyw ˛
a jak ˛
a mo˙zna na sferze poł ˛
aczy´c
punkty
A
i
X
. Linie geodezyjne, (zwane te˙z odległo´sciami sferycznymi) pełni ˛
a na sferze rol˛e
analogiczn ˛
a jak linie proste w geometrii euklidesowej.
Poniewa˙z promie´n sfery
r
=
1
, długo´s´c łuku koła wielkiego równa jest k ˛
atowi w radianach
2
jaki ten łuk rozpina wzgl˛edem ´srodka sfery.
Dwuk ˛
at sferyczny.
W wyniku przeci˛ecia tej samej sfery dwoma kołami wielkimi wyznaczone zostan ˛
a cztery obszary
(powierzchnie) zwane dwuk ˛
atami sferycznymi (rysunek 2.2-a). Dwuk ˛
aty przeciwległe s ˛
a parami
przystaj ˛
ace. Dwuk ˛
at sferyczny okre´slony jest k ˛
atem sferycznym, np. k ˛
atem
A
na rysunku 2.2-a.
K ˛
at ten jest równy k ˛
atowi liniowemu okre´slonemu przez płaszczyzny, na których le˙z ˛
a koła wielkie
tworz ˛
ace dany dwuk ˛
at. Półokr˛egi tych kół wielkich nazywamy bokami dwuk ˛
ata. Na danej sferze,
boki wszystkich dwuk ˛
atów s ˛
a równe i maj ˛
a długo´s´c
r
, gdzie
r
jest promieniem sfery.
3
Pole
2
Oprócz radianów mo˙zna oczywi´scie u˙zywa´c innych miar k ˛
ata.
3
Ze wzgl˛edu na do´s´c cz˛este nieporozumienia, warto zapami˛eta´c, ˙ze dwuk ˛
at sferyczny (podobie trójk ˛
at sferyczny) jest
fragmentem powierzchni sfery, a nie k ˛
atem.
12
Poło˙zenia ciał na sferze
A
a)
b)
A
O
A’
Rysunek 2.2: Dwuk ˛
aty sferyczne: a) cztery dwuk ˛
aty powstaj ˛
a w wyniku przeci˛ecia sfery dwoma
kołami wielkimi; b) dwuk ˛
aty sferyczne s ˛
a w pełni opisane przez promie´n sfery i k ˛
at sferyczny
A
.
pwierzchni dwuk ˛
ata sferycznego mo˙zna obliczy´c ze wzoru
S
=
2r
2
A
gdzie
A
jest k ˛
atem dwuk ˛
ata wyra˙zonym w radianach.
K ˛
at pomi˛edzy płaszczyznami kół wielkich cz˛esto jest nazywany k ˛
atem sferycznym. O K ˛
acie
sferycznym mo˙zna te˙z powiedzie´c, ˙ze jest to k ˛
at pomi˛edzy stycznymi wystawionymi w punkcie
wzajemnego przeci˛ecia si˛e kół wielkich (patrz rysunek 2.2-b).
Trójk ˛
at sferyczny.
Trzy koła wielkie nie przecinaj ˛
ace si˛e w jednej parze punktów diametralnych tworz ˛
a na sferze
osiem obszarów zwanych trójk ˛
atami sferycznymi (rysunek 2.3-a). Znaj ˛
ac elementy jednego z nich
(czyli trzy boki i trzy k ˛
aty sferyczne, np. łuki
a;
b;
kół wielkich i k ˛
aty wewn˛etrzne
A;
B
;
C
),
łatwo wyznaczy´c elementy wszystkich pozostałych trójk ˛
atów. Dlatego zwykle rozpartuje si˛e
zale˙zno´sci pomi˛edzy elementami tylko jednego trójk ˛
ata, którego wszystkie boki s ˛
a krótsze od
połowy obwodu koła wielkiego. Taki trójk ˛
at sferyczny nosi nazw˛e trójk ˛
ata paralaktycznego lub
trójk ˛
ata eulerowskiego.
Boki trójk ˛
ata sferycznego tradycyjnie oznaczane s ˛
a małymi literami, a ich długo´sci mierzone
s ˛
a za pomoc ˛
a płaskich k ˛
atów k ˛
ata trój´sciennego
O
AB
C
, np. na rysunku 2.3
=
℄AO
B
.
Wewn˛etrzne k ˛
aty przy wierzchołkach trójk ˛
ata oznaczane du˙zymi literami
A;
B
;
C
mierzone s ˛
a
k ˛
atami dwu´sciennymi (k ˛
atami sferycznymi) tego samego k ˛
ata trój´sciennego.
Trójk ˛
aty sferyczne eulerowskie maj ˛
a pewne własno´sci wspólne z trójk ˛
atami płaskimi, np.
dowolny bok trójk ˛
ata jest mniejszy od sumy, ale wi˛ekszy od ró˙znicy dwóch boków pozostałych.
Jednak mamy te˙z mi˛edzy nimi istotne ró˙znice, np. w trójk ˛
acie sferycznym suma k ˛
atów wewn˛etrz-
nych nie jest stała, bowiem, ˙ze suma nale˙zy do przedziału
<
s
=
A
+
B
+
C
<
3
Trójk ˛
at płaski mo˙ze mie´c tylko jeden k ˛
at prosty, trójk ˛
at sferyczny niekoniecznie, mo˙ze mie´c ich
dwa a nawet trzy. Ró˙znica
s
=
"
gdzie nazywanajest nadmiarem sferycznym.
2.2 Elementy geometryczne na sferze
13
A
C
B
a)
b)
B
C
A
b
O
a
c
Rysunek 2.3: Trójk ˛
aty sferyczne: a) osiem trójk ˛
atów sferycznych mo˙zna otrzyma´c z przeci˛ecia
trzech kół wielkich; b) Trójk ˛
at paralaktyczny A,B,C rozpi˛ety na trzech wektorach jednostkowych.
Elementy trójk ˛
ata to k ˛
aty wierzchołkowe (k ˛
aty sferyczne)
A;
B
;
C
i naprzeciw nich poło˙zone boki
a;
b;
. Zarówno k ˛
aty jak i boki mierzone s ˛
a w jednostkach k ˛
atowych.
O
P
E
A
B
D
C
Q
F
r
M
S
θ
ψ
kolo male
θ
ψ
Rysunek 2.4: Elementy sfery: koło małe
AE
B
. Jest ono odległe o k ˛
at
od swego bieguna
P
.
Wycinek
AE
koła małego opiera si˛e ramionach
AS
i
E
S
rozwartych o k ˛
at
. Długo´s´c wycinka
wynosi
AE
=
sin
.
Pole powierzchni trójk ˛
ata sferycznego dana jest z pomoc ˛
a formuły
S
=
r
2
"
gdzie
r
jest promieniem sfery wyra˙zonym w radianach.
Koło małe.
´Slad przeci˛ecia sfery płaszczyzn ˛a nie przechodz ˛ac ˛a przez ´srodek sfery jest okr˛egiem, tradycyjnie
zwanym kołem małym. Jego biegunami s ˛
a punkty skrajne ´srednicy sfery wystawionej prostopadle
do płaszczyzny koła małego.
Promie´n koła małego jest zawsze mniejszy od promienia sfery. Na rysunku 2.4 widzimy koło
małe
AE
B
, równoległe do ´n koło wielkie
C
F
D
oraz ich wspólne w tym przypadku bieguny
P
i
Q
.
Wyprowadzimy formuł˛e na długo´s´c łuku koła małego, cz˛esto stosowan ˛
a w dalszej cz˛e´sci wykładu.
Niech
r
M
b˛edzie promieniem koła małego, łuk
AP
=
jest miar ˛
a odległo´sci koła małego od
14
Poło˙zenia ciał na sferze
O
Z
A
A’
x
y
z
O
S
R
A
Z
X
Y
θ
r
ψ
r
r
r
kolo male
Z’
ψ
B
X
Y
θ
r
ψ
Z’
Rysunek 2.5: Sferyczne współrz˛edne biegunowe
(
;
)
punktu
A
. Współrz˛edna polarna
okre´sla
k ˛
atow ˛
a odległo´s´c punktu
A
od bieguna
Z
układu współrzednych. Współrz˛edna azymutalna
jest
k ˛
atem dwu´sciennym. Ustala ona k ˛
atow ˛
a odległo´s´c płaszczyzny południka, w którym le˙zy punkt
A
, od płaszczyzny południka
Z
X
. Południk
Z
X
wybrano jako pocz ˛
atek rachuby współrz˛ednej
azymutalnej.
bieguna
P
. Z trójk ˛
ata płaskiego
AO
S
wynika
AS
=
AO
sin
AO
S
r
M
=
sin
(2.7)
Punkt
E
koła małego poł ˛
aczmy z biegunem
P
łukiem koła wielkiego i powstały łuk przedłu˙zmy
do przeci˛ecia z kołem
C
D
w punkcie
F
. Je´sli k ˛
at sferyczny
AP
E
oznaczymy przez
, to mamy
tak˙ze, ˙ze
C
O
F
=
. A poniewa˙z odcinki
AS
i
C
O
s ˛
a do siebie równoległe, podobnie
E
S
i
F
O
mamy jeszcze, ˙ze k ˛
at
AS
E
=
.
Zatem długo´s´c łuku
AE
koła małego wynosi
AE
=
r
M
AS
E
=
sin
(2.8)
2.3
Sferyczne współrz˛edne biegunowe
W celu ustalenia poło˙zenia punktu na sferze, mo˙zna wykorzysta´c ró˙zne układy współrz˛ednych.
Przykładowo, we´zmy prawoskr˛etny prostok ˛
atny zbiór osi kartezja´nskich
O
xy
z
, okre´slonych trójk ˛
a
wersorów
i;
j;
k
o pocz ˛
atkach w ´srodku
O
sfery jednostkowej. Dodatnie kierunki tych osi przeci-
naj ˛
a sfer˛e w punktach
X
;
Y
;
Z
, natomiast koła wielkie
X
Y
;
Z
X
poło˙zone s ˛
a w płaszczyznach
xO
y
i
z
O
x
, odpowiednio, patrz (rysunek 2.5).
Obierzmy na sferze punkt
A(xy
z
)
, wówczas. dla sfery jednostkowej prawdziwy jest zwi ˛
azek
x
2
+
y
2
+
z
2
=
1
(2.9)
Czyli jedna z prostok ˛
atnych współrz˛ednych
x;
y
;
z
punktu
A
jest zbyteczna, co oznacza, ˙ze w celu
ustalenia poło˙zenia ciała niebieskiego, tzn. kierunku do tego ciała, wystarczy posłu˙zy´c si˛e dwoma
liczbami.
Obok współrz˛ednych prostok ˛
atnych, praktyce astronomicznej wygodnie jest stosowa´c współ-
rz˛edne biegunowe
(r
;
;
)
jako bli˙zsze naszemu intuicyjnemu wyczuciu kierunku. Zgodnie z
2.4 Współrz˛edne prostok ˛
atne punktów na sferze
15
tradycyjn ˛
a definicj ˛
a (patrz rysunek 2.5):
r
=
O
A
jest współrz˛edn ˛
a radialn ˛
a punktu
A
,
jest współrz˛edn ˛
a polarn ˛
a punktu
A
, jest ona identyczna z k ˛
atem
Z
O
A
,
jest współrz˛edn ˛
a azymutaln ˛
a punktu
A
równ ˛
a k ˛
atowi dwu´sciennemu pomi˛edzy płasz-
czyzn ˛
a
Z
O
A
i płaszczyzn ˛
a
Z
O
X
.
Poniewa˙z dla sfery jednostkowej promie´n sfery
r
=
1
, zatem dwie współrz˛edne k ˛
atowe
(
;
)
w
pełni okre´slaj ˛
a poło˙zenie ciała na sferze: k ˛
at
jest długo´sci ˛
a łuku
Z
A
, natomiast
jest k ˛
atem
sferycznym
X
Z
A
.
W celu ustalenia poło˙zenia punktów na całej sferze, wystarczy je´sli współrz˛edne
(
;
)
przyjm ˛
a
warto´sci nale˙z ˛
ace do dziedziny
0
0
2
(2.10)
Z układem współrz˛ednych sferycznych niekiedy wi ˛
a˙ze si˛e siatk˛e współrz˛ednych, która defin-
iowana jest nast˛epuj ˛
aco:
dla
=
onst
, krzywe siatki s ˛
a małymi kołami o biegunach w
Z ;
Z
0
,
dla
=
onst
, krzywe siatki s ˛
a półkolami wielkimi przecinaj ˛
acymi si˛e w biegunach
Z ;
Z
0
.
Podsumujmy: w celu zdefiniowania jakiegokolwiek układu współrz˛ednych sferycznych musimy
(patrz rysunek 2.5):
dokona´c wyboru bieguna
Z
układu, wzgl˛edem którego mierzona b˛edzie współrz˛edna — k ˛
at
polarny
,
dokona´c wyboru koła wielkiego
Z
X
pełni ˛
acego rol˛e płaszczyzny odniesienia, wzgl˛edem
której mierzony b˛edzie dwu´scienny k ˛
at azymutalny
.
ustali´c skr˛etno´s´c układu,
poda´c jednostki miary i dziedzin˛e warto´sci k ˛
atów
i
.
Wszystkie sferyczne astronomiczne układy współrz˛ednych s ˛
a konstruowane w taki sposób. Ró˙zni ˛
a
si˛e doborem bieguna i koła odniesienia, mo˙zna ´sród nich napotka´c zarówno układy lewoskr˛etne
jak i prawoskr˛etne, cz˛esto zamiast k ˛
ata polarnego brane jest jego dopełnienie
(
=2
)
, natomiast
warto´sci k ˛
atów podawane s ˛
a w stopniach albo w jednostkach czasu.
Poza tymi ró˙znicami astronomiczne układy współrz˛ednych zawsze stanowi ˛
a realizacje sfer-
ycznych współrz˛ednych biegunowych
(
;
)
omówionych powy˙zej.
2.4
Współrz˛edne prostok ˛
atne punktów na sferze
Pomimo nieodł ˛
acznego nadmiaru (patrz równanie (2.9)), w astronomii sferycznej warto stosowa´c
współrz˛edne kartezja´nskie. Uj˛ete w postaci uporz ˛
adkowanych trójek współrz˛edne te, pozwalaj ˛
a
nada´c równaniom eleganck ˛
a i uniwersaln ˛
a wektorow ˛
a form˛e.
16
Poło˙zenia ciał na sferze
B
A
C
b
c
r
C
r
B
D
r
D
O
90
r
A
a
Rysunek 2.6: Trójk ˛
at paralaktyczny ABC rozpi˛ety na trójce wektorów jednostkowych
r
A
;
r
B
;
r
C
.
Wektor
r
D
jest prostopadły do wektorów
r
B
;
r
C
.
Je´sli
i;
j;
k
s ˛
a jednostkowymi wektorami o własno´sciach okre´slonych równaniami (2.1), je´sli
wzdłu˙z tych wersorów zorientowano dodatnie kierunki osi
x;
y
;
z
, wówczas zgodnie z równaniem
(2.5) poło˙zenie punktu
A(x;
y
;
z
)
na sferze (rysunek 2.5) okre´slone jest za pomoc ˛
a wektora poło˙ze-
nia
r
A
r
A
=
x
i
+
y
j
+
z
k
(2.11)
gdzie składowe
(x;
y
;
z
)
s ˛
a cosinusami kierunkowymi odcinka
O
A
, obliczonymi odpowiednio,
wzdłu˙z osi
X
;
Y
;
Z
, patrz rysunek 2.5
x
=
os
X
A
y
=
os
Y
A
z
=
os
Z
A
(2.12)
Pomi˛edzy współrz˛ednymi
x;
y
;
z
oraz współrz˛ednymi sferycznymi
;
tego samego punktu
A
mamy znane zwi ˛
azki
x
=
sin
os
y
=
sin
sin
z
=
os
(2.13)
Dowolny problem w astronomii sferycznej mo˙zna rozwi ˛
aza´c za pomoc ˛
a metod trygonometrii sfer-
yczej i współrz˛ednych sferycznych. W czasach przedkomputerowych pozwalało to na uzyski-
wanie rozwi ˛
aza´n oszcz˛ednych pod wzgl˛edem obliczeniowym. Obecnie dawne i nowe problemy
rozwi ˛
azujemy stosuj ˛
ac bardziej ogólne podej´scie wektorowe. Jednak poniewa˙z astronomia jest
nauk ˛
a, w której dane obserwacyjne maj ˛
a inn ˛
a rang˛e ni˙z to ma miejsce np. w fizyce, w pewnych
wypadkach znajomo´s´c metod redukcyjnych jakimi kiedy´s posługiwali si˛e astronomowie jest niezb˛edna,
by przykładowo, obserwacje komety pochodz ˛
ace z odległych epok wykorzysta´c razem z ob-
serwacjami współczesnymi.
2.5
Podstawowe wzory trygonometrii sferycznej
Wyprowadzimy cz˛esto wykorzystywane przez astronomów zwi ˛
azki pomi˛edzy elementami trójk ˛
ata
paralaktycznego. Nasze podej´scie b˛edzie bardziej współczesne: skorzystamy z zale˙zno´sci wek-
2.5 Podstawowe wzory trygonometrii sferycznej
17
torowych, przy czym składowe wektorów wyrazimy poprzez współrz˛edne sferyczne za pomoc ˛
a
formuł (2.13).
Niech
AB
C
b˛edzie trójk ˛
atem sferycznym (rysunek 2.6). Jak wida´c jest on rozpi˛ety na trzech
wektorach jednostkowych
r
A
;
r
B
;
r
C
. Wybierzmy układ współrz˛ednych
(
;
)
o biegunie w
punkcie
A
, łuk AB obierzmy za koło wielkie odniesienia miary współrz˛ednej azymutalnej
.
Po tych ustaleniach: poło˙zenie punktu
B
okre´slone jest przez
(
=
,
=
0)
, poło˙zenie punktu
C
przez
(
=
b;
=
A)
. A zgodnie z równaniami (2.13) składowe wektorów poło˙ze´n punktów
B
i
C
, wynosz ˛
a
r
B
=
(sin
;
0;
os
)
r
C
=
(sin
b
os
A;
sin
b
sin
A;
os
b)
(2.14)
K ˛
at mi˛edzy
r
A
;
r
B
jest równy długo´sci boku
B
C
trójk ˛
ata sferycznego
A;
B
;
C
, a poniewa˙z s ˛
a to
wektory jednostkowe, ich iloczyn skalarny wynosi
r
B
r
C
=
os
a
Podstawiaj ˛
ac za
r
B
i
r
C
prawe strony równa´n (2.14) otrzymamy
os
a
=
os
b
os
+
sin
b
sin
os
A
(2.15)
Jest to jedna z najbardziej podstawowych formuł trygonometrii sferycznej nazywana wzorem cos-
inusów.
4
Za pomoc ˛
a tej formuły, drog ˛
a odpowiednich przekształce´n, mo˙zna otrzyma´c dalsze wzory.
Jednak bardziej bezpo´srednio, dwa z nich uzyskamy badaj ˛
ac wektor
r
C
r
B
. Poniewa˙z k ˛
at
mi˛edzy wersorami
r
B
;
r
C
równy jest łukowi
B
C
, ich iloczyn wektorowy jest wektorem o długo´sci
sin
a
skierowanym ku punktowi
D
poło˙zonemu o
90
Æ
zarówno od
C
i
B
. Jak widzimy na rysunku
2.6, punkt
D
jest biegunem boku
B
C
trójk ˛
ata sferycznego
A;
B
;
C
. my´sl tego co powiedziano,
mamy wi˛ec
r
C
r
B
=
sin
a
r
D
(2.16)
gdzie
r
D
jest wektorem jednostkowym punktu
D
, Z pomoc ˛
a równa´n (2.14) lew ˛
a stron˛e równania
wektorowego (2.16) mo˙zemy napisa´c w postaci
r
C
r
B
=
(sin
b
os
sin
A;
os
b
sin
sin
b
os
os
A;
sin
b
sin
sin
A)
(2.17)
Jak widzimy na rysunku 2.6, sferyczne współrz˛edne punktu
D
wynosz ˛
a
(
=
B
AD
;
=
AD
)
,
zatem korzystaj ˛
ac z formuł (2.13) praw ˛
a stron˛e równania (2.16) mo˙zemy wyrazi´c jako
sin
a
r
D
=
sin
a(sin
AD
os
B
AD
;
sin
AD
sin
B
AD
;
os
AD
)
(2.18)
Mogliby´smy teraz porówna´c odpowiednie składowe w równaniach (2.17) i (2.18), jednak warto
przedtem pozby´c si˛e sinusów i cosinusów k ˛
atów
AD
i
B
AD
.
Dokonamy tego za pomoc ˛
a zwi ˛
azków mi˛edzy elementami trójk ˛
ata
AB
C
. W tym celu popa-
trzmy na trójk ˛
at sferyczny
B
AD
. Skoro
D
jest biegunem koła wielkiego
B
C
, to łuk
B
D
=
90
Æ
i
jest prostopadły do koła
B
C
. St ˛
ad k ˛
at sferyczny
AB
D
=
90
Æ
+
B
i w trójk ˛
acie sferycznym
B
AD
ze wzóru cosinusów mamy
os
AD
=
os
90
Æ
os
+
sin
90
Æ
sin
os
(90
Æ
+
B
)
os
AD
=
sin
sin
B
4
Komplet wzorów postaci (2.15) daje si˛e otrzyma´c np. poprzez cykliczn ˛
a permutacj˛e symboli abcABC.
18
Poło˙zenia ciał na sferze
Podstawiaj ˛
ac ten rezultat do składowej
z
-towej w równaniu (2.18), porównuj ˛
ac j ˛
a ze składow ˛
a
z
-tow ˛
a z równania (2.17) otrzymamy
sin
A
sin
a
=
sin
B
sin
b
Na mocy symetrii, równanie to mo˙zna uzupełni´c o dodatkowy człon, i w tej pełnej postaci nosi
ono nazw˛e wzoru sinusów
sin
A
sin
a
=
sin
B
sin
b
=
sin
C
sin
(2.19)
Równo´sci (2.19) wykorzystamy do wyprowadzenia kolejnych wzorów trygonometrii sferycznej.
W trójk ˛
acie sferycznym
B
AD
z rysunku 2.6, dla boku
AD
i k ˛
ata wierzchołkowego
AD
B
na
mocy wzoru sinusów b˛edzie
sin
(90
Æ
+
B
)
sin
AD
=
sin
B
AD
sin
90
Æ
sin
AD
sin
B
AD
=
sin(90
Æ
+
B
)
=
os
B
Kład ˛
ac ten rezultat do y-kowej składowej równania (2.18), przyrównuj ˛
ac j ˛
a ze składow ˛
a y-kow ˛
a
równania (2.17) otrzymujemy wa˙zny wzór zwany wzorem pi˛ecioelementowym
sin
a
os
B
=
os
b
sin
sin
b
os
os
A
(2.20)
Pozostałe pi˛e´c wzorów typu (2.20) otrzymamy dzi˛eki odpowiednim permutacjom symboli w
trójk ˛
acie sferycznym
AB
C
. Np. zmieniaj ˛
ac w ci ˛
agu symboli
aB
b b A
rol˛e
B
z
C
oraz
b
z
dostaniemy nast˛epny wzór pi˛ecioelementowy
sin
a
os
C
=
os
sin
b
sin
os
b
os
A
(2.21)
Z wektorowego równania (2.16) nie otrzymamy ju˙z ˙zadnych nowych formuł. Ostatni wa˙zny wzór
trygonometrii sterycznej tzw. wzór cotangensowy (czterocz˛e´sciowy) mo˙zna wydedukowa´c ze
wzorów cosinusów i sinusów. W tym celu stosujmy wzór cosinusów do boków
b
i
trójk ˛
ata
AB
C
(rysunek 2.6), mamy
os
b
=
os
a
os
+
sin
a
sin
os
B
os
=
os
a
os
b
+
sin
a
sin
b
os
C
Eliminuj ˛
ac
os
w pierwszym równaniu za pomoc ˛
a prawej strony drugiego równania, podstawia-
j ˛
ac za
sin
odpowiednie wyra˙zenie ze wzoru sinusów dostaniemy
os
b
=
os
a( os
a
os
b
+
sin
a
sin
b
os
C
)
+
sin
a
sin
b
sin
C
sin
B
os
B
a po podzieleniu obu stron przez
sin
b
, b˛edzie
ot
b
=
os
2
a
ot
b
+
os
a
sin
a
os
C
+
sin
a
sin
C
ot
B
sin
2
a
ot
b
=
sin
a( os
a
os
C
+
sin
C
ot
B
)
Dziel ˛
ac w ostatnim równaniu obie strony przez
sin
a
otrzymujemy ostatecznie
os
a
os
C
=
sin
a
ot
b
sin
C
ot
B
(2.22)
Istnieje komplet sze´sciu takich wzorów, mo˙zna go wypisa´c odpowiednio permutuj ˛
ac symbole w
trójk ˛
acie sferycznym
AB
C
.
2.6 Małe przesuni˛ecie na sferze niebieskiej
19
O
U
P
X
X’
γ
C
90−δ0
90−δ
α−αο
θ
χ
s
x
o
s
ds
s
o
C
s
L
X
O
Rysunek 2.7: Małe przesuni˛ecie na sferze niebieskiej. a) Ciało niebieskie uległo drobnemu prze-
suni˛eciu z punktu
X
do
X
0
, wzdłu˙z koła wielkiego
O
X
. b) wersja wektorowa małego przesuni˛e-
cia.
2.6
Małe przesuni˛ecie na sferze niebieskiej
W wielu problemach astronomii sferycznej mamy do czynienia z niewielkimi zmianami poło˙ze´n
ciał niebieskich (tzw. małe przesuni˛ecie). Przyczyny zmian bywaj ˛
a ró˙zne, natomiast same zmi-
any niemal zawsze przebiegaj ˛
a w taki sam sposób, dlatego warto zapozna´c si˛e ze standardowym
opisem małego przesuni˛ecia.
Wielko´s´c przesuni˛ecia mo˙ze by´c ró˙zna, m.in. jest zale˙zna od poło˙zenia obiektu, ale zawsze
przesuni˛ecie odbywa si˛e po kole wielkim ł ˛
acz ˛
acym dany obiekt z jakim´s ustalonym punktem
sfery, wspólnym dla wszystkich obiektów. Np. małe przesuni˛ecie zwane paralaks ˛
a roczn ˛
a, ma
miejsce zawsze wzdłu˙z koła wielkiego zawieraj ˛
acego kierunki ku Sło ´nca i do danego ciała nie-
bieskiego, przesuni˛ecie zwane aberacj ˛
a dobow ˛
a, przebiega po kole wielkim rozpi˛etym na kierunku
do danego ciała i na kierunku do punktu wschodu horyzontu miejsca obserwacji. Wszystkie tego
typu przesuni˛ecia mo˙zna traktowa´c jako szczególne przypadki ogólniejszego małego przesuni˛ecia
opisanego poni˙zej.
Przyjmijmy, ˙ze np. kierunek do gwiazdy
X
(;
Æ
)
z rysunku 2.7a, uległ niewielkiemu prze-
suni˛eciu do punktu
X
0
, oraz ˙ze odbyło si˛e to wzdłu˙z koła wielkiego ł ˛
acz ˛
acego
X
z punktem
O
(
0
;
Æ
0
)
. Oznaczmy łuk
O
X
przez
a łuk
X
X
0
przez
d
. Załó˙zmy, ˙ze
d
jest małym k ˛
atem
dodatnim. Niech dalej b˛edzie, ˙ze przesuni˛ecie
d
da si˛e wyrazi´c za pomoc ˛
a formułki
X
X
0
=
d
=
k
sin
(2.23)
gdzie
k
jest stał ˛
a dodatni ˛
a lub ujemn ˛
a niezale˙zn ˛
a od obranej gwiazdy, co zreszt ˛
a nie jest a˙z tak
wa˙zne dla dalszego wywodu. W ten sposób umówili´smy si˛e, ˙ze interesujemy si˛e przesuni˛eciami,
których wielko´s´c dla danego
k
, zale˙zy jedynie od odległo´sci obiektu od pewnego punktu wspól-
nego
O
.
Niech b˛edzie, ˙ze punkt
X
0
ma współrz˛edne
(
+
d;
Æ
+
dÆ
)
. Poszukamy wyra˙ze´n pozwalaj ˛
a-
cych na obliczenie przyrostów
d;
dÆ
spowodowanych małym przesunieciem
d
. W celu znalezie-
nia takich zwi ˛
azków, konstruujemy na sferze obiekt geometryczny, najlepiej taki, który ma znane
własno´sci matematyczne (np. trójk ˛
at), i którego elementy b˛ed ˛
a miały zwi ˛
azek z wyst˛epuj ˛
acymi w
naszym problemie wielko´sciami:
d;
dÆ
;
d
:
:
:
.
20
Poło˙zenia ciał na sferze
Poprowad´zmy koło małe o biegunie w
P
, przechodz ˛
ace przez
X
0
, przecinaj ˛
ace
P
X
w punkcie
U
. Mo˙zemy zatem utworzy´c trójk ˛
at
U
X
X
0
, dla którego poszukamy wyra˙ze´n na niektóre jego
elementy.
Poniewa˙z
P
X
=
,
P
X
0
=
+
d
mamy, ˙ze k ˛
at
U
P
X
=
d
. Skoro
P
X
0
=
P
U
=
90
Æ
(Æ
+
dÆ
)
, st ˛
ad z równania (2.8), z dokładno´sci ˛
a do wyrazów pierwszego rz˛edu bok
U
X
0
wynosi
U
X
0
=
d
sin(90
Æ
(Æ
+
dÆ
))
=
d
os
(Æ
+
dÆ
)
d
os
Æ
Drugi bok
U
X
=
dÆ
, bowiem
P
X
=
90
Æ
Æ
.
Poszukamy teraz wyra˙ze´n na długo´sci tych boków, przy czym chcemy by wyst˛epowały w
nich wielko´sci mo˙zliwie bezpo´srednio zwi ˛
azane z parametrami opisuj ˛
acymi małe przesuni˛ecie.
Oznaczmy k ˛
at sferyczny
O
X
P
przez
, wówczas
U
X
X
0
=
180
Æ
. Ze wzgl˛edu na bardzo
małe rozmiary w stosunku do promienia sfery, trójk ˛
at
U
X
X
0
5
w przybli˙zeniu mo˙zna traktowa´c
jako trójk ˛
at płaski, o k ˛
acie prostym w wierzchołku
U
. Dla takiego przybli˙zenia, mamy
U
X
=
X
X
0
os
(180
Æ
)
=
X
X
0
os
U
X
0
=
X
X
0
sin(180
Æ
)
=
X
X
0
sin
Kład ˛
ac za
U
X
0
i
U
X
rezultaty uzyskane wcze´sniej, bior ˛
ac równie˙z pod uwag˛e równanie (2.23)
mamy, ˙ze
os
Æ
d
=
k
sin
sin
dÆ
=
k
sin
os
(2.24)
Pozostaje nam jeszcze wyeliminowanie k ˛
atów
i
z pomoc ˛
a wyra˙ze´n, w których wyst˛epuj ˛
a
wył ˛
acznie wielko´sci znane tzn.
k
;
0
;
Æ
0
. W tym celu popatrzmy na rysunek 2.7a, na trójk ˛
at
sferyczny
O
P
X
. Mamy tu, ˙ze
P
O
=
0
co poci ˛
aga
O
P
X
=
0
. Dalej mamy
P
X
=
90
Æ
Æ
,
P
O
=
90
Æ
Æ
0
,
O
X
=
oraz
O
X
P
=
. Stosuj ˛
ac do trójk ˛
ata
O
P
X
wzor sinusów
(2.19) i wzór pi˛ecioelementowy (2.21) otrzymamy
sin
sin
=
sin(90
Æ
Æ
0
)
sin(
0
)
sin
os
=
os
(90
Æ
Æ
0
)
sin(90
Æ
Æ
)
sin(90
Æ
Æ
0
)
os
(90
Æ
Æ
)
os
(
0
)
Podstawiaj ˛
ac prawe strony tych wyra˙ze´n do równa´n (2.24) ostatecznie mamy
d
=
k
se
Æ
os
Æ
0
sin(
0
)
dÆ
=
k
(sin
Æ
os
Æ
0
os
(
0
)
os
Æ
sin
Æ
0
)
(2.25)
Aby te równania zastosowa´c w jakim´s konkretnym przypadku np. do opisu skutków zjawiska
refrakcji, wystarczy poło˙zy´c odpowiedni ˛
a warto´s´c
k
oraz współrz˛edne
(
0
;
Æ
0
)
punktu
O
.
Opis małego przesuni˛ecia dany równaniami (2.25) mo˙zna równie˙z otrzyma´c w formie wek-
torowej. Niech
s
b˛edzie wersorem poło˙zenia punktu
X
a
s
0
b˛edzie wersorem poło˙zenia punktu
O
(rysunek 2.7b). Poniewa˙z s ˛
a to wektory jednostkowe st ˛
ad wektor
s
s
0
ma długo´s´c
sin
, i
skierowany jest do punktu
L
na sferze, odległego o
90
Æ
zarówno od
O
jak i od
X
. Punkt
L
jest
wi˛ec biegunem koła wielkiego
O
X
.
Niech kierunek do punktu
X
0
opisuje wektor
s
+
ds
. By okre´sli´c poło˙zenie punktu
X
0
musimy
policzy´c składowe wektora
ds
i w tym celu szukamy zale˙zno´sci, w których ten wektor wyst˛epuje.
6
Poniewa˙z iloczyn skalarny
s
s
=
1
, ró˙zniczkuj ˛
ac to wyra˙zenie dostaniemy
s
ds
=
0
(2.26)
5
Uwaga, nie jest to trójk ˛
at sferyczny!
6
Wymaga to wyczucia, czyli nosa i dlatego najcz˛e´sciej stosowan ˛
a metod ˛
a rozwi ˛
azywania takich problemów jest metoda
prób i bł˛edów.
2.7 Wektorowe transformacje współrz˛ednych sferycznych
21
Wynika st ˛
ad, ˙ze wektor
ds
jest prostopadły do
s
, a skoro przesuni˛ecie odbywa si˛e wzdłu˙z łuku
O
X
, to wektor
ds
jest tak˙ze prostopadły do wektora
s
s
0
. Korzystuj ˛
ac z reguły prawej dłoni
mo˙zna przekona´c si˛e, ˙ze
ds
skierowany jest zgodnie z kierunkiem wektora
s
(s
s
0
)
. Mo˙zna te˙z
pokaza´c, ˙ze długo´s´c tego iloczynu wektorowego równa si˛e
sin
. St ˛
ad wektorowy odpowiednik
równania (2.23) ma posta´c
ds
=
k
s
(s
s
0
)
(2.27)
Równanie to jest bardziej ogólne ni˙z równanie (2.25), bowiem dotyczy ono dowolnego układu
współrz˛ednych, nie tylko układu równikowego.
2.7
Wektorowe transformacje współrz˛ednych sferycznych
Macierze obrotów i lustrzanych odbi´c
Prostok ˛
atne układy współrz˛ednych definiowane s ˛
a za pomoc ˛
a trójek wzajemnie ortogonalnych
jednostkowych wektorów np.
i;
j;
k
, patrz rysunek 2.8a. Tak ˛
a trójk˛e mo˙zna uj ˛
a´c w formie macierzy
R
zwanej triad ˛
a, o rozmiarach
3
3
R
=
[i;
j;
k℄
(2.28)
W jej kolumnach znajduj ˛
a si˛e kosinusy kierunkowe wersorów
i;
j;
k
. Wzajemna ortogonalno´s´c
tych wersorów poci ˛
aga ortogonalno´s´c triady
R
, co łatwo sprawdzi´c z pomoc ˛
a warunku
R
T
R
=
I
(2.29)
gdzie
I
jest macierz ˛
a jednostkow ˛
a.
Współrz˛edne wektora
a
mog ˛
a by´c podane wzgl˛edem dowolnych układów współrz˛ednych,
czyli innymi słowy wzgl˛edem dowolnych triad, np.
R
i
P
. Je´sli
[a
1
;
a
2
;
a
3
℄
s ˛
a współrz˛ednymi
wektora
a
wzgl˛edem triady
R
, mo˙zemy to samo wyrazi´c w postaci
a
=
R
2
4
a
1
a
2
a
3
3
5
(2.30)
Je´sli
[
1
;
2
;
3
℄
s ˛
a współrz˛ednymi wektora
a
wzgl˛edem triady
P
, to piszemy
a
=
P
2
4
1
2
3
3
5
(2.31)
Mno˙z ˛
ac lewostronnie równania (2.30) i (2.31) przez
R
T
widzimy, ˙ze
R
T
a
=
2
4
a
1
a
2
a
3
3
5
=
R
T
P
2
4
1
2
3
3
5
(2.32)
Z tego zapisu wynika, ˙ze transformacja składowych
[
1
;
2
;
3
℄
w składowe
[a
1
;
a
2
;
a
3
℄
mo˙ze by´c
dokonana za po´srednictwem iloczynu macierzy
R
T
P
. Transformacja odwrotna mo˙ze by´c doko-
nana za pomoc ˛
a iloczynu
P
T
R
, b˛ed ˛
acego transpozycj ˛
a poprzedniego iloczynu.
22
Poło˙zenia ciał na sferze
k
i
j
a
j
1
1
k
i
1
a)
j
k
θ
j
1
k
1
θ
cos
sin
θ
θ
b)
Rysunek 2.8: a) Wektor
a
mo˙ze mie´c składowe okre´slone wzgledem wielu trójek wersorów
i;
j;
k
.
b) Zwi ˛
azki pomi˛edzy elementami triad ró˙zni ˛
acych si˛e obrotem wokół jednej osi pokrywaj ˛
acej si˛e
z wersorem
i
. Wersor
j
1
mo˙zna otrzyma´c za po´srednictwem wersorów
j;
k
odkładaj ˛
ac w kierunku
j
odcinek o długo´sci
os
po czym w kierunku
k
odcinek
sin
.
7
Wyprowadzimy teraz formuły umo˙zliwiaj ˛
ace pewn ˛
a szczególn ˛
a transformacj˛e współrz˛ednych.
Mianowicie, interesuje nas przeliczenie współrz˛ednych z układu zdefiniowanego prostok ˛
atn ˛
a tri-
ad ˛
a
R
=
[i;
j;
k℄
do układu współrz˛ednych danego triad ˛
a
R
1
=
[i;
j
1
;
k
1
℄
, przy czym transforma-
cja
R
T
1
R
ma posta´c
i
1
=
i
j
1
=
( os
)j
+
(sin
)k
k
1
=
(
sin
)j
+
( os
)k
co oznacza, ˙ze oba układy ró˙zni ˛
a si˛e jedynie o dodatni obrót o k ˛
at
wokół osi
i
. A w postaci
dogodniejszej do umacierzowienia zapisu b˛edzie
i
1
=
1i
+0j
+0k
j
1
=
0i
+( os
)j
+(sin
)k
k
1
=
0i
+(
sin
)j
+( os
)k
Zatem, zgodnie z definicj ˛
a triady, triad˛e
R
1
mo˙zemy okre´sli´c jako iloczyn macierzowy
R
1
=
[i
1
;
j
1
;
k
1
℄
=
R
2
4
1
0
0
0
os
sin
0
sin
os
3
5
a wobec tego co pokazano wcze´sniej, transformacji składowych wektora wyznaczonych wzgl˛edem
R
do składowych podanych wzgl˛edem triady
R
1
da si˛e dokona´c za pomoc ˛
a macierzy
R
T
1
R
, czyli
R
1
T
R
=
2
4
R
2
4
1
0
0
0
os
sin
0
sin
os
3
5
3
5
T
R
=
2
4
1
0
0
0
os
sin
0
sin
os
3
5
Mamy zatem nast˛epuj ˛
acy wniosek. Transformacja współrz˛ednych (składowych) wektora
[x;
y
;
z
℄
T
wyznaczonych wzgl˛edem jednego układu, do współrz˛ednych
[x
1
;
y
1
;
z
1
℄
T
tego˙z wektora wzgl˛e-
dem układu powstałego przez obrót
p(
)
wokół osi
x
o dodatni k ˛
at
, mo˙ze by´c dokonana za
7
Jest to konsekwencj ˛
a ortogonalno´sci obu macierzy.
2.7 Wektorowe transformacje współrz˛ednych sferycznych
23
pomoc ˛
a formuły
2
4
x
1
y
1
z
1
3
5
=
p(
)
2
4
x
y
y
3
5
(2.33)
gdzie
p(
)
=
2
4
1
0
0
0
os
sin
0
sin
os
3
5
(2.34)
Analogiczne formuły daj ˛
a si˛e wyprowadzi´c dla układów ró˙zni ˛
acych si˛e obrotami o dodatni k ˛
at
wokół osi
y
i
z
. Odpowiednie macierze transformacyjne oznaczone jako
q(
)
i
r(
)
maj ˛
a posta´c
q(
)
=
2
4
os
0
sin
0
1
0
sin
0
os
3
5
(2.35)
r(
)
=
2
4
os
sin
0
sin
os
0
0
0
1
3
5
(2.36)
W przypadku gdy układy współrz˛ednych ró˙zni ˛
a si˛e dwoma lub trzema obrotami, macierz transfor-
macji współrz˛ednych wyznaczona jest jako iloczyn dwóch lub trzech macierzy odpowiadaj ˛
aych
pojedynczym obrotom.
Macierz transformacyjna dla układów ró˙zni ˛
acych si˛e skr˛etno´sci ˛
a (przypadek układów lewo i
prawoskr˛etnych) jest macierz ˛
a modyfikuj ˛
ac ˛
a jedynie współrz˛edn ˛
a y-ow ˛
a (lustrzane odbicie wzgl˛e-
dem płaszczyzny
x
z
), ma ona posta´c
M
y
=
2
4
1
0
0
0
1
0
0
0
1
3
5
(2.37)
W praktyce mo˙zemy napotka´c przypadek takich układów, dla których korzystnym b˛edzie zas-
tosowanie macierzy
M
x
, pozwalaj ˛
acej na lustrzane odbicie wzgl˛edem płaszczyzny
y
z
czyli na
zmian˛e znaku współrz˛ednej x-owej.
M
x
=
2
4
1
0
0
0
1
0
0
0
1
3
5
(2.38)
W celu praktycznego zastosowania formuł podanych wy˙zej musimy dysponowa´c wzorami umo˙zli-
wiaj ˛
acymi transformacj˛e współrz˛ednych sferycznych do współrz˛ednych prostok ˛
atnych i odwrot-
nie. Je´sli dane s ˛
a współrz˛edne sferyczne
(u;
v
)
poło˙zenia ciała niebieskiego,
u
jest współrz˛edn ˛
a
azymutaln ˛
a a
v
jest dopełnieniem do
90
Æ
odległo´sci biegunowej. Odpowiadaj ˛
ace im prostok ˛
atne
składowe wersora poło˙zenia tego ciała wyliczamy ze pomoc ˛
a wzorów
x
=
os
u
os
v
y
=
sin
u
os
v
z
=
sin
v
(2.39)
24
Poło˙zenia ciał na sferze
x
1
z
1
x
z
y
N
O
φ
ψ
θ
Rysunek 2.9: Ilustracja k ˛
atów Eulera. Układy współrz˛ednych o wspólnym pocz ˛
atku, mo˙zna trans-
formowa´c jeden w drugi za pomoc ˛
a trzech obrotów o k ˛
aty Eulera.
Zale˙zno´sci odwrotne maj ˛
a posta´c
v
=
ar sin
z
u
=
ar tan
y
x
(2.40)
przy czym w celu ustalenia wła´sciwej ´cwiartki k ˛
ata
u
musimy zastosowa´c stosown ˛
a procedur˛e
normuj ˛
ac ˛
a.
Transformacja współrz˛ednych z wykorzystaniem k ˛
atów Eulera.
Za pomoc ˛
a macierzy podanych wy˙zej mo˙zemy z powodzeniem przelicza´c współrz˛edne pomi˛e-
dzy dowolnymi układami. W zale˙zno´sci od potrzeby, formułu transformacyjne b˛ed ˛
a zło˙zeniami
obrotów w okół osi
X
;
Y
;
Z
. Jednak przy takim podej´sciu dla ka˙zdego przypadku, kolejno´s´c
obrotów musimy odgadn ˛
a´c sami.
Ale mo˙zliwe jest inne podej´scie, w którym wzajemna orientacja dwóch układów współrz˛e-
dnych okre´slona jest za pomoc ˛
a tzw. k ˛
atów Eulera, natomiast transformacja składa si˛e z trzech
obrotów wzgl˛edem osi w ustalonej, zawsze takiej samej kolejno´sci.
Na rysunku 2.9 widzimy dwa prostok ˛
atne układy współrz˛ednych zorientowane wzgl˛edem
siebie tak, ˙ze ˙zadna para osi tych układów nie jest do siebie wzajemnie równoległa, płaszczyzny
X
Y
tych układów przecinaj ˛
a si˛e wzdłó˙z kierunku
O
N
. Jest to najbardziej ogólny przypadek
je´sli chodzi o orientacj˛e układów.
Na rysunku 2.9 zaznaczono trzy k ˛
aty wykorzystywane do transformacji współrz˛ednych wyznac-
zonych wzgl˛edem obu układów, s ˛
a to:
k ˛
at
zawarty pomi˛edzy osiami
Z
i
Z
1
obu układów,
k ˛
at
zawarty pomi˛edzy osi ˛
a
X
i lini ˛
a
O
N
przeci˛ecia płaszczyzn
X
Y
obu układów,
k ˛
at
pomi˛edzy osi ˛
a
X
1
i lini ˛
a
O
N
, liczony jako dodatni od linii
O
N
do osi
X
1
.
Transformacja współrz˛ednych
[x;
y
;
z
℄
T
we współrz˛edne
[x
1
;
y
1
;
z
1
℄
T
, co łatwo odczyta´c z ry-
2.8 Dygresja na temat miary małych k ˛
atów
25
sunku 2.9, jest zło˙zeniem trzech obrotów
2
4
x
1
y
1
z
1
3
5
=
r(
)p(
)r()
2
4
x
y
z
3
5
(2.41)
A zatem je´sli tylko mamy do dyspozycji k ˛
aty Eulera okre´slaj ˛
ace wzajemn ˛
a orientacj˛e dwóch
dowolnych układów współrz˛ednych kartezja´nskich, to transformacja pomi˛edzy współrz˛ednymi
wzgl˛edem tych układów zawsze b˛edzie miała posta´c równania (2.41).
2.8
Dygresja na temat miary małych k ˛
atów
Analizuj ˛
ac teoretyczne problemy na sferze, wygodnie jest wyra˙za´c k ˛
aty w radianach. Kiedy jed-
nak wykorzystujemy rezultaty analiz, opłaca si˛e stosowa´c inne miary k ˛
atów, np. małe k ˛
aty takie
jak k ˛
at paralaksy najcz˛e´sciej podane s ˛
a w sekundach łuku.
Na podstawie znanej definicji mamy w przybli˙zeniu
1
r
ad
=
57
Æ
17
0
45
00
=
206265
00
(2.42)
Radiany wykorzystywane s ˛
a w przybli˙zeniach małych k ˛
atów, tzn. je´sli warto´sci małych k ˛
atów s ˛
a
podane w radianach, wówczas dopuszczalne s ˛
a nas˛epuj ˛
ace przybli˙zenia niektórych funkcji try-
gonometrycznych
sin
os
1
tan
(2.43)
Z równa´n (2.42), (2.43) wynika, ˙ze
sin
1
00
=
1
206265
(2.44)
Wyra˙zenie to znakomicie nadaje si˛e do zamiamy radianów na sekundy łuku. Np. je´sli zapis
00
oznacza liczb˛e sekund w małym k ˛
acie
, to
sin
00
206265
=
00
sin
1
00
Jak pami˛etamy, równania (2.25) na zmiany współrz˛ednych w rezultacie małego przesuni˛ecia s ˛
a
dokładne tylko do rz˛edu pierwszego, dlatego stosuj ˛
ac je nale˙zy liczy´c si˛e z bł˛edem
"
wynosz ˛
acym
w radianach
"
=
O
(k
2
)
=
O
((k
00
sin
1
00
)
2
)
W sekundach łuku
"
wynosi
8
"
00
=
O
(k
002
sin
1
00
)
Daje to po˙zyteczn ˛
a formuł˛e na oszacowanie dokładno´sci. Np. dla przemieszcze´n na sferze o
warto´sci
1
00
bł ˛
ad formuły pierwszego przybli˙zenia jest rz˛edu
5
10
6
sekundy łuku, co jest do
8
Kwadrat sinusa znikn ˛
ał gdy˙z trzeba było podzieli´c poprzednie wyra˙zenie przez
sin
1
00
.
26
Poło˙zenia ciał na sferze
zaniedbania w ka˙zdym przypadku. Dla przesuni˛e´c o wielko´sci
15
00
bł ˛
ad ten wynosi około
0:001
00
co jest jeszcze poni˙zej precyzji najdokładniejszych teleskopów astrometrycznych. Ale dla prze-
suni˛e´c rz˛edu jednej minuty łuku bł˛edy wynosz ˛
a około
0:02
00
co mie´sci si˛e ju˙z w zakresie precyzji
współczesnych obserwacji optycznych i radiowych. Formuły pierwszego rz˛edu s ˛
a bardzo przy-
datne, ale trzeba je stosowa´c z rozwag ˛
a i to ilo´sciow ˛
a.
Przypu´s´cmy, ˙ze wykorzystuj ˛
ac równania (2.25) stosowano sekundy łuku. Dla jasno´sci umaw-
iamy si˛e, ˙ze
k
jest w radianach a zapis
k
00
oznacza t ˛
a sam ˛
a wielko´s´c podan ˛
a w sekundach łuku.
Poniewa˙z stosujemy ten sam współczynnik zamiany jednostek do obu stron równa´n (2.25),
d
i
dÆ
otrzymamy od razu w sekundach zast˛epuj ˛
ac
k
przez
k
00
. St ˛
ad je´sli wyrazimy parametr
k
w
sekundach łuku, równania (2.25) w jednostkach praktycznych przyjm ˛
a posta´c
d
s
=
1
15
k
00
se
Æ
os
Æ
0
sin(
0
)
dÆ
00
=
k
00
(sin
Æ
os
Æ
0
os
(
0
)
os
Æ
sin
Æ
0
)
(2.45)
2.9
Zadanka
1. Podaj definicje nast˛epuj ˛
acych poj˛e´c: sfera niebieska, koło wielkie, koło małe, dwuk ˛
at sfer-
yczny, trójk ˛
at sferyczny, nadmiar sferyczny, kat sferyczny, k ˛
at dwu´scienny, k ˛
at trój´scienny.
2. Wykonaj pełne wyprowadzenie równa´n (2.17) i (2.22).
3. Drog ˛
a przestawie´n symboli
ab AB
C
wypisz pozostałe wzory cosinusów, wzory pi˛ecioele-
mentowe oraz wzory cotangensowe, tak by otrzyma´c komplet podstawowych formuł try-
gonometrii sferycznej.
4. Wykonaj pełne wyprowadzenie równa´n (2.25).
5. Poka˙z, ˙ze długo´s´c wektora
s
(s
s
0
)
wynosi
sin
.
6. Podaj lepsze przybli˙zenie zale˙zno´sci (2.42).
7. Udowodnij równo´s´c (2.44).
8. Poka˙z, ˙ze warto´s´c k ˛
ata sferycznego na sferze jednostkowej równa si˛e odległo´sci na powierz-
chni sfery pomi˛edzy biegunami kół wielkich tworz ˛
acych ten k ˛
at.
9. Dla ka˙zdego trójk ˛
ata sferycznego
AB
C
mo˙zna zdefiniowa´c tzw. trójk ˛
at biegunowy
A
0
B
0
C
0
.
Mianowicie:
A
0
jest biegunem boku
B
C
tak, ˙ze
AA
0
<
90
Æ
.
B
0
i
C
0
definiuje si˛e podobnie.
Poka˙z, ˙ze boki i k ˛
aty obydwu trójk ˛
atów zwi ˛
azane s ˛
a wzorami
A
0
=
180
Æ
a
a
0
=
180
Æ
A
B
0
=
180
Æ
b
b
0
=
180
Æ
B
C
0
=
180
Æ
0
=
180
Æ
C