Освоение численного интегрирования: правило Симпсона и не только

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

  1. Правило Симпсона.
    Правило Симпсона аппроксимирует интеграл функции путем деления интервала на подинтервалы и использования квадратичных полиномов для оценки площади под кривой. Формула правила Симпсона выглядит следующим образом:
def simpsons_rule(f, a, b, n):
    h = (b - a) / n
    x = a
    result = f(x)

    for i in range(1, n):
        x += h
        if i % 2 == 0:
            result += 2 * f(x)
        else:
            result += 4 * f(x)

    result += f(b)
    result *= h / 3

    return result
  1. Правило трапеций.
    Правило трапеций — еще один широко используемый метод численного интегрирования. Он аппроксимирует интеграл, заменяя кривую серией трапеций. Формула правила трапеций:
def trapezoidal_rule(f, a, b, n):
    h = (b - a) / n
    result = (f(a) + f(b)) / 2

    for i in range(1, n):
        x = a + i * h
        result += f(x)

    result *= h

    return result
  1. Правило средней точки.
    Правило средней точки — это простой метод численного интегрирования, который аппроксимирует интеграл с помощью прямоугольников одинаковой ширины. Формула правила средней точки:
def midpoint_rule(f, a, b, n):
    h = (b - a) / n
    result = 0.0

    for i in range(n):
        x = a + (i + 0.5) * h
        result += f(x)

    result *= h

    return result
  1. Другие методы.
    Помимо правила Симпсона, правила трапеций и правила средней точки, существует несколько других методов численного интегрирования, которые стоит изучить, например квадратура Гаусса, интегрирование Ромберга и методы адаптивного интегрирования, такие как адаптивное правило Симпсона. или адаптивное правило Гаусса-Кронрода.

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