-
Data: 2009-07-07 11:53:13
Temat: Re: Pocedura całkowania
Od: bartekltg <b...@g...com> szukaj wiadomości tego autora
[ pokaż wszystkie nagłówki ]On 7 Lip, 07:34, "slawek" <s...@h...pl> wrote:
> y[1],...,y[n]. Oczywiście metodę Romberga można zastosować do funkcji
> interpolującej - ale dla takiej funkcji łatwiej i szybciej policzyć całkę
> analitycznie (np. gdy interpolacja wielomianem).
Jeszcze inaczej:
No wlasnie! na tym polegaja nasze kwadratury interpolacyjne.
Bierzemy wielomian/funkcje sklajana, dopasowujemy do wartosci
w wybranych wezlach i liczymy scisle calke.
Kwadratura parabol bierze trzy punkty, (0,1,2) dopasowyje tam
parabole,
po czym z jej wspolczynnikow okresla 'pole pod parabolą'. Nastepnie
bierze kolejne dwa przedziały/trzy punkty (2,3,4)..
Jadnak pracujac chwile na kartce papieru mozna zrobic to samo
nie wyliczajac wprost wspolczynnikow paraboli, ale od razu napisac
wzorek wiazacy trzy punkty z _calką pod parabolą interpolujaca te
punkty_.
Podobnie simpson interpoluje na kazdych polejnych trzech przedialach
(4 punkty) wielomian 3 stopnia i wylica calke tego wialomianu.
Jesli odleglosci meidzy wezlami nie sa stale, romberg rzeczywiscei
odpada,
ale kadratury interpolacyjne (wialomianowe, o splajnach juz wiesz jak
to zrobic)
nadal dzialaja. Bierzesz jakas ladna baze, najlepiej lagrange'a
('kanoniczna')
i nia interpolujesz.
Przypomne, ze jesli mamy funkcja f, L_i to baza lagrangea oparta na
wezlach x_j,
to wielomianem interpolacyjnym funkcji f jest po prostu
ff= \sum f(x_i) *L_i
Calke takiego wielomianu mozemy zapisac \int ff= \int \sum f(x_i)
*L_i =
= \sum f(x_i) \int L_i
L_i(x) = (x-x_1)..(x-x_{i-1})*(x-x_{i+1})(x-x_n) /[ (x_i-x_1)..(x_i-x_
{i-1})*(x_i-x_{i+1})(x_i-x_n) ]
calke z tego wyrazenia w granicach x1-xn jestes w stanie policzyc
analitycznie i wstawic
w praogram (jaklo funkcje x_1..x_n). Dla przyzwoitych wezlow powinne
byc wszystkie
dodatnie (ma to pewne znaczenie dla stabilnosci numerycznej). n rzedu
2-10.
Bedzie to na pewno szybsze niz splajny dla duzej liczby punktow.
Mowiles, ze punktow
przybywa. Tutaj poprawiasz jedynie ostatni fragment, dla splajnow
musisz
poprawci wszystkie m(10 000) punktow.
Tylko uwaga, aproksymujac wielomianem stopnia n-1 potrzebujesz n
punktow, czyli
n-1 nowych punktow. Teraz juz Twoj wybor, czy uaktialniac calke tylko,
gdy ich tyle
przybedzie, czy ostatni fragment przyblizac inaczej (wielomianem
mniejszego stopnia
[lepsze wlasnosci] czy wrecz kilkoma - nawet trapezami[najgorsze, dla
niewielkeigo m
i duzego n mozesz wrecz widziac na wykresie poprawke]).
pozdr
bartekltg
Następne wpisy z tego wątku
- 07.07.09 12:20 slawek
- 07.07.09 12:25 slawek
- 07.07.09 12:40 slawek
- 07.07.09 12:53 Mateusz Ludwin
- 07.07.09 13:02 A.L.
- 07.07.09 13:23 slawek
- 07.07.09 13:54 bartekltg
- 07.07.09 13:58 slawek
- 07.07.09 14:12 slawek
- 07.07.09 14:22 slawek
- 07.07.09 14:24 bartekltg
- 07.07.09 16:15 slawek
- 07.07.09 16:31 slawek
- 07.07.09 16:36 Wit Jakuczun
- 07.07.09 16:53 slawek
Najnowsze wątki z tej grupy
- 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
- Młodzi programiści i tajna policja
- Ada 2022 Language Reference Manual to be Published by Springer
Najnowsze wątki
- 2024-10-04 Warszawa => Senior PHP Laravel Developer (e-commerce) <=
- 2024-10-04 Warszawa => Data Scientist / Data Engineer (predictive modelling) <=
- 2024-10-03 Nieparzyste dmuchanie
- 2024-10-03 Prognozowanie zużycia energii przez PGE?
- 2024-10-03 Re: Drugi ekran na Androidzie
- 2024-10-03 sprawiedliwosc nierychliwa
- 2024-10-03 zloto
- 2024-10-03 Odkurzacz mnie bije :(
- 2024-10-03 Gdańsk => Technical Lead ( (Java Background)) <=
- 2024-10-03 Warszawa => Mid IT Recruiter <=
- 2024-10-03 Olsztyn => Sales Specialist <=
- 2024-10-03 Leszczyna nie zna prawa?
- 2024-10-03 Warszawa => OpenText ECM Specialist <=
- 2024-10-03 Blokowanie informacji - test
- 2024-10-02 Warszawa => Fullstack Developer <=