-
Path: news-archive.icm.edu.pl!newsfeed.gazeta.pl!newsfeed.neostrada.pl!nemesis.news.n
eostrada.pl!atlantis.news.neostrada.pl!news.neostrada.pl!not-for-mail
From: "slawek" <s...@h...pl>
Newsgroups: pl.comp.programming
Subject: Re: Całkowanie numeryczne - reaktywacja
Date: Sun, 26 Jul 2009 16:48:42 +0200
Organization: TP - http://www.tp.pl/
Lines: 32
Message-ID: <h4hqv0$ft7$1@nemesis.news.neostrada.pl>
References: <h4eu8v$e00$1@nemesis.news.neostrada.pl> <h4fl6n$fm4$1@inews.gazeta.pl>
<h4gu5f$h3e$1@atlantis.news.neostrada.pl> <h4h98b$q6j$1@inews.gazeta.pl>
NNTP-Posting-Host: 62.69.219.25
Mime-Version: 1.0
Content-Type: text/plain; format=flowed; charset="iso-8859-2"; reply-type=original
Content-Transfer-Encoding: 8bit
X-Trace: nemesis.news.neostrada.pl 1248620320 16295 62.69.219.25 (26 Jul 2009
14:58:40 GMT)
X-Complaints-To: u...@n...neostrada.pl
NNTP-Posting-Date: Sun, 26 Jul 2009 14:58:40 +0000 (UTC)
In-Reply-To: <h4h98b$q6j$1@inews.gazeta.pl>
X-Priority: 3
X-MSMail-Priority: Normal
Importance: Normal
X-Newsreader: Microsoft Windows Live Mail 14.0.8064.206
X-MimeOLE: Produced By Microsoft MimeOLE V14.0.8064.206
Xref: news-archive.icm.edu.pl pl.comp.programming:182795
[ ukryj nagłówki ]
Użytkownik "Mariusz Marszałkowski" <b...@N...gazeta.pl> napisał w
wiadomości grup dyskusyjnych:h4h98b$q6j$...@i...gazeta.pl...
> Szczerze to też doznałem szoku. Ale gdy zrozumiałem że tyle samo traci
> "pod górkę" co zyskuje "z górki" to szok minął. :)
Nie, to nie tak. Po prostu algorytm z trapezami (I) i algorytm ze splinami
(II) dają dla danego przedziału (a,b) liczby różne od wyniku analitycznie
otrzymanego U o odpowiednio -epsilon i +epsilon, gdzie epsilon jest funkcją
kroku h. Oczywiście, że dla (0,pi/2) jest trochę np. nadmiar a dla (pi/2,pi)
jest niedomiar w danej metodzie, więc się znosi - to też jest, ale nie o to
chodzi. Chodzi o to, że dwa zupełnie różne algorytmy dają wartość odchylenia
od wyniku takie same co do wartości bezwzględej (z różnicą do około 1%) dla
różnych kroków itd. - tylko znaki są różne. Czyli średnia z obu metod jest
dużo dokładniejsza i błąd jej malej znacznie szybciej niż O(h^2). Patrz
http://drop.io/aqtdw2y - dla małych h, czyli dużych n wynik jest 10 tysięcy
razy dokładniejszy! Błąd jest rzędu O(h^3) - wykres to ładnie prezentuje.
> Nie wiem, dlaczego wtedy jest źle?
Niekoniecznie źle być musi. Po prostu "cały wielki program" liczy y[j+1]
jako funkcję poprzednich y[1],y[2],...,y[j] - całki i takie tam. Możliwe
jest, że y[j] leci do nieskończoności. Biorąc dużo punktów opóźnia się
moment wykrycia takiej rozbieżności - procedura wygładza zachowanie funkcji,
która w pobliżu osobliwości mało przypomina grzeczny wielomian. To trochę
bardziej skomplikowane - bo to i układ równań jest, i jeszcze y jest
zespolone i parę jeszcze zagwozdek po drodze.
slawek
Następne wpisy z tego wątku
- 27.07.09 18:10 Wit Jakuczun
- 27.07.09 18:54 Mariusz Marszałkowski
- 27.07.09 19:01 slawek
- 27.07.09 19:09 Wit Jakuczun
- 27.07.09 19:11 Wit Jakuczun
- 27.07.09 19:15 Wit Jakuczun
- 27.07.09 19:15 slawek
- 27.07.09 19:21 slawek
- 27.07.09 19:27 Mariusz Marszałkowski
- 27.07.09 19:32 slawek
- 27.07.09 19:36 bartekltg
- 27.07.09 19:38 Wit Jakuczun
- 28.07.09 07:44 slawek
- 28.07.09 08:18 slawek
- 28.07.09 15:31 bartekltg
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-10 wyobrazcie sobie ze
- 2024-06-10 malowanie samochodu
- 2024-06-10 News from Poland
- 2024-06-10 Czy na pewno będą CŁA na chińskie samochody?
- 2024-06-09 Dlaczego w Polsce sie nic nie udaje, na przykładzie niebieskiego lasera a teraz perskowitów
- 2024-06-09 Dlaczego w Polsce sie nic nie udaje, na przykładzie niebieskiego lasera a teraz perskowitów
- 2024-06-09 Wykrywanie przerwy w długim przewodzie zakopanym w ziemi.
- 2024-06-09 Czemu news.chmurka.nwt jest taki wolny?
- 2024-06-11 Funbox 3.0 zakres adresów DHCP
- 2024-06-11 Re: Funbox 3.0 zakres adresów DHCP
- 2024-06-09 Miernik szybkości netu
- 2024-06-11 Panele PV w pionie (prawie).
- 2024-06-11 czy ta grupa żyje?
- 2024-06-11 Warszawa => Senior React Native Developer <=
- 2024-06-11 Gdańsk => Kierownik Działu Spedycji Międzynarodowej <=