-
Data: 2011-12-10 11:06:27
Temat: Re: pętla fizyki
Od: bartekltg <b...@g...com> szukaj wiadomości tego autora
[ pokaż wszystkie nagłówki ]W dniu 2011-12-10 11:29, Roman W pisze:
> On Dec 9, 11:12 pm, bartekltg<b...@g...com> wrote:
>>> A jak to mozna ulepszyc?
>>
>> Ojej, tego jest na pęczki.http://en.wikipedia.org/wiki/Numerical_ordina
ry_differential_equations
>> I oglądaj te 'wysokiego rzędu'.
>>
>> Na dzień dobry poszukaj Runge-Kutta czwartego rzędu albo
>> adaptatywnej 'RK45'.http://en.wikipedia.org/wiki/Runge%E2%80%93Ku
tta%E2%80%93Fehlberg_method
>> Do zabaw najcześćiej całkowicie wystarczają.
>>
>> Po coś lepszego możesz przeszukać ten wątek:
>> Newsgroups: pl.sci.kosmos,pl.sci.fizyka
>> Subject: Domowe symulacje Układu Słonecznego [numerki]
>> Date: Tue, 15 Nov 2011 17:42:11 +0100
>> Message-ID:<j9u4p5$7m8$1@node2.news.atman.pl>
>>
>> Głownie jechaliśmy tam na gotowcach (np nieźle sprawdzała
>> się wielokrotowa metoda z matalba), ale m.in A.L podawał
>> trochę linków do bardziej zaawansowanych metod.
>>
>> Ale do piania skryptów na www zmęczyłbym po prostu tą RK45,
>> i to też raczej poszukał gotowca;)
>
> O ile dobrze pamietam, to do symulacji ukladow oddzialujacych
> grawitacyjnie RK nie wystarcza, bo nie gwarantuje zachowania energii.
Ale takich cudownych metod niestety nie ma:)
Są schematy geometryczne / symplektyczne (i do nich
właśnie linka podawał A.L.) które są temu dość bliskie
(dokładniej, zachowują formę symplektyczną i 'lekko
zmodyfikowany hamiltonian' - zachowyją enerigę, ale
nieco innego układu).
Z drugiej strony nie ma co popadać w przesadyzm. Porządna
metoda ogólna z odpowiednim krokiem/błędem lokalnym
też da radę (We wspominanym wątku zwykłe metody
dawały zadowalające wyniki przy 1000 latach symulowania
Układu Słonecznego).
> Wez jakas ksiazke o symulowaniu dynamiki molekularnej, byly nawet
> polskie tlumaczenia wydane X lat temu przez PWN (pewnie pojawily sie
> nowe, bo to popularny temat).
Dynamika molekularna nie kojarzy mi sie z precyzyjnymi obliczeniami.
Ma nieco co innego za zadanie policzyć. Symulujemy nawet Eulerem
i dodajemy 'termostat'. Nikt nie oczekuje ani zachowania energii.
pzdr
bartekltg
Następne wpisy z tego wątku
- 11.12.11 16:33 M.M.
- 14.12.11 02:24 bartekltg
- 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
- 7. Raport Totaliztyczny: Sprawa Qt Group wer. 424
- TCL - problem z escape ostatniego \ w nawiasach {}
- Nauka i Praca Programisty C++ w III Rzeczy (pospolitej)
- testy-wyd-sort - Podsumowanie
- Tworzenie Programów Nieuprzywilejowanych Opartych Na Wtyczkach
- 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
Najnowsze wątki
- 2024-12-02 Akumulatorki Ni-MH AA i AAA Green Cell
- 2024-12-02 Usiłowanie zabójstwa
- 2024-12-01 Rambo 2024. Co z radio-stopem
- 2024-12-01 Pijani kierowcy
- 2024-12-01 "Chciałem zamówić kurs tym"
- 2024-11-30 Windykatorzy ścigają spadkobierców z mandat nieboszczyka za przekroczenie prędkości???
- 2024-11-30 Łódź => Technical Artist <=
- 2024-11-30 Lublin => Inżynier Serwisu Sprzętu Medycznego <=
- 2024-11-30 Warszawa => Microsoft Dynamics 365 Business Central Developer <=
- 2024-11-30 Bieruń => Team Lead / Tribe Lead FrontEnd <=
- 2024-11-30 Zielona Góra => Senior PHP Symfony Developer <=
- 2024-11-30 Gdańsk => Specjalista ds. Sprzedaży <=
- 2024-11-30 Lublin => Spedytor międzynarodowy <=
- 2024-11-30 Warszawa => Mid IT Recruiter <=
- 2024-11-30 Warszawa => Fullstack Developer <=