Метод прямоугольников. Методы численного интегрирования

Методы численного интегрирования

Идея численного интегрирования предельно проста и вытекает из геометрического смысла определенного интеграла – значение определенного интеграла численно равно площади криволинейной трапеции, ограниченной графиком функции y=f(x), осью абсцисс и прямыми х=а, х=b. Находя приближенно площадь криволинейной трапеции, мы получаем значение интеграла. Формально процедура численного интегрирования заключается в том, что отрезок [а, b] разбивается на n частичных отрезков, а затем подынтегральная функция заменяется на нем легко интегрируемой функцией, по определенной зависимости интерполирующей значения подынтегральной функции в точках разбиения. Рассмотрим теперь простейшие из численных методов интегрирования.

Итак, функция у=f(x) интегрируема на сегменте [a,b] и требуется вычислить ее интеграл . Составим интегральную сумму для f(x) на сегменте [a,b]. Для этого разобьем сегмент [a,b] на n равных между собой частей с помощью точек: x₁, x₂, …, x_k, …, x_n-1.

Если длину каждой части мы обозначим через Dх, так что , то для каждой точки x_k будем иметь: x_k = a + kDx (k=0, 1, 2, … n).

Обозначим теперь через y_k значение подынтегральной функции f(x) при x = x_k = a + kDx, то есть положим (k=0, 1, … n).

Тогда суммы и будут интегральными для функции f(x) на отрезке [a,b]. (При составлении первой суммы мы рассматриваем значения функции y=f(x) в точках, являющихся левыми концами частичных сегментов, а при составлении второй суммы – в точках, являющихся правыми концами этих сегментов.)

По определению интеграла имеем:

и .

Поэтому в качестве приближенного значения естественно

взять интегральную сумму и , т.е. положить:

,

а также

,

т.е.

(1)

и

(2)

Эти приближенные равенства называются формулами прямоугольников.

В том случае, когда f(x)³0, формулы (1) и (2) с геометрической точки зрения означают, что площадь криволинейной трапеции aABb, ограниченной дугой кривой y=f(x), осью Ох и прямыми х=а и х=b, принимается приближенно равной площади ступенчатой фигуры, образованной из n прямоугольников с основаниями и высотами: y₀, y₁, y₂, … y_n-1 – в случае формулы (1) (рис. 8) и y₁, y₂, y₃, … y_n – в случае формулы (2) (рис. 9).

Рис. 8 Рис. 9

Исходя из приведенного выше геометрического смысла формул (1) и (2) способ приближенного вычисления определенного интеграла по этим формулам принято называть методом прямоугольников.

Всякое приближенное вычисление имеет определенную ценность лишь тогда, когда оно сопровождается оценкой допущенной при этом погрешности. Поэтому формулы прямоугольников будут практически пригодны для приближенного вычисления интегралов лишь в том случае, если будет существовать удобный способ оценки получающейся при этом погрешности (при заданном n), позволяющий к тому же находить и число частей n разбиения сегмента, гарантирующее требуемую степень точности приближенного вычисления.

Будем предполагать, что функция f(x) имеет ограниченную производную на сегменте [a, b], так что существует такое число М>0, что для всех значений х из [a, b] выполняется неравенство |f'(x)| £ M. Качественный смысл этого неравенства заключается в том, что скорость изменения значения функции ограничена. В реальных природных системах это требование практически всегда выполнено. В этих условиях абсолютная величина погрешности R_n, которую мы допускаем, вычисляя интеграл по формуле прямоугольников, может быть оценена по формуле [27]:

| R_n| £ M(b-a)²/2n. (3)

При неограниченном возрастании n выражение M(b-a)²/2n, а следовательно, и абсолютная величина погрешности R_n будет стремиться к нулю, т.е. точность приближения будет тем больше, чем на большее число равных частей будет разделен сегмент [a, b]. Абсолютная погрешность результата будет заведомо меньше заданного числа e>0, если взять

n > M(b-a)²/2e.

Следовательно, для вычисления интеграла с указанной степенью точности достаточно сегмент [a, b] разбить на число частей, большее числа M(b-a)²/2e. [27].

Метод прямоугольников – это наиболее простой и вместе с тем наиболее грубый метод приближенного интегрирования. Заметно меньшую погрешность дает другой метод – метод трапеций.

Дата добавления: 2014-12-16