Задача численного интегрирования функции


Скачать 92.05 Kb.
НазваниеЗадача численного интегрирования функции
Дата публикации17.09.2013
Размер92.05 Kb.
ТипЗадача
referatdb.ru > Математика > Задача

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


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

, (1) где .

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

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

Задача численного интегрирования функции f(x) заключается в вычислении значения определенного интеграла на основании ряда значений подынтегральной функции yk = f(xk).

Идея численного интегрирования вытекает из геометрического смысла определенного интеграла – значение определенного интеграла численно равно площади криволинейной трапеции, ограниченной графиком функции y=f(x), осью абсцисс и прямыми х=а, х=b. (рис 1)
y рис 1

f(x)

f(x
x

a b a b
Формально процедура численного интегрирования заключается в том, что отрезок [а, b] разбивается на n частичных отрезков, а затем подинтегральная функция заменяется легко интегрируемой функцией, интерполирующей значения подинтегральной функции в точках разбиения.

Наиболее часто приближенное значение интеграла ищут в виде линейной комбинации значений функции на отрезке :

. (2)

Приближенное равенство (2) называют квадратурной формулой, определяемой узлами и коэффициентами .

Выражение в правой части (2) называют квадратурной суммой,
а разность остаточным членом, или остатком, квадратурной формулы.
Рассмотрим общие вид квадратурных формул.

^

7.4. Квадратурные формулы Ньютона — Котеса


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

где — некоторые постоянные коэффициенты.

Выведем явные выражения для коэффициентов .

Для этого вспомним некоторые формулы из предыдущих лекций.

Введя переменную , полином Лагранжа можно переписать в виде:

(4)

Заменяя в (3) функцию полиномом Лагранжа в виде (4), учитывая сделанные выше выкладки, имеем: ,

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

, .

Так как , то обычно полагают ,

Где , , (5)

постоянные, называемые коэффициентами Ньютона — Котеса.

Квадратурная формула (3) при этом принимает вид:

, (6)
где , , .

Квадратурные формулы Ньютона — Котеса (6) дают на одном участке интегрирования различные представления для различного числа .

^

Формула трапеций и ее остаточный член.


Из формулы (5) при имеем: , ,

Отсюда . Мы получили известную формулу трапеций.

Остаточный член квадратурной формулы (3) равен . Предполагая, что , выведем простую формулу для остаточного члена. Будем рассматривать как функцию шага , тогда можно положить:

.

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



и ,

причем и .
Интегрируя по и используя теорему о среднем, последовательно выводим:

, где , и

, где .

Таким образом, окончательно имеем: , (7)

где .

Отсюда, в частности следует, если , то формула трапеций дает значение интеграла с избытком, если же — то с недостатком.

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



или . (8)
Геометрически формула (8) получается в результате замены кривой ломаной линией.
Если , то остаточный член квадратурной формулы (8) в силу (7) равен: , (9)

где .

Рассмотрим среднее арифметическое: . (10)

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

Из формул (9) и (10) имеем: , (11)

где .

^

Формула Симпсона и ее остаточный член.


Из формулы (5) при получаем: ,

, .

Так как , имеем: .

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

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




Остаточный член формулы Симпсона равен: .

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

.

Дифференцируя функцию по последовательно три раза, получим:



,



,



,

где .

Кроме того, очевидно, что , , . Последовательно интегрируя и используя теорему о среднем, находим:

,

где .

,

где .

,

где .
Таким образом, остаточный член формулы Симпсона равен:

, .

Очевидно, что формула Симпсона является точной для полиномов не только второй, но и третьей степени.

Пусть есть четное число и — значения функции для равноотстоящих точек с шагом .

Применяя формулу Симпсона к каждому удвоенному промежутку длины будем иметь:

.

Отсюда получаем общую формулу Симпсона

.

Если , то ошибка формулы Симпсона на каждом промежутке , , дается формулой: .

Суммируя все эти ошибки, получим остаточный член общей формулы Симпсона в виде: .

Так как непрерывна на , то найдется точка такая, что

.

Поэтому будем иметь: , (12)

где .




^

Об оценке точности квадратурных формул


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

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

Учитывая (12), получим: , где h и h шаги интегрирования. Тогда R=16R

Если I истинное значение интеграла (1), то I+R или I+R.

Из этого следует, что R=.

Эта формула удобна для практической оценки погрешности метода Симпсона.

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

Исследуем формулы с равноотстоящими узлами. К их числу относятся формулы трапеций, Симпсона, Ньютона — Котеса. Точность квадратурной формулы в основном характе-ризуется порядком остаточного члена R = O(h), где h= — шаг и mнатуральн. число.

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

Например, для f(x) = –8 + 45x  25x имеем:

.

При = 1 формула трапеций дает точное решение:

,

тогда как формула Симпсона при = 1 не обеспечивает даже знака интеграла:

.

Точность квадратурной формулы при фиксированном числе узлов существенно зависит от расположения этих узлов. При неудачном расположении узлов квадратурная формула может дать сильно искаженные результаты. При наличии значительного числа нулей подынтегральной функции или при большом числе ее экстремумов (т. е. когда много нулей производной ) точность квадратурных формул сильно понижается. Поэтому шаг h следует выбирать так, чтобы он был намного меньше расстояний между соседними нулями функции f(x) и ее производной . Для этого рекомендуется, например, основной отрезок интегрирования [a,b] разбить на частичные отрезки [α,β], внутри каждого из которых, функции f(x) и сохраняют постоянный знак (если это возможно!), и производить вычисление интеграла по частям, выбирая для каждого отрезка, вообще говоря, свой шаг.

Похожие рефераты:

Задача численного интегрирования функции
Если функция непрерывна на отрезке и известна ее первообразная, то определенный интеграл от этой функции в пределах от до может быть...
§ 9 Обзор методов интегрирования
Сводится к интегрированию произведения с помощью формулы кратного интегрирования по частям
Метод интегрирования по частям применяется, например, к интегралам вида
Сводится к интегрированию произведения с помощью формулы кратного интегрирования по частям
Тема: Гамма-функция и функция Бесселя
Определение неполной гамма-функции, построение её асимптотики с помощью метода интегрирования по частям
1 Неопределенный интеграл
Основной задачей дифференциального исчисления является нахождение производной или дифференциала функции. В интегральном исчислении...
Формулы интегрирования Основные формулы интегрирования Интегралы...

Неопределённый интеграл
Основной задачей дифференциального исчисления является задача нахождения производной данной функции. Многие задачи физики, механики,...
Работа 9 квантование энергии и волновые функции электрона в атоме водорода
Цель работы Нахождение уровней энергии и радиальных волновых функций атома водорода путём численного решения стационарного уравнения...
Задача Обеспечение развития экономики в пределах прогнозируемых параметров...
Запланированные мероприятия по реализации стратегического направления и цели государственного органа
Вопросы к экзамену по курсу «Экономико-математические методы»
Использование моделей оптимизации в промышленности (задача планирования производства, задача об оптимальном составе смеси, общая...

Вы можете разместить ссылку на наш сайт:
Школьные материалы


При копировании материала укажите ссылку © 2013
контакты
referatdb.ru
referatdb.ru
Рефераты ДатаБаза