-
Data: 2009-07-26 14:48:42
Temat: Re: Całkowanie numeryczne - reaktywacja
Od: "slawek" <s...@h...pl> szukaj wiadomości tego autora
[ pokaż wszystkie 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 <=