Večkrat smo omenili, da enačba y = y(x) opisuje ravninsko krivuljo, če sta spremenljivki x in y dolžinski koordinati. Enačba z = z(x, y) pa na podoben način opisuje ploskev v prostoru. Čas je, da se opisa krivulj in ploskev lotimo sistematično (Descartes, Gauss).
Koordinate
Osnova za opis krivulj in ploskev z enačbami je "poimenovanje" vsake prostorske točke z njenimi koordinatami (x, y, z) v poljubno izbranem koordinatnem sistemu, katerega osi so med seboj pravokotne in umerjene v enakih dolžinskih enotah. Rečemo, da so to kartezične koordinate. Razdalja med dvema točkama potem znaša, po hipotenuznem izreku,
(16.1)
s2 = (x2 − x1)2 + (y2 − y1)2 + (z2 − z1)2 . |
Pametno je sistem izbrati tako, da bo enačba krivulje ali ploskve v njem čim bolj preprosta.
Enačba premice
Najpreprostejša "krivulja" je premica. Uteleša jo, na primer, brazda ladje, ki pluje po morju v stalni smeri φ glede na sever. Pot ne sme biti predolga, da se ne pokaže zakrivljenost morja. Koordinatni sistem postavimo v začetno pristanišče, ordinatno os y usmerimo proti severu in abscisno os x proti vzhodu. Enačba brazde-premice se potem glasi
(16.2)
y = kx |
k = tan φ. |
Slika 16.1
Premica. Najkrajša pot med dvema točkama v
prostoru.
Smerni koeficient k ima nazoren pomen: to je prirast ordinatne razdalje na prirast abscisne razdalje. Če je koeficient pozitiven, premica narašča, sicer upada.
Premico, ki ne gre skozi izhodišče, ampak seka ordinatno os v y0, opišemo kot (y − y0) = kx. Če seka abscisno os v točki x0, velja y = k(x − x0). Ako pa gre skozi točko (x0, y0), se enačba premice glasi (y − y0) = k(x − x0).
Parametrični zapis
Pri ladji, ki pluje z enakomerno hitrostjo, sta njeni koordinati enolično določeni s pretečenim časom t:
(16.3)
x = At |
y = Bt . |
Ladja zariše isto premico ne glede na to, kako hitro pluje oziroma kako hitro teče čas (to je ura, ki jo imamo). Zato bomo opustili časovne enote in uporabljali kar brezdimenzijska števila. Takšen "čas", ki zavzema vrednosti na intervalu (−∞, +∞), bomo poimenovali parameterski čas oziroma parameter in ga označevali kar s t. Vsaki vrednosti parametra ustreza natanko ena vrednost koordinat. Primerjava parametričnega in eksplicitnega zapisa pove k = B/A.
Enačba krožnice
Iz sive davnine je poznana krožnica: krivulja, katere vsaka točka je enako oddaljena od izbrane točke, središča. Že stara ljudstva so jo risala s količkom in vrvico pri gradnji kolib in obzornih krogov. Mi bomo postavili koordinatni sistem v središče kroga. Potem pove hipotenuzni izrek
(16.4)
x2 + y2 = r2 . |
V translatorno zamaknjenem koordinatnem sistemu pa ima središče kroga koordinati (x0, y0). Tedaj očitno velja (x − x0)2 + (y − y0)2 = r2.
Slika 16.2
Krožnica. Vsaka njena točka je enako oddaljena
od izbrane točke, središča.
Parametrični zapis
Tudi krožnico lahko opišemo parametrično. Spomnimo se enakomernega kroženja nihala (iz fizike), pa takoj uvidimo
(16.5)
x = r cos t |
y = r sin t , |
pri čemer leži parameter t na intervalu [0,2π].
Prečno presekano bambusovo steblo ima rob v obliki krožnice. Če ga presekamo poševno, pa je rob "raztegnjena" krožnica, elipsa. Kako bi tako elipso narisali na tleh? Prej ali slej – morda kot kakšen kraljevi vrtnar – odkrijemo postopek: središče kroga "raztegnemo" v dve središči, nanju privežemo vrv, jo nategnemo z risalnim količkom in začrtamo željeno krivuljo. Elipsa je s tem definirana kot množica točk, pri katerih je vsota razdalj do dveh izbranih točk, gorišč, konstantna.
Točko na polovici zveznice med obema goriščema poimenujemo središče elipse. Skozi središče potekata dva odlikovana premera: dolga os 2a in kratka os 2b. Razdaljo med središčem in (katerimkoli) goriščem poimenujemo ekscentričnost e. Ko je risalni količek v temenu velike osi, vidimo, da velja r1 + r2 = 2a. Ko je v temenu male osi, pa hipotenuzni izrek pove b2 + e2 = a2.
Slika 16.3
Elipsa. Vsota razdalj iz dveh izbranih točk, gorišč, je do
vsake njene točke enaka.
Enačba elipse
Koordinatni križ postavimo v središče elipse in ga zavrtimo tako, da njegove osi sovpadajo z veliko in malo osjo. Levo gorišče ima potem koordinato (−e, 0) in desno (+e, 0). Razdalji od gorišč do izbrane točke na elipsi znašata r12 = (x + e)2 + y2 in r22 = (x − e)2 + y2. Njuna vsota mora biti r1 + r2 = 2a in iz tega pogoja sledi, z nekaj računanja, enačba
(16.6)
x2 a2 |
+ | y2 b2 |
= 1 . |
Elipso v premaknjenem koordinatnem sistemu (oziroma premaknjeno elipso v obstoječem sistemu) pa opišemo z zamenjavo x → x − x0 in y → y − y0.
Parametrični zapis
Pri a = b preide elipsa v krog, kakor je tudi prav. Parametrični opis zato kar uganemo:
(16.7)
x = a cos t |
y = b sin t . |
Parameter t leži na intervalu [0, 2π]. Da je to res pravi opis, preverimo z vstavitvijo v implicitno enačbo.
Krogelno zrcalo (katerega presek je krožni lok) zbira vzporeden snop žarkov v goriščno točko, vendar samo tedaj, kadar je snop ozek. Bolj oddaljeni žarki se po odboju sekajo v gorišču, ki je bliže temenu. Morda obstaja kakšna krivulja, ki bi vse vzporedne žarke združevala v isti točki? Drugače povedano: tako krivuljo – parabolo – bi morale sestavljati točke, ki so enako oddaljene od premice in goriščne točke.
Slika 16.4
Parabola. Vsaka njena točka je enako oddaljena od izbrane
točke, gorišča, in od vodilne premice.
Enačba parabole
Postavimo koordinatni sistem tako, da bo premica "vodilja" vodoravna pri koordinati (0, −p/2). Gorišče je potem v točki (0, +p/2). Razdalja poljubne točke na iskani krivulji od gorišča je r12 = (y − p/2)2 + x2 in razdalja te točke od premice je r2 = |y + p/2|. Iz pogoja r1 = r2 sledi, z nekaj računanja,
(16.8)
2py = x2 . |
Parametrični zapis
Enačba ima obliko y ∝ x2. Spomnimo se, da prav takšna enačba opisuje tir kamna pri vodoravnem metu (v fiziki). Tam narašča vodoravna koordinata s časom in navpična s kvadratom časa, kar nas navede na naslednji parametrično zapis parabole z navpično simetrijsko osjo:
(16.9)
x = At |
y = Bt2 . |
Vstavitev polarnih enačb v implicitno enačbo pove 2p = A2/B.
Hodograf vektorja
Namesto s koordinatami lahko delamo z ustreznimi vektorji lege: r = (x, y). Razdaljo med dvema točkama potem zapišemo kot absolutno vrednost razlike dveh vektorjev: s = |r2 − r1|.
Parametrski zapis krivulje pove, kako se vsaka koordinata spreminja s časom: x = x(t) in y = y(t). To zapišemo v vektorski obliki kot
(16.10)
r(t) = (x(t), y(t)) . |
S časom se vektor spreminja – obrača, daljša in krajša – in s svojo konico zarisuje hodograf – krivuljo. Naraščajoči parameter t definira pozitivno smer gibanja po krivulji.
Odvodi po parametru
Kako se odvod ene koordinate po drugi izraža z odvodoma koordinat po parametru? Verižno pravilo pove dy/dt = (dy/dx) · (dx/dt), torej
(16.11)
dy dx |
= | y' x' |
. |
Odvod po parametru smo označili s črtico. Drugi odvod pa računamo takole. Posredno odvajamo (d/dx)(dy/dx) = (d/dt)(dy/dx) · dt/dx. Ker dy/dx = y'/x' in dt/dx = 1/x', velja
(16.12)
d2y dx2 |
= | x'y" − y'x" x'3 |
. |
Kako parametrizirati
Kako za funkcijo y = y(x) določiti parametrično obliko? Izberemo (skoraj) poljubno funkcijo x = x(t) in nato izračunamo y = y(x(t)). Očitno je možnosti za izbiro neskončno. Poiščemo takšno, da je rezultat najbolj preprost. Posebno zanimiva izbira je kar x = t. Tedaj velja r(x) = (x, y(x)). Parabolo, na primer, zapišemo kot r(t) = (At, Bt2) ali kot r(x) = (x, x2/2p). Očitno je parametrični zapis krivulje zelo nazoren in vsestranski.
Kako pa iz parametrične oblike x = x(t), y = y(t) določiti eksplicitno oziroma implicitno obliko funkcije? Iz prve in druge enačbe izrazimo t, ju izenačimo in dobimo iskano enačbo, ki jo po potrebi še preoblikujemo v lepšo obliko.
Ločni element
Prirast parametra za dt se odraža kot sprememba vektorja dr = (dx, dy) oziroma kot kratek kos krivulje, ločni element ds2 = dx2 + dy2.
Slika 16.5
Ločni element krivulje. Njegova dolžina je
limitno enaka spremembi vektorja lege.
Velja ds = |dr|. Enačbo delimo na obeh straneh z dt, pa dobimo
(16.13)
ds = |dr| = |r'| dt = √(x'2 + y'2) dt . |
Dolžina krivulje
Dolžina poti, ki jo zariše vektor med začetno in končno lego, znaša
(16.14)
s = ∫ √(x'2 + y'2) dt . |
Če je parameter koordinata t = x, pomeni odvajanje na parameter kar odvajanje na koordinato: x' = dx/dx = 1 in y' = dy/dx, torej ds = √(1 + y'2) dx.
Naravna parameterizacija
Dolžina krivulje od izbrane začetne točke naprej in nazaj je odličen parameter za opis krivulje. Krivulja je tedaj kot cesta, na kateri so v enakih dolžinskih presledkih postavljeni mejniki. Vsak tak mejnik ima svoje koordinate in krivuljo opišemo kot r(s) = (x(s), y(s)). Parameter je sedaj vezan zgolj na krivuljo in nič na okolico. Pri takšni parametrizaciji seveda velja x'2 + y'2 = 1 (črtica označuje odvod po parametru s).
Kako dolžinsko parametrizirati krivuljo, ki je podana s splošnim parametrom t? — Izračunamo dolžino vzdolž krivulje kot funkcijo časa s(t). — Izračunamo obratno funkcijo t(s). — Vstavimo jo v prvotno enačbo r(t(s)). Za krog, na primer, dobimo x = r cos (s/r) in y = r sin (s/r).
Tangenta
Smer krivulje v izbrani točki je podana z normaliziranim premikom
(16.15)
τ = | dr ds |
. |
Števec in imanovalec ulomka delimo z dt in dobimo enotni tangentni vektor r'/|r'|, to je
(16.16)
τ = | (x', y') √(x'2 + y'2) |
. |
Tangenta, na kateri leži enotni tangentni vektor, ima smerni koeficient k = y'/x'. Če se dve krivulji sekata, je kot med njunima tangentnima vektorja določen s skalarnim produktom τ1 .τ2 = cos φ.
Normala
S tangentnim vektorjem je definiran normalni vektor, ki stoji nanj pravokotno:
(16.17)
n = k × τ , |
pri čemer je k enotni vektor v smeri osi z. Normalni vektor dobimo s križnim množenjem vektorskega produkta (ali z množenjem z rotacijsko matriko za 90°):
(16.18)
n = | (−y',x') √(x'2 + y'2) |
. |
Normala, na kateri leži normalni vektor, ima smerni koeficient k = −x'/y'. To je negativna recipročna vrednost smernega koeficienta tangente.
Ukrivljenost
Koliko se zasuče enotni vektor preko dolžinskega elementa, je mera za lokalno ukrivljenost krivulje
(16.19)
K = | | dτ ds |
| . |
Izračunamo jo takole. — Vektor r(t) odvajamo po času posredno: r' = (dr/ds) · (ds/dt) in dobimo τv. — Vektor r' odvajamo po času posredno: r" = (d/ds)(τv) · (ds/dt), upoštevamo pravilo za odvod produkta in dτ/ds = Kn ter dobimo Kv2n + τdv/dt. — Izračunamo produkt r' × r" = Kv3τ × n. — Iz slednjega izrazimo K, pri čemer upoštevamo τ × n = k, in dobimo K = (r' × r")k/v3, torej:
(16.20)
K = | x'y" − y'x" (x'2 + y'2)3/2 |
. |
Enačbe za tangento, normalo in ukrivljenost se ustrezno poeneostavijo, če vzamemo t = x ali t = s. Ukrivljenost se, na primer, izrazi kot K = y" / (1 + y'2)3/2 oziroma kot K = √(x"2 + y"2) .
Krivinski radij
Ko izračunamo ukrivljenost krožnice z radijem R, dobimo v vsaki točki vrednost
(16.21)
K = | 1 R |
. |
Če je ukrivljenost krivulje K, zato rečemo, da je njen lokalni krivinski radij R = 1/K. Krivulja je lokalno "nerazločljiva" od takega "pritisnjenega" kroga. Pritisnjeni krog je lokalno enak krivulji v tem smislu, da imata enak "ničti", prvi in drugi odvod.
Slika 16.6
Krivinski radij krivulje. To je radij kroga,
ki se najtesneje prilega krivulji.
Invariante krivulj
Nekatere značilnosti krivulje so odvisne od njene lege v izbranem koordinatnem sistemu. Primer so nagibi tangent ali normal glede na abscisno ali ordinatno os. Pri vrtenju koordinatnega sistema se takšni nagibi ne ohranjajo. Po drugi strani pa je ukrivljenost v izbrani točki krivulje neodvisna od izbire koordinatnega sistema. Rečemo, da je to invariantna lastnost krivulje oziroma njena invarianta. Invariante se ne izražajo s koordinatami, marveč le z njihovimi diferenciali.
Ravnina
Ravnina, ki gre skozi izhodišče koordinatnega sistema, zareže v ravnini xz enotni vektor r1 = (cos θ1, 0, sin θ1). V ravnini yz zareže vektor r2 = (0, cos θ2, sin θ2). Poljubna linearna kombinacija teh dveh vektorjev r = Ar1 + Br2 je krajevni vektor do ustrezajoče točke na preučevani ravnini. Zapišimo to kombinacijo v komponentah. Iz prve enačbe x = A cos θ1 izrazimo A, iz druge y = B cos θ2 izrazimo B in oboje vstavimo v tretjo enačbo z = A sin θ1 + B sin θ2. Tako dobimo eksplicitno enačbo ravnine
(16.22)
z = k1x + k2y , |
pri čemer sta k1 in k2 smerna koeficienta, torej tangensa obeh naklonskih kotov θ1 in θ2.
Valj
Vodoravno krožnico x2 + y2 = r2 premikamo v navpični smeri. Pri tem zariše plašč valja. Enačba zanj je kar enaka enačbi krožnice:
(16.23)
x2 + y2 = r2 . |
Stožec
Premico z = kx zavrtimo okrog navpične osi z. Nobeni točki se pri tem koordinata z ne spreminja, njena koordinata x pa prehaja v koordinate ρ = √(x2 + y2). Enačbo z = kρ kvadriramo in dobimo enačbo stožca
(16.24)
z2 k2 |
= x2 + y2 . |
Krogla
Krožnico x2 + z2 = r2 zavrtimo okoli navpične osi z. Transformacija x2 → x2 + y2 da enačbo krogle
(16.25)
x2 + y2 + z2 = r2 . |
Rotacijski elipsoid
Elipso x2/a2 + z2/c2 = 1 zavrtimo okrog navpične osi z. Dobimo rotacijski elipsoid
(16.26)
x2 a2 |
+ | y2 a2 |
+ | z2 c2 |
= 1 . |
Rotacijski paraboloid
Parabolo 2pz = x2 zavrtimo okrog navpične osi z. Nastane rotacijski paraboloid
(16.27)
2pz = x2 + y2 . |
Vse zapisane enačbe veljajo v posebno skrbno izbranih sistemih. Tako so tudi enačbe preproste. Seveda pa lahko koordinatni sistem translatorno premaknemo, kar je isto, kot da premaknemo ploskev v nasprotni smeri. Premik vzdolž osi z, na primer, je ekvivalenten transformaciji spremenljivke z → z − z0. Enačba se temu ustrezno "pogrša". Še hujše lepotne spremembe dosežemo z rotacijo.
Izbrane ploskve smo zapisali implicitno ali eksplicitno. Pojavi se vprašanje, ali (in kako) jih lahko zapišemo parametrično oziroma vektorsko. Poskusimo z najpomembnejšo ploskvijo, kroglo.
Točka na krogli radija R je enolično določena z vektorjem lege r = (x, y, z). Komponente vektorja izrazimo, kot že znamo, z azimutnim kotom φ in s polarnim kotom θ (14.2):
(16.28)
x = R sin θ cos φ |
y = R sin θ sin φ |
z = R cos θ . |
Parametrska ravnina
Vsaki dvojici kotov torej pripada ustrezna trojica koordinat. Na podoben način se lotimo tudi drugih ploskev. Valj in stožec, na primer, parametriziramo z azimutnim kotom in višino. Ne predivje ploskve nasploh opišemo z dvema parametroma:
(16.29)
r(u, v) = (x(u,v), y(u,v), z(u,v)) . |
Parametra sta lahko karkoli. V posebnem primeru izberemo kar dve koordinati: r = (x, y, z(x,y)). Tedaj preide parametrični opis v eksplicitnega. Hkrati nam ponudi še naslednjo nazorno sliko: dva splošna parametra tvorita posebno "parametrično" ravnino. Točke te ravnine se preslikajo v točke na aktualni ploskvi.
Krivulja (u(t), v(t)) v parametrični ravnini se preslika v ustrezno krivuljo na ploskvi. Poseben primer je preslikava, ko je eden izmed parametrov konstanten, recimo v = const. Tedaj se na ploskvi zariše ena izmed izo-parametričnih krivulj. Pri različnih vrednostih konstante se nariše množica takih krivulj – krivočrtnih koordinat na ploskvi. Tako se na krogli, na primer, zarišejo poldnevniki φ = const in vzporedniki θ = const.
Slika 16.7
Parcialna premika na ploskvi. To sta
prirastka vektorja lege vzdolž krivočrtnih koordinat na
ploskvi.
Parametrski kot
Vektorja ru in rv ležita v lokalni tangentni ravnini vzdolž obeh krivočrtnih koordinat. Kakšen je sekalni kot teh koordinat, α, pove skalarni produkt:
(16.30)
cos α = | ru · rv |ru∥rv| |
. |
Pri lepo izbranih parametrizacijah je kot v vsaki točki (morda s kakšno izjemo) enak 90°. Tedaj so krivočrtne koordinate med seboj pravokotne. Takšni so poldnevniki in vzporedniki na krogli.
Dolžinski element
V tangentni ravnini leži tudi totalni diferencial – "poševni" premik dr = rudu + rvdv. S kvadratom tega premika je določena njegova dolžina ds2 = dr · dr, torej:
(16.31)
ds2 = ru2du2 + 2rurvdudv + rv2dv2 = |
g11du2 + 2g12dudv + g22dv2 . |
V komponentah zapišemo
(16.32)
g11 = xu2 + yu2 + zu2 |
g12 = xuxv + yuyv + zuzv |
g22 = xv2 + yv2 + zv2 . |
Koeficienti g11, g12 in g22 so realna števila. Vsaka točka na ploskvi ima svojo trojico teh števil. Rečemo, da so to metrični koeficienti ploskve. Njihov pomen je, da diferenciale parametrov "povežejo" z diferenciali dolžin. Če izberemo drugačno parametrizacijo ploskve, se metrični koeficienti seveda spremenijo. V pravokotni koordinatni mreži je koeficient g12 = 0. Za kroglo v standardni parametrizaciji izračunamo g11 = r2 sin2 θ in g22 = r2. Za valj pa g11 = 1 in g22 = r2.
Dolžina krivulje na ploskvi je limitna vsota vseh dolžinskih diferencialov, torej (če označimo odvod po času s črtico)
(16.33)
s = ∫ √(g11u'2 + 2g12u'v' + g22v'2)dt . |
Geodetke
Med dvema oddaljenim točkama A in B na ploskvi poteka neskončno mnogo krivulj. Ena od njih je najkrajša. Rečemo, da je to geodetka. Na krogli je geodetka glavni krog, to je tak, ki ima središče v središču Zemlje. Nazorno si geodetko predstavljamo kot elastično nit, napeto med obema točkama: elastičnost jo skrči na najkrajšo dolžino.
Ploščinski element
Dolžinska elementa vzdolž krivočrtnih koordinat, pravokotnih ali ne, sta (ds)u = √g11du in (ds)v = √g22dv. Ploščina paralelograma, ki ga zamejujeta, pa znaša
(16.34)
dS = (ds)u (ds)v sin α = √(g11g22 − g122)dudv . |
Ploščina ploskve je limitna vsota ploščinskih elementov, torej dvojni integral
(16.35)
S = ∫∫ √(g11g22 − g122) dudv . |
Za parametra, ki sta kar koordinati, se enačba poenostvi v obliko
(16.36)
S = ∫∫ √(1 + zx2 + zy2) dxdy . |
Normala
Tangentna vektorja ležita v tangentni ravnini. Njun vektorski produkt je pravokoten nanjo. Če ga normiramo, dobimo normalo
(16.37)
n = | ru × rv |ru × rv | |
= | ru × rv √(g11g22 − g122) |
. |
Za parametra, ki sta kar koordinati, se enačba zapiše v obliki
(16.38)
n = | (−zx, −zy, 1) √(1 + zx2 + zy2) |
. |
Odmik tangentne ravnine
Vektor iz opazovane točke v bližnjo okolišnjo točko na ploskvi, torej vektor r(u + du, v + dv) − r(u,v), aproksimiramo s potenčno vrsto z linearnim členom (rudu + rvdv) in s kvadratnim členom 1/2 · (ruudu2 + 2ruvdudv + rvvdv2). Prvi člen je pomik po tangentni ravnini. Drugi člen je pomik do pritisnjenega kroga v smeri pravokotno na krog. Če ga pomnožimo z normalo, dobimo pravokotno razdaljo od tangentne ravnine:
(16.39)
2dh = L11du2 + 2L12dudv + L22dv2 |
L11 = ruu · n |
L12 = ruv · n |
L22 = rvv · n . |
Za kroglo v standardni parametrizaciji izračunamo L11 = r sin2 θ in L22 = r. Za valj pa velja L11 = 0 in L22 = r.
Slika 16.8
Odmik ploskve od tangentne ravnine. Limitno je enak odmiku
pritisnjene paraboloidne ploskve.
Ukrivljenost
Ukrivljenost ploskve je enaka ukrivljenosti pritisnjenega kroga: dh = ds2/2R, torej 1/R = 2dh/ds2, zato:
(16.40)
K = | L11du2 + 2L12dudv + L22dv2 g11du2 + 2g12dudv + g22dv2 |
. |
To je ukrivljenost ploskve v smeri, ki jo določata du in dv. Skozi izbrano točko potekajoče krivulje imajo večjo ali manjšo ukrivljenost. Izmed njih ima ena maksimalno ukrivljenost Kmax = 1/Rmin in druga minimalno Kmin = 1/Rmax. Najdemo ju kot ekstremalne vrednosti po vseh smereh. V to se ne bomo spuščali. Ko takšni vrednosti najdemo, se lahko igramo z njunima vrednostima: tvorimo, na primer, "povprečno" ukrivljenost K = (Kmax + Kmin)/2 ali "metrično" ukrivljenost K = Kmax · Kmin ter poskušamo najti, kako se izražata s koeficienti g11 … L22. Tudi to zahtevno zabavo prepustimo drugim, ki jih to zanima.
Poglejmo še nekaj zgledov. Ravnina ima v vsaki točki vse ukrivljenosti nič. Zato ji tudi rečemo ravnina. Na krogli so poldnevniški krivinski radiji večji od vzporedniških. Vsi glavni krogi skozi vsako točko pa imajo enak radij, ki je enak poldnevniškemu. Najmanjši krivinski radij na valju je enak polmeru valja in največji je neskončen. Podobno je pri stožcu. Vidimo, da se da marsikaj dognati tudi brez računanja.
Na majhnih razdaljah je zemeljska površina ravna in koti, premice in trikotniki na njej se pokoravajo že spoznanim pravilom, recimo pravilu o vsoti notranjih kotov v trikotniku ali hipotenuznemu pravilu o razdalji med dvema točkama. Na večjih razdaljah pa je treba upoštevati zemljino zakrivljenost. "Ravne" premice na njej postanejo glavni krogi. Vsota notranjih kotov trikotnika postane večja od 180°, kar se lepo vidi na primeru trikotnika z bazo na ekvatorju in vrhom na polu. Hipotenuzni, kosinusni in sinusni izrek za trikotnike pa bo treba na novo premisliti.
Dolžina geodetke
Za lažje preučevanje bomo vse dolžine na krogli merili z radijem kot enoto. S tem postane radij brezdimenzijska količina z velikostjo 1, dolžinski odsek vsakega glavnega kroga pa številsko enak središčnemu kotu, v radianih, na katerega je napet. Prvo vprašanje, ki si ga zastavimo, je: kolikšna je dolžina glavnega kroga med dvema točkama?
Na točki naj kažeta vektorja r1(θ1, φ1) in r2(θ2, φ2) iz središča krogle. Njuna velikost je enaka ena. Kot med njima, torej brezdimenzijska dolžina glavnega kroga, je določen s skalarnim produktom r1 · r2 = cos α. Zmnožimo komponente, upoštevamo kosinus razlike in dobimo
(16.41)
cos α = sin θ1 sin θ2 cos (φ2 − φ1) + cos θ1 cos θ2 . |
Razdalja, v dolžinskih enotah, je potem d = R α. Poseben primer φ1 = φ2 pove dolžino poldnevnika: α = |θ2 − θ1|, kakor tudi mora biti.
Hipotenuzni izrek
Pravokotni trikotnik na krogli določajo trije enotni vektorji iz njenega izhodišča do trikotnikovih oglišč. Vseeno je, kako je koordinatni sistem postavljen. Izberemo ga tako, da kaže vektor r1 vzdolž osi x, vektor r2 leži v ekvatorski ravnini xy pod dolžinskim kotom a in vektor r3 leži v poldnevniški ravnini pod širinskim kotom h. Vektorji so torej naslednji: r1 = (1, 0, 0), r2 = (cos a, sin a, 0) in r3 = (cos a cos h, sin a sin h, sin h). Kot d med r1 in r3 je hipotenuza trikotnika in je določen s skalarnim produktom cos d = r1 · r2. Pomnožimo komponente in dobimo hipotenuzni izrek
(16.42)
cos d = cos a cos h. |
Pri kratkih stranicah aproksimiramo cos x ≈ 1 − x2/2, zanemarimo visoke potence in izrek preide v ravninskega.
Kosinusni izrek
Podobno se lotimo poševnega trikotnika na krogli. Omejimo se na "prave" trikotnike, katerih koti so manjši od π in katerih stranice so tudi manjše od π.
Slika 16.9
Poševni trikotnik na krogli.
(Mercator, 2013)
Na tri oglišča trikotnika kažejo vektorji OA, OB in OC. Koordinatni sistem usmerimo, kot kaže slika. V njem velja OA = (0, 0, 1) in OB = (sin c, 0, cos c). Vektor OC se projicira v ON pod kotom A, torej OC = (sin b cos A, sin b sin A, cos b). Skalarni produkt OB · OC = cos a. Zmnožimo komponente in dobimo:
(16.43)
cos a = cos b cos c + sin b sin c cos A . |
Stranica a je podana z drugima dvema stranicama in kotom med njima. To je iskani kosinusni izrek. Velja seveda za vsakršno permutacijo zapisanih količin. Opazimo tudi, da je kosinusni izrek povsem enak izrazu za dolžino geodetke (16.41). To pa ni nič čudnega, saj je slednji le poseben primer prvega: za glavne kroge uporablja poldnevnike in ekvator.
Pri majhnih razdaljah aproksimiramo sin x ≈ x in cos x ≈ 1 − x2/2, zanemarimo visoke potence in izrek preide v ravninskega. V posebnem primeru, ko A = 90°, je trikotnik pravokoten in izrek se reducira v hipotenuzni izrek.
Sinusni izrek
Ideniteta sin2 A = 1 − cos2 A nas navede na misel, da vanjo vstavimo cos A iz kosinusnega izreka in upamo, da se bo izcimil sinusni izrek. Res pridelamo izraz sin A/sin a = f(a, b, c). Desna stran izraza je invariantna glede na ciklično permutacijo stranic, kar pomeni, da mora veljati
(16.44)
sin A sin a |
= | sin B sin b |
= | sin C sin c |
. |
To je sinusni izrek. Pri majhnih razdaljah preide v že znano ravninsko obliko.
Hipotenuzni, kosinusni in sinusni izrek nam pomagajo pri računanju kotov in stranic na krogli točno na tak način, kot to počnemo v ravnini. Ko delamo s kroglo polmera R namesto 1, moramo vse stranice trikotnika, podane v dolžinskih enotah, deliti z R. Drugače rečeno: namesto brezdimenzijske stranice a moramo povsod pisati a/R in podobno za druge stranice.
Projekcija krogle
Točke na zemeljski površini so enolično določene s svojimi zemljepisnimi koordinatami: širino δ (oziroma polarnim kotom θ = π/2 − δ) in dolžino φ. Zemljo verodostojno predstavimo s pomanjšanim krogelnim modelom. Takšen globus pa je, žal, neprimeren za prenašanje in tudi ni dovolj velik za podroben prikaz manjših območij. Naravno je torej pomisliti, kako bi ga preslikali – v celoti ali deloma – na ravno ploskev, zemljevid. Iščemo torej primerne preslikave
(16.45)
(x, y) ← (θ, φ). |
Rečemo jim zemljepisne projekcije.
Slika 16.10
Preslikava krogle na ravnino s svetlobnimi žarki.
Oblika sence je zanimiva tudi za slikarje.
(Rubens, 1613)
Napake projekcij
Vsaka preslikava, ki jo vpeljemo, preslika Zemljine poldnevnike in vzporednike v dve družini ravninskih krivulj. Dva bližnja vzporednika in poldnevnika na Zemlji oblikujeta ploščinski element, približni pravokotnik. Ko se takšen pravokotnik preslika, pričakujemo naslednje nevšečnosti: kot med stičnima stranicama se spremeni; razmerje med tema stranicama se spremeni; enaki pravokotniki na različnih lokacijah se preslikajo neenako, bodisi po dolžini, širini ali ploščini. Seveda hočemo najti take preslikave, ki bodo obremenjene s čim manj nevšečnostmi. Posebej pomembno je, da se ohranjajo lokalni koti, to je lokalna razmerja stranic. Tedaj se oblika in orientacija drobnih likov pri preslikavi ohranja. Drobni krogi se, na primer, preslikajo kot krogi. Takim preslikavam rečemo konformne.
Preslikava z žarki
Preslikajmo severno poloblo na tangentno ravnino na severnem polu! Preslikujemo lahko z žarki, ki izhajajo is središča krogle, iz njenega južnega pola ali iz južne neskončnosti. V vsakem primeru se Zemljini poldnevniki preslikajo v radialne premice, vzporedniki pa v koncentrične kroge. Razdalja med krogi je odvisna od izbire žarkov. Središčni žarki "preveč" raztegnejo ekvatorske predele, neskončni pa jih "preveč" stisnejo. Osredotočimo se torej na južni pol kot izvor žarkov. To je polarna stereografska projekcija.
Slika 16.11 Polarna stereografska projekcija. Projekcija je primerna za prikaz polarnih dežel, pa tudi za prikaz zvezdnega neba.
Polarni izvor žarkov
Slika pokaže, da se točka P(θ) preslika v točko P'(ρ). Ker je obodni kot enak polovici središčnega, razberemo iz pravokotnega trikotnika SNP' povezavo
(16.46)
ρ = 2R tan | θ 2 |
. |
Za radij Zemlje izberemo primerno pomanjšano vrednost: R = M · RE, na primer M = 1:107. Namesto polarnega kota θ lahko uporabimo tudi zemljepisni kot δ = 90° − θ. V tangentni ravnini vpeljemo koordinatni sistem z izhodiščem v polu; os y kaže vzdolž poljubnega poldnevnika φ0. Potem velja
(16.47)
x = ρ sin (φ − φ0) |
y = −ρ cos (φ − φ0) . |
S tem je preslikava zaključena. Seveda ni treba projicirati celotne hemisfere, ampak le kakšen njen del. Tedaj na tangentni ravnini vpeljemo lokalni koordinatni sistem, ki je glede na polarnega ustrezno translatorno zamaknjen.
Konformnost
Je projekcija morda konformna? Ploščinski element na krogli je približno pravokotnik z vzporedniško stranico R sin θ dφ in s poldnevniško stranico R dθ. Ustrezajoči ploščinski element v tangentni ravnini je tudi približno pravokotnik s stranicama ρdφ in dρ. Z razmerjem istoležnih stranic sta podana raztezna faktorja H = ρdφ / R sin θ dφ in K = dρ / R dθ. Če je preslikava konformna, mora veljati H = K. Izračunamo odvod dρ/dθ in ga vstavimo v enačbo. Pokaže se, da je kvocient razteznih faktorjev enak 1. Preslikava je povsod konformna.
Polarna stereografska projekcija je primerna za prikaz dežel v visokih zemljepisnih širinah, pa tudi za prikaz zvezdnega neba.
Morska navigacija
Ko mora ladja pluti iz kraja A v oddaljeni kraj B, ima na voljo neomejeno mnogo poti. Če odmislimo tokove, vetrove in neurja, je najboljša pot tista, ki je najkrajša, torej geodetka, to je glavni krog na krogli. Takšna geodetka je na polarni stereografski projekciji v splošnem krivulja, ki seka poldnevnike pod različnimi koti. Določiti in zarisati jo brez obsežnega računanja ni možno. Pa tudi sledenje taki črti bi zahtevalo, da krmar stalno spreminja magnetni kurz ladje.
Druga možnost je krivulja, ki seka vse poldnevnike pod istim kotom – loksodroma. Je sicer daljša od geodetke, vendar je za krmarjenje mnogo bolj primerna. Seveda je tudi loksodroma kriva črta na polarni stereografski projekciji (razen če pluje ladja po poldnevniku). Kaj ne bi bilo čudovito, če bi imel navigator na mizi zemljepisno projekcijo, na kateri bi bila loksodroma povsod ravna črta? Med krajema A in B bi potegnil ravno črto in s tem določil kurz ladje. Bolj preprosto ne gre. Poizkusimo, kot navigatorji, najti tako projekcijo!
Slika 16.12 Ekvatorska valjna konformna projekcija. Projekcija je primerna za prikaz ekvatorskih dežel in za pomorsko navigacijo.
Valjna projekcija
Da bo loksodroma ravna, morajo očitno biti poldnevniki ekvidistantne premice, vzporedniki pa nanje pravokotne premice v takšnih medsebojnih razmakih, da je mreža povsod lokalno konformna. To pomeni, da moramo projicirati kroglo na valj, ovit okoli njenega ekvatorja. Valj se seveda da razviti v ravnino. Na valju postavimo koordinatni sistem z osjo x vzdolž ekvatorja in y vzdolž poljubnega poldnevnika φ0. Točke s poldnevnika φ se vse preslikajo v
(16.48)
x = R (φ − φ0) . |
Vpeljava konformnosti
Pri tem se točke iz različnih širin θ preslikajo v ustrezne y, kakor določa zahteva po konformnosti. Ravnamo tako kot pri polarni stereografski projekciji. Izenačimo raztezna faktorja H = dx/R sin θ dφ in K = dy/R dθ. V dobljeni enačbi sta vsebovana dva odvoda. Prvega dx/dφ zlahka izračunamo in s tem je določen drugi: dy/dθ = R / cos θ. Ločitev spremenljivk in integracija pove
(16.49)
y = R ln tan | θ 2 |
. |
Razmiki med vzporedniki torej naraščajo z oddaljenostjo od ekvatorja. To je tudi pričakovati, saj projekcija na silo širi in paralelizira krogelne poldnevnike. Seveda ni treba projicirati celotne krogle, ampak le kakšen njen del. Tam postavimo lokalni koordinatni sistem, ki je ustrezno translatorno premaknjen.
Ekvatorska valjna konformna projekcija je odlična za pomorsko navigacijo in primerna za prikaz dežel v nizkih zemljepisnih širinah.
Razrast industrializacije, širjenje železniškega in cestnega omrežja ter nenehna vojskovanja zahtevajo natančne zemljevide velikih držav. Pokaže se potreba po ustrezni projekciji za srednje zemljepisne širine. Smer raziskave je hitro pri roki: zemeljsko kroglo je treba projicirati na plašč stožca, ki se je dotika v izbranem vzporedniku. Poldnevniki so tedaj radialne premice, vzporednike – koncentrične kroge – pa želimo razmestiti tako, da bo projekcija konformna. Tako kot valj lahko tudi stožec nato razvijemo v ravnino.
Slika 16.13 Stožčna konformna projekcija. Projekcija je primerna za prikaz dežel v zmernem pasu.
Razvoj stožca v ravnino
Naj se stožec dotika vzporednika δ0 = π/2 − θ0, ki je za ρ0 oddaljen od vrha stožca. Vrhnji polkot stožca je potem tudi enak δ0. Obseg stožca po tem vzporedniku znaša L1 = 2π ρ0 sin δ0. Ko plašč stožca razvijemo v ravnino, nastane izsekan krog, katerega celotni obseg je L2 = 2π ρ0. Razmerje teh dveh obsegov L1 / L2 = k = sin δ0. (Spomnimo se na stožčaste šotore, tipije, prerijskih severnoameriških domorodcev! Plašč tipija je točno polovica kroga: k = 1/2. To pomeni, da ima vrhnji polkot δ0 = 30°.)
V izsekani krog vpeljimo ravninski koordinatni sistem z vrhom v presečišču tangentnega vzporednika in poljubnega poldnevnika φ0. Os x je usmerjena vzdolž vzporednika in os y vzdolž poldnevnika. Krogelni poldnevnik φ postane na razvitem plašču stolpca poldnevnik kφ.
Vpeljava konformnosti
Ploskovni element na razvitem plašču stožca ima vzporedniško stranico ρ kdφ in poldnevniško stranico dρ, s čimer sta določena raztezna faktorja glede na ploskovni element na krogli. Izenačitev razteznih faktorjev vodi do enačbe dρ/ρ = kdθ/sin θ. Integriranje obeh strani da rešitev
(16.50)
ρ = C tank | θ 2 |
k = sin (π/2 − θ0). |
Konstanto C določimo iz razteznega pogoja: ρ(θ0) = R tan θ0. S tem sta določeni tudi koordinati
(16.51)
x = ρ cos k(φ − φ0) |
y = ρ0 − ρ sin k(φ − φ0) . |
Vzdolž tangentnega vzporednika so razdalje točne. Če za tangentni vzporednik izberemo pol, preide stožec v tangentno ravnino in projekcija v polarno stereografsko. Če za tangentni vzporednik izberemo ekvator, pa preide stožec v valj in projekcija v ekvatorsko valjno konformno.
Stožčna konformna projekcija je dobra za prikaz dežel na srednjih zemljepisnih širinah.
Različice projekcij
Vsaka izmed obravnavanih tipov projekcij – ravninska, valjna in stožčna – ima več različic. Če, na primer, razvrstimo vzporednike na enake medsebojne razdalje, dobimo ekvidistantne projekcije. Razdalje vzdolž poldnevnikov so tedaj pravilne. Spet drugače izbrana razvrstitev poldnevnikov pa zagotovi, da so pravilne ploščine. To so ekvivalentne projekcije. Jasno je, da spremenjene projekcije niso več konformne.
Zemlja kot rotacijski elipsoid
Zemlja je krogla le v prvem, čeravno zelo dobrem približku. Tisti, ki želijo večjo natančnost, jo aproksimirajo z rotacijskim elipsoidom s kratko polosjo med poloma. Projekcijske enačbe se močno zapletejo in vprašanje je, kdaj jih je sploh smiselno uporabljati. Sploščenost Zemlje je namreč zelo majhna: (a − b) / a ≈ 1 / 300.
Globalne projekcije
Nobena izmed naštetih projekcij ni primerna za prikaz celotne zemeljske oble. Obliko velikih in "oddaljenih" kontinentov namreč močno popačijo. So pa ljudje iznašli mnogo kar sprejemljivih globalnih projekcij. Žal to, da je teh projekcij mnogo, pove, da nobena ni povsem zadovoljujoča. Ena izmed boljših je eliptična projekcija z naslednjimi značilnostmi. Slika sveta je elipsa z razmerjem polosi 1:2. Ekvator in vzporedniki so vzporedne daljice z enakomernim presledkom. Centralni poldnevnik je daljica. Vsi drugi so polelipse, simetrične glede na ekvator in na centralni poldnevnik. Polelipsi skozi ±90° tvorita krog. Projekcijski obrazci so ustrezno zamotani in jih ne bomo izpeljevali. □