Методы парабол (Симпсона) и более высоких степеней (Ньютона - Котеса)

Материал из MachineLearning.

(Различия между версиями)
Перейти к: навигация, поиск
(Постановка математической задачи)
(Список литературы)
Строка 41: Строка 41:
== Список литературы ==
== Список литературы ==
 +
* ''А.А.Самарский, А.В.Гулин'' Численные методы М.: Наука, 1989.
 +
* ''А.А.Самарский'' Введение в численные методы М.: Наука, 1982.

Версия 12:41, 26 сентября 2008

Содержание

Введение

Постановка математической задачи

Задача численного интегрирования состоит в нахождении приближенного значения интеграла

J[f]=\int_a^b{f(x)dx}

где f(x) - заданная на отрезке [a,b] функция. На отрезке вводится сетка \omega=\{x_i:x_0=a<x_1<\ldots<x_i<\ldots<x_N=b\} и в качестве приближенного значения интеграла рассматривается число

J_N[f]=\sum_{i=0}^N {c_i f(x_i)}

где f(x_i) - значения функции f(x) в узлах x=x_i , где c_i - весовые множители, зависящие только от узлов, но не зависящие от выбора f(x). Данная формула называется квадратурной формулой.

Задача численного интегрирования при помощи квадратур состоит в отыскании таких узлов \{x_i\} и таких весов \{c_i\}, чтобы погрешность квадратурной формулы

D[f]=\sum_{i=0}^N{c_i f(x_i)} - \int_a^b{f(x)dx} = J_N[f] - J[f]

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

\int_\alpha^\beta{f(x)dx}

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

L[f]=\int_0^1{f(s)ds}

с помощью замены

x=\alpha+(\beta-\alpha)s,

f(x)=f(\alpha+(\beta-\alpha)s)

Изложение метода

Анализ метода

Числовой пример

Рекомендации программисту

Заключение

Список литературы

  • А.А.Самарский, А.В.Гулин Численные методы М.: Наука, 1989.
  • А.А.Самарский Введение в численные методы М.: Наука, 1982.
Личные инструменты