Gaussovy kvadratury Gaussian Quadrature
Gaussovy kvadratury Gaussian Quadrature
Gaussovy kvadratury Gaussian Quadrature
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
4 INTERPOLAČNÍ POLYNOMY 20<br />
Je nasnadě ověřit, jak bude tato funkce vypadat v případě, že místo proměnné t dosadíme<br />
hodnotyx 1 ,x 2 ,...,x n . Dosazením dostaneme<br />
F(x j ) = (x−x 1 )(x−x 2 )...(x−x n )R(x j ), j ∈ {1,2,...,n}.<br />
Jak ale bylo řečeno, R(x j ) = 0, takže F(x j ) = 0 pro každé j ∈ {1,2,...,n}.<br />
Stejný postup provedeme také pro pevně zvolený bodx, tedy<br />
F(x) = (x−x 1 )(x−x 2 )...(x−x n )R(x)−(x−x 1 )(x−x 2 )...(x−x n )R(x) = 0.<br />
Z toho vyplývá, že funkceF je nulová pro n+1 hodnot intervaluI.<br />
V původním intervalu I budeme uvažovat n subintervalů s krajními body x 1 ,x 2 , ...,x n ,x,<br />
ve kterých má funkce F hodnotu 0. Z předpokladů výše uvedené věty vyplývá, že jsou splněny<br />
i předpoklady Rolleovy věty. Víme, že v každém z těchto subintervalů existuje bod, v němž má<br />
funkce F nulovou první derivaci. První derivace funkce F tedy má na intervalu I alespoň n nulových<br />
bodů. Opět aplikujeme Rolleovu větu a získámen−1 bodů, v nichž má funkceF nulovou<br />
druhou derivaci.<br />
Stejný postup můžeme opakovat. Zjistili jsme, že funkce F má n-tou vlastní derivaci, jejíž<br />
hodnota je nulová alespoň pro jeden bodξ ∈ I. Vztah pro n-tou derivaci bude<br />
Pro každé t ∈ I však platí, že<br />
F (n) (t) = (x−x 1 )(x−x 2 )...(x−x n )R (n) (t)−n!R(x).<br />
R (n) (t) = f (n) (t), nebot’P (n) (t) = 0.<br />
Vztah můžeme přepsat do tvaru<br />
F (n) (t) = (x−x 1 )(x−x 2 )...(x−x n )f (n) (t)−n!R(x).<br />
Za t dosadímeξ. Víme, že v tomto bodě je n-tá derivace funkceF nulová.<br />
Po úpravě pak chybu Lagrangeova polynomu můžeme vyjádřit jako<br />
kde ξ ∈ I.<br />
R(x) = (x−x 1)(x−x 2 )...(x−x n )<br />
f (n) (ξ), (9)<br />
n!<br />
Zbývá vyšetřit situaci, kdyxnení různé odx 1 ,x 2 ,...,x n , tznx = x j pro některé<br />
j ∈ {1,2,...,n}. Po dosazení do pravé strany vztahu (9) bude rozdíl(x−x j ) nulový. Vzhledem<br />
k (8) platí rovněž R(x j ) = 0. Vztah (9) tedy platí pro libovolné ξ ∈ I.