ВУЗ:
Рубрика:
26 ИНТЕГРАЛЬНОЕ ИСЧИСЛЕНИЕ
Значит, как это следует из примера 25, 0 < K 6
π
2
. С помощью очень непростых искусственных
приёмов, которые выходят за рамки нашего курса, можно показать, что
Z
+∞
0
e
−x
2
dx =
√
π
2
. (76)
6. Приближённые вычисления определённого интеграла
Поскольку далеко не все (а точнее, большинство) определённые интегралы можно вычислить,
используя известные (элементарные) функции, для практических целей используют приближён-
ные методы.
Итак, пусть требуется вычислить интеграл
Z
b
a
f(x) dx; (77)
мы будем предполагать, что (77 ) — собственный интеграл.
Метод трапеций. Этот метод основан на разбиении отрезка интегрирования [a, b] на n равных
частей,
a = x
0
, x
1
, . . . , x
n
= b, x
i+1
− x
i
= h =
b − a
n
,
и замене каждой «маленькой» криволинейной трапеции на участке [x
i
, x
i+1
] «настоящей» прямо-
линейной. В силу определения 2, площадь получаемой фигуры приблизительно равна интегра-
лу (77), так что мы имеем
Z
b
a
f(x) dx ∼
b − a
n
·
f(a) + f (b)
2
+ f(x
1
) + ··· + f (x
n−1
)
. (78)
Формула (78) называется формулой т рапеций.
Поскольку формула трапеций даёт лишь оценку, приближённое значение искомого интеграла,
важно оценить, какова её точность. Пусть ε — разность между левой и правой частями в соот-
ношении (78). Тогда можно показать, что если подынтегральная функция имеет непрерывные
первую и вторую производные на отрезке [a, b], то справедливо неравенство
|ε| 6
|b − a|
3
12n
2
· max
[a,b]
|f
′′
(x)|. (79)
Это и есть оценка точности вычисления интеграла (77) по формуле трапеций.
Метод Симпсона. Метод Симпсона также основан на разбиении отрезка интегрирования на n
равных частей, а затем каждого из полученных отрезков — пополам (в итоге весь отрезок разбива-
ется на 2n частей). Интегрируемая функция заменяется в этом случае дугой параболы, проходя-
щей через три соседних точки с номерами 2i, 2i+1 и 2i+2. Получаемое при этом приблизительное
значение интеграла имеет вид
Z
b
a
f(x) dx ∼
b − a
6n
·
(f(a) + f(b)) + +2(f (x
2
) + ··· + f (x
2n−2
)) + 4(f(x
1
) + ··· + f (x
2n−1
)
. (80)
Это — так называемая параболическая формула, или формула Симпсона.
Если подынтегральная функция обладает первыми четырьмя непрерывными производными,
то точность вычислений оценивается по формуле
|ε| 6
|b − a|
5
180 · (2n)
4
· max
[a,b]
|f
(4)
(x)|. (81)
26 ИНТЕГРАЛЬНОЕ ИСЧИСЛЕНИЕ Значит, как это следует из примера 25, 0 < K 6 π2 . С помощью очень непростых искусственных приёмов, которые выходят за рамки нашего курса, можно показать, что Z +∞ √ −x2 π e dx = . (76) 0 2 6. Приближённые вычисления определённого интеграла Поскольку далеко не все (а точнее, большинство) определённые интегралы можно вычислить, используя известные (элементарные) функции, для практических целей используют приближён- ные методы. Итак, пусть требуется вычислить интеграл Z b f (x) dx; (77) a мы будем предполагать, что (77) — собственный интеграл. Метод трапеций. Этот метод основан на разбиении отрезка интегрирования [a, b] на n равных частей, b−a a = x0 , x1 , . . . , xn = b, xi+1 − xi = h = , n и замене каждой «маленькой» криволинейной трапеции на участке [xi , xi+1 ] «настоящей» прямо- линейной. В силу определения 2, площадь получаемой фигуры приблизительно равна интегра- лу (77), так что мы имеем Z b b − a f (a) + f (b) f (x) dx ∼ · + f (x1 ) + · · · + f (xn−1 ) . (78) a n 2 Формула (78) называется формулой трапеций. Поскольку формула трапеций даёт лишь оценку, приближённое значение искомого интеграла, важно оценить, какова её точность. Пусть ε — разность между левой и правой частями в соот- ношении (78). Тогда можно показать, что если подынтегральная функция имеет непрерывные первую и вторую производные на отрезке [a, b], то справедливо неравенство |b − a|3 |ε| 6 · max |f ′′ (x)|. (79) 12n2 [a,b] Это и есть оценка точности вычисления интеграла (77) по формуле трапеций. Метод Симпсона. Метод Симпсона также основан на разбиении отрезка интегрирования на n равных частей, а затем каждого из полученных отрезков — пополам (в итоге весь отрезок разбива- ется на 2n частей). Интегрируемая функция заменяется в этом случае дугой параболы, проходя- щей через три соседних точки с номерами 2i, 2i+1 и 2i+2. Получаемое при этом приблизительное значение интеграла имеет вид Z b b−a f (x) dx ∼ · (f (a) + f (b)) + +2(f (x2 ) + · · · + f (x2n−2 )) + 4(f (x1 ) + · · · + f (x2n−1 ) . (80) a 6n Это — так называемая параболическая формула, или формула Симпсона. Если подынтегральная функция обладает первыми четырьмя непрерывными производными, то точность вычислений оценивается по формуле |b − a|5 |ε| 6 · max |f (4) (x)|. (81) 180 · (2n)4 [a,b]