Екатеринбург

Вычисление определенного интеграла

Введение

Задача численного интегрирования функций заключается в вычислении приближенного значения определенного интеграла:

на основе ряда значений подынтегральной функции.{ f(x) |x=x k = f(x k) = y k }.

Формулы численного вычисления однократного интеграла называются квадратурными формулами, двойного и более кратного – кубатурными.

Обычный прием построения квадратурных формул состоит в замене подынтегральной функции f(x) на отрезке интерполирующей или аппроксимирующей функцией g(x) сравнительно простого вида, например, полиномом, с последующим аналитическим интегрированием. Это приводит к представлению

В пренебрежении остаточным членом R[f] получаем приближенную формулу

Обозначим через y i = f(x i) значение подинтегральной функции в различных точках

В качестве приближенной функции g(x) рассмотрим интерполяционный полином на

Формула (1) дает

В формуле (2) величины {

Подведем итог.

2. В квадратурных формулах интерполяционного типа остаточный член R n [f] может быть представлен в виде значения конкретного дифференциального оператора на функции f(x). Для

3. Для полиномов до порядка n включительно квадратурная формула (2) точна, т.е.

Рассмотрим частные случаи формул (2) и (3): метод прямоугольников, трапеций, парабол (метод Симпсона). Названия этих методов обусловлены геометрической интерпретацией соответствующих формул.

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

Определенный интеграл функции от функции f(x):

Метод, представленный формулой (4), называется методом левых прямоугольников, а метод, представленный формулой(5) – методом правых прямоугольников:

Для нахождения определенного интеграла методом средних прямоугольников площадь, ограниченная прямыми a и b, разбивается на n прямоугольников с одинаковыми основаниями h, высотами прямоугольников будут точки пересечения функции f(x) с серединами прямоугольников (h/2). Интеграл будет численно равен сумме площадей n прямоугольников (рисунок 3).

n – количество разбиений отрезка .

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

Для нахождения определенного интеграла методом трапеций площадь криволинейной трапеции также разбивается на n прямоугольных трапеций с высотами h и основаниями у 1 , у 2 , у 3 ,..у n , где n - номер прямоугольной трапеции. Интеграл будет численно равен сумме площадей прямоугольных трапеций (рисунок 4).

n – количество разбиений

Погрешность формулы трапеций оценивается числом

Погрешность формулы трапеций с ростом

Формула Симпсона

Если для каждой пары отрезков

5.3 Метод трапеций

Выведем формулу трапеций так же, как и формулу прямоугольников, из геометрических соображений. Заменим график функции y = f(x) (рис.5.1) ломаной линией (рис.5.7), полученной следующим образом. Из точек a = x 0 , x 1 , x 2 ,…, x n = b проведем ординаты до пересечения с кривой y = f(x). Концы ординат соединим прямолинейными отрезками.

Тогда площадь криволинейной трапеции приближенно можно считать равной площади фигуры, составленной из трапеций. Так как площадь трапеции, построенной на отрезке длины h = , равна h , то, пользуясь этой формулой для i = 0, 2, … , n – 1, получим квадратурную формулу трапеций:

I=»I тр =h= (5.7)

Оценка погрешности. Для оценки погрешности формулы трапеций воспользуемся следующей теоремой.

Теорема 5.2. Пусть функция f дважды непрерывно дифференцируема на отрезке . Тогда для формулы трапеций справедлива следующая оценка погрешности:

| I – I тр | £ h 2 , (5.8)

где M 2 = |f "(x)|.

Пример 5.2.

Вычислим значение интеграла по формуле трапеций (5.7) и сравним полученный результат с результатом примера 5.1.

Используя таблицу значений функции eиз примера 5.1 и производя вычисления по формуле трапеций (5.7), получим: I тр = 0.74621079.

Оценим погрешность полученного значения. В примере (5.1) получили оценку: | f "(x)| £ M 2 = 2. Поэтому по формуле (5.8)

I – I тр | £ (0.1) 2 » 1.7× 10 -3 .

Сравнивая результаты примеров 5.1 и 5.2, видим, что метод средних прямоугольников имеет меньшую погрешность, т.е. он более точный.

5.4 Метод Симпсона (метод парабол)

Заменим график функции y = f(x) на отрезке , i = 0, 2, … , n – 1, параболой, проведенной через точки (x i , f(x i)), (x,f(x)), (x i+ 1 , f(x i+ 1)), где x - середина отрезка . Эта парабола есть интерполяционный многочлен второй степени L 2 (x) с узлами x i , x, x i+ 1 . Нетрудно убедиться, что уравнение этой параболы имеет вид:

f(x) + (x – x) + (x - x) 2 , (5.9)

Проинтегрировав функцию (5.9) на отрезке , получим

I i = » = (f(x i) + 4f(x) + f(x i+ 1)). (5.10)

Суммируя выражение (5.10) по i = 0, 1, 2, … , n – 1, получим квадратурную формулу Симпсона (или формулу парабол):

I =» I С = (f(x 0) + f(x n) + 4 + 2). (5.11)

Оценка погрешности. Для оценки погрешности формулы Симпсона воспользуемся следующей теоремой.

Теорема 5.2. Пусть функция f имеет на отрезке непрерывную производную четвертого порядка f (4) (x). Тогда для формулы Симпсона (5.9) справедлива следующая оценка погрешности:

| I – I С | £ h 4 , (5.12)

где M 4 = | f (4) (x)|.

Замечание. Если число элементарных отрезков, на которые делится отрезок , четно, т.е. n = 2m, то параболы можно проводить через узлы с целыми индексами, и вместо элементарного отрезка длины h рассматривать отрезок длины 2h. Тогда формула Симпсона примет вид:

I » (f(x 0) + f(x 2m) + 4 + 2), (5.13)

а вместо оценки (5.10) будет справедлива следующая оценка погрешности:

| I – I С | £ h 4 , (5.14)

Пример 5.3.

Вычислим значение интеграла по формуле Симпсона (5.11) и сравним полученный результат с результатами примеров 5.1 и 5.2.

Используя таблицу значений функции eиз примера 5.1 и производя вычисления по формуле Симпсона (5.11) , получим:

I С = 0.74682418.

Оценим погрешность полученного значения. Вычислим четвертую производную f (4) (x).

f (4) (x) = (16x 4 – 48x 2 + 12) e, | f (4) (x)| £ 12.

| I – I С | £ (0.1) 4 » 0.42 × 10 -6 .

Сравнивая результаты примеров 5.1, 5.2 и 5.3, видим, что метод Симпсона имеет меньшую погрешность, чем метод средних прямоугольников и метод трапеций.

Учебно-воспитательные задачи:

• Дидактическая цель. Познакомить учащихся с методами приближённого вычисления определённого интеграла.
• Воспитательная цель. Тема данного занятия имеет большое практическое и воспитательное значение. Наиболее просто к идее численного интегрирования можно подойти, опираясь на определение определённого интеграла как предела интегральных сумм. Например, если взять какое-либо достаточно мелкое разбиение отрезка [a ; b ] и построить для него интегральную сумму, то её значение можно приближённо принять за значение соответствующего интеграла. При этом важно быстро и правильно производить вычисления с привлечением вычислительной техники.

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

Обеспечение занятия

• Раздаточный материал. Карточки-задания для самостоятельной работы.
• ТСО. Мультипроектор, ПК, ноутбуки.
• Оснащение ТСО. Презентации: “Геометрический смысл производной”, “Метод прямоугольников”, “Метод трапеций”. (Презентации можно взять у автора).
• Вычислительные средства: ПК, микрокалькуляторы.
• Методические рекомендации

Вид занятия. Интегрированное практическое.

Мотивация познавательной деятельности учащихся. Очень часто приходится вычислять определённые интегралы, для которых невозможно найти первообразную. В этом случае применяют приближённые методы вычисления определённых интегралов. Иногда приближённый метод применяют и для “берущихся” интегралов, если вычисление по формуле Ньютона-Лейбница не рационально. Идея приближённого вычисления интеграла заключается в том, что кривая заменяется новой, достаточно “близкой” к ней кривой. В зависимости от выбора новой кривой можно использовать ту или иную приближённую формулу интегрирования.

Последовательность занятия.

1. Формула прямоугольников.
2. Формула трапеций.
3. Решение упражнений.

План занятия

1. Повторение опорных знаний учащихся.

Повторить с учащимися: основные формулы интегрирования, сущность изученных методов интегрирования, геометрический смысл определённого интеграла.

1. Выполнение практической работы.

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

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

Простейшим приближённым методом является метод прямоугольников. Геометрически идея способа вычисления определённого интеграла по формуле прямоугольников состоит в том, что площадь криволинейной трапеции АВСD заменяется суммой площадей прямоугольников, одна сторона которых равна , а друга - .

Если суммировать площади прямоугольников, которые показывают площадь криволинейной трапеции с недостатком [Рисунок1], то получим формулу:

[Рисунок1]

то получим формулу:

Если с избытком

[Рисунок2],

то

Значения у 0 , у 1 ,..., у n находят из равенств , к = 0, 1..., n .Эти формулы называются формулами прямоугольников и дают приближённый результат. С увеличением n результат становится более точным.

Итак, чтобы найти приближённое значение интеграла , нужно:

Для того, чтобы найти погрешность вычислений, надо воспользоваться формулами:

Пример 1. Вычислить по формуле прямоугольников . Найти абсолютную и относительную погрешности вычислений.

Разобьём отрезок [a, b ] на несколько (например, на 6) равных частей. Тогда а = 0, b = 3 ,

х k = a + k х
х 0 = 2 + 0 = 2
х 1 = 2 + 1 = 2,5
х 2 = 2 + 2 =3
х 3 = 2 + 3 = 3
х 4 = 2 + 4 = 4
х 5 = 2 + 5 = 4,5

f (x 0) = 2 2 = 4
f (x 1) = 2 ,5 2 = 6,25
f (x 2) = 3 2 = 9
f (x 3) = 3,5 2 = 12,25
f (x 4) = 4 2 = 16
f (x 5) = 4,5 2 = 20,25.

По формуле (1):

Для того, чтобы вычислить относительную погрешность вычислений, надо найти точное значение интеграла:

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

Для нахождения определённого интеграла методом прямоугольников необходимо ввести значения подынтегральной функции f(x) в рабочую таблицу Excel в диапазоне х с заданным шагом х = 0,1.

1. Составляем таблицу данных (х и f(x)). х f(x). Аргумент , а в ячейку В1 – слово Функция 2 2,1 ). Затем, выделив блок ячеек А2:А3, автозаполнением получаем все значения аргумента (за правый нижний угол блока протягиваем до ячейки А32, до значения х=5 ).
2. Далее вводим значения подынтегральной функции. В ячейку В2 необходимо записать её уравнение. Для этого табличный курсор необходимо установить в ячейку В2 и с клавиатуры ввести формулу =А2^2 (при английской раскладке клавиатуры). Нажимаем клавишу Enter . В ячейке В2 появляется 4 . Теперь необходимо скопировать функцию из ячейки В2. Автозаполнением копируем эту формулу в диапазон В2:В32.
   В результате должна быть получена таблица данных для нахождения интеграла.
3. Теперь в ячейке В33 может быть найдено приближённое значение интеграла. Для этого в ячейку В33 вводим формулу = 0,1*, затем вызываем Мастер функций (нажатием на панели инструментов кнопки Вставка функции (f(x)) . В появившемся диалоговом окне Мастер функции-шаг 1 из 2 слева в поле Категория выбираем Математические. Справа в поле Функция - функцию Сумм. Нажимаем кнопку ОК. Появляется диалоговое окно Сумм. В рабочее поле мышью вводим диапазон суммирования В2:В31. Нажимаем кнопку ОК. В ячейке В33 появляется приближённое значение искомого интеграла с недостатком (37,955 ) .

Сравнивая полученное приближённое значение с истинным значением интеграла (39 ), можно видеть, что ошибка приближения метода прямоугольников в данном случае равна

= |39 - 37 , 955| = 1 ,045

Пример 2. Используя метод прямоугольников, вычислить с заданным шагом х = 0,05.

Сравнивая полученное приближённое значение с истинным значением интеграла , можно видеть, что ошибка приближения метода прямоугольников в данном случае равна

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

[Рисунок3]

Пример 3. Методом трапеций найти с шагом х = 0,1.

1. Открываем чистый рабочий лист.
2. Составляем таблицу данных (х и f(x)). Пусть первый столбец будет значениями х , а второй соответствующими показателями f(x). Для этого в ячейку А1 вводим слово Аргумент , а в ячейку В1 – слово Функция . В ячейку А2 вводится первое значение аргумента – левая граница диапазона (0 ). В ячейку А3 вводится второе значение аргумента – левая граница диапазона плюс шаг построения (0,1 ). Затем, выделив блок ячеек А2:А3, автозаполнением получаем все значения аргумента (за правый нижний угол блока протягиваем до ячейки А33, до значения х=3,1 ).
3. Далее вводим значения подынтегральной функции. В ячейку В2 необходимо записать её уравнение (в примере синуса). Для этого табличный курсор необходимо установить в ячейку В2. Здесь должно оказаться значение синуса, соответствующее значению аргумента в ячейке А2. Для получения значения синуса воспользуемся специальной функцией: нажимаем на панели инструментов кнопку Вставка функции f(x) . В появившемся диалоговом окне Мастер функции-шаг 1 из 2 слева в поле Категория выбираем Математические. Справа в поле Функция - функцию SIN . Нажимаем кнопку ОК. Появляется диалоговое окно SIN . Наведя указатель мыши на серое поле окна, при нажатой левой кнопке сдвигаем поле вправо, чтобы открыть столбец данных (А ). Указываем значение аргумента синуса щелчком мыши на ячейке А2. Нажимаем кнопку ОК. В ячейке В2 появляется 0. Теперь необходимо скопировать функцию из ячейки В2. Автозаполнением копируем эту формулу в диапазон В2:В33. В результате должна быть получена таблица данных для нахождения интеграла.
4. Теперь в ячейке В34 может быть найдено приближённое значение интеграла по методу трапеций. Для этого в ячейку В34 вводим формулу = 0,1*((В2+В33)/2+, затем вызываем Мастер функций (нажатием на панели инструментов кнопки Вставка функции (f(x)) . В появившемся диалоговом окне Мастер функции-шаг 1 из 2 слева в поле Категория выбираем Математические. Справа в поле Функция - функцию Сумм. Нажимаем кнопку ОК. Появляется диалоговое окно Сумм. В рабочее поле мышью вводим диапазон суммирования В3:В32. Нажимаем кнопку ОК и ещё раз ОК. В ячейке В34 появляется приближённое значение искомого интеграла с недостатком (1,997 ) .

Сравнивая полученное приближённое значение с истинным значением интеграла можно видеть, что ошибка приближения метода прямоугольников в данном случае вполне приемлемая для практики.

1. Решение упражнений.

Как вычислить определенный интеграл
по формуле трапеций и методом Симпсона?

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

Что нужно знать, чтобы освоить данные методы? Прозвучит забавно, но можно вообще не уметь брать интегралы. И даже вообще не понимать, что такое интегралы. Из технических средств потребуется микрокалькулятор. Да-да, нас ждут рутинные школьные расчёты. А еще лучше – закачайте мой калькулятор-полуавтомат для метода трапеций и метода Симпсона . Калькулятор написан в Экселе и позволит в десятки раз уменьшить время решения и оформления задач. Для экселевских чайников прилагается видеомануал! К слову, первая видеозапись с моим голосом.

Сначала зададимся вопросом, а зачем вообще нужны приближенные вычисления? Вроде бы можно найти первообразную функции и использовать формулу Ньютона-Лейбница, вычислив точное значение определенного интеграла. В качестве ответа на вопрос сразу рассмотрим демонстрационный пример с рисунком.

Вычислить определенный интеграл

Всё было бы хорошо, но в данном примере интеграл не берётся – перед вами неберущийся, так называемый интегральный логарифм . А существует ли вообще этот интеграл? Изобразим на чертеже график подынтегральной функции :

Всё нормально. Подынтегральная функция непрерывна на отрезке и определенный интеграл численно равен заштрихованной площади. Да вот только одна загвоздка – интеграл не берётся. И в подобных случаях на помощь как раз приходят численные методы. При этом задача встречается в двух формулировках:

1) Вычислить определенный интеграл приближенно, округляя результат до определённого знака после запятой . Например, до двух знаков после запятой, до трёх знаков после запятой и т.д. Предположим, получился приближенный ответ 5,347. На самом деле он может быть не совсем верным (в действительности, скажем, более точный ответ 5,343). Нашазадача состоит лишь в том , чтобы округлить результат до трёх знаков после запятой.

2) Вычислить определенный интеграл приближенно, с определённой точностью . Например, вычислить определённый интеграл приближенно с точностью до 0,001. Что это значит? Это значит, мы должны отыскать такое приближенное значение, которое по модулю (в ту или другую сторону) отличается от истины не более чем на 0,001.

Существуют несколько основных методов приближенного вычисления определенного интеграла, который встречается в задачах:

Отрезок интегрирования разбивается на несколько частей и строится ступенчатая фигура, которая по площади близка к искомой площади:

Не судите строго за чертежи, точность не идеальна – они лишь помогают понять суть методов.

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

Таким образом, наша площадь (синяя штриховка) приближается суммой площадей трапеций (красный цвет). Отсюда и название метода. Легко заметить, что метод трапеций даёт значительно лучшее приближение, чем метод прямоугольников (при одинаковом количестве отрезков разбиения). И, естественно, чем больше более мелких промежуточных отрезков мы рассмотрим, тем будет выше точность. Метод трапеций время от времени встречается в практических заданиях, и в данной статье будет разобрано несколько примеров.

Метод Симпсона (метод парабол) . Это более совершенный способ – график подынтегральной функции приближается не ломаной линией, а маленькими параболками. Сколько промежуточных отрезков – столько и маленьких парабол. Если взять те же три отрезка, то метод Симпсона даст ещё более точное приближение, чем метод прямоугольников или метод трапеций.

Чертеж строить не вижу смысла, поскольку визуально приближение будет накладываться на график функции (ломаная линия предыдущего пункта – и то практически совпала).

Задача на вычисление определенного интеграла по формуле Симпсона – самая популярное задание на практике. И методу парабол будет уделено значительное внимание.

Как вычислить определенный интеграл методом трапеций?

Сначала формула в общем виде. Возможно, она будет не всем и не сразу понятна… да Карлссон с вами – практические примеры всё прояснят! Спокойствие. Только спокойствие.

Рассмотрим определенный интеграл , где – функция, непрерывная на отрезке . Проведём разбиение отрезка на равных отрезков:
. При этом, очевидно: (нижний предел интегрирования) и (верхний предел интегрирования). Точки также называют узлами .

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

Пример 1

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

а) Разбив отрезок интегрирования на 3 части.
б) Разбив отрезок интегрирования на 5 частей.

Решение:
а) Специально для чайников я привязал первый пункт к чертежу, который наглядно демонстрировал принцип метода. Если будет трудно, посматривайте на чертёж по ходу комментариев, вот его кусок:

По условию отрезок интегрирования нужно разделить на 3 части, то есть .
Вычислим длину каждого отрезка разбиения: . Параметр , напоминаю, также называют шагом .

Сколько будет точек (узлов разбиения)? Их будет на одну больше , чем количество отрезков:

Ну а общая формула трапеций сокращается до приятных размеров:

Для расчетов можно использовать обычный микрокалькулятор:

Обратите внимание, что, в соответствии с условием задачи, все вычисления следует округлять до 3-го знака после запятой .

Окончательно:

С геометрической точки зрения мы вычислили сумму площадей трёх трапеций (см. рис. выше) .

б) Разобьём отрезок интегрирования на 5 равных частей, то есть . Зачем это нужно? Чтобы Фобос-Грунт не падал в океан – увеличивая количество отрезков, мы увеличиваем точность вычислений.

Если , то формула трапеций принимает следующий вид:

Найдем шаг разбиения:
, то есть, длина каждого промежуточного отрезка равна 0,6.

При чистовом оформлении задачи все вычисления удобно оформлять расчетной таблицей:

В первой строке записываем «счётчик»

Как формируется вторая строка, думаю, всем видно – сначала записываем нижний предел интегрирования , остальные значения получаем, последовательно приплюсовывая шаг .

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

В результате:

Ну что же, уточнение, и серьёзное, действительно есть! Если для 3 отрезков разбиения приближённое значение составило, то для 5 отрезков . Таким образом, с большой долей уверенности можно утверждать, что, по крайне мере .

Пример 2

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

Решение: Почти та же задача, но немного в другой формулировке. Принципиальное отличие от Примера 1 состоит в том, что мы не знаем , НА СКОЛЬКО отрезков разбивать отрезок интегрирования, чтобы получить два верных знака после запятой. Иными словами, мы не знаем значение .

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

Сначала отрезок интегрирования разбивается на несколько больших отрезков, как правило, на 2-3-4-5. Разобьем отрезок интегрирования, например, на те же 5 частей. Формула уже знакома:

И шаг, естественно, тоже известен:

Но возникает еще один вопрос, до какого разряда округлять результаты ? В условии же ничего не сказано о том, сколько оставлять знаков после запятой. Общая рекомендация такова: к требуемой точности нужно прибавить 2-3 разряда . В данном случае необходимая точность 0,01. Согласно рекомендации, после запятой для верности оставим пять знаков (можно было и четыре):

В результате:
, обозначим приближение через .

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

Для формула трапеций приобретает следующий вид:

В бумажной версии запись можно спокойно перенести на следующую строчку.

Вычислим шаг разбиения:

Результаты расчётов сведём в таблицу:


При чистовом оформлении в тетрадь длинную таблицу выгодно превратить в двухэтажную.

В результате:

Теперь вычислим расхождение между приближениями:

Здесь используем знак модуля, поскольку нас интересует абсолютная разность , а не какой результат больше, а какой – меньше.

Что касается дальнейших действий, то лично мне на практике встречалось 2 пути решения:

1) Первый способ – это «сравнение в лоб». Поскольку полученная оценка погрешности больше , чем требуемая точность:, то необходимо ещё раз удвоить количество отрезков разбиения до и вычислить уже . С помощью экселевского калькулятора готовый результат можно получить в считанные секунды: . Теперь снова оцениваем погрешность: . Полученная оценка меньше , чем требуемая точность: , следовательно, вычисления закончены. Осталось округлить последний (наиболее точный) результат до двух знаков после запятой и дать ответ.

2) Другой, более эффективный способ основан на применении так называемого правила Рунге , согласно которому мы ошибаемся в оценке определённого интеграла на самом деле не более чем на . В нашей задаче: , таким образом, надобность в вычислении отпадает. Однако за скорость решения в данном случае пришлось расплатиться точностью: . Тем не менее, такой результат приемлем, поскольку наш «лимит на ошибку» как раз и составляет одну сотую.

Что выбрать? Ориентируйтесь на вашу методичку или предпочтения преподавателя.

Ответ: с точностью до 0,01 ( при использовании правила Рунге) .

Пример 3

Вычислить приближенно определенный интеграл по формуле трапеций с точностью до 0,001.

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

Как вычислить определенный интеграл по формуле Симпсона?

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

И снова, начнём с общей формулы
Рассмотрим определенный интеграл , где – функция, непрерывная на отрезке . Проведём разбиение отрезка на чётное количество равных отрезков. Чётное количество отрезков обозначают через .

На практике отрезков может быть:
два :
четыре :
восемь :
десять :
двадцать :
Другие варианты не припоминаю.

Внимание! Число понимается как ЕДИНОЕ ЧИСЛО. То есть, НЕЛЬЗЯ сокращать, например, на два, получая . Запись лишь обозначает , что количество отрезков чётно . И ни о каких сокращениях речи не идёт

Итак, наше разбиение имеет следующий вид:

Термины аналогичны терминам метода трапеций:
Точки называют узлами .

Формула Симпсона для приближенного вычисления определенного интеграла имеет следующий вид:
, где:
– длина каждого из маленьких отрезков или шаг ;
– значения подынтегральной функции в точках .

Детализируя это нагромождение, разберу формулу подробнее:
– сумма первого и последнего значения подынтегральной функции;
– сумма членов с чётными индексами умножается на 2;
– сумма членов с нечётными индексами умножается на 4.

Пример 4

Вычислить приближенно определенный интеграл по формуле Симпсона с точностью до 0,001. Разбиение начать с двух отрезков

Интеграл, кстати, опять неберущийся.

Решение: Сразу обращаю внимание на тип задания – необходимо вычислить определенный интеграл с определенной точностью . Что это значит, уже комментировалось в начале статьи, а также на конкретных примерах предыдущего параграфа. Как и для метода трапеций, существует формула, которая сразу позволит определить нужное количество отрезков (значение «эн») чтобы гарантированно достичь требуемой точности. Правда, придётся находить четвертую производную и решать экстремальную задачу. Кто понял, о чём я, и оценил объем работы, тот улыбнулся. Однако здесь не до смеха, находить четвертую производную от такой подынтегральной функции будет уже не мегаботан, а клинический психопат. Поэтому на практике практически всегда используется упрощенный метод оценки погрешности.

Начинаем решать. Если у нас два отрезка разбиения , то узлов будет на один больше : . И формула Симпсона принимает весьма компактный вид:

Вычислим шаг разбиения:

Заполним расчетную таблицу:

Еще раз комментирую, как заполняется таблица:

В верхнюю строку записываем «счётчик» индексов

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

В третью строку заносим значения подынтегральной функции. Например, если , то . Сколько оставлять знаков после запятой? Действительно, в условии опять об этом ничего не сказано. Принцип тот же, что и в методе трапеций, смотрим на требуемую точность: 0,001. И прибавляем дополнительно 2-3 разряда. То есть, округлять нужно до 5-6 знаков после запятой.

В результате:

Первичный результат получен. Теперь удваиваем количество отрезков до четырёх: . Формула Симпсона для данного разбиения принимает следующий вид:

Вычислим шаг разбиения:

Заполним расчетную таблицу:


Таким образом:

Найдём абсолютное значение разности между приближениями:

Правило Рунге для метода Симпсона очень вкусное. Если при использовании метода средних прямоугольников и метода трапеций нам даётся «поблажка» в одну треть, то сейчас – аж в одну пятнадцатую:
, и точность здесь уже не страдает:

Но для полноты картины я приведу и «простецкое» решение, где придётся сделать дополнительный шаг: так как больше требуемой точности: , то необходимо еще раз удвоить количество отрезков: .

Формула Симпсона растёт, как на дрожжах:

Вычислим шаг:

И снова заполним расчетную таблицу:

Таким образом:

Заметьте, что здесь вычисления желательно уже расписать более подробно, поскольку формула Симпсона достаточно громоздка, и если сразу бУхнуть:
, то выглядеть сиё бухло будет как халтура. А при более детальной записи у преподавателя сложится благостное впечатление, что вы добросовестно стирали клавиши микрокалькулятора в течение доброго часа. Детальные вычисления для «тяжелых» случаев присутствуют в моём калькуляторе.

Оцениваем погрешность:

Погрешность меньше требуемой точности: . Осталось взять наиболее точное приближение , округлить его до трёх знаков после запятой и записать:

Ответ : с точностью до 0,001

Пример 5

Вычислить приближенно определенный интеграл по формуле Симпсона с точностью до 0,0001. Разбиение начать с двух отрезков

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

В заключительной части урока рассмотрим еще пару распространенных примеров

Пример 6

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

Метод трапеций является одним из методов численного интегрирования. Он позволяет вычислять определенные интегралы с заранее заданной степенью точности.

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

Навигация по странице.

Суть метода трапеций.

Поставим перед собой следующую задачу: пусть нам требуется приближенно вычислить определенный интеграл , где подынтегральная функция y=f(x) непрерывна на отрезке .

Разобьем отрезок на n равных интервалов длины h точками . В этом случае шаг разбиения находим как и узлы определяем из равенства .

Рассмотрим подынтегральную функцию на элементарных отрезках .

Возможны четыре случая (на рисунке показаны простейшие из них, к которым все сводится при бесконечном увеличении n ):

На каждом отрезке заменим функцию y=f(x) отрезком прямой, проходящей через точки с координатами и . Изобразим их на рисунке синими линиями:

В качестве приближенного значения интеграла возьмем выражение , то есть, примем .

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

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

Таким образом, мы подошли к сути метода трапеций , которая состоит в представлении определенного интеграла в виде суммы интегралов вида на каждом элементарном отрезке и в последующей приближенной замене .

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

В силу пятого свойства определенного интеграла .

Если вместо интегралов подставить их приближенные значения, то получится :

Оценка абсолютной погрешности метода трапеций.

Абсолютная погрешность метода трапеций оценивается как
.

Графическая иллюстрация метода трапеций.

Приведем графическую иллюстрацию метода трапеций :

Примеры приближенного вычисления определенных интегралов методом трапеций.

Разберем на примерах применение метода трапеций при приближенном вычислении определенных интегралов.

В основном встречаются две разновидности заданий:

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

Следует заметить, что при заданном n промежуточные вычисления следует проводить с достаточной степенью точности, причем, чем больше n , тем выше должна быть точность вычислений.

Если требуется вычислить определенный интеграл с заданной точностью, к примеру, до 0.01 , то промежуточные вычисления рекомендуем проводить на два-три порядка точнее, то есть, до 0.0001 - 0.00001 . Если указанная точность достигается при больших n , то промежуточные вычисления следует проводить с еще более высокой точностью.

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

Итак, .

Пример.

Вычислить определенный интеграл методом трапеций для n = 10 .

Решение.

Формула метода трапеций имеет вид . То есть, для ее применения нам достаточно вычислить шаг h по формуле , определить узлы и вычислить соответствующие значения подынтегральной функции .

Вычислим шаг разбиения: .

Определяем узлы и вычисляем значения подынтегральной функции в них (будем брать четыре знака после запятой):

Результаты вычислений для удобства представляем в виде таблицы:

Подставляем их в формулу метода трапеций:

Полученное значение совпадает до сотых со значением, вычисленным по формуле Ньютона-Лейбница.

Пример.

Вычислите определенный интеграл методом трапеций с точностью до 0.01 .

Решение.

Что мы имеем из условия: a = 1; b = 2 ; .

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

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

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

Вернемся к нашему неравенству и подставим в него полученное значение:

Так как n – число натуральное (n - количество элементарных интервалов, на которые разбивается отрезок интегрирования), то можно брать n = 6, 7, 8, ... Возьмем n = 6 . Это позволит нам достичь требуемой точности метода трапеций при минимуме расчетов (хотя для нашего случая при n = 10 производить вычисления вручную удобнее).

Итак, n найдено, теперь действуем как в предыдущем примере.

Вычисляем шаг: .

Находим узлы сетки и значения подынтегральной функции в них:

Занесем в таблицу результаты расчетов:

Подставляем полученные результаты в формулу трапеций:

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

Следовательно, требуемая точность достигнута.

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

Приближенное значение определенного интеграла, полученное по методу трапеций для n узлов, будем обозначать .

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