ПРИБЛИЖЁННОЕ ВЫЧИСЛЕНИЕ ИНТЕГРАЛОВ

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

4.1. Концепция численного интегрирования

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

(4.1)

Чем меньше интервалы, на которых производят замену, тем точнее вы­чис­ляется интеграл. Поэтому исходный отрезок [а, b] для повышения точ­нос­ти делят на несколько равных или неравных интервалов, на каждом из кото­рых применяют формулу интегрирования, а затем складывают результаты.

Все методы различаются значениями ординат x_i, и весов _i.

В большинстве случаев погрешность численного интегрирования опре­де­­ляется путем двойного интегрирования: с исходным шагом (шаг опреде­ляется путем равномерного деления отрезка (b–а) на число отрезков n:
h = (b–а)/n)и с шагом, увеличенным в 2 раза. Разница вычисленных зна­чений интегралов определяет погрешность.

Сравнение эффективности различных методов проводится по степени полинома, который данным методом интегрируется точно, без ошибки. Чем выше степень такого полинома, тем выше точность метода, тем он эффективнее [4].

4.2. Метод прямоугольников

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

В этом случае подынтегральная функция заменяется горизонтальной прямой (у = с₀) со значением ординаты, т. е. значения функции соответст­венно слева или справа участка.

Вычисление определенного интеграла (геометрическая интерпретация опре­­деленного интеграла) – это вычисление площади криволинейной трапеции.

Существуют методы: левых прямоугольников, правых прямоугольников и средних прямоугольников [4].

Разбиваем отрезок [а, b] на n частей с равномерным шагом h.

Метод левых прямоугольников

Формулы интегрирования:

• для одного участка интегрирования (рис. 4.1, a):

; (4.2)

(4.3)

 для n участков интегрирования (рис. 4.1, б):

; (4.4)

(4.5)

а	б
Рис. 4.1. Иллюстрация метода левых прямоугольников: а – для одного участка интегрирования; б – для n участков интегрирования

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

Формулы интегрирования:

• для одного участка интегрирования (рис. 4.2, a):

; (4.6)

(4.7)

• для n участков интегрирования (рис. 4.2, б):

; (4.8)

(4.9)

аб

	y		y

Рис. 4.2. Иллюстрация метода правых прямоугольников:
a – для одного участка интегрирования; б – для n участков интегрирования

Метод средних прямоугольников

Формулы интегрирования:

• для одного участка интегрирования (рис. 4.3, a):

; (4.10)

• для n участков интегрирования (рис. 4.3, б):

; (4.11)

где x₀ = a, x₁ = b, h = b–a.

	a

Нетрудно заметить, что в методе прямоугольников интеграл вычисляется абсолютно точно только при f(х) = с (соnst).

На рис. 4.4 для сравнения приведены примеры прямоугольников при различном числе участков. Наглядно видно, что площадь всех прямоугольников на правом рисунке отличается меньше от площади под кривой f(х), чем на левом.

Метод прямоугольников не находит практического применения в силу значительных погрешностей, что подтверждается рис. 4.4.

a б

Рис. 4.4. Иллюстрация метода левых прямоугольников:
а – с 3 участками разбиения отрезка интегрирования [а, b];
б – с 6 участками разбиения отрезка интегрирования [а, b]

4.3. Метод трапеций

Метод трапеций так же, как и метод прямоугольников, относится к простейшим методам. В этом случае подынтегральная функция заменяется наклонной прямой (у = с₁х + с_о) [2].

На рис. 4.5 приведен пример вычисления интеграла методом трапеций. По сравнению с методом прямоугольников метод трапеций более точен, так как трапеция точнее заменяет соответствующую криволинейную трапецию, чем прямоугольник.

Нетрудно заметить, что в методе трапеций интеграл вычисляется абсолютно точно только при f(х) линейной или кусочно-линейной.

Формулы интегрирования:

• для одного участка интегрирования:

; (4.12)

(4.13)

где x₀ = a, x₁ = b, h = b–a;

• для n участков интегрирования (рис. 4.5):

_(4.14)

; (4.15)

(4.16)

Погрешность R вычисления интеграла методом трапеций при использовании двойного просчета на практике может быть определена из следующего соотношения:

(4.17)

где I_n и I_n/2 – соответственно значения интеграла при числе разбиений n и n/2.

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

Пример.

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

Решение.