Численное интегрирование

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

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

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

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

  1. Метод прямоугольников
  2. Метод трапеций
  3. Метод парабол (Симпсона)

PLANETCALC, Интеграл численным методом по формулам Ньютона-Котеса

Интеграл численным методом по формулам Ньютона-Котеса

Знаков после запятой: 6
Формула
 
Значение определенного интеграла
 
Квадратурная функция
 
Погрешность метода
 
Интервал
 
Геометрический вид интеграла
Файл очень большой, при загрузке и создании может наблюдаться торможение браузера.
Источник формулы
 
Файл очень большой, при загрузке и создании может наблюдаться торможение браузера.

Численное интегрирование с использованием функций Ньютона Котеса

При использовании функций Ньютона-Котеса отрезок интегрирования разбивается на несколько равных отрезков точками x1,x2,x3..xn.
Подинтегральную функцию заменяют интерполяционным многочленом Лагранжа различной степени, интегрируя который, получают формулу численного интегрирования различного порядка точности.

В итоге, приближенное значение определенного интеграла вычисляется, как сумма значений подинтегральной функции в узлах, помноженных на некоторые константы Wi (веса):
I\approx \sum _{{i=1}}^{{n}}{W_i}f(x_{i})+R_n

  • Rn - остаток или погрешность.
  • n - общее количество точек.
  • Сумма в формуле - квадратурное правило (метод).

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

Границы отрезка интегрирования

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

Открытые правила, (правила, в которых граничные точки не включаются в расчет) удобно использовать в том случае, если подинтегральная функция не определена в некоторых точках.
Например, используя метод прямоугольников мы сможем вычислим приблизительное значение интеграла функции ln(x) на отрезке (0,1), несмотря на то, что ln(0) не существует.

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

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

Погрешность вычисления

В целом с увеличением количества узлов в правиле (при повышении степени интерполирующего полинома) возрастает точность вычисления интеграла. Однако для некоторых функций это может и не быть справедливо.
Впервые анализ этой особенности опубликовал Карл Рунге, немецкий математик, занимавшийся исследованием численных методов.
Он заметил, интерполирующий полином с равномерным разбиением отрезка для функции \frac{1}{1+25x^2} перестает сходиться в диапазоне значений 0.726.. ≤ |x| <1 при увеличении степени полинома.
В выражении для вычисления погрешности участвует интервал h, факториал от количества разбиений, которые при увеличении степени полинома уменьшают значение погрешности, но для некоторых функций значения производной, также участвующие в выражении погрешности, растут быстрее с увеличением ее порядка.

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

Принимая во внимание эти особенности, правила с полиномами степеней >10 применять не рекомендуется.

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

Для исследования работы с заданной функцией новых, основанных на формулах Ньютона-Котеса правил, можно воспользоваться базовым калькулятором, в котором веса задаются в явном виде:

PLANETCALC, Численное интегрирование с заданными весами Ньютона-Котеса

Численное интегрирование с заданными весами Ньютона-Котеса

Перечислите веса через запятую, в самом начале укажите общий множитель. Можно указывать коэффициенты в виде простой дроби, например, так: 3/4. Пример весов для метода Симпсона: 1/3,1,4,1.
Знаков после запятой: 6
Значение определенного интеграла
 
Формула
 
Квадратурная функция
 
Геометрический вид интеграла
Файл очень большой, при загрузке и создании может наблюдаться торможение браузера.
Файл очень большой, при загрузке и создании может наблюдаться торможение браузера.



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

Например, веса: 3/8,1,3,3,1 определяют Метод Симпсона 3/8

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


Литература:

  1. Н.С.Бахвалов Численные методы, 2012
  2. У.Г.Пирумов Численные методы, 2006
  3. Д.Каханер, К.Моулер, С.Нэш Численные методы и программное обеспечение, 1989
  4. Р.В. Хемминг Численные методы для научных работников и инженеров, 1972
  5. M. Abramovitz и I. Stegun Handbook of Mathematical Functions With Formulas, Graphs and Mathematical Tables, 1973
Ссылка скопирована в буфер обмена
PLANETCALC, Численное интегрирование

Комментарии