Grafični prikaz gibanja satelita okrog Zemlje

V poljubno verzijo Pythona dodamo katero od grafičnih knjižnic. Prav preprosta je graphic.py . Namestimo jo na disk tako, da postane Pythonu vidna, preberemo še navodila na začetku knjižnice in lahko začnemo s programiranjem grafike. Tako recimo koda

from math import *
from graphics import *
def delay(m):
    for i in range(1000*m):
            continue
def main():
    mx=600    # širina in višina okna
    my=400
    win=GraphWin("Moj Krog",mx,my,autoflush=False)
    p=Rectangle(Point(0,0),Point(mx,my))
    p.setFill("white")
    c=Circle(Point(mx/2,my/2),10)      #Zemlja je modri krogec
    c.setFill("blue")
    p.draw(win)
    c.draw(win)
    dt=10                  #interval med računi leg
    x=mx/3                 #začetna lega
    y=0
    vx=0                   # začetna hitrost
    vy=0.04                # to komponento malo spremeni 
    while True:
        r=sqrt(x*x+y*y)   # račun razdalje satelit-Zemlja
        ax=-x/(r*r*r)     # pospešek satelita sledi iz 
        ay=-y/(r*r*r)     # gravitacijskeg azakona
        vx=vx+ax*dt       # račun nove hitrosti
        vy=vy+ay*dt       
        x=x+vx*dt         #račun nove lege satelita 
        y=y+vy*dt
        t=Point(mx/2+x,my/2-y)     # risanje satelita
        t.setFill("red")
        t.draw(win)
        delay(100)
        #t.setFill("white")
        #t.draw(win)
        update()
    win.getMouse()
    win.close()
main()

spravimo v gibanje satelit okrog Zemlje.

PTR v srednji šoli(9)

Lorenzove transformacije lahko zapišemo v kompaktnejši matrični obliki:

\[
\begin{bmatrix}ct\\x\end{bmatrix}=\begin{bmatrix}\gamma&\gamma\beta\\\gamma\beta&\gamma\end{bmatrix}\begin{bmatrix}ct^\prime\\x^\prime\end{bmatrix}
\]

V njej nastopa Lorenzova matrika

\[
\begin{bmatrix}\frac{1}{\sqrt{1-\beta^2}}&\frac{\beta}{\sqrt{1-\beta^2}}\\\frac{\beta}{\sqrt{1-\beta^2}}&\frac{1}{\sqrt{1-\beta^2}}
\end{bmatrix}\quad\quad\quad(1)
\]

Prvo koordinato v levem vektorju enačbe (1) dobimo tako, da skalarno pomnožimo 1. vrstico matrike z desnim vektorjem in podobno tudi 2. koordinato. Pred matriko je relativistični faktor.

Opazimo, da se s svetlobno hitrostjo c pomnoženi čas v zapisu obnaša tako kot koordinata x. Če pišemo še koordinati y in z, ki sta prečni na smer gibanja, dobimo

\[
\begin{bmatrix}
ct\\x\\y\\z
\end{bmatrix}=
\begin{bmatrix}
\gamma& \gamma\beta&0&0\\\gamma\beta&\gamma&0&0\\0&0&1&0\\0&0&0&1
\end{bmatrix}\begin{bmatrix}
ct^\prime\\x^\prime\\y^\prime\\z^\prime\end{bmatrix}
\]

Seveda tudi tu velja, da dobimo i-to komponento levega vektorja tako, da skalarno pomnožimo i-to vrstico matrike z desnim vektorjem. Še obratna Lorenzova transformacija:

\[
\begin{bmatrix}ct^\prime\\x^\prime\\y^\prime\\z^\prime
\end{bmatrix}=
\begin{bmatrix}
\gamma& -\gamma\beta&0&0\\-\gamma\beta&\gamma&0&0\\0&0&1&0\\0&0&0&1
\end{bmatrix}\begin{bmatrix}ct\\x\\y\\z
\end{bmatrix}
\]

Relacije nam ponujajo odgovor na vprašanje, kaj je čas. Čas je pač ena od koordinat štirirazsežnega prostora-časa. Vektorju

\[
\begin{bmatrix}
ct\\x\\y\\z
\end{bmatrix}\]

pravimo dogodek  v prostoru- času. Lorentzove transformacije nam pomagajo preračunavati dogodke iz enega v drug inercialni sistem v prostoru-času.

PTR v srednji šoli(4)

NOVE TRANSFORMACIJE

Iščemo torej linearno transformacijo, ki prevede točko [math](u,v)[/math] v točko [math](u^\prime,v^\prime)[/math] tako, da  velja zveza
[math]u^{\prime 2}-v^{\prime 2}=u^2-v^2.[/math]
Ker je transformacija linearna,  jo iščemo v obliki
[math]u=Au^\prime +Bv^\prime \qquad v=Cu^\prime +Dv^\prime ,[/math]
pri čemer so A, B,C in D konstante, ki jih je treba določiti. Vstavimo zato te transformacije v zgornjo enačbo, pa dobimo
[math](Au^\prime+Bv^\prime )^2-(Cu^\prime+Dv^\prime )^2=u^{\prime2}-v^{\prime 2}.[/math]
Po kvadriranju in primerjanju koeficientov dobimo naslednje enačbe
[math]A^2-C^2=1,\quad AB=CD, \quad C^2-D^2=-1.[/math]
Imamo torej tri enačbe in štriri neznanke. Zato uvedemo parameter
[math]\beta=\frac{C}{A}=\frac{B}{D}[/math]
ter z njim izrazimo vse koeficiente. Dobimo
[math]A=D=\frac{1}{\sqrt{1-\beta^2}},\quad B=C=\frac{\beta}{\sqrt{1-\beta^2}}[/math]
Iskane transformacije so torej
[math]u=\frac{u^\prime+\beta v^\prime}{\sqrt{1-\beta^2}},\quad v=\frac{\beta u^\prime+v^\prime}{\sqrt{1-\beta^2}},[/math]
obratne transformacije pa
[math]u^\prime=\frac{u-\beta v}{\sqrt{1-\beta^2}},\quad v^\prime=\frac{-\beta u+v}{\sqrt{1-\beta^2}}.[/math]
Naslednjič pa jim bomo dali fizikalni pomen.

 

PTR v srednji šoli(3)

Oba inercialna sistema so povezovale Galilejeve transformacije, a te smo v zadnjem poglavju razglasi za neveljavne. Potrebujemo torej nove transformacije.

Vprašamo se, ali je kaj, kar opišeta enako oba,  sprevodnik in postajenačelnik.  Odgovor nam spet prinese naslednji miselni poskus:

Mislimo si, da se vlak s sprevodnikom približuje postajenačelniku. Koordinatna sistema obeh imata vzporedne osi, vlak in z njim gibajoći se koordinatni sistem pa se giblje  v smeri  postajenačelnikove x-osi.  V trenutku, so oba koordinatna sistema pokrijeta, na vlaku zasveti okrogla luč..  Sprevodnik  opiše svetlobo kot krogelni val, ki se širi s hitrostjo svetlobe c na vse strani, torej

\[x^{\prime 2}+y^{\prime 2}+z^{\prime 2}=c^2t^{\prime 2}.\].

Postajenačelnik pa val opiše podobno, a v svojih koordinatah

\[x^{2}+y^{2}+z^{ 2}=c^2t^{2}.\]

Vidimo torej, da sta opisa enaka. Velja torej

\[x^{\prime 2}+y^{\prime 2}+z^{\prime 2}-c^2t^{\prime 2}=x^{2}+y^{2}+z^{ 2}-c^2t^{2}.\].

Ker se vlak giblje prečno na koordinate \[y,y^\prime, z, z^\prime,\] je

\[y=y^\prime,\quad z=z^\prime,\]

pa se zgornja zveza še poenostavi v

\[x^{\prime 2}-c^2t^{\prime 2}=x^{2}-c^2t^{2}.\].

Iščemo torej linearno transformacijo, ki bo zadostila zgornji enačbi.

PTR v srednji šoli (2)

Morda ste uganili, kaj je bilo treba prečrtati – Galilejeve transformacije. Na prvi pogled je nenavadno, da ne veljajo enačbe, ki so se do takrat izkazale za dobro preizkušene. A vendar imamo sedaj razmere, ki so posebne – zelo velike hitrosti. Hitrosti, ki niso majhne v primeri s hitrostjo svetlobe.  Galilejeve transformacije dobro veljajo pri majhnih hitrostih, pri ekstremnih hitrostih pa odpovejo.

Če Galilejeve transformacije ne veljajo več, potem čas ni več absoluten, temveč relativen – odvisen od opazovalnega sistema.  To pa ima v primerjavi z dosedašnjim gledanjem na svet nenavadne posledice. Oglejmo si eno od njih v naslednjem (miselnem) poskusu:

SOČASNOST JE RELATIVNA

Imejmo vlak, ki je zelo dolg in se zelo hitro giblje. Njegova dolžina naj bo 300.000km, njegova hitrost pa  4c/5. Vlak ima dvoje vrat, ki se odpirata takrat, ko žarek svetlobe iz žarnice na sredini vlaka posveti na fotocelico na vratih. Sprevodnik stoji na sredini vlaka in s stikalom odpira vrata. Ker vrata  glede na sprevodnika mirujejo, le-ta opiše odpiranje vrat takole: Od trenutka, ko posveti žarnica potuje svetloba proti zadnjim vratom in po

\[t_1=\frac{150000km\cdot s}{300000km}=0,5s\]

se le-ta odprejo. Ravno tako ponovi račun za odpiranje prvih vrat

\[t_2=\frac{150000km\cdot s}{300000km}=0,5s.\]

Oba časa sta enaka, za sprevodnika se oboja vrata torej odprejo sočasno.

Dogajanje na drvečem vlaku pa opazuje tudi postajenačelnik, ki vidi dogajanje  nekoliko drugače : Medtem, ko potuje svetloba s hitrostjo c proti zadnjim vratom, se ji le-ta približujejo s hitrostjo 4c/5. Zato se zadnja vrata odprejo po času

\[t_1=\frac{150000km\cdot s}{300000km+240000km}=\frac{5}{18}s,\]

prva vrata pa se žarku odmikajo, zato je

\[t_2=\frac{150000km\cdot s}{300000km-240000km}=2,5s.\]

Postajenačelniku se torej vrata ne odprejo istočasno, zadnja vrata se odprejo prej kot prva. Dogodka, ki sta sočasna v enem opazovalnem sistemu, nista sočasna v drugem. Pravimo, da je sočasnost dogodkov relativna, torej odvisna od opazovalnega sistema.