----------------------------------- kepler 30 Mai 2008 09:43 despre un punct material ce se misca liber sub influentsa c. ----------------------------------- va rog sa-mi scrie cine poate, sa-mi dee macar un indiciu unde pot sa gasesc ecuatia unui pct material ce se misca liber sub influentsa cimpului gravitational creat de o sursa pct-forma. e relativ important pentru mine. cine poate discuta pe tema asta va rog sa imi scrie cit mai repede... ----------------------------------- Doru Dragan 30 Mai 2008 19:43 ----------------------------------- http://scienceworld.wolfram.com/physics/Two-BodyProblem.html http://en.wikipedia.org/wiki/Three_body_problem Gasesti de asemenea informatiile astea in cartea UNIVERSUL a lui V. Ureche aparuta cam demult. Din pacate e vorba de matematica destul de serioasa :( Am un program BASIC care calculeaza efemeridele planetelor si care include formulele astea. ----------------------------------- kepler 30 Mai 2008 19:51 ----------------------------------- va sunt foarte recunoscator p-u ajutorul acordat sper ca voi gasi ceia ce ma intereseaza in matematika "serioasa" care mi-o propunetsi, defapt la facultate sunt profesori care explika dak rogi frumos. chestia e in faptul ca am facut si eu nuste studii amanuntite (dupa mine) care mi-au dat ca rezultat o functie graficul cariua vb cam straniu de aia vreau sa ma consult cu diferite materiale, inca odata multsumesc. ----------------------------------- Andi 30 Mai 2008 20:04 ----------------------------------- http://scienceworld.wolfram.com/physics/Two-BodyProblem.html http://en.wikipedia.org/wiki/Three_body_problem Gasesti de asemenea informatiile astea in cartea UNIVERSUL a lui V. Ureche aparuta cam demult. Din pacate e vorba de matematica destul de serioasa :( Am un program BASIC care calculeaza efemeridele planetelor si care include formulele astea. Pe mine mă interesează programul in BASIC. Dacă se poate lista ??? Mii de mulțumiri !!! ----------------------------------- Doru Dragan 30 Mai 2008 20:37 ----------------------------------- Programul e adaptat dupa o fotocopie pe care am avut-o a unui articol din revista Ciel et Espace. Liniile nu sunt numerotate. Merge in GWBASIC si TurboBasic. Pe margine sunt si niste comentarii lamuritoare. Daca nu te descurci poti sa mai intrebi. :) Daca gasesti erori sau faci imbunatatiri, m-as bucura sa-mi comunici. DIM V#(9,9):DEF FNF(X)=60*(X-INT(X)): PI#=4*ATN(1):GR#=PI#/180 DEF FNS(X)=ATN(X/SQR(1-X*X)):DEF FNC(X)=PI#/2-FNS(X) DEF FNG$(X)=STR$(INT(X))+STR$(INT(FNF(X)))+STR$(INT(FNF(FNF(X)))) DATA MERCUR,.24085,231.2973,77.1442128,.00004253 DATA .2056306,.3870986,7.0043579,48.0941733,.00003244 DATA VENUS,.61521,355.73352,131.289579,.00003757 DATA .0067826,.7233316,.394435,76.4997524,.00002502 DATA TERRA,1.00004,98.83354,102.5964032,0,.016718,1,0,0,0 DATA MARTE,1.88089,126.30783,335.690817,.00005038 DATA .0933865,1.5236883,1.8498011,49.4032001,.00002127 DATA JUPITER,11.86224,146.966365,14.0095493,.00004408 DATA .0484658,5.202561,1.3041819,100.252018,.00002767 DATA SATURN,29.45771,165.322242,92.6653974,.00005362 DATA .0556155,9.554747,2.4893741,113.4888341,.00002391 DATA URANUS,84.01247,228.070855,172.736329,.00004064 DATA .0463232,19.21814,.7729895,73.8768642,.00001365 DATA NEPTUN,164.79558,260.357899,47.8672148,.000039 DATA .0090021,30.10957,1.7716017,131.560649,.00003009 DATA PLUTO,247.686,209.439,222.972,0,.25387,39.78459,17.137,109.941,0 98 FOR P = 1 TO 9: READ P$(P): FOR I = 1 TO 9: READ V#(P, I): NEXT I: NEXT P 100 INPUT "ZI,LUNA,AN:"; Z, L, A: IF L < 3 THEN A = A - 1: L = L + 12 130 J# = INT(365.25 * A) + INT(30.6001 * (L + 1)) + Z + 1720994.5# 140 IF A > 1582 THEN J# = J# - INT(A / 100) + INT(A / 400) + 2 150 NJ# = J# - 2444238.5# ' Epoca 1980 200 T# = V#(3, 1): L# = V#(3, 2): S# = V#(3, 3): E# = V#(3, 5): A# = 1: IT = 0: OT = 0 210 GOSUB 700: TL# = I#: RT# = R# ' Calcul pt. Terra 300 FOR P = 1 TO 9: IF P = 3 THEN 460 302 T# = V#(P, 1): L# = V#(P, 2): S# = V#(P, 3) + NJ# * V#(P, 4): E# = V#(P, 5): A# = V#(P, 6) 304 IP# = GR# * V#(P, 7): OP# = V#(P, 8 ) + NJ# * V#(P, 9) 306 GOSUB 700: PL# = I#: RP# = R#: X# = GR# * (PL# - OP#)' Calcul pt. planeta 310 HP# = FNS(SIN((PL# - OP#) * GR#) * SIN(IP#)) ' Latitudinea heliocentrica 320 PL# = ATN(TAN(X#) * COS(IP#)) 330 IF COS(X#) < 0 THEN PL# = PL# + PI# 335 IF PL# < 0 THEN PL# = PL# + 2 * PI# 340 PL# = PL# * 180 / PI# + OP#: IF PL# > 360 THEN PL# = PL# - 360 ' Longitudinea heliocentrica 345 RP# = RP# * COS(HP#) ' Raza vectoare 347 REM Coordonate ecliptice geocentrice - Longitudine 350 IF P < 3 THEN LE# = ATN(RP# * SIN((TL# - PL#) * GR#) / (RT# - RP# * COS((TL# - PL#) * GR#))) + PI# + TL# * GR# 360 IF P > 3 THEN LE# = ATN(RT# * SIN((PL# - TL#) * GR#) / (RP# - RT# * COS((PL# - TL#) * GR#))) + PL# * GR# 365 REM Latitudine ecliptica geocentrica 370 LA# = ATN(RP# * TAN(HP#) * SIN((LE# * 180 / PI# - PL#) * GR#) / (RT# * SIN((PL# - TL#) * GR#))) 400 REM Conversie coordonate ecliptice => coordonate ecuatoriale 410 Q# = 23.44133493#: DE# = FNS(SIN(LA#) * COS(Q# * GR#) + COS(LA#) * SIN(Q# * GR#) * SIN(LE#)) 420 ZS# = ATN(COS(Q# * GR#) * TAN(LE#) - SIN(Q# * GR#) * TAN(LA#) / COS(LE#)) 430 IF COS(LE#) < 0 THEN ZS# = ZS# + PI# 440 IF ZS# < 0 THEN ZS# = ZS# + 2 * PI# 445 DE# = DE# * 180 / PI#: ZS# = ZS# * 12 / PI# 450 PRINT P$(P),"Asc.:";FNG$(ZS#);" ";"Dec.:";FNG$(DE#) 460 NEXT P 470 END 700 REM Calculul anomaliei medii si adevarate 710 N#=(360*NJ#)/(365.2422*T#):M#=N#+L#-S#:IF M#<0 THEN M#=M#+360 720 IF M#>360 THEN M#=M#-(360*INT(M#/360)) 721 GOTO 740 722 I#=N#+L#+(360*E#*SIN(M#))/PI#:IF I#>360 THEN I#=I#-(360*INT(I#/360)) 723 IF V#<0 THEN V#=V#+360 724 V#=I#-S#:GOTO 770 740 M#=M#*GR#:K#=M# 750 FOR W=1 TO 5:K#=M#+E#*SIN(K#):NEXT W 760 V#=SQR((1+E#)/(1-E#))*TAN(K#/2):V#=2*ATN(V#) 770 R#=(A#*(1-E#*E#))/(1+E#*COS(V#)):V#=V#*180/PI#:I#=V#+S# 780 RETURN ----------------------------------- Andi 30 Mai 2008 20:46 ----------------------------------- Este interesant. Îl pun la treabă. Mulțumesc. ----------------------------------- kepler 30 Mai 2008 22:03 ----------------------------------- Domnul Dragan dak stiti solutia ecuatiei diferentiale din pagina web care mi-ati dat-o va rog sa mi-o transmitetsi ----------------------------------- Doru Dragan 31 Mai 2008 07:59 ----------------------------------- Ei mai baiete, daca as sti ... n-as mai fi aici cu voi pe acest forum :) Am facut facultatea de matematica acum ... 30 de ani. Si atunci eram lenes :( si cam chiulangiu. Iar acum daca ma pui sa ridic pixul de pe masa mi se pare un efort deosebit :D Ti-am dat un punct de pornire. Am impresia ca tu esti mai tanar decat mine. Hai, fa un efort :lol: ----------------------------------- kepler 31 Mai 2008 08:51 ----------------------------------- Eu defapt creadeam ca si la chestia asta aveti un link:) de aceia si v-am rugat, ei dar referitor la rezolvare, nustiu, am rugat doi profesori sa o rezolve unul de matematica si unul de fizica, cel de matematica a rezolvat-o dar a obtinut o solutie generala, dar ce rost din ea dak trebuie problema Coshi de rezolvat, ei shi, de treaba asta s-a apukat prof-ul de fizica. a spus ca pina pe 9 o rezolva, asa ca eu am vrut mai repede:) dak gasesc ceva interesant in ecuatie numaidecit voi scrie aici :) ----------------------------------- Doru Dragan 31 Mai 2008 09:37 ----------------------------------- Ia-o usurel :) Asemenea probleme nu pot fi rezolvate cum bati din palme. Parerea mea este sa nu te repezi prea tare si sa nu astepti cine stie ce. Ecuatiile astea au fost intoarse si pe-o fata si pe alta, studiate de-a fir a par (cum s-ar zice). La ora actuala toata teoria este cuprinsa in asanumitul model VSOP (Variations Seculaires des Orbites Planetaires) pe care-l poti consulta aici: http://adsabs.harvard.edu/abs/1988A&A...202..309B ----------------------------------- kepler 31 Mai 2008 10:30 ----------------------------------- ok ms orikum nu e ceia ce cautam eu vorbeam de un caz general, nasa vad ca rezova probleme in cadrul sistemului solar, fiind strins in timp de sesiunile ce urmeaza sa le dau nu prea am posibilitate sa analizez aparatul matematik din link. Dupa vedem ce facem. Eu vreau sa va intreb dak dumneavoastra cunoastetsi ecuatia asta r''*r^2-C=0 eu ma gindesc ca ar trebui sa fie o constanta sau cel putin o relatie care ar determina o marime constanta dependenta de masa pm. Ei asta este un vis din clasa a 10 asa ca poate sa fie solutia undeva pe net sau prin carti. la noi in Republika MOldova putsini se okupa de chestii din astea de aia nu se gaseste ink nimeni sa raspunda orikum am sa am in vedere de informatia care mi-ati dat-o. PS am impresia ca activati in domeniul astronomiei, nu-i asa? ----------------------------------- Doru Dragan 31 Mai 2008 12:23 ----------------------------------- Uite pentru ca mi-esti simpatic si pentru ca cine stie de unde sare iepurele am mai cautat doua adrese interesante pentru tine. Prima este un referat la disciplina fizica asa cum e moda acum. Se ia referatul de-a gata de pe internet :) : http://www.referatele.com/referate/fizica/online9/Atractia-universala-Camp-gravitational-Proprietati-ale-atractiei-universale-Newton-referatele-com.php A doua este teza de doctorat a lui Ovidiu Vaduvescu: http://ovidiuvh.tripod.com/Teza-Romania/ Dupa ce o citesti cu atentie poti sa incerci sa iei legatura cu el. Iti dau adresa lui de pe internet ca el se misca mult si repede si nu stiu altfel cu sa-l gasesti: http://www.geocities.com/ovidiuv/ Da, activez in domeniul astronomiei ... de amatori :) ----------------------------------- kepler 31 Mai 2008 23:05 ----------------------------------- mdea, dupa cum vad voi avea de afacere cu o persoana destul de serioasa, va trebui ca prealabil discutiei noastre sa ma pregatesc destul de bine :) am sa ma starui sa analizez lukrarea care mi-ati propus-o (cel putsin ce ma intereseaza din ea). Si dupa aia sa ma apuk sa-l caut pe domnul Vaduvescu. sau cel putsin sa ma starui sa-l gasesc dupa cum ati mentionat, caci dak ma apuk aku sa vb cu el apoi nici nu voi avea cuvinte in fatsa lui:) mai in scurt trebuie sa ma pregatesc:) multsumesc pentru ajutor poate ceva se primeste :) ----------------------------------- Doru Dragan 01 Iun 2008 00:21 ----------------------------------- Cine se ocupa de mecanica cereasca nu poate sa fie decat o persoana serioasa :) Si eu, pe vremuri cand am inceput sa ma preocupe mecanica cereasca (problema celor doua corpuri, problema celor N corpuri, teoria perturbatiilor, teoria stabilitatii sistemelor etc.) eram o persoana foarte serioasa :D Dar intre timp mi-a trecut :lol: ----------------------------------- kepler 01 Iun 2008 03:19 ----------------------------------- Eu cind am inceput sa-mi pun imtrebarea asta nici nu banuiam ca face parte din mecanica cereasca, defapt evident lucru numai formula care urmeaza sa fie analizata multe spune. Am vrut sa vin la studii in Romania dar aveam un bal super de jos si nu facea concurentsa semenilor :( si acum sunt nevoit sa-mi fac studiile la facultatea de energetika nu rau dar nu ceia ce imi place si... nu merge bine :) ----------------------------------- mobutu 17 Iun 2008 10:59 ----------------------------------- :arrow: http://www.phys.spbu.ru/content/File/Library/studentlectures/schlippe/twobdy.pdf ----------------------------------- kepler 20 Iun 2008 13:09 ----------------------------------- Iti multumesc mobutu, am sa analizez referatul si am sa vad daca este acolo ceia ce am vazut eu "ca prin ceata" :) in two body problem, da intr-adevar aici e chestia de functii armonice care au proprietati dracoase si deaceia trebuie de posedat aparatul matematic necesar, dupa cum am previzualizat referatul 70% este pe inteles resul ramine de cunoscut inca o data multumesc