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.

6 GAUSSOVA-LEGENDROVA KVADRATURNÍ FORMULE 48<br />

a získáváme<br />

|R(f)| ≤<br />

k∑<br />

j=1<br />

Ω<br />

(2n)!<br />

b−a<br />

∫k<br />

0<br />

(<br />

1 +1<br />

2<br />

u− x〈−1,1〉<br />

(<br />

u− x〈−1,1〉<br />

) 2<br />

b−a<br />

k<br />

2 +1<br />

2<br />

) 2<br />

b−a<br />

...(<br />

k<br />

n +1<br />

2<br />

u−<br />

x〈−1,1〉<br />

) 2<br />

b−a<br />

du.<br />

k<br />

Ze všech n závorek (kořeny Legendrova polynomu mají indexy 1,2,...,n) můžeme vytknout<br />

zlomek b−a<br />

2k a získáváme<br />

|R(f)| ≤ (b−a)2n<br />

2 2n k 2n k∑<br />

j=1<br />

Opět zavedeme substituci<br />

Ω<br />

(2n)!<br />

b−a<br />

∫k<br />

0<br />

( ) 2ku 2<br />

b−a −x〈−1,1〉 1 −1<br />

( 2ku 2 ( 2ku 2<br />

b−a −x〈−1,1〉 2 −1)<br />

...<br />

b−a −x〈−1,1〉 n −1)<br />

du.<br />

b−a<br />

2k<br />

v = 2ku<br />

b−a −1<br />

dv =du<br />

a dostaneme<br />

b−a<br />

k<br />

0 ↦→ 0−1 = −1<br />

↦→ 2k(b−a)<br />

k(b−a)<br />

−1 = 2−1 = 1<br />

|R(f)| ≤ (b−a)2n+1<br />

2 2n+1 k 2n+1 k∑<br />

j=1<br />

Ω<br />

(2n)!<br />

∫ 1<br />

−1<br />

(<br />

v −x 〈−1,1〉<br />

1<br />

) 2 (<br />

v −x 〈−1,1〉<br />

2<br />

) 2... (<br />

v −x 〈−1,1〉<br />

n<br />

) 2<br />

dv.<br />

V souladu s dokazovanou větou označíme proměnnou v jako x a uzly x 〈−1,1〉<br />

1 ,x 〈−1,1〉<br />

2 , ...,<br />

x n 〈−1,1〉 jakox 1 ,x 2 ,...,x n .<br />

Všimněme si, že hodnota integrálu je nezávislá naj, takže ji můžeme vytknout. Totéž platí<br />

a Ω, takže obdržíme vztah<br />

ipro<br />

1<br />

(2n)!<br />

|R(f)| ≤ Ω (b−a) 2n+1 ∫1<br />

(2n)! 2 2n+1 k 2n (x−x 1 ) 2 (x−x 2 ) 2 ...(x−x n ) 2 dx.<br />

−1

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

Saved successfully!

Ooh no, something went wrong!