-
Data: 2012-11-08 14:19:56
Temat: Re: Błędny epsilon - this is not a bug, this is ?
Od: bartekltg <b...@g...com> szukaj wiadomości tego autora
[ pokaż wszystkie nagłówki ]W dniu 2012-11-08 01:51, AK pisze:
> Użytkownik "bartekltg" <b...@g...com> napisał:
>
>> Nie no, nie sprowadzał _wszystkich_ metod całkowania numerycznego
>> do trapezów, a jedynie kwadratury Newtona-Cotesa (interpolacja
>> wielomianem na równo odległych węzłach, sprowadza się do sumy
>> sum_i w_i*f(x_i)).
>
> Hm.. Jakby tu rzec.. Co mi "swita we lbie", ze przy pewnym zalozeniu
> to slawek moze miec sporo racji..
Tak, dla całkiem sporej rodziny funkcji wyniki będą te same.
Np liniowych;) czy sinusa po pełnym okresie. Tylko co z tego.
Dla każdej funkcji ciaglej znajdziesz kwadrature jednopunktową,
dająca dokładny wynik. Liczy się ogolna sprawność metody
dla funkcji danej klasy.
> O ile pamietam (wybacz, 25lat to naprawde dosc duzo aby zapomniec)
> to trapezy maja sie tak:
> 1/2*(f[0] + 2*SUM(i=1,n-1: f[i] + f[n])
> a "klasyczny" Simpson:
> 1/3*(f[0] + 4*SUM(i=1,n-1, 2: f[n]) + 2*SUM(i=2,n-2, 2: f[n]) + f[n])
> ale.. gdyby tak w ty Simpsonie brac tylko/sumowac pierwszy przedzial
> i przesuwac parabolke o jeden krok a nie dwa, albo tym szesciennym (3/8)
> Simpsonie brac srodkowy przedzial, tez przesuwajac co jeden krok
> to kto wie czy nie sprowadzi sie to (poza punktami skrajnymi przedzialu,
> ale t mozna pomonac/zaniedbac) do tych trapezow ?
Podobnie argumentował slawek. I jest to właśnie nasza bzdura;)
Wyobraź sobie to tak: wydziel z simpsona sumę dającą trapezy,
pozostanie poprawka (o dodatnich i ujemnych wagach!).
Ta poprawka odpowiada za scałkowanie 'krzywizny' - tego,
co wystaje ponad linię prostą z dokładnością do paraboli.
Po uśrednieniu ten wkład się kasuje. Tak, problem występuje
na brzegach. Jeśli jest to funkcja okresowa, można simpsona
przesunąć 'o pół', dostać równie dobry wynik. Ich średnia
też bedzie tak samo dobra => trapezy są tak samo dobre.
Ale co w tym dziwnego, skoro dla takich funkcji całki po okresie
z kolejnych pochodnych wynoszą 0:)
Zresztą, ostatnio chciałeś zabaw programistycznych. Jest
okazja. Porównanie kilku takich metod to chwila, a zobaczysz
na własne oczy różnicę, nie będziesz musiał wierzyć na słowo,
że taka drobna zmiana daje dużo.
Złośliwie można zaproponować jako funkcję próbną np x^3:)
Bardziej uczciwie, np sin(x)^2/x na [0,pi/2]. To ma znaną
analitycznie wartość.
Dywagacje dywagacjami, a metody wyższego rzędu działają lepiej;)
[oczywiście dla odpowiednio gładkich funkcji, simpson da ciała
na nieciągłości pochodnej, i z uwagą, że N-S przy rzędzie chyba
coś koło 8 przestaje być stabilny - pojawiają się się ujemne wagi].
BTW, kwadratura romberga oparta ejst na tych samych węzłach,
a zbiega jak szalona (dla funkcji gładkich). A to przecież
tylko średnia ważona kilku trapezów.
pzdr
bartekltg
Następne wpisy z tego wątku
- 08.11.12 14:51 Baranosiu
- 08.11.12 20:36 AK
- 08.11.12 20:49 AK
- 08.11.12 21:07 Michoo
- 08.11.12 21:36 Baranosiu
- 09.11.12 00:04 e...@g...com
- 09.11.12 01:12 Andrzej Jarzabek
- 09.11.12 01:21 Andrzej Jarzabek
- 08.11.12 09:19 slawek
- 08.11.12 09:23 slawek
- 08.11.12 09:34 slawek
- 08.11.12 09:31 slawek
- 08.11.12 09:38 slawek
- 10.11.12 19:25 Roman W
- 14.11.12 23:26 Michoo
Najnowsze wątki z tej grupy
- Arch. Prog. Nieuprzywilejowanych w pełnej wer. na nowej s. WWW energokod.pl
- 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
Najnowsze wątki
- 2025-01-04 gasik
- 2025-01-04 13. Raport Totaliztyczny: Powszechna Deklaracja Praw Człowieka Nie Chroni Przed Wyzyskiem Ani Przed Eksploatacją
- 2025-01-04 Zbieranie danych przez www
- 2025-01-04 reverse engineering i dodawanie elementów do istniejących zamkniętych produktów- legalne?
- 2025-01-04 w Nowym Roku 2025r
- 2025-01-04 Warszawa => Specjalista ds. IT - II Linia Wsparcia <=
- 2025-01-04 Warszawa => Java Developer <=
- 2025-01-04 Warszawa => Spedytor Międzynarodowy <=
- 2025-01-04 Warszawa => System Architect (Java background) <=
- 2025-01-04 Wrocław => Application Security Engineer <=
- 2025-01-04 Chrzanów => Specjalista ds. public relations <=
- 2025-01-04 Katowice => Key Account Manager (ERP) <=
- 2025-01-03 Problem z odczytem karty CF
- 2025-01-03 Jazda z Warszawy do Krakowa teslą
- 2025-01-03 Wrocław => Konsultant Wdrożeniowy Comarch XL/Optima (Księgowość i