Разобьем отрезок интегрирования [а , 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.

Формула

Формулой Симпсона называется интеграл от интерполяционного многочлена второй степени на отрезке :

где , и - значения функции в соответствующих точках (на концах отрезка и в его середине).

Погрешность

При условии, что у функции на отрезке существует четвёртая производная, погрешность , согласно найденной Джузеппе Пеано формуле равна:

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

Представление в виде метода Рунге-Кутты

Формулу Симпсона можно представить в виде таблицы метода Рунге-Кутты следующим образом:

Составная формула (формула Котеса)

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

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

Также формулу можно записать используя только известные значения функции, то есть значения в узлах:

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

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

.

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

.

Примечания

Литература

• Костомаров Д. П., Фаворский А. П. «Вводные лекции по численным методам»
• Петров И. Б., Лобанов А. И. Лекции по вычислительной математике

Wikimedia Foundation . 2010 .

• Western Union
• Патагонский попугай

Смотреть что такое "Формула Симпсона" в других словарях:

СИМПСОНА ФОРМУЛА - (формула парабол) формула для приближенного вычисления определенных интегралов (квадратурная формула), Названа по имени Т. Симпсона (1743) … Большой Энциклопедический словарь

СИМПСОНА ФОРМУЛА - (формула парабол), формула для приближённого вычисления определ. интегралов (квадратурная формула), имеющая вид где А = (b а)/2n, fk = f(a + kh), k = 0, 1, 2, ..., 2n. Названа по имени Т. Симпсона (1743) …

Симпсона формула - формула для приближённого вычисления определённых интегралов, имеющая вид: , где h = (b а)/2n; fi, = f (a + ih), i = 0, 1, 2,..., 2n. С. ф. называют иногда формулой парабол, т. к. вывод этой формулы основан на… … Большая советская энциклопедия

Симпсона формула - формула парабол, формула для приближённого вычисления определённых интегралов (квадратурная формула), имеющая вид, где h = (b–a)/2n, fk = f(а + kh), k = 0, 1, 2, ..., 2n. Названа по имени Т. Симпсона (1743). * * * СИМПСОНА ФОРМУЛА СИМПСОНА… … Энциклопедический словарь

Формула прямоугольников

Формула трапеций - Определённый интеграл как площадь фигуры Численное интегрирование (историческое название: квадратура) вычисление значения определённого интеграла (как правило, приближённое), основанное на том, что величина интеграла численно равна площади… … Википедия

СИМПСОНА ФОРМУЛА - частный случай Ньютона Котеса квадратурной формулы, в к рой берутся три узла: Пусть промежуток [а, b]разбит на пчастичных промежутков , i=0, 1, 2, ..., n 1, длины h=(b а)/п, при этом n считается четным числом, и для вычисления интеграла … Математическая энциклопедия

Симпсона формула - … Википедия

Метод Симпсона - Формула Симпсона относится к приёмам численного интегрирования. Получила название в честь британского математика Томаса Симпсона (1710 1761). Рассмотрим отрезок . Пусть известны значения вещественной функции f(x) в точках a, (a+b)/2, b.… … Википедия

КВАДРАТУРНАЯ ФОРМУЛА - формула, служа щая для приближённого вычисления определ. интегралов по значениям подынтегральной функции в конечном числе точек. Примеры К. ф. прямоугольников формула, трапеций формула, Симпсона формула … Естествознание. Энциклопедический словарь

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

Рис. 4

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

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

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

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

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

В методе Симпсона для вычисления определенного интеграла весь интервал интегрирования разбивается на подинтервалы равной длины h=(b-a)/n. Число отрезков разбиения является четным числом. Затем на каждой паре соседних подинтервалов подинтегральная функция f(x) заменяется многочленом Лагранжа второй степени (рисунок 5).

Рис. 5 Функция y=f(x) на отрезке заменяется многочленом 2-го порядка

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

Проинтегрируем на отрезке.:

Введем замену переменных:

Учитывая формулы замены,

Выполнив интегрирование, получим формулу Симпсона:

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

В формуле параболы значение функции f(x) в нечетных точках разбиения х 1 , х 3 , ..., х 2n-1 имеет коэффициент 4, в четных точках х 2 , х 4 , ..., х 2n-2 - коэффициент 2 и в двух граничных точках х 0 =а, х n =b - коэффициент 1.

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

Если функция f(x) имеет на непрерывную производную четвертого порядка, то абсолютная величина погрешности формулы Симпсона не больше чем

где М - наибольшее значение на отрезке . Так как n 4 растет быстрее, чем n 2 , то погрешность формулы Симпсона с ростом n уменьшается значительно быстрее, чем погрешность формулы трапеций.

Вычислим интеграл

Этот интеграл легко вычисляется:

Возьмем n равным 10, h=0.1, рассчитаем значения подынтегральной функции в точках разбиения, а также полуцелых точках.

По формуле средних прямоугольников получим I прям =0.785606 (погрешность равна 0.027%), по формуле трапеций I трап =0.784981 (погрешность около 0,054. При использовании метода правых и левых прямоугольников погрешность составляет более 3%.

Для сравнения точности приближенных формул вычислим еще раз интеграл

но теперь по формуле Симпсона при n=4. Разобьем отрезок на четыре равные части точками х 0 =0, х 1 =1/4, х 2 =1/2, х 3 =3/4, х 4 =1 и вычислим приближенно значения функции f(x)=1/(1+x) в этих точках: у 0 =1,0000, у 1 =0,8000, у 2 =0,6667, у 3 =0,5714, у 4 =0,5000.

По формуле Симпсона получаем

Оценим погрешность полученного результата. Для подынтегральной функции f(x)=1/(1+x) имеем: f (4) (x)=24/(1+x) 5 , откуда следует, что на отрезке . Следовательно, можно взять М=24, и погрешность результата не превосходит величины 24/(2880 4 4)=0.0004. Сравнивая приближенное значение с точным, заключаем, что абсолютная ошибка результата, полученного по формуле Симпсона, меньше 0,00011. Это находится в соответствии с данной выше оценкой погрешности и, кроме того, свидетельствует, что формула Симпсона значительно точнее формулы трапеций. Поэтому формулу Симпсона для приближенного вычисления определенных интегралов используют чаще, чем формулу трапеций.

Остаточный член квадратурной формулы Симпсона равен , где ξ∈(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 равных частей. Вычисления производить с округлением до четвертого десятичного знака.

Для построения формулы Симпсона предварительно рассмотрим такую задачу: вычислить площадь S криволинейной трапеции, ограниченной сверху графиком параболы y = Ax 2 + Bx + C, слева прямой х = - h, справа прямой x = h и снизу отрезком [-h; h]. Пусть парабола проходит через три точки (рис.8): D(-h; y 0) E(0; y 1) и F(h; y 2), причем х 2 - х 1 = х 1 - х 0 = h. Следовательно,

x 1 = x 0 + h = 0; x 2 = x 0 + 2h.

Тогда площадь S равна интегралу:

Выразим эту площадь через h, y 0 , y 1 и y 2 . Для этого вычислим коэффициенты параболы А, В, С. Из условия, что парабола проходит через точки D, E и F, имеем:

Решая эту систему, получаем: C = y 1 ; A =

Подставляя эти значения А и С в (3), получаем искомую площадь

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

Для этого отрезок интегрирования разобьем на 2n равных частей длиной

В точках деления (рис.4).а = х 0 , х 1 , х 2 , ...,х 2n-2 , x 2n-1 , x 2n = b,

Вчисляем значения подынтегральной функции f: y 0 , y 1 , y 2 , ...,y 2n-2 , y 2n-1 , y 2n , де y i = f(x i), x i = a + ih (i = 0, 1, 2,...,2n).

На отрезке подынтегральную функцию заменяем параболой, проходящей через точки (x 0 ; y 0), (x 1 ; y 1) и (x 2 ; y 2), и для вычисления приближенного значения интеграла от х 0 до х 2 воспользуемся формулой (4). Тогда (на рис. 4 заштрихованная площадь):

Аналогично находим:

................................................

Сложив полученные равенства, имеем:

Формула (5) называется обобщенной формулой Симпсона или формулой парабол , так как при ее выводе график подынтегральной функции на частичном отрезке длины 2h заменяется дугой параболы.

Задание на работу:

1. По указанию преподавателя или в соответствии с вариантом из Таблицы 4 заданий (см. Приложение) взять условия – подынтегральную функцию, пределы интегрирования.

2. Составить блок-схему программы и программу, которая должна:

Запросить точность вычисления определенного интеграла, нижний и верхний пределы интегрирования;

Вычислить заданный интеграл методами: для вариантов 1,4,7, 10… - правых, для вариантов 2,5,8,… - средних; для вариантов 2,5,8,… - левых прямоугольников. Вывести количество разбиений диапазона интегрирования, при котором достигнута заданная точность вычисления;

Вычислить заданный интеграл методом трапеций (для четных вариантов) и методом Симпсона (для нечетных вариантов).

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

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

Контрольные вопросы

1. Что такое определенный интеграл?

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

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

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

5. Как вычислить интеграл любым методом с заданной точностью?