nauka.102markov,
HI!!!
Trenutno radim na jednom projekticu i potrebna mi je pomoc.
Tema rada mi je: Sistemi 2 i 3 tela (naravno gravitacioni).
Koristim kompis da mi iterativnom metodom racuna pozicije
tela, ali je do problema doslo kada je trebalo naci jednacine koje
bi korak po korak opisivale pozicije u sistemu 3 tela koja
medjusobno reaguju.
Svaka pomoc je veoma dobrodosla!!!!!!
Unapred hvala ,
Marko
nauka.103kriss,
-> #101, dr.grba˙˙ Gde šta? Može li to malo jasnije? Ili je greška u tvom execu?
˙˙ (;
˙˙
˙˙ BTW, eno u EXTRA:zvezde opet pričaju o brzini svetlosti i
˙˙ problemima kada se zdravim razumom pokušava da shvati
˙˙ relativistička mehanika. <───┐
│
To sam sve ja započeo ;)) ───────┘
žovek je pitao (a i ja pa niko nije odgovorio :( ) gde može da se nađe
elektronsko ;) izdanje Galaksije - koliko se sećam bilo je reči da je
na SEZAM-u?
nauka.104ilazarevic,
-> #102, markov> tela, ali je do problema doslo kada je trebalo naci jednacine koje
> bi korak po korak opisivale pozicije u sistemu 3 tela koja
> medjusobno reaguju.
Da li radiš kretanje po Keplerovoj jednačini, ili po masama?
nauka.105dejanr,
-> #103, kriss>> žovek je pitao (a i ja pa niko nije odgovorio :( ) gde može da se nađe
>> elektronsko ;) izdanje Galaksije - koliko se sećam bilo je reči da je
>> na SEZAM-u?
Bilo je na Sezamu, dok je bilo reči o tome da će "Galaksija" redovno
izlaziti u tom obliku. Pošto od toga ne bi ništa, taj podsistem je
ukinut.
nauka.106markov,
-> #104, ilazarevicKretanje racunam preko Njutnovog zakona grav. , sto ce reci
preko masa.
Imam problema sa izvlacenjem jednacina koje bi mi dale potrebne
podatke za racunanje pozicija kada su tri tela u pitanju i kada
ona reaguju medju sobom (ili bar kada dva "satelita" reaguju medju sobom).
Pomoc vrlo dobrodosla.
Marko
nauka.107ilazarevic,
-> #106, markov> Kretanje racunam preko Njutnovog zakona grav. , sto ce reci
> preko masa.
> Imam problema sa izvlacenjem jednacina koje bi mi dale potrebne
> podatke za racunanje pozicija kada su tri tela u pitanju i kada
> ona reaguju medju sobom (ili bar kada dva "satelita" reaguju medju
> sobom).
Ja sam za maturski radio keplerovsko kretanje, bez učešća masa. Problem ti je
dosta širok - ne znam tačno gde si zaglavio u tim jednačinama. Pretpostavljam
da ti treba rezultujuća sila, pa da onda iz nje računaš položaje kroz vreme.
Da li u tvom radu postoji mogućnost zadavanja bilo kojih masa i rastojanja
tela, da bi se videlo šta će se desiti, ili radiš samo sa stabilnim
sistemima?
nauka.108dikla,
-> #99, kriss> ˙˙ Ima ga na Sezamu !
>
> Gde?
Pokusao sam da ga nadjem ali mi nije uspelo, mozda je skinut ..:(
nauka.109dikla,
-> #105, dejanr> Bilo je na Sezamu, dok je bilo reci o tome da ce "Galaksija" redovno
> izlaziti u tom obliku. Posto od toga ne bi nista, taj podsistem je
> ukinut.
A ja bio ubedjen da je bio okacen uz neku poruku ..:)
nauka.110srdjan.j,
-> #105, dejanr==> Bilo je na Sezamu, dok je bilo reči o tome da će "Galaksija"
==> redovno izlaziti u tom obliku. Pošto od toga ne bi ništa, taj
==> podsistem je ukinut.
Može li da se to izdanje vrati na SEZAM, bar kao datoteka, kada je ono
izašlo, ja nisam bio član SEZAM-a :(.
Srđan.
nauka.111maksa,
-> #106, markov>> Imam problema sa izvlacenjem jednacina koje bi mi dale
>> potrebne podatke za racunanje pozicija kada su tri tela u
>> pitanju i kada
Može malo detalja ?
nauka.112mirijevac,
Evo, ovako.....
Daje se na znanje, da ce Vladimir Sukalovic, aka Johnny ili Bandit
BRANITI SVOJ DIPLOMSKI RAD na temu:
INTERAKCIJE FOSFOLIPIDA I HEMOGLOBINA U HIPERINSULINIZMU
u PETAK, 1.7.94 god. u Svecanoj Sali, Hemijskog Fakulteta...
U 10:30 h <<<<<<<<<<<< ulaz slobodan...
mirijevac
nauka.113dr.grba,
-> #112, mirijevac>> Daje se na znanje, da ce Vladimir Sukalovic, aka Johnny ili Bandit
Hoće li doći baroness? <:
nauka.114mirijevac,
-> #113, dr.grba
Pa sada, pise ulaz slobodan....
nisam samo siguran da bi Johnny bio cool da mu se tamo
pojavljuju bivse...
ma, dodji pa cemo da ga vidimo, cuo sam da ce biti koktel posle
odbrane... ;))
Mirijevac
nauka.115ognjen,
-> #112, mirijevac)-> Daje se na znanje, da ce Vladimir Sukalovic, aka Johnny ili
)-> Bandit BRANITI SVOJ DIPLOMSKI RAD na temu:
Da li je taj Bandit neki rod Karadjordjevica kada u tri
konferencije saljes jednu te istu poruku?
nauka.116m.hristodulo,
-> #102, markov>> Koristim kompis da mi iterativnom metodom racuna
>> pozicije tela, ali je do problema doslo kada je
>> trebalo naci jednacine koje bi korak po korak
>> opisivale pozicije u sistemu 3 tela koja
>> medjusobno reaguju.
Ne mogu da se setim tačno, ali mi se čini da je pokojni Mika na
jednom predavanju uzgred rekao da ne postoji u zatvorenom i
konačnom obliku jednačina rezultujućih međusobnih gravitacionih
sila po masi, kada više od dva tela slobodno interaguju
međusobno...
Moraćeš da upotrebiš neku približnu metodu. Mislim da je pokojni
(bokte, sve pokojni, :-( ) dr. Slavić imao nekojiko članaka u
Računarima o približnim metodama u astronomiji.
nauka.117mirijevac,
-> #115, ognjen)-> Daje se na znanje, da ce Vladimir Sukalovic, aka Johnny ili
)-> Bandit BRANITI SVOJ DIPLOMSKI RAD na temu:
Da li je taj Bandit neki rod Karadjordjevica kada u tri
konferencije saljes jednu te istu poruku?
U cetiri konferencije....
i ako mene pitas, hteo sam da je dam kao info kada se
logujes...
Mirijevac
nauka.118markov,
-> #107, ilazarevicHI!
Vec si rekao na koji nacin pronalazim pozicije tela i iz cega
one proizilaze.
Posto ovaj rad treba da ima najvise 4000 reci , sto ce reci
10-15 strana , nisam siguran da li uopste treba da se petljam u
problem tri tela. Uzeo sam neke knjige i kada sam video na kom
je to nivou zamal' sto nisam capnuo. Ovaj rad je reda velicine
naseg maturskog , ili malo jace.
Sve u svemu , ja preko jednacina koje imam izvodim racunicu za
silu , pa preko sile za ubrzanje , pa preko ubrzanja za promenu
polozaja. Na kraju to sve ispadne lepo i krasno :))
Ovde prilazem programcic koji to sve lepo izvodi u delo (u C-u).
Inace ja zadajem pocetne brzine, a odnos (G*M*m)/m je uzet da
bude jedan , pa jednacine mogu da budu malo zbunjujuce :)
Meni u stvari treba predlog sta sve da radim sa tim sistemom
od dva (ili eventualno tri tela). Ja mogu da proverim
Keplerove zakone kroz graficku analizu , mogu da isprobavam
razne uslove sistema i kako ce se on ponasati (granicne brzine) ,
mogu da napravim grafik energije sistema .
Da li ima predloga sta bi tu jos moglo da se ubaci i gde bi
mogla da se nabavi literatura iz koje bi ja mogao , zajedno sa
svojim znanjem, da ubacim "meso" (tj. 4000 reci , bre :) ) .
Ako ima nekih pitanja i predloga veoma sam TU.
Marko
Sledi program:
#include<stdio.h>
#include<math.h>
#include<graphics.h>
#include<conio.h>
#define faktor 10
#define tacke 0.1
int promena;
int ko,GraphDriver,GraphMode,ErrorCode;
char *stri;
main()
{
int i;
double u,u2,k=0,x2=1,x=.5,y2=0,y=0,Vx2=0,Vx=0,Vy=1.7,Vy2=1.2,Ax2,Ay2,r1,r2,Ax;
double r,Ay,
clrscr();
GraphDriver = DETECT;
initgraph( &GraphDriver, &GraphMode, "c:\tc" );
line(getmaxx()/2,0,getmaxx()/2,getmaxy());
line(0,getmaxy()/2+2,getmaxx(),getmaxy()/2+2);
Vx=Vx+(-4)*(tacke/2);
for (i=0;i<150;i++){
r1=sqrt(((x2-x)*(x2-x))+((y2-y)*(y2-y)));
if (i<56) {
x=x+Vx*k;
y=y+Vy*k;
r=sqrt(x*x+y*y);
u=1/(r*r*r);
Ax=-(x*u);
Ay=-(y*u);
Vx=Vx+Ax*k;
Vy=Vy+Ay*k;
k=tacke;
{
int c,b;
c=((floor(x*1000)+ceil(x*1000))/(2*faktor)+getmaxx()/2);
b=(-(floor(y*1000)+ceil(y*1000))/(2*faktor)+getmaxy()/2);
putpixel(c,b,12);
}}
getch();
}
P.S. Ako ti nije problem , da li bi mogao nekako da mi doturis taj tvoj
maturski da vidim ima li ista sto bi mi pomoglo.
nauka.119ilazarevic,
-> #118, markov> P.S. Ako ti nije problem , da li bi mogao nekako da mi doturis taj tvoj
> maturski da vidim ima li ista sto bi mi pomoglo.
Doturiću ti ga, ali sumnjam da će ti išta koristiti. Ja sam radio samo preko
Keplerove jednačine, a ona je cela smeštena u jednoj proceduri. Ostali,
nesređeni deo programa, čine meniji, grafika... Imaš ga na mail.