Возникает задача о численном вычислении определенного интеграла, решаемая с помощью формул, носящих название квад­ратурных.

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

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

Величина

называется интервалом интегрирования или шагом. Отметим, что в практике -вычислений число я выбирают небольшим, обычно оно не больше 10-20.На частичном интервале

подынтегральную функцию заменяют интерполяционным много­членом

который на рассматриваемом интервале приближенно представ­ляет функцию f (х).

а) Удержим в интерполяционном многочлене только один первый член, тогда

Полученная квадратная формула

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

б) Удержим в интерполяционном многочлене два первых члена, тогда

(2)

Формула (2) называется формулой трапеций.

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

Полученная квадратурная формула называется формулой Симп­сона

(3)

Формулы (1), (2) и (3) имеют простой геометрический смысл. В формуле прямоугольников подынтегральная функция f(х) на интервале
заменяется отрезком прямой у = ук, параллельной оси абсцисс, а в формуле трапеций - отрезком прямой
и вычисляется соответственно площадь прямо­угольника и прямолинейной трапеции, которые затем сумми­руются. В формуле Симпсона функция f(х) на интервале
длиной 2h заменяется квадратным трехчленом - параболой
вычисляется площадь криволинейной параболической трапеции, затем площади суммируются.

ЗАКЛЮЧЕНИЕ

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

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

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

Таким образом, интегрирование всякой рациональной дроби сводится к интегрированию многочлена и нескольких простейших дробей.

Интеграл от любой рациональной функции может быть выражен через элементарные функции в конечном виде, а именно:

через логарифмы- в случаях простейших дробей 1 типа;

через рациональные функции- в случае простейших дробей 2 типа

через логарифмы и арктангенсы- в случае простейших дробей 3 типа

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

Формула Ньютона – Лейбница представляет собой общий подход к нахождению определенных интегралов.

Что касается приемов вычисления определенных интегралов, то они практически ничем не отличаются от всех тех приемов и методов.

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

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

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

Если необходимо получить наиболее точный результат, идеально подходит метод Симпсона .

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

Остаточный член квадратурной формулы Симпсона равен , где ξ∈(x 0 ,x 2) или

Назначение сервиса . Сервис предназначен для вычисления определенного интеграла по формуле Симпсона в онлайн режиме.

Инструкция . Введите подынтегральную функцию f(x) , нажмите Решить. Полученное решение сохраняется в файле Word . Также создается шаблон решения в Excel .

Правила ввода функции

Примеры правильного написания F(x):
1) 10 x e 2x ≡ 10*x*exp(2*x)
2) x e -x +cos(3x) ≡ x*exp(-x)+cos(3*x)
3) x 3 -x 2 +3 ≡ x^3-x^2+3

Вывод формулы Симпсона

Из формулы
при n = 2 получаем

Т.к. x 2 -x 0 = 2h, то имеем . (10)
Это формула Симпсона . Геометрически это означает, что кривую y=f(x) мы заменяем параболой y=L 2 (x), проходящей через три точки: M 0 (x 0 ,y 0), M 1 (x 1 ,y 1), M 2 (x 2 ,y 2).

Остаточный член формулы Симпсона равен

Предположим, что y∈C (4) . Получим явное выражение для R . Фиксируя среднюю точку x 1 и рассматривая R=R(h) как функцию h, будем иметь:
.
Отсюда дифференцируя последовательно три раза по h , получим






Окончательно имеем
,
где ξ 3 ∈(x 1 -h,x 1 +h). Кроме того, имеем: R(0) = 0, R"(0)=0. R""(0)=0. Теперь, последовательно интегрируя R"""(h), используя теорему о среднем, получим


Таким образом, остаточный член квадратурной формулы Симпсона равен
, где ξ∈(x 0 ,x 2). (11)
Следовательно, формула Симпсона является точной для полиномов не только второй, но и третьей степени.
Получим теперь формулу Симпсона для произвольного интервала [a ,b ]. Пусть n = 2m есть четное число узлов сетки {x i }, x i =a+i·h, i=0,...,n, и y i =f(x i). Применяя формулу Симпсона (10) к каждому удвоенному промежутку , ,..., длины 2h , будем иметь


Отсюда получаем общую формулу Симпсона
.(12)
Ошибка для каждого удвоенного промежутка (k=1,...,m) дается формулой (11).

Т.к. число удвоенных промежутков равно m , то

С учетом непрерывности y IV на [a ,b ], можно найти точку ε, такую, что .
Поэтому будем иметь
. (13)
Если задана предельно допустимая погрешность ε, то, обозначив , получим для определения шага h
.
На практике вычисление R по формуле (13) бывает затруднительным. В этом случае можно поступить следующим образом. Вычисляем интеграл I(h)=I 1 с шагом h , I(2h)=I 2 с шагом 2h и т.д. и вычисляем погрешность Δ:
Δ = |I k -I k-1 | ≤ ε. (14)
Если неравенство (14) выполняется (ε - заданная погрешность), то за оценку интеграла берут I k = I(k·h).
Замечание. Если сетка неравномерная, то формула Симпсона приобретает следующий вид (получить самостоятельно)
.
Пусть число узлов n = 2m (четное). Тогда

где h i =x i -x i-1 .

Пример №1 . С помощью формулы Симпсона вычислить интеграл , приняв n = 10.
Решение: Имеем 2m = 10. Отсюда . Результаты вычислений даны в таблице:

i	x i	y 2i-1	y 2i
0	0		y 0 = 1.00000
1	0.1	0.90909
2	0.2		0.83333
3	0.3	0.76923
4	0.4		0.71429
5	0.5	0.66667
6	0.6		0.62500
7	0.7	0.58824
8	0.8		0.55556
9	0.9	0.52632
10	1.0		y n =0.50000
∑		σ 1	σ 2

По формуле (12) получим .
Рассчитаем погрешность R=R 2 . Т.к. , то .
Отсюда max|y IV |=24 при 0≤x≤1 и, следовательно . Таким образом, I = 0.69315 ± 0.00001.

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

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

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

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

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

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

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

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

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

Пример 4

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

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

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

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

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

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

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

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

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

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

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

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

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

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


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

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

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

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

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

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

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

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

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

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

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

Пример 5

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

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

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

Пример 6

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

Этот интеграл берётся, правда, новичку взломать его не так-то просто, соответствующий метод решения рассмотрен в примере 5 урока Сложные интегралы . В задачах на приближенное вычисление интеграл не обязан быть непременно неберущимся! Любознательные студенты могут вычислить его точно и оценить погрешность относительно приближенного значения.

Решение: Обратите внимание на формулировку задания: «Точность вычислений 0,001». Смысловой нюанс данной формулировки предполагает, что результаты нужно только округлить до третьего знака после запятой, а не достигнуть такой точности. Таким образом, когда вам предлагается для решения задача на метод трапеций, метод Симпсона, всегдавнимательно вникайте в условие ! Спешка, как известно, нужна при охоте на блох.

Используем формулу Симпсона:

При десяти отрезках разбиения шаг составляет

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

Таблицу рациональнее сделать двухэтажной, чтобы не пришлось «мельчить» и всё разборчиво вместилось на тетрадный лист.

Вычисления, не ленимся, расписываем подробнее:

Ответ:

И еще раз подчеркну, что о точности здесь речи не идет. На самом деле, ответ может быть не , а, условно говоря, . В этой связи в ответе не нужно машинально приписывать «дежурную» концовку: «с точностью до 0,001»

Пример 7

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

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

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

А сейчас настала пора раскрыть страшную тайну интегрального исчисления. Я создал уже больше десятка уроков по интегралам, и это, так скажем, теория и классика темы. На практике же, в частности, при инженерных расчетах – приблизить объекты реального мира стандартными математическими функциями практически невозможно. Невозможно идеально точно рассчитать, площадь, объем, плотность, к примеру, асфальтового покрытия.Погрешность , пусть с десятого, пусть с сотого знака после запятой – но она всё равно будет . Именно поэтому по приближенным методам вычисления написаны сотни увесистых кирпичей и создано серьёзное программное обеспечение для приближенных вычислений. Классическая же теория интегрального исчисления в действительности применяется заметно реже. Но, кстати, без неё – тоже никуда!

Данный урок не рекорден по объему, но на его создание у меня ушло необычно много времени. Я правил материал и переделывал структуру статьи несколько раз, поскольку постоянно прорисовывались новые нюансы и тонкости. Надеюсь, труды были не напрасны, и получилось вполне логично и доступно.

Всего вам доброго!

Решения и ответы:

Пример 3: Решение: Разбиваем отрезок интегрирования на 4 части:
Тогда формула трапеций принимает следующий вид:

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

В этом методе предлагается подынтегральную функцию на частичном отрезке аппроксимировать параболой, проходящей через точки
(x j , f (x j )), где j = i -1; i -0.5; i , то есть подынтегральную функцию аппроксимируем интерполяционным многочленом Лагранжа второй степени:

Проведя интегрирование, получим:

Это и есть формула Симпсона или формула парабол. На отрезке
[a, b ] формула Симпсона примет вид

Графическое представление метода Симпсона показано на рис. 2.4.

Рис. 10.4. Метод Симпсона

Избавимся в выражении (2.16) от дробных индексов, переобозначив переменные:

Тогда формула Симпсона примет вид

Погрешность формулы (2.18) оценивается следующим выражением:

где h·n = b - a , . Таким образом, погрешность формулы Симпсона пропорциональна O (h 4 ).

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

10.5. Вычисление определенных интегралов методами
Монте–Карло

Рассматриваемые ранее методы называются детерминированными , то есть лишенными элемента случайности.

Методы Монте–Карло (ММК) – это численные методы решения математических задач с помощью моделирования случайных величин. ММК позволяют успешно решать математические задачи, обусловленные вероятностными процессами. Более того, при решении задач, не связанных с какими-либо вероятностями, можно искусственно придумать вероятностную модель (и даже не одну), позволяющую решать эти задачи. Рассмотрим вычисление определенного интеграла

При вычислении этого интеграла по формуле прямоугольников интервал [a, b ] разбиваем на N одинаковых интервалов, в серединах которых вычислялись значения подынтегральной функции. Вычисляя значения функции в случайных узлах, можно получить более точный результат:

Здесь γ i - случайное число, равномерно распределенное на интервале
. Погрешность вычисления интеграла ММК ~ , что значительно больше, чем у ранее изученных детерминированных методов.

На рис. 2.5 представлена графическая реализация метода Монте-Карло вычисления однократного интеграла со случайными узлами (2.21) и (2.22).

(2.23)

Рис. 10.6. Интегрирование методом Монте-Карло (2-й случай)

Как видно на рис. 2.6, интегральная кривая лежит в единичном квадрате, и если мы сумеем получать пары случайных чисел, равномерно распределенных на интервале , то полученные значения (γ 1, γ 2) можно интерпретировать как координаты точки в единичном квадрате. Тогда, если этих пар чисел получено достаточно много, можно приблизительно считать, что
. Здесь S – число пар точек, попавших под кривую, а N – общее число пар чисел.

Пример 2.1. Вычислить следующий интеграл:

Поставленная задача была решена различными методами. Полученные результаты сведены в табл. 2.1.

Таблица 2.1

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

11 ПРИБЛИЖЕННОЕ РЕШЕНИЕ НЕЛИНЕЙНЫХ
И ТРАНСЦЕНДЕНТНЫХ УРАВНЕНИЙ

Разобьем отрезок интегрирования [а , b ] на четное число n равных частей с шагом h . На каждом отрезке [х 0, х 2], [х 2, х 4],..., [x i-1, x i+1],..., [x n-2, x n] подынтегральную функцию f (х ) заменим интерполяционным многочленом второй степени:

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

Сумму элементарных площадей и (рис. 3.3) можно вычислить с помощью определенного интеграла. Учитывая равенства получаем

-

Рис. 3.3. Иллюстрация к методу Симпсона

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

Данное выражение для S принимается в качестве значения определенного интеграла:

(3.35)

Полученное соотношение называется формулой Симпсона или формулой парабол .

Эту формулу можно получить и другими способами, например двукратным применением метода трапеций при разбиениях отрезка [а , b ] на части с шагами h и 2h или комбинированием формул прямоугольников и трапеций (см. разд. 3.2.6).

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

(3.36)

Легко видеть, что формула (3.36) совпадет с (3.35), если формулу (3.35) применить для числа отрезков разбиения 2n и шага h /2.

Пример . Вычислить по методу Симпсона интеграл

Значения функции при n = 10, h = 0.1 приведены в табл. 3.3. Применяя формулу (3.35), находим

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

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

Рис. 3.4. Алгоритм метода Симпсона

Первоначально отрезок разбивается на две части с шагом h =(b - a)/2. Вычисляется значение интеграла I 1. Потом число шагов удваивается, вычисляется значение I 2 с шагом h /2. Условие окончание счета принимается в виде . Если это условие не выполнено, происходит новое деление шага пополам и т.д.

Отметим, что представленный на рис. 3.4 алгоритм не является оптимальным: при вычислении каждого приближения I 2 не используются значения функции f (x ), уже найденные на предыдущем этапе. Более экономичные алгоритмы будут рассмотрены в разд. 3.2.7.