background image

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

background image

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

.

background image

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.

background image

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.

background image

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

background image

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

background image

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 ˛

punktu

A

,





jest współrz˛edn ˛

a polarn ˛

punktu

A

, jest ona identyczna z k ˛

atem

Z

O

A

,



 

jest współrz˛edn ˛

a azymutaln ˛

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.

background image

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-

background image

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.

background image

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

.

background image

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 ;

Æ

+

)

. Poszukamy wyra˙ze´n pozwalaj ˛

a-

cych na obliczenie przyrostów

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

:

:

:

.

background image

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

Æ

+

)

, 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

os

+

)



d

os

Æ

Drugi bok

U

X

=

, 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



=

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

)

=

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.

background image

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.

background image

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.

background image

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)

background image

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-

background image

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

.

background image

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

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

)

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

background image