Re: pętla fizyki - Grupy dyskusyjne w eGospodarka.pl

eGospodarka.pl › Grupy › pl.comp.programming › pętla fizyki › Re: pętla fizyki

Path: news-archive.icm.edu.pl!news.icm.edu.pl!news.nask.pl!news.nask.org.pl!newsfeed2
.atman.pl!newsfeed.atman.pl!.POSTED!not-for-mail
From: bartekltg <b...@g...com>
Newsgroups: pl.comp.programming
Subject: Re: pętla fizyki
Date: Thu, 15 Dec 2011 23:43:28 +0100
Organization: ATMAN - ATM S.A.
Lines: 62
Message-ID: <jcdt6g$r7d$1@node2.news.atman.pl>
References: <jbnfs8$ihp$1@inews.gazeta.pl> <jbsrmu$svg$1@inews.gazeta.pl>
<jbto5f$tb3$1@node2.news.atman.pl> <jbtpqm$1cu$1@inews.gazeta.pl>
<jbu4kt$b5m$1@node2.news.atman.pl> <jc2lvv$ig4$1@inews.gazeta.pl>
<jc91dk$nfh$1@node2.news.atman.pl> <jccumu$468$1@inews.gazeta.pl>
NNTP-Posting-Host: 144-mi3-6.acn.waw.pl
Mime-Version: 1.0
Content-Type: text/plain; charset=UTF-8; format=flowed
Content-Transfer-Encoding: 8bit
X-Trace: node2.news.atman.pl 1323989008 27885 85.222.69.144 (15 Dec 2011 22:43:28
GMT)
X-Complaints-To: u...@a...pl
NNTP-Posting-Date: Thu, 15 Dec 2011 22:43:28 +0000 (UTC)
User-Agent: Mozilla/5.0 (Windows NT 6.1; WOW64; rv:8.0) Gecko/20111105
Thunderbird/8.0
In-Reply-To: <jccumu$468$1@inews.gazeta.pl>
Xref: news-archive.icm.edu.pl pl.comp.programming:194115
[ ukryj nagłówki ]
W dniu 2011-12-15 15:03, M.M. pisze:
> bartekltg<b...@g...com> napisał(a):
>
>>> Zarowno zblizanie jak i oddalanie odbywa sie z pewnym krokiem.
>> Nie mamy 'ciÄgĹych' komputerĂłw.
> Zadziwiajace jak szybko kazdy temat kieruje sie ku filozofii :)
>
> Przejrzalem ta prace:
> http://www.twins.pk.edu.pl/~awosatko/pocz_w11s.pdf
> W koncu cos prostego i w miare zrozumialego :)
>
> Na wykresach widac, ze metoda Lsode jest jednak duzo dokladniejsza
> niz RK.

Nigdy nie mówiłem, że RK jest najlepsza;)
No i co to jest lsode. Bo równie dobrze za pomocą lsode
można odpalić eulera, jak i metode adamsa 12 rzędu;)
Jeszcze się okaże, że domyślnie (bez opcji) używa
RK45;) Porównujesz wynik z metoda RK trzeciego rzędu.

> Przy ugieciu belki mozna zalozyc ze oddzialujace sily sa niezmienne.
> Belka po ugieciu jest odrobine blizej osrodka ktory ja przyciaga i
> sila grawitacyjna jest minimalnie wieksza, ale mozna ja pominac.

> Nie widze jeszcze jak uwzglednic owe przesuniecie.

Olać.
Inaczej, policz najpierw zmienne tej siły przez zbliżenie
do jądra ziemi, a potem olej, jako nieistotną.

>> 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.
> Hmmm, nie rozumiem. Zeby zastosowac RK np. drugiego rzedu to musze

A przeczytałeś to, co było po 'nie'.

> policzyc k1 i k2. Jesli RR jest drugiego rzedu to o ile sie nie myle
> musze 3 razy liczyc k1 i 3 razy k2? 2 razy dla drugiej pochodnej i jeden
> raz dla pierwszej? Czyli gdy uzyje RK 4go rzedu to 20 "obliczen" w jednym
> kroku?

Nie przeczytałeś:(
Prościej nie umiem napisać, musisz sięgnąć do materiałów:)

Nie mamy żadnych prędkości i przyspieszeń. Do 'solvera' wkładamy
już układ pierwszego rzędu. I dalej masz wszytko wektorowo,
a k1,k2 są skalarmi.

K1, k2... liczysz raz (w kroku). W RK rzędu n wywołujesz
funkcję n razy. W adaptatywnych ciut gorzej, RK45 wywołuje
6 razy.

pzdr
bartekltg