страница 1
7. Численное интегрирование
Если функция  непрерывна на отрезке  и известна ее первообразная  , то определенный интеграл от этой функции в пределах от  до  может быть вычислен по формуле Ньютона — Лейбница:
 , (1) где  .
Однако во многих случаях функция не может быть найдена с помощью элементарных средств или является слишком сложной, вследствие этого вычисление определенного интеграла по формуле (1) может быть практически невыполнимым. Кроме того, на практике подынтегральная функция часто задается таблично и тогда само понятие первообразной теряет смысл.
Поэтому особое значение приобретают численные (приближенные) методы вычисления определенных интегралов.
Постановка задачи.
Задача численного интегрирования функции f(x) заключается в вычислении значения определенного интеграла на основании ряда значений подынтегральной функции yk = f(xk).
Идея численного интегрирования вытекает из геометрического смысла определенного интеграла – значение определенного интеграла численно равно площади криволинейной трапеции, ограниченной графиком функции y=f(x), осью абсцисс и прямыми х=а, х=b. (рис 1)
y рис 1
f (x)
f(x
x
a b a b
Формально процедура численного интегрирования заключается в том, что отрезок [а, b] разбивается на n частичных отрезков, а затем подинтегральная функция заменяется легко интегрируемой функцией, интерполирующей значения подинтегральной функции в точках разбиения.
Наиболее часто приближенное значение интеграла ищут в виде линейной комбинации значений функции на отрезке :
 . (2)
Приближенное равенство (2) называют квадратурной формулой, определяемой узлами  и коэффициентами  .
Выражение в правой части (2) называют квадратурной суммой,
а разность  остаточным членом, или остатком, квадратурной формулы.
Рассмотрим простейшие примеры численных методов интегрирования.
7.1. Метод прямоугольников
Составим интегральную сумму для функции f(x) на сегменте [a,b] .
Для этого разобьем [a,b] на n равных частей точкками: x1, x2, … , xk, … , xn-1.
Длина каждой части равна
О бозначим значение подынтегральной функции f(x)
в узлах
Составим суммы (3)
При составлении первой суммы мы рассматриваем значения функции y=f(x) в точках, являющихся левыми концами частичных сегментов (рис 2а), а при составлении второй суммы – в точках, являющихся правыми концами (рис 2б).
y
yk-1
y1 B
y0 yn-1
A
0 a x1 xk-1 xn-1 b x
   y
  y1 B
  yk yn-1
    
  А
0 a x1 xk xn-1 b x
П  о определению интеграла:
и
Поэтому в качестве приближенного значения интеграла естественно  взять суммы
 
и (4)
т.е (5)
и (5')
Эти приближенные равенства называются формулами прямоугольников.
В случае, когда f(x)0, они означают, что площадь криволинейной трапеции aABb, ограниченной дугой кривой y=f(x), осью Ох и прямыми х=а и х=b, принимается приближенно равной площади ступенчатой фигуры, образованной
из n прямоугольников с основаниями и высотами yк,
Всякое приближенное вычисление имеет определенную ценность лишь тогда, когда оно сопровождается оценкой допущенной при этом погрешности. Поэтому формулы прямоугольников будут практически пригодны для приближенного вычисления интегралов лишь в том случае, если будет существовать удобный способ оценки получающейся при этом погрешности (при заданном n), позволяющий к тому же находить и число частей n разбиения сегмента, гарантирующее требуемую степень точности приближенного вычисления.
Будем предполагать, что функция f(x) имеет ограниченную производную на сегменте [a, b], так что существует такое число М>0, что для всех значений х из [a, b] выполняется неравенство |f'(x)| M. Качественный смысл этого неравенства заключается в том, что скорость изменения значения функции ограничена.
В реальных физических системах это требование практически всегда выполнено!
Тогда погрешность Rn, вычисления интеграла по формулам прямоугольников, может быть оценена по формуле:
|Rn| M(b-a)2/2n (6)
При неограниченном возрастании n выражение M(b-a)2/2n, а следовательно, и погрешность Rn стремится к нулю, т.е. точность приближения будет тем больше, чем на большее число частей будет разделен сегмент [a, b].
Следовательно, для вычисления интеграла с заданной точностью достаточно сегмент [a, b] разбить на число частей n = M(b-a)2/2 .
Метод прямоугольников – это наиболее простой и вместе с тем наиболее грубый метод приближенного интегрирования.
7.2. Метод трапеций
Очевидно, что чем больше число отрезков разбиения n , тем более точный результат дадут формулы (4 , 5). Однако увеличение числа разбиений интервала интегрирования не всегда возможно. Поэтому большой интерес представляют методы, дающие более точные результаты при том же числе точек разбиения.
П ростейшая из таких формул получается как среднее арифметическое правых частей формул (4):
(7)
Л егко усмотреть геометрический смысл этой формулы. Если на каждом отрезке разбиения дугу графика подинтегральной функции y=f(x) заменить стягивающей ее хордой (линейная интерполяция), то мы получим трапецию,
площадь которой равна
С
у0 у1 у2 уn-1 yn В
А
a x1 x2 xn-1 b
          ледовательно, формула (7) представляет собой площадь фигуры, состоящей из таких трапеций (рис. 3) .
И з геометрических соображений понятно, что площадь такой фигуры будет, вообще говоря, более точно выражать площадь криволинейной трапеции, нежели площадь ступенчатой фигуры, рассматриваемая в методе прямоугольников.
Приведя подобные, окончательно получим
(8)
Формулу (5) называют формулой трапеций.
Формулой трапеций часто пользуются для практических вычислений.
Для оценки погрешности Rn, возникающей при замене левой части (7) правой, доказывается, что абсолютная величина ее удовлетворяет неравенству:
(9)
где М2 – максимум модуля второй производной функции на отрезке [a,b], т.е.
Погрешность Rn будет меньше заданного числа 0, если взять
7.3. Метод парабол (Симпсона)
З  начительное повышение точности приближенных формул может быть достигнуто за счет повышения порядка интерполяции. Одним из таких методов приближенного интегрирования является метод парабол.
Идея метода исходит из того, что на частичном промежутке дуга некоторой параболы в общем случае теснее прилегает к кривой y=f(x), чем хорда, соединяющая концы дуги этой кривой, и поэтому значения площадей соответствующих элементарных трапеций, ограниченных «сверху» дугами парабол, являются более близкими к значениям площадей соответствующих частичных криволинейных трапеций, ограниченных сверху дугой кривой y=f(x), чем значения площадей соответствующих прямолинейных трапеций.
x2k-2 x2k-1 x2k
Сущность метода состоит в следующем.
Отрезок [a,b] делится на 2n
р
авных частей.
Пусть точки деления
х0=а, x1, x2, …x2n-2, x2n-1, x2n=b,
y0, y1, …y2n – соответствующие значения подинтегральной функции.
Произведем квадратичную интерполяцию данной подинтегральной функции на каждом из отрезков разбиения, т.е. заменим дугу графика подинтегральной функции дугой параболы (рис.4).
Приведем без вывода формулу парабол в окончательном виде:
(10)
Если подинтегральная функция имеет на отрезке [a,b] непрерывную четвертую производную, то для поправочного члена формулы (10) оценка
(11)
где М4- максимум модуля четвертой производной функции на отрезке [a,b].
Cравнивая между собой оценки, замечаем, что с увеличением n погрешность формулы прямоугольников уменьшается пропорционально величине n,
  формулы трапеций , а формулы парабол – пропорционально величине
Метод парабол сходится значительно быстрее метода трапеций, тогда как с точки зрения техники вычислений оба метода одинаковы.
страница 1
|