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.

3 ORTOGONÁLNÍ SYSTÉMY POLYNOMŮ 14<br />

3.4 Hermitovy polynomy<br />

Ačkoli tomu název napovídá, jedná se o systém ortogonálních polynomů, nikoli o Hermitovy<br />

interpolační polynomy. Těmi se budeme zabývat v jedné z následujících kapitol.<br />

Necht’ I = (−∞,∞) a w(x) = e −x2 . Pak prostor F bude prostorem všech funkcí splňujících<br />

podmínku (2) a skalární součin (3) bude vypadat následovně<br />

〈u,v〉 =<br />

∫<br />

+∞<br />

−∞<br />

e −x2 u(x)v(x)dx.<br />

Ortogonalizací posloupnosti 1,x,x 2 ,... vzhledem k výše uvedenému skalárnímu součinu<br />

získáváme takzvané Hermitovy polynomy.<br />

Hermitovy polynomy získané ortogonalizací mají tvar<br />

M (−∞,+∞)<br />

0 (x) = 1,<br />

M (−∞,+∞)<br />

1 (x) = x,<br />

M (−∞,+∞)<br />

2 (x) = x 2 − 1 2 ,<br />

M (−∞,+∞)<br />

3 (x) = x 3 − 3 2 x.<br />

...<br />

3.5 Vlastnosti kořenů ortogonálních polynomů<br />

Věta 3.1. Necht’ n ∈ N. Pak polynom p I n stupně n má právě n různých reálných kořenů, které<br />

navíc všechny leží v intervaluI.<br />

Důkaz. Necht’ α 1 < α 2 < ··· < α k , kde k ∈ N 0 , k ≤ n, jsou všechny kořeny polynomu p I n<br />

s lichou násobností, které leží v intervaluI. Pak můžeme polynomp I n přepsat do tvaru<br />

p I n(x) = (x−α 1 ) l 1<br />

(x−α 2 ) l 2<br />

...(x−α k ) l k<br />

s(x),<br />

kde l 1 ,l 2 ,...,l k jsou lichá čísla udávající násobnost jednotlivých kořenů α 1 ,α 2 , ...,α k .<br />

Polynomsmá přitom kořeny sudé násobnosti, případně kořeny, které se nachází mimo interval<br />

I, z čehož vyplývá, že na intervalu I nemění znaménko. Pro polynom s na intervalu I tedy platí<br />

jedna z variants(x) ≤ 0 nebos(x) ≥ 0. Nejprve předpokládejme, žes(x) ≥ 0 prox ∈ I.<br />

Dále definujme polynom<br />

r(x) = (x−α 1 )(x−α 2 )(x−α 3 )...(x−α k ).

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

Saved successfully!

Ooh no, something went wrong!