marcinka all, 20030107


0x01 graphic

0x01 graphic
(wielomian interpolacyjny funkcji sklejanej)0x01 graphic
kwadratury

TWIERDZENIE

Kwadratury N-C oparte na n+ 1 węzłach są rzędu:

0x01 graphic
dla n parzystych

0x01 graphic
dla n nieparzystych

(Sprawdzenie z przykładu 3 - poprzedni wykład)

0x01 graphic
nieparzyste

0x01 graphic
- dla wielomianu stopnia 2 z resztą 0x01 graphic
gdyż zależą od 0x01 graphic

Przykład 4

0x01 graphic
węzły: 0x01 graphic
, 0x01 graphic
, 0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic
(7)

Wzór paraboliczny, wzór Simpsona

Kwadratura N-C dla 0x01 graphic

Reszta: 0x01 graphic
; 0x01 graphic
0x01 graphic

Zbieżność ciągu kwadratur

0x01 graphic
ciąg kwadratur

0x08 graphic
0x01 graphic
(8)

numer ciągu

przybliżających całkę 0x01 graphic

Załóżmy, że dane są skończone macierze trójkątne węzłów 0x01 graphic
i współczynników 0x01 graphic

0x01 graphic
0x01 graphic

... ...

definiujące ciąg kwadratur 0x01 graphic

TWIERDZENIE

Ciąg kwadratur (8) jest rozbieżny dla dowolnych funkcji ciągłych na 0x01 graphic
czyli:

0x01 graphic

wtedy i tylko wtedy gdy:

  1. ciąg (8) jest zbieżny dla dowolnego wielomianu

oraz

  1. istnieje taka stała K, że dla n, (0x01 graphic
    ) zachodzi nierówność:

0x01 graphic

Istnieją funkcje, dla których ciąg kwadratur N-C nie jest zbieżny.

Akumulacja błędów.

Nie stosuje się kwadratur wyższych rzędów.

Złożone kwadratury Newtona-Cotesa

Przedział dzielimy na N części:

0x01 graphic
0x01 graphic
0x01 graphic

w każdym z nich kwadratura N-C niskiego rzędu.

Przykład 5 - złożona kwadratura

W każdym podprzedziale 0x01 graphic
stosujemy wzór trapezów (kwadratura N-C 2 rzędu z przykładu 1)

0x01 graphic

0x08 graphic

0x01 graphic
0x01 graphic

--> 0x01 graphic
[Author:JDz] - złożony wzór trapezów

0x01 graphic

jeśli 0x01 graphic
dla 0x01 graphic
to

0x01 graphic
0x01 graphic

dla 0x01 graphic

0x01 graphic
0x01 graphic

Przykład 6 (złożona kwadratura Simpsona)

Założenie: 0x01 graphic

0x08 graphic

0x01 graphic

0x08 graphic

0x01 graphic

0x08 graphic
nieparzyste

błąd

0x01 graphic

0x01 graphic

Ogólnie złożone kwadratury N-C są postaci:

0x01 graphic

0x01 graphic
; 0x01 graphic

TWIERDZENIE

Dla dowolnego n ciąg złożonych kwadratur N-C 0x01 graphic
jest przy 0x01 graphic
(czyli 0x01 graphic
) zbieżny dla wszystkich funkcji 0x01 graphic
.

0x01 graphic

Praktycznie

Szukam przybliżenia 0x01 graphic
mając żądaną dokładność:

0x01 graphic

oszacowanie reszty często trudne lub niemożliwe. Jak wyznaczyć N?

Obliczyć kolejne kwadratury 0x01 graphic
, 0x01 graphic
, ..., 0x01 graphic

Ilość przedziałów

Zwykle 0x01 graphic
, potem każdy przedział dzielimy na m części, m dobieramy tak, aby w następnej kwadraturze korzystać z poprzednio obliczonych wartości f

0x01 graphic

np. dla złożonej kwadratury trapezów należy przyjąć m = 2, gdyż (z przykładu 5) dla 0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

gdzie 0x01 graphic

Obliczenia kontynuujemy aż:

0x01 graphic
0x01 graphic

Skąd się wzięło to przejście?

0x01 graphic



Wyszukiwarka

Podobne podstrony:
marcinka all, 20030121
marcinka all, 20021203, Ciąg dalszy:
marcinka all, 20021119
marcinka all, 20021008
marcinka all, 20021015, SZUKANIE ZER W FUNKCJACH NIELINIOWYCH
marcinka all, 20021126, (RYSUNEK)
marcinka all, 20021112, INTERPOLACJA FUNKCJAMI SKLEJANYMI:
IO ALL
ZLL ALL
All Flesh Must Be Eaten Two Rotted Thumbs Up
Jim Hall at All About Jazz
all
PDH, Broadband ISDN, ATM and all that
marcinstolp pro
mo all
Twarde dyski, Informatyka -all, INFORMATYKA-all
farmacja 12czerwca2007, Receptura, Farma - pytania, testy egzaminacyjne-all
Opis programu komputerowego Twierdzenie Pitagorasa-dowód i z, wrzut na chomika listopad, Informatyka

więcej podobnych podstron