02 Położenia ciał na sferze

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 ˛

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.

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

Wyszukiwarka

Podobne podstrony:
02 OGÓLNY POGLĄD NA ZDROWIE I CHOROBĘid 3432 ppt
Sciaga, Odległość sferyczną między dwoma punktami leżącymi na sferze nazywamy kąt środkowy oparty na
Kopia Sciaga, Odległość sferyczną między dwoma punktami leżącymi na sferze nazywamy kąt środkowy opa
TPL WYK 13 02 25 Płyn na odciski
Wpływ położenia budynku na zapotrzebowanie energii
TPL PRAC 13 02 25 Płyn na odciski
02 Wykonywanie obliczeń na liczbach przybliżonych, proste zadania(1)id 3492
02 06 Roboty na dachach v1 1
Jednostka wobec położenia narodu na przykładzie Konrada Wallenroda i Kordiana
02.Instrukcja higieny na st. pracy, Haccp-Dokumentacja-przykład
NaTropie 8'02 pazdziernik, ZHP, Na tropie, Rocznik 2002
NaTropie 3'02 maj, ZHP, Na tropie, Rocznik 2002
NaTropie 7'02 wrzesien, ZHP, Na tropie, Rocznik 2002
Ruch roczny Słońca na sferze niebieskiej, StUdiA
2014-15-TRD-T1 02 Pomiar-ruchu-na-skrzyżowaniu formularz, Budownictwo UTP, III rok, VI semestr - DUL
NaTropie 5'02 lipiec, ZHP, Na tropie, Rocznik 2002
02, Wykonywanie obliczeń na liczbach przybliżonych, proste zadania(1)
02 07 Roboty na duzych wysokosc Nieznany (2)

więcej podobnych podstron