Файл: Методические указания по выполнению курсовой работы по дисциплине (модулю) Введение в математический анализ.docx
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 30.11.2023
Просмотров: 141
Скачиваний: 1
ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.
2 x1
2. Приближенные вычисления определенных интегралов
-
Некоторые теоретические сведения
Если функция
y fx
непрерывна на отрезке a,b
и известна ее
первообразная F(x) , то определенный интеграл удается вычислить
непосредственно с помощью формулы Ньютона-Лейбница:
b
fxdx Fb Fa,
a
где
Fb и
Fa
-
значения первообразной
Fx
функции
fx
на концах
отрезка интегрирования.
Однако во многих случаях в реальных исследовательских задачах первообразная функции F(x) не может быть выражена через элементарные функции или является слишком сложной; вследствие этого вычисление определенного интеграла по формуле Ньютона-Лейбница может быть затруднительным или невыполнимым. Кроме того, на практике подынтегральная функция f(x) часто задается таблично и тогда само понятие первообразной теряет смысл. Поэтому большое значение имеют приближенные и в первую очередь численные методы вычисления определенных интегралов.
Назначение большинства приближенных методов вычисления
определенных интегралов состоит в замене подынтегральной функции
fx
аппроксимирующей функцией
x, для которой можно легко записать
первообразную в элементарных функциях, то есть
b b
I fxdx xdx R,
a a
где R– погрешность вычисления интеграла.
Чаще всего функцию
fx
заменяют интерполяционным полиномом
(под интерполяцией понимают приближенное вычисление неизвестных зна- чений функции по известным ее значениям в заданных точках), для построе-
ния которого используются значения функции в узлах
-
, i 0, n:
где
rx – остаточный член.
fx fxi rx,
n
i0
Подставляя полученное выражение в определенный интеграл вместо подынтегральной функции, получим общую формулу численного интегрирования
b n
I f xdx Aif xi R,
a i0
где
fxi
-
значения подынтегральной функции в узловых точках
xi ,
Ai –
весовые коэффициенты, а R– погрешность или остаточный член формулы.
С целью уменьшения погрешности, связанной с заменой
подынтегральной функции, отрезок интегрирования a,b
разбивают на n
отрезков и на каждом из полученных (частичных) отрезков xi1, xi, заменяют подынтегральную функцию аппроксимирующей функцией
i 1, n,
ix.
Тогда приближенное значение интеграла определяется суммой частичных
интегралов от функций ix, взятых в пределах от
xi1 до xi
для i 1, n:
b n xi
I fxdx ixdx. (3.1)
a i1 xi1
Методы численного интегрирования можно классифицировать в зависимости
от способа аппроксимации подынтегральной функции функций ix, i 1, n.
fx
с помощью
Методы Ньютона-Котеса основаны на полиномиальной
аппроксимации подынтегральной функции. Методы данного класса отличаются друг от друга степенью используемого полинома, от которой
зависит количество узлов, в которых необходимо вычислять функцию
fx.
В методах Ньютона-Котеса отрезок интегрирования разбивается, как правило, на отрезки равной длины, величина которых определяется как
h b a
n
и называется шагом интегрирования. Алгоритмы данных методов
просты и легко поддаются программной реализации.
В процессе численного интегрирования необходимо вычислить приближенное значение интеграла и оценить погрешность R . Погрешность, возникающая при численном интегрировании (также как и при численном дифференцировании), имеет два основных источника. Первым источником погрешности является замена подынтегральной функции аппроксимирующей функцией – погрешность аппроксимации. Как будет показано далее, погрешность аппроксимации уменьшается с увеличением количества nотрезков разбиения исходного отрезка интегрирования за счет более точной аппроксимации подынтегральной функции. Второй источник погрешности – неточности в вычислении подынтегральной функции в узловых точках и ошибки округления. Данная погрешность возрастает с ростом nи с
некоторого значения
n* начинает преобладать над погрешностью
аппроксимации. Это обстоятельство должно предостеречь от выбора чрезмерно большого числа n.
Способ получения формул для вычисления приближенного значения
интеграла в методах Ньютона-Котеса состоит в следующем. Разобьем
отрезок интегрирования a,b
на n частичных (как правило, равных по
длине) отрезков, точки разбиения
x0 , x1,..., xn,
a x0 ,
b xn
будем называть
узлами интегрирования, а расстояния между узлами
hi xi xi1 ,
i 1, n, –
шагамиинтегрирования. В частном случае шаг интегрирования может быть
постоянным ( h b a). На каждом из частичных отрезков интегрирования
n
xi1, xi,
i 1, n, будем аппроксимировать подынтегральную функцию
полиномом некоторой степени. В результате вычисление частичных интегралов на отрезках xi1, xi, i 1, n, по формуле (3.1) не составит труда.
-
Методы численного интегрирования
-
Методы прямоугольников
-
Рассмотрим простейшие методы из класса методов Ньютона-Котеса, в
которых подынтегральную функцию
fx
на отрезках интегрирования
xi1, xi,
i 1, n, заменяют полиномом нулевой степени (константой):
fx ix ci. Подставляя интегрирование, получаем
ix
в формулу (3.1) и выполняя
n
i
n
n
b
n
x
x
I f
a
xdx ci
i1
dxcix
x i1
i
xi1
ci
i1
xi xi1
cihi. (3.2)
i1
i1
Таким образом, в геометрической интерпретации приближенное значение интеграла определяется суммой площадей прямоугольников, одна из сторон которых соответствует отрезкам интегрирования длиной
hi xi xi1 , а другая – аппроксимирующим константам. Отсюда происходит
и название методов.
Далее будем использовать обозначения:
fi1
fxi1 ,
fi
fxi,