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.
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