Файл: Методические указания по выполнению курсовой работы по дисциплине (модулю) Введение в математический анализ.docx

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 30.11.2023

Просмотров: 141

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

2 x1

2. Приближенные вычисления определенных интегралов

    1. Некоторые теоретические сведения


Если функция

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

заменяют интерполяционным полиномом


(под интерполяцией понимают приближенное вычисление неизвестных зна- чений функции по известным ее значениям в заданных точках), для построе-

ния которого используются значения функции в узлах


  1. , i 0, n:


где
rx остаточный член.

fx fxi rx,


n
i0

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

b n

I f xdx Aif xi R,

a i0

где

fxi

  • значения подынтегральной функции в узловых точках

xi ,

Ai

весовые коэффициенты, а R погрешность или остаточный член формулы.

С целью уменьшения погрешности, связанной с заменой

подынтегральной функции, отрезок интегрирования a,b

разбивают на n

отрезков и на каждом из полученных (частичных) отрезков xi1, xi, заменяют подынтегральную функцию аппроксимирующей функцией

i 1, n,

ix.


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

интегралов от функций ix, взятых в пределах от

xi1 до xi

для i 1, n:

b n xi

I fxdx ixdx. (3.1)

a i1 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) не составит труда.
    1. Методы численного интегрирования

      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

i1

dxcix

x i1

i

xi1

ci

i1

xi xi1

cihi. (3.2)

i1

i1

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

hi xi xi1 , а другая аппроксимирующим константам. Отсюда происходит

и название методов.

Далее будем использовать обозначения:
fi1
fxi1 ,
fi
fxi,