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

Суть метода, формула, оценка погрешности.

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

Разобьем отрезок на n элементарных

отрезков [;], i = 1., n длины 2*h = (b-a)/ n точками

a = < < < < = b. Пусть точки, i = 1., n являются серединами отрезков [;], i = 1., n соответственно. В этом случае все «узлы» определяются из равенства = a + i*h, i = 0,1., 2*n.

На каждом интервале [;], i = 1,2., n подынтегральная функция

приближается квадратичной параболой y = a* + b*x + c, проходящей через точки (; f ()), (; f ()), (; f ()). Отсюда и название метода - метод парабол.

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

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

Для получения формулы метода парабол (Симпсона) нам осталось вычислить

Покажем, что через точки (; f ()), (; f ()), (; f ()) проходит только одна квадратичная парабола y = a* + b*x + c. Другими словами, докажем, что коэффициенты, определяются единственным образом.

Так как (; f ()), (; f ()), (; f ()) - точки параболы, то справедливо каждое из уравнений системы

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

Перейдем к нахождению интеграла.

Очевидно:

f () = f(0) = + + =

f () = f(h) = + +

f () = f (2*h) = + +

Используем эти равенства, чтобы осуществить последний переход в следующей цепочке равенств:

= = (++) = h/3*(f ()+4*f ()+f ())

Таким образом, можно получить формулу метода парабол:

Пример метода Симпсона.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Формула Симпсона становится больше:

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

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

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

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

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

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

Формула

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

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

Погрешность

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

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

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

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

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

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

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

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

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

Примечания

Литература

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

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.… … Википедия

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

В этом методе предлагается подынтегральную функцию на частичном отрезке аппроксимировать параболой, проходящей через точки
(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).

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

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

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

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

Таблица 2.1

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

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

Кафедра «Высшей математики»

Выполнил: Матвеев Ф.И.

Проверила: Бурлова Л.В.

Улан-Удэ.2002

1.Численные методы интегрирования

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

3.Геометрическая иллюстрация

4.Выбор шага интегрирования

5.Примеры

1. Численные методы интегрирования

Задача численного интегрирования заключается в вычислении интеграла

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

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

Численные методы условно можно сгруппировать по способу аппроксимации подынтегральной функции.

Методы Ньютона-Котеса основаны на аппроксимации функции

Методы сплайн-интегрирования базируются на аппроксимации функции

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

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

погрешность округления

Независимо от выбранного метода в процессе численного интегрирования необходимо вычислить приближенное значение интеграла и оценить погрешность. Погрешность уменьшается при увеличении n-количества

разбиений отрезка

за счет суммирования значений интегралов, вычисленных на частичных отрезках.

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

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

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

Проинтегрируем

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

Полученное для интеграла

Оценим теперь погрешность интегрирования по формуле Симпсона. Будем считать, что у

К каждому из этих двух интегралов уже можно применить теорему о среднем, поскольку

(мы воспользовались теоремой о среднем, поскольку

Дифференцируя

Из обеих оценок для

Если отрезок

Запишем формулу Симпсона в общем виде.

Кафедра «Высшей математики»

Выполнил: Матвеев Ф.И.

Проверила: Бурлова Л.В.

Улан-Удэ.2002

1.Численные методы интегрирования

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

3.Геометрическая иллюстрация

4.Выбор шага интегрирования

5.Примеры

1. Численные методы интегрирования

Задача численного интегрирования заключается в вычислении интеграла

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

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

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

Численные методы условно можно сгруппировать по способу аппроксимации подынтегральной функции.

Методы Ньютона-Котеса основаны на аппроксимации функции полиномом степени . Алгоритм этого класса отличается только степенью полинома. Как правило, узлы аппроксимирующего полинома – равноотносящие.

Методы сплайн-интегрирования базируются на аппроксимации функции сплайном-кусочным полиномом.

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

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

суммарная погрешность

погрешность усечения

погрешность округления

Независимо от выбранного метода в процессе численного интегрирования необходимо вычислить приближенное значение интеграла и оценить погрешность. Погрешность уменьшается при увеличении n-количества

разбиений отрезка . Однако при этом возрастает погрешность округления

за счет суммирования значений интегралов, вычисленных на частичных отрезках.

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

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

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

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

Проинтегрируем :

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

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

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

К каждому из этих двух интегралов уже можно применить теорему о среднем, поскольку непрерывна на и функция неотрицательна на первом интервале интегрирования и неположительна на втором (то есть не меняет знака на каждом из этих интервалов). Поэтому:

(мы воспользовались теоремой о среднем, поскольку - непрерывная функция; ).

Дифференцируя дважды и применяя затем теорему о среднем, получим для другое выражение:

, где

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

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

Запишем формулу Симпсона в общем виде:

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

, (3)

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

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

3. Геометрическая иллюстрация

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

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

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

(4)

Это формула Симпсона «трех восьмых».

Для произвольного отрезка интегрирования формула (4) может быть «продолжена»; при этом число частичных отрезков должно быть кратно трем ( точек).

, m=2,3,... (5)

Целая часть

Можно получить формулы Ньютона-Котеса старших порядков:

(6)

Количество отрезков разбиения;

Степень используемого полинома;

Производная -го порядка в точке ;

Шаг разбиения.

В таблице 1 выписаны коэффициенты . Каждая строка соответствует одному набору промежутков узлами для построения многочлена k-ой степени. Чтобы воспользоваться этой схемой для большего количества наборов (например, при k=2 и n=6), нужно «продолжить» коэффициенты, а затем сложить их.

Таблица 1:

Алгоритм оценки погрешности формул трапеции и Симпсона можно записать в виде: (7),

где - коэффициент, зависящий от метода интегрирования и свойств подынтегральной функции;

h - шаг интегрирования;

p - порядок метода.

Правило Рунге применяют для вычисления погрешности путем двойного просчета интеграла с шагами h и kh.

(8) - апостериорная оценка. Тогда Iуточн.= +Ro (9), уточненное значение интеграла .

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

из системы трех уравнений:

с неизвестными I,А и p получаем:

Из (10) следует (11)

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

Аналогично, для интегралов, вычисленных по формуле с шагами и , справедливы соотношения:

,

(13)

4. Выбор шага интегрирования

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

Если ê ê, то ê ê.

По заданной точности e метода интегрирования из последнего неравенства определяем подходящий шаг.

, .

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

Разберем один из таких приемов. Пусть

,

где - приближенное значение интеграла с шагом . Уменьшим шаг в два раза, разбив отрезок на две равные части и ().

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

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

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

При выводе правила Рунге вы существенно пользовались предположением, что . Если имеется только таблица значений , то проверку «на постоянство» можно сделать непосредственно по таблице Дальнейшее развитие приведенных алгоритмов позволяет перейти к адаптивным алгоритмам, в которых за счет выбора различного шага интегрирования в разных частях отрезка интегрирования в зависимости от свойств уменьшается количество вычислений подынтегральной функции.

Другая схема уточнения значений интеграла - процесс Эйтнена. Производится вычисление интеграла с шагами, причем . Вычисление значений . Тогда (14).

За меру точности метода Симпсона принимают величину:

5. Примеры

Пример 1. Вычислить интеграл по формуле Симпсона, если задана таблицей. Оценить погрешность.

Таблица 3.

Решение: Вычислим по формуле (1) при и интеграл .

По правилу Рунге получаем Принимаем .

Пример 2. Вычислить интеграл .

Решение: Имеем . Отсюда h==0.1. Результаты вычислений приведены в таблице 4.

Таблица 4.

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