3/8-Regel y × × × × × × f1 · L1 (t) × 1 3 2 3 x0 a t 1 f1 · L1 (x) × x1 x2 x3 b Die Regel hat die Form (4 Stützstellen): Z b f (x) dx ≈ (b − a)[ w0 f0 + w1 f1 + w2 f2 + w3 f3 ] a Um die Gewichte wi zu ermitteln, approximieren wir f durch eine Summe von Lagrange-Polynomen. f (x) ≈ f0 · L0 (x) + f1 · L1 (x) + f2 · L2 (x) + f3 · L3 (x) mit L0 = (x − x1 )(x − x2 )(x − x3 ) , (x0 − x1 )(x0 − x2 )(x0 − x3 ) L1 , L2 , L3 Das ergibt: Z b Z b f (x) dx ≈ [f0 · L0 (x) + f1 · L1 (x) + f2 · L2 (x) + f3 · L3 (x) ] dx a a Z 1 = (b − a) [f0 · L0 (t) + f1 · L1 (t) + f2 · L2 (t) + f3 · L3 (t) ] dt 0 wk = Z 1 Lk (t) dt = 0 L0 = (Substitution, d. h. Verschiebung, Streckung) Z 1 0 1 2 3 1 2 3 j=3 Y 3t−j k−j dt es wurde mit 3 erweitert j=0 j 6= k (t − 3 )(t − 3 )(t − 3 ) (0 − 3 )(0 − 3 )(0 − 3 ) (3t − 1)(3t − 2)(3t − 3) = , (−1)(−2)(−3) Z 0 1 L0 (t) dt = . . . = 18 1 3/8-Regel Zb f (x) dx ≈ b −8 a [f (a) + 3f (a + b −3 a ) + 3f (a + 2 b −3 a ) + f (b) ] a Zb f (x) dx ≈ b −8 a [f0 + 3f1 + 3f2 + f3 ] a w3 = Z 1 L3 (t) dt = 0 w2 = Z Z 1 dt = 3−j 0 1 L2 (t) dt = 0 Z j=3 Y 3t−j j=0 j 6= 3 1 j=3 Y 3t−j dt = 2−j 0 j=0 j 6= 2 Z Z 1 0 1 0 3t·(3t−1)·(3t−2) dt 3·2·1 = 18 3t·(3t−1)·(3t−3) dt 2·1·(−1) = 38 allgemein für Polynome vom Grad N (N + 1 Stützstellen) wk = Z 1 Lk (t) dt = 0 Z 0 1 j=N Y N t−j k−j dt j=0 j 6= k Grad N = 3: N t (Nt − 1)(N t − 2)(N t − 3) − 2)(k − 3) k (k − 1)(k k = 1, L1 = k = 2, N t(N t − 1) (Nt − 2)(N t − 3) L2 = − 3) k(k − 1) (k −2)(k Nächste Newton-Cotes-Regel (5 Stützstellen), Milne-Regel Zb −a f (x) dx ≈ b 90 [7f (a) + 32f (a + b −4 a ) + 12f (a + 2 b −4 a ) + 32f (a + 3 b −4 a ) + 7f (b) ] a Zb −a f (x) dx ≈ b 90 [7f0 + 32f1 + 12f2 + 32f3 + 7f4 ] a w1 = Z 0 1 L1 (t) dt = Z 0 1 j=4 Y 4t−j 1−j j=0 j 6= 1 dt = Z 1 0 2 4t·(4t−2)·(4t−3)·(4t−4) dt 1·(−1)·(−2)·(−3) = 32 90 Siehe auch: Numerische Integration Lagrange-Interpolation Startseite 3
© Copyright 2024 ExpyDoc