Формула трапеций.

Приближенное вычисление определенных интегралов

 

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

 

Формула трапеций.

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

Где f(x_k-1) и f(x_k) - соответственно основания трапеций; x_k - x_k-1 = (b-a)/n - их высоты.

Таким образом, получена приближенная формула

 

которая и называется формулой трапеций. Эта формула тем точнее, чем больше n.

Рассмотрим в качестве примера интеграл . Точное значение этого интеграла находится просто:

Вычислим теперь по формуле трапеций его приближенное значение. Пусть n=5. Тогда имеем: a=x₀=0, x₁=0,2, x₂=0,4, x₃=0,6, x₄=0,8, x₅=1=b и соответственно f(x₀)=0, f(x₁)=0,04, f(x₂)=0,16, f(x₃)=0,36, f(x₄)=0,64, f(x₅)=1. Следовательно,

Точное значение интеграла равно 0,3333...., поэтому абсолютная ошибка меньше 0,007. Во многих технических задач эта точность достаточна.

Если увеличить число n, то точность будет большей. Так, например, при n=10

т.е. абсолютная ошибка меньше 0,002.

В более полных курсах высшей математики доказывается, что если функция f(x) имеет на [a, b] непрерывную вторую производную, то абсолютная величина погрешности формулы трапеций не больше, чем

где k -наибольшее значение на отрезке [a, b].

Следует отметить, что с увеличением n увеличивается не только точность вычисления определенного интеграла, но и объем вычислительной работы. Однако здесь на помощь приходят ЭВМ.

Вычислим по формуле трапеции интеграл при n=10. Разобьем отрезок [0, 1] на 10 равных частей точками х₀=0, х₁=0,1, ..., х₉=0,9, х₁₀=1. Вычислим приближенно значения функции f(x)= в этих точках: f(0)=1,0000, f(0,1)=0.9091, f(0,2)=0,8333, f(0,3)=0.7692, f(0,4)=0,7143, f(0,5)=0,6667, f(0,6)=0,6250, f(0,7)=0,5882, f(0,8)= 0,5556, f(0,9)=0,5263, f(1)=0,5000.

По формуле трапеций получаем

 

Оценим погрешность полученного результата. Так как f(x)=1/(1+x), то На отрезке [0, 1] имеем . Поэтому погрешность полученного результата не превосходит величины

Вычислим точное значение данного интеграла по формуле Ньютона-Лейбница:

Абсолютная ошибка результата, полученного по формуле трапеций, меньше 0,0007. Это находится в соответствии с данной выше оценкой погрешности.

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