-
Data: 2011-12-14 02:24:49
Temat: Re: pętla fizyki
Od: bartekltg <b...@g...com> szukaj wiadomości tego autora
[ pokaż wszystkie nagłówki ]W dniu 2011-12-11 17:33, M.M. pisze:
>
> Planety najpierw zblizaja sie do siebie, a potem oddalaja.
Tak mają;)
> Zarowno zblizanie jak i oddalanie odbywa sie z pewnym krokiem.
Nie mamy 'ciągłych' komputerów.
> Idealnie byloby, gdyby odleglosci od planet we wszystkich odpowiadajacych
> sobie krokach przyblizajacych i oddalajacych byly rowne, ale taka
> sytuacja jest nieprawdopodobna. Przy rownoodleglych krokach planety by
> tak samo przyspieszaly jak i spowalnialy i animacja bylaby stablina nawet
> po zastosowaniu Eulera.
le przestań myśleć eulerem. Na tym polega dobry schemat, że nie
ważne gdzie są kroki, wynik ma wyjść odpowiednio dokładny.
Taki schemat działa tak, że przybliża sobie naszą funkcję
czymś. Np RK robi coś w rodzaju rozwinięcia taylora
(oczywiście teoretyczine, w algorytmie spakowane jest
to do czystego wyliczania funkcji i artmetyki).
Taki euler korzysta z prostej i tyle. Przylizenie styczne.
A schematy wyższego rzędu - z wyższych przybliżeń;)
Takie kawałki 1/r^2 dość dobrze przybliżaką się
wielomianami 4 rzędu;)
W ciągu 1000 lat moj US zgubił/zyskał nieco energii
delta E/E = 8.58401867654852e-014
Dość dokładnie;)
Trzeci wykres to odchyłka energii w czasie:
http://students.mimuw.edu.pl/~bs219608/US/stale.png
>
> Wydaje sie ze wynik zaburzaja obliczenia w kroku ktory przypada na najblizsza
> odleglosc pomiedzy planetami. Jesli ten krok wypadnie na etap
> przyspieszania, to przyspieszy tak mocno planety, ze nastepne kroki tego
> przyspieszenia nie wyhamuja. Wydaje sie to sporym problemem, czy na
> pewno RK45 podola?
Działa znacznie lepiej niż euler;)
Problem, który opisujesz wiąże się ze stabilnością
schematów i cierpią na niego wszytskie schematy otwarte
(w tym RK).
A problem o którym mówisz byłby poważniejszy jakby planeta wpadała
w 'osobliwość', czyli b. blisko Słońca.
> Pozatym RR jest rzedu drugiego. Jesli dobrze rozumiem, to najpierw musze
> rozwiazac jedno w kilku punktach aby miec dane dla drugiego.
Nie.
1) te schematy sa dla układów równań!
wektor wartości Y [y1;y2;y3...yn]
Y' = F (Y,t) gdzie F jest funkcją wektorową.
2)
Układ X'' = F(X,X',t)
jest równoważny ukladowi po podstawieniu V=X'
V' = F(X,X',t)
X' = V
Mając więc ukłąd fizyczny o n stopniach swobody
wsadzamy do maszynki wektor Y = [V,X]
i Ff(Y) = [F(X,V,t), V]
> obliczen szybko sie naloza na siebie. Jestes pewien ze ta metoda starcza
> na 1000 lat symulacji ukladu slonecznego?
Dostałeś message ID, grupę i tytuł wątku.
Przeczytaj pierwszego posta, dowiesz się
wszytkiego;)
Spoiler: starcza;) Ale liczyło się długo.
No i nie było to RK45. Tym też przeliczyłem,
'na 100 lat', wyszło to samo, ale liczył się
znacznie dłużej.
pzdr
bartekltg
Następne wpisy z tego wątku
- 14.12.11 02:49 A.L.
- 14.12.11 06:12 bartekltg
- 14.12.11 16:47 A.L.
- 15.12.11 14:03 M.M.
- 15.12.11 22:43 bartekltg
- 15.12.11 23:31 bartekltg
- 16.12.11 00:04 A.L.
- 16.12.11 00:16 A.L.
- 16.12.11 06:07 M.M.
Najnowsze wątki z tej grupy
- Do czego nadaje się QDockWidget z bibl. Qt?
- Bibl. Qt jest sztucznie ograniczona - jest nieprzydatna do celów komercyjnych
- Co sciaga kretynow
- AEiC 2024 - Ada-Europe conference - Deadlines Approaching
- Jakie są dobre zasady programowania programów opartych na wtyczkach?
- sprawdzanie słów kluczowych dot. zła
- Re: W czym sie teraz pisze programy??
- Re: (PDF) Surgical Pathology of Non-neoplastic Gastrointestinal Diseases by Lizhi Zhang
- CfC 28th Ada-Europe Int. Conf. Reliable Software Technologies
- Młodzi programiści i tajna policja
- Ada 2022 Language Reference Manual to be Published by Springer
- Press Release - AEiC 2023, Ada-Europe Reliable Softw. Technol.
- Ada-Europe - AEiC 2023 early registration deadline approaching
- Ada-Europe Int.Conf. Reliable Software Technologies, AEiC 2023
- Ile cykli zajmuje mnożenie liczb 64-bitowych?
Najnowsze wątki
- 2024-06-05 [ot] spec od renowacji/reperacji kurtek skorzanych
- 2024-06-05 Koszt przywrócenia wychodnego numerowi w Plusie
- 2024-06-06 korki prawie takie same
- 2024-06-05 Takie elektryki mają sens ale czy z Francuską MARŻĄ?
- 2024-06-05 Warta S.A. - przyjęta odpowiedzialność?
- 2024-06-04 nie zna życia ten
- 2024-06-06 A jednak nie kondensatory
- 2024-06-06 Re: A jednak nie kondensatory
- 2024-06-06 Wymiana SIM Aero2
- 2024-06-06 Gdańsk => Programista Full Stack .Net <=
- 2024-06-06 Warszawa => Senior React Native Developer <=
- 2024-06-06 Gdańsk => Head of International Freight Forwarding Department <=
- 2024-06-06 Warszawa => Kierownik Działu Spedycji Międzynarodowej <=
- 2024-06-05 Olsztyn => Sales Specialist <=
- 2024-06-05 Ulm => Integration & Test Engineer <=