21.01.2015 Views

Gaussovy kvadratury Gaussian Quadrature

Gaussovy kvadratury Gaussian Quadrature

Gaussovy kvadratury Gaussian Quadrature

SHOW MORE
SHOW LESS

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.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!