eGospodarka.pl
eGospodarka.pl poleca

eGospodarka.plGrupypl.comp.programmingpętla fizyki › Re: p?tla fizyki
  • Path: news-archive.icm.edu.pl!news.icm.edu.pl!plix.pl!newsfeed1.plix.pl!goblin1!gobli
    n3!goblin.stu.neva.ru!news.ripco.com!news.glorb.com!news-in-01.newsfeed.easynew
    s.com!easynews!core-easynews-01!easynews.com!en-nntp-11.dc1.easynews.com.POSTED
    !not-for-mail
    From: A.L. <l...@a...com>
    Newsgroups: pl.comp.programming
    Subject: Re: p?tla fizyki
    Message-ID: <o...@4...com>
    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>
    <9...@4...com>
    <jc9ep7$3qe$1@node2.news.atman.pl>
    X-Newsreader: Forte Agent 4.2/32.1118
    MIME-Version: 1.0
    Content-Type: text/plain; charset=ISO-8859-2
    Content-Transfer-Encoding: 8bit
    Lines: 54
    X-Complaints-To: a...@e...com
    Organization: Forte Inc. http://www.forteinc.com/apn/
    X-Complaints-Info: Please be sure to forward a copy of ALL headers otherwise we will
    be unable to process your complaint properly.
    Date: Wed, 14 Dec 2011 10:47:34 -0600
    Xref: news-archive.icm.edu.pl pl.comp.programming:194077
    [ ukryj nagłówki ]

    On Wed, 14 Dec 2011 07:12:52 +0100, bartekltg <b...@g...com>
    wrote:

    >W dniu 2011-12-14 03:49, A.L. pisze:
    >> On Wed, 14 Dec 2011 03:24:49 +0100, bartekltg<b...@g...com>
    >> wrote:

    >>
    >> Co to znaczy "lepiej"?...
    >
    >Osiągamy dany wynik przy mniejszej długości kroku,
    >czyli mniejszym kosztem numerycznem.

    Niezupelnie jest to oczywiste, bo RK wymaga wielu ewaluacji prawych
    stron rownania do wykonania jednego kroku, a Euler tylko jednego. Wiec
    nalezaloby sie zastanowic jak wypada proownanie.

    >Wszystkie wspominane tu schematy są zbieżne,
    >więc gdy h-> 0 nasze rozwiązanie numeryczne
    >zbiegnie do prawdziwego (wszytko to pod założeniami,
    >żę rozwiazanie istnieje, jest jednoznaczne,
    >i raczej w artmetyce dokładnej).
    >

    Zbieznosc jest wlasnoscia "trywialna"; ewidentnie nei stosujemy
    schematow ktore nie sa zbiezne. Natomiast kazdy schemat charakteryzuje
    sie dwoma parametrami: rzedem aproksymacji i obszarem stabilnosci. I
    tylko w tych kategoriach schematy mozna porownywac, dodajac
    ewentualnie zlozonosc obliczeniowa rozumiana jako ilosc ewaluacji
    prawych stron rownania na krok.

    Przy okazji dobrze wspomniec ze "Runge Kutta" nei oznacza jednej
    metody, a klase metod z ktorych kazda ma inne wlasciwosci w sensie
    podanych wyzej, to znaczy obszar stabilnosci, rzad aproksymacji i
    zlozonosc obliczeniowa. RK45 jest nazwa jednej z konkretnych
    implementacji.

    >O dwóch porządnych schamatach tak bym nie mówił,
    >bo czy RK, czy wielokrotowy, czy jeszcze inny
    >'jest lepszy' będzie zależało od zagadnienia.
    >
    >Ale porównując adaptatywne RK45 i otwartego eulera,
    >ciężko mi znaleść zalet tego drugiego. Więc
    >pierwszy lepszy;)

    Obszary stabilnosci popularnych metod RK i Eulera sa bardzo zblizone,
    wiec z punktu widzenia maksymalnego dopuszczalnego kroku niewiele sie
    obie metody roznia. Jezeli zas jestesmy w obszarze stabilnosci, RM
    wysokiego rzedu aproksymacji da rozwiazanie dokladniejsze niz Euler,
    ale bedzie to okupione wieksza zlozonoscia obliczeniowa. Co sie oplaca
    a co nie - pozostaje problemem otwartym., a odpowiedz zalezy miedzy
    innymi od tego jak kosztowne jest obliczanie prawych stron roznania.
    >
    A.L.

Podziel się

Poleć ten post znajomemu poleć

Wydrukuj ten post drukuj


Następne wpisy z tego wątku

Najnowsze wątki z tej grupy


Najnowsze wątki

Szukaj w grupach

Eksperci egospodarka.pl

1 1 1

Wpisz nazwę miasta, dla którego chcesz znaleźć jednostkę ZUS.

Wzory dokumentów

Bezpłatne wzory dokumentów i formularzy.
Wyszukaj i pobierz za darmo: