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



Скачать 33.03 Kb.
Дата16.01.2013
Размер33.03 Kb.
ТипДокументы


- -


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


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

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

,

(4.11)

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

Любую формулу численного интегрирования можно записать в виде суммы всех ординат c весовыми коэффициентами , 1, 2, ..., ,

,

(4.12)
где - это приближенное значение . Значения весовых коэффициентов определяются шагом дискретизации и методом интерполяции. Самые известные методы численного интегрирования  - метод трапеций и метод Симпсона - используют линейную и параболическую интерполяцию соответственно.

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

png" name="graphics1" align=bottom width=245 height=140 border=0>

Рис.4.2.

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

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

Парабола в методе Симпсона проводится через каждые три точки, т.е. значение   должно быть нечетным. Разобьем отрезок интегрирования на четное число отрезков , .

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

.В результате суммирования получается формула Симпсона:



(4.14)

В (4.14) все четные ординаты имеют множитель 4, а все нечетные, кроме первой и последней, множитель 2.

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

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

Видим, что полученное значение не зависит от шага . Зависимость суммарной погрешности от шага дискретизация показана на рис. 4.3. 



Рис 4.3. Зависимость суммарной погрешности численного интегрирования от шага дискретизации .


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

Контроль точности вычисления интеграла любым численным методом аналогичен численному дифференцированию.

В системе MathCAD для непрерывных функций применяется сложный метод вычисления интеграла — метод Ромберга, основанный на методе трапеций. Метод трапеций с шагами и используется в нем для контроля точности полученного значения интеграла, автоматического выбора шага и для уменьшения погрешности полученного значения . Отметим, что при этом величина начального шага для вычисления площади трапеции по (4.13) равна . Для обращения к стандартному оператору численного интегрирования следует ввести символ &.

Похожие:

4 Численное интегрирование iconКонтрольные вопросы к части 6 Часть 2 (дополнительная). Численное интегрирование методами Симпсона и Гаусса 7 Метод Симпсона 7
Обязательная. Численное интегрирование методами прямоугольников и трапеций 1
4 Численное интегрирование iconЧисленное интегрирование. Физические задачи, приводящие к интегрированию
Интегрирование функций является составной частью многих научных и технических задач. Поскольку аналитическое интегрирование не всегда...
4 Численное интегрирование iconОтчет по лабораторной работе №4 по теме «Численное интегрирование»
При этом значения в столбцах t и m означают: 1 –интегрирование методом средних прямоугольников, 2 – методом трапеций, 3 – методом...
4 Численное интегрирование iconЧисленное интегрирование
Функция y=f(x) заменяется интерполяционным многочленом P(x), который в точках xi = значению функции
4 Численное интегрирование iconВопросы по теме «Численное интегрирование»
Симпсона, непосредственно вычисляя интеграл разности подынтегральной функции и формулы
4 Численное интегрирование iconЛабораторная работа №6. Численное интегрирование
Рассмотрим вопрос о применении некоторых классов квадратурных формул к вычислению интегралов вида
4 Численное интегрирование iconПрограмма составлена кандидатом физ мат наук Барановым В. Н
Среднеквадратическая непрерывность случайного процесса. Интегрирование случайных процессов. Среднеквадратические интегралы с переменными...
4 Численное интегрирование iconЛабораторная работа №2 численное интегрирование на отрезке [ a, b
На отрезке [a, b] задана функция, будем считать, что функция достаточно гладкая. Требуется вычислить
4 Численное интегрирование iconЧисленное интегрирование
В данном интеграле пределы интегрирования а=-1, b Возьмем число разбиений n=100. Тогда номера точек i. 100. Шаг численного интегрирования...
4 Численное интегрирование iconМетоды вычислений Лабораторная работа №2 1
Напишите программу, реализующую численное интегрирование функции одной переменной методом (одним из трех по вариантам) левых прямоугольников,...
Разместите кнопку на своём сайте:
ru.convdocs.org


База данных защищена авторским правом ©ru.convdocs.org 2016
обратиться к администрации
ru.convdocs.org