Programiranje
Objavljeno: 22.8.2007 14:50
ahoj!
med programiranjem sem naletel na precej zanimivo obnasanje.
Zanimalo me je cist poreprosto kako se giblje satelit v sferno simetričnem potencialu
\(V = \frac{GMm}{r}.\)
Vzel sem preprosto Eulerjevo integracijsko metodo
\(v_{i}^{j} = v_{i}^{j-1} + a_{i} dt\)
\(x_{i}^{j} = x_{i}^{j-1} + v_{i}^{j} dt\)
kjer indeks \(i\) oznacuje koordinato (recimo i=0, za x, ipd.) in indeks \(j\) zaporedno vrednost koraka.
Najprej delal sem v sfericnih koordinatah in za orbiti okoli ekvatorja in pola je vse vredu. Ko pa sem nagnil tirnico glede na ekvator za 30 stopinj pa sem dobil za orbito cudne vzorce in ne kroznice kot pricakovano.
Zatem se enak poskus ponovil v karteznicnih koordinatah, kjer sem pa dobil pricakovani razultat,
Zanima me ali je kdo ze naletel na kaj podobnega in ce mogoce pozna vzrok?
Ce koga zanima lahko pripnem tukaj se kodo programa (je v C jeziku).
lp vid
med programiranjem sem naletel na precej zanimivo obnasanje.
Zanimalo me je cist poreprosto kako se giblje satelit v sferno simetričnem potencialu
\(V = \frac{GMm}{r}.\)
Vzel sem preprosto Eulerjevo integracijsko metodo
\(v_{i}^{j} = v_{i}^{j-1} + a_{i} dt\)
\(x_{i}^{j} = x_{i}^{j-1} + v_{i}^{j} dt\)
kjer indeks \(i\) oznacuje koordinato (recimo i=0, za x, ipd.) in indeks \(j\) zaporedno vrednost koraka.
Najprej delal sem v sfericnih koordinatah in za orbiti okoli ekvatorja in pola je vse vredu. Ko pa sem nagnil tirnico glede na ekvator za 30 stopinj pa sem dobil za orbito cudne vzorce in ne kroznice kot pricakovano.
Zatem se enak poskus ponovil v karteznicnih koordinatah, kjer sem pa dobil pricakovani razultat,
Zanima me ali je kdo ze naletel na kaj podobnega in ce mogoce pozna vzrok?
Ce koga zanima lahko pripnem tukaj se kodo programa (je v C jeziku).
lp vid