Образовательный портал - Varnavinschool

Метод трапеций. Оценка точности вычисления «неберущихся» интегралов Примеры приближенного вычисления определенных интегралов методом парабол

При вычислении определенного интеграла не всегда получаем точное решение. Не всегда удается представление в виде элементарной функции. Формула Ньютона-Лейбница не подходит для вычисления, поэтому необходимо использовать методы численного интегрирования. Такой метод позволяет получать данные с высокой точностью. Метод Симпсона является таковым.

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

Метод парабол – суть, формула, оценка, погрешности, иллюстрации

Задана функция вида y = f (x) , имеющая непрерывность на интервале [ a ; b ] , необходимо произвести вычисление определенного интеграла ∫ a b f (x) d x

Необходимо разбить отрезок [ a ; b ] на n отрезков вида x 2 i - 2 ; x 2 i , i = 1 , 2 , . . . , n с длиной 2 h = b - a n и точками a = x 0 < x 2 < x 4 < . . . < x 2 π - 2 < x 2 π = b . Тогда точки x 2 i - 1 , i = 1 , 2 , . . . , n считаются серединами отрезков x 2 i - 2 ; x 2 i , i = 1 , 2 , . . . , n . Данный случай показывает, что определение узлов производится через x i = a + i · h , i = 0 , 1 , . . . , 2 n .

Каждый интервал x 2 i - 2 ; x 2 i , i = 1 , 2 , . . . , n подынтегральной функции приближен при помощи параболы, заданной y = a i x 2 + b i x + c i , проходящей через точки с координатами x 2 i - 2 ; f (x 2 i - 2) , x 2 i - 1 ; x 2 i - 1 , x 2 i ; f (x 2 i) . Поэтому метод и имеет такое название.

Данные действия выполняются для того, чтобы интеграл ∫ x 2 i - 2 x 2 i a i x 2 + b i x + c i d x взять в качестве приближенного значения ∫ x 2 i - 2 x 2 i f (x) d x . Можем вычислить при помощи формулы Ньютона-Лейбница. Это и есть суть метода парабол.Рассмотрим рисунок, приведенный ниже.

Графическая иллюстрация метода парабол (Симпсона)

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

Исходя из пятого свойства определенного интеграла получаем ∫ a b f (x) d x = ∑ i = 1 n ∫ x 2 i - 2 x 2 i f (x) d x ≈ ∑ i = 1 n ∫ x 2 i - 2 x 2 i (a i x 2 + b i x + c i) d x

Для того чтобы получить формулу методом парабол, необходимо произвести вычисление:

∫ x 2 i - 2 x 2 i (a i x 2 + b i x + c i) d x

Пусть x 2 i - 2 = 0 . Рассмотрим рисунок, приведенный ниже.

Изобразим, что через точки с координатами x 2 i - 2 ; f (x 2 i - 2) , x 2 i - 1 ; x 2 i - 1 , x 2 i ; f (x 2 i) может проходить одна квадратичная парабола вида y = a i x 2 + b i x + c i . Иначе говоря, необходимо доказать, что коэффициенты могут определяться только единственным способом.

Имеем, что x 2 i - 2 ; f (x 2 i - 2) , x 2 i - 1 ; x 2 i - 1 , x 2 i ; f (x 2 i) являются точками параболы, тогда каждое из представленных уравнений является справедливым. Получаем, что

a i (x 2 i - 2) 2 + b i · x 2 i - 2 + c i = f (x 2 i - 2) a i (x 2 i - 1) 2 + b i · x 2 i - 1 + c i = f (x 2 i - 1) a i (x 2 i) 2 + b i · x 2 i + c i = f (x 2 i)

Полученная система разрешается относительно a i , b i , c i , где необходимо искать определитель матрицы по Вандермонду. Получаем, что

(x 2 i - 2) 2 x 2 i - 2 1 x 2 i - 1) 2 x 2 i - 1 1 (x 2 i) 2 x 2 i 1 , причем он считается отличным от нуля и не совпадает с точками x 2 i - 2 , x 2 i - 1 , x 2 i . Это признак того, что уравнение имеет только одно решение, тогда и выбранные коэффициенты a i ; b i ; c i могут определяться только единственным образом, тогда через точки x 2 i - 2 ; f (x 2 i - 2) , x 2 i - 1 ; x 2 i - 1 , x 2 i ; f (x 2 i) может проходить только одна парабола.

Можно переходить к нахождению интеграла ∫ x 2 i - 2 x 2 i (a i x 2 + b i x + c i) d x .

Видно, что

f (x 2 i - 2) = f (0) = a i · 0 2 + b i · 0 + c i = c i f (x 2 i - 1) = f (h) = a i · h 2 + b i · h + c i f (x 2 i) = f (0) = 4 a i · h 2 + 2 b i · h + c i

Для осуществления последнего перехода необходимо использовать неравенство вида

∫ x 2 i - 2 x 2 i (a i x 2 + b i x + c i) d x = ∫ 0 2 h (a i x 2 + b i x + c i) d x = = a i x 3 3 + b i x 2 2 + c i x 0 2 h = 8 a i h 3 3 + 2 b i h 2 + 2 c i h = = h 3 8 a i h 2 + 6 b i h + 6 c i = h 3 f x 2 i - 2 + 4 f 2 2 i - 1 + f x 2 i

Значит, получаем формулу, используя метод парабол:

∫ a b f (x) d x ≈ ∑ i = 1 n ∫ x 2 i - 2 x 2 i a i x 2 + b i x + c i d x = = ∑ i = 1 n h 3 (f (x 2 i - 2) + 4 f (x 2 i - 1) + f (x 2 i)) = = h 3 f (x 0) + 4 f (x 1) + f (x 2) + f (x 2) + 4 f (x 3) + f (x 4) + . . . + + f (x 2 n - 2) + 4 f (x 2 n - 1) + f (x 2 n) = = h 3 f (x 0) + 4 ∑ i = 1 n f (x 2 i - 1) + 2 ∑ i = 1 n - 1 f (x 2 i) + f (x 2 n)

Определение 1

Формула метода Симпсона имеет вид ∫ a b f (x) d x ≈ h 3 f (x 0) + 4 ∑ i = 1 n f (x 2 i - 1) + 2 ∑ i = 1 n - 1 f (x 2 i) + f (x 2 n) .

Формула оценки абсолютной погрешности имеет вид δ n ≤ m a x [ a ; b ] f (4) (x) · (b - a) 5 2880 n 4 .

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

Метод Симпсона предполагает приближенное вычисление определенных интегралов. Чаще всего имеются два типа задач, для которых применим данный метод:

  • при приближенном вычислении определенного интеграла;
  • при нахождении приближенного значения с точностью δ n .

На точность вычисления влияет значение n , чем выше n , тем точнее промежуточные значения.

Пример 1

Вычислить определенный интеграл ∫ 0 5 x d x x 4 + 4 при помощи метода Симпсона, разбивая отрезок интегрирования на 5 частей.

Решение

По условию известно, что a = 0 ; b = 5 ; n = 5 , f (x) = x x 4 + 4 .

Тогда запишем формулу Симпсона в виде

∫ a b f (x) d x ≈ h 3 f (x 0) + 4 ∑ i = 1 n f (x 2 i - 1) + 2 ∑ i = 1 n - 1 f (x 2 i) + f (x 2 n)

Чтобы применить ее в полной мере, необходимо рассчитать шаг по формуле h = b - a 2 n , определить точки x i = a + i · h , i = 0 , 1 , . . . , 2 n и найти значения подынтегральной функции f (x i) , i = 0 , 1 , . . . , 2 n .

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

h = b - a 2 n = 5 - 0 2 · 5 = 0 . 5

Найдем значение функции в точках

i = 0: x i = x 0 = a + i · h = 0 + 0 · 0 . 5 = 0 ⇒ f (x 0) = f (0) = 0 0 4 + 4 = 0 i = 1: x i = x 1 = a + i · h = 0 + 1 · 0 . 5 = 0 . 5 ⇒ f (x 1) = f (0 . 5) = 0 . 5 0 . 5 4 + 4 ≈ 0 . 12308 . . . i = 10: x i = x 10 = a + i · h = 0 + 10 · 0 . 5 = 5 ⇒ f (x 10) = f (5) = 5 5 4 + 4 ≈ 0 . 00795

Наглядность и удобство оформляется в таблице, приведенной ниже

i 0 1 2 3 4 5
x i 0 0 . 5 1 1 . 5 2 2 . 5
f x i 0 0 . 12308 0 . 2 0 . 16552 0 . 1 0 . 05806
i 6 7 8 9 10
x i 3 3 . 5 4 4 . 5 5
f x i 0 . 03529 0 . 02272 0 . 01538 0 . 01087 0 . 00795

Необходимо подставить результаты в формулу метода парабол:

∫ 0 5 x d x x 4 + 4 ≈ h 3 f (x 0) + 4 ∑ i = 1 n f (x 2 i - 1) + 2 ∑ i = 1 n - 1 f (x 2 i) + f (x 2 n) = = 0 . 5 3 0 + 4 · 0 . 12308 + 0 . 16552 + 0 . 05806 + + 0 . 02272 + 0 . 01087 + 2 · 0 . 2 + 0 . 1 + + 0 . 03529 + 0 . 01538 + 0 . 00795 ≈ ≈ 0 . 37171

Для вычисления мы выбрали определенный интеграл, который можно вычислить по Ньютону-Лейбницу. Получим:

∫ 0 5 x d x x 4 + 4 = 1 2 ∫ 0 5 d (x 2) x 2 2 + 4 = 1 4 a r c t g x 2 2 0 5 = 1 4 a r c t g 25 2 ≈ 0 . 37274

Ответ: Результаты совпадают до сотых.

Пример 2

Вычислить неопределенный интеграл ∫ 0 π sin 3 x 2 + 1 2 d x при помощи метода Симпсона с точностью до 0 , 001 .

Решение

По условию имеем, что а = 0 , b = π , f (x) = sin 3 x 2 + 1 2 , δ n ≤ 0 . 001 . Необходимо определить значение n . Для этого используется формула оценки абсолютной погрешности метода Симпсона вида δ n ≤ m a x [ a ; b ] f (4) (x) · (b - a) 5 2880 n 4 ≤ 0 . 001

Когда найдем значение n , то неравенство m a x [ a ; b ] f (4) (x) · (b - a) 5 2880 n 4 ≤ 0 . 001 будет выполняться. Тогда, применив метод парабол, погрешность при вычислении не превысит 0 . 001 . Последнее неравенство примет вид

n 4 ≥ m a x [ a ; b ] f (4) (x) · (b - a) 5 2 . 88

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

f " (x) = sin 3 x 2 + 1 2 " = 3 2 cos 3 x 2 ⇒ f "" (x) = 3 2 cos 3 x 2 " = - 9 4 sin 3 x 2 ⇒ f " " " (x) = - 9 4 sin 3 x 2 " = - 27 8 cos 3 x 2 ⇒ f (4) (x) = - 27 8 cos 3 x 2 " = 81 16 sin 3 x 2

Область определения f (4) (x) = 81 16 sin 3 x 2 принадлежит интервалу - 81 16 ; 81 16 , а сам отрезок интегрирования [ 0 ; π) имеет точку экстремума, из этого следует, что m a x [ 0 ; π ] f (4) (x) = 81 16 .

Производим подстановку:

n 4 ≥ m a x [ a ; b ] f (4) (x) · (b - a) 5 2 . 88 ⇔ n 4 ≥ 81 16 · π - 0 5 2 . 88 ⇔ ⇔ n 4 > 537 . 9252 ⇔ n > 4 . 8159

Получили, что n – натуральное число, тогда его значение может быть равным n = 5 , 6 , 7 … для начала необходимо взять значение n = 5 .

Действия производить аналогично предыдущему примеру. Необходимо вычислить шаг. Для этого

h = b - a 2 n = π - 0 2 · 5 = π 10

Найдем узлы x i = a + i · h , i = 0 , 1 , . . . , 2 n , тогда значение подынтегральной функции будет иметь вид

i = 0: x i = x 0 = a + i · h = 0 + 0 · π 10 = 0 ⇒ f (x 0) = f (0) = sin 3 · 0 2 + 1 2 = 0 . 5 i = 1: x i = x 1 = a + i · h = 0 + 1 · π 10 = π 10 ⇒ f (x 1) = f (π 10) = sin 3 · π 10 2 + 1 2 ≈ 0 . 953990 . . . i = 10: x i = x 10 = a + i · h = 0 + 10 · π 10 = π ⇒ f (x 10) = f (π) = sin 3 · π 2 + 1 2 ≈ - 0 . 5 7 π 10

4 π 5 9 π 10 π f (x i) 1 . 207107 0 . 809017 0 . 343566 - 0 . 087785 - 0 . 391007 - 0 . 5

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

∫ 0 π sin 3 x 2 + 1 2 ≈ h 3 f (x 0) + 4 ∑ i = 1 n f (x 2 i - 1) + 2 ∑ i = 1 n - 1 f (x 2 i) + f (x 2 n) = = π 30 · 0 , 5 + 4 · 0 . 953990 + 1 . 487688 + 1 . 207107 + + 0 . 343566 - 0 . 391007 + 2 · 1 . 309017 + 1 . 451056 + + 0 . 809017 - 0 . 87785 - 0 . 5 = = 2 . 237650

Метод Симпсона позволяет нам получать приближенное значение определенного интеграла ∫ 0 π sin 3 x 2 + 1 2 d x ≈ 2 . 237 с точностью до 0 , 001 .

При вычислении формулой Ньютона-Лейбница получим в результате

∫ 0 π sin 3 x 2 + 1 2 d x = - 2 3 cos 3 x 2 + 1 2 x 0 π = = - 3 2 cos 3 π 2 + π 2 - - 2 3 cos 0 + 1 2 · 0 = π 2 + 2 3 ≈ 2 . 237463

Ответ: ∫ 0 π sin 3 x 2 + 1 2 d x ≈ 2 . 237

Замечание

В большинстве случаях нахождение m a x [ a ; b ] f (4) (x) проблематично. Поэтому применяется альтернатива – метод парабол. Его принцип подробно разъясняется в разделе метода трапеций. Метод парабол считается предпочтительным способом для разрешения интеграла. Вычислительная погрешность влияет на результат n . Чем меньше его значение, тем точнее приближенное искомое число.

Если вы заметили ошибку в тексте, пожалуйста, выделите её и нажмите Ctrl+Enter

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

Рассмотрим произвольный интеграл

Воспользуемся заменой переменной таким образом, чтобы границы отрезка интегрирования вместо стали [-1,1], для этого введем переменную z:

Тогда и

Рассмотрим задачу интерполирования полиномом второй степени (параболой) подынтегральной функции, используя в качестве узлов три равноудаленные узловые точки – z = -1, z = 0, z = +1 (шаг равен 1, длина отрезка интегрирования равна 2). Обозначим соответствующие значения подынтегральной функции в узлах интерполяции

Система уравнений для нахождения коэффициентов полинома

Проходящего через три точки , и

примет вид

или

Коэффициенты легко могут быть получены

Вычислим теперь значение интеграла от интерполяционного многочлена

Путем обратной замены переменной вернемся к исходному интегралу. Учтем, что

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

При необходимости, исходный отрезок интегрирования может быть разбит на N сдвоенных отрезков, к каждому из которых применяется формула Симпсона. Шаг интерполирования при этом составит

Для первого отрезка интегрирования узлами интерполирования будут являться точки a, a+h, a+2h, для второго – a+2h, a+3h, a+4h, третьего a+4h, a+5h, a+6h и т.д. Приближенное значение интеграла получается суммированием N площадей:

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

Что эквивалентно

Так как

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

Увеличение точности

Здесь мы рассмотрим так называемый процесс Эйткена. Он дает возможность оценить погрешность метода и указывает алгоритм уточнения результатов. Расчет проводится последовательно три раза при различных шагах разбиения h 1 , h 2 , h 3 , причем их отношения постоянны: h 2 / h 1 = h 3 / h 2 = q (например, при делении шага пополам q=0.5). Пусть в результате численного интегрирования получены значения интеграла I 1 , I 2 , I 3 . Тогда уточненное значение интеграла вычисляется по формуле

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

.

Уточнение значения интеграла можно также проводить методом Рунге-Ромберга.

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

На практике нередко встречаются случаи, когда подынтегральная функция меняется по-разному на отдельных участках отрезка интегрирования. Это обстоятельство требует такой организации экономичных численных алгоритмов, при которой они автоматически приспосабливались бы к характеру изменения функции. Такие алгоритмы называются адаптивными (приспосабливающимися). Они позволяют вводить разные значения шага интегрирования на отдельных участках отрезка интегрирования. Это дает возможность уменьшить машинное время без потери точности результатов расчета. Подчеркнем, что этот подход используется обычно при задании подынтегральной функции y=f(x) в виде формулы, а не в табличном виде.

Рассмотрим принцип работы адаптивного алгоритма. Первоначально отрезок разбиваем на n частей. В дальнейшем каждый такой элементарный отрезок делим последовательно пополам. Окончательное число шагов, их расположение и размеры зависят от подынтегральной функции и допустимой погрешности e .

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

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

Процесс деления отрезка пополам и вычисления уточненных значений продолжается до тех пор, пока их разность станет не больше некоторой заданной величины d i, зависящей от e и h:

.

Аналогичная процедура проводится для всех n элементарных отрезков. Величина принимается в качестве искомого значения интеграла. Условия и соответствующий выбор величин d i обеспечивают выполнение условия

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

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

Приравнивая его нулю, получим квадратное уравнение

Тот из двух корней квадратного уравнения (36), который меньше по модулю, определяет новое приближение

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

Метод парабол построен по образцу методов третьего порядка. Однако замена производных разделенными разностями приводит к существенному уменьшению скорости сходимости. Рассуждениями, аналогичными рассуждениям в п. 7, можно показать, что вблизи простого корня выполняется соотношение

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

В методе парабол «разболтка» счета вблизи корня сказывается еще сильней, чем в методе секущих, ибо в расчете участвуют вторые разности. Тем не менее корни можно найти, с хорошей точностью; для определения оптимального числа итераций удобно пользоваться приемом Гарвика, описанном в п. 7.

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

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

Для многочлена частное есть тоже многочлен; поэтому последовательно удаляя найденные корни, можно найти все корни исходного многочлена.

Замечание 1. Если - многочлен высокой степени, то возникают дополнительные трудности. Многочлен быстро возрастает при увеличении аргумента, поэтому в программе для ЭВМ должна быть страховка от переполнения. Обычно вводят масштабные множители, величина которых связана с диапазоном изменения аргумента.

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

являются последовательные целые числа слегка измененный многочлен имеет такие корни:

(здесь приведен только один знак после запятой). Кратные или близкие корни могут быть слабо устойчивыми даже при меньших степенях многочлена.

Замечание 3. Для удаления вычисленных корней надо... делить многочлен. Это вносит погрешность округления в коэффициенты и влияет на точность нахождения следующих корней. На практике отмечено, что если сначала удалять меньшие по модулю корни, точность падает мало, но если надать удаление с больших корней, точность может упасть катастрофически. Поэтому за начальное приближение берут тогда итерации обычно сходятся к наименьшему по модулю корню. Его удаляют и по такому же начальному приближению ищут следующий корень и т. д. При такой организации вычислений потеря точности будет небольшой.


Возникает задача о численном вычислении определенного интеграла, решаемая с помощью формул, носящих название квад­ратурных.

Напомним простейшие формулы численного интегрирования.

Вычислим приближенное численное значение . Интервал интегрирования [а, b] разобьем на п равных частей точками деле­ния
, называемыми узлами квадра­турной формулы. Пусть в узлах известны значения
:


Величина

называется интервалом интегрирования или шагом. Отметим, что в практике -вычислений число я выбирают небольшим, обычно оно не больше 10-20.На частичном интервале

подынтегральную функцию заменяют интерполяционным много­членом


который на рассматриваемом интервале приближенно представ­ляет функцию f (х).

а) Удержим в интерполяционном многочлене только один первый член, тогда


Полученная квадратная формула

называется формулой прямоугольников.

б) Удержим в интерполяционном многочлене два первых члена, тогда

(2)

Формула (2) называется формулой трапеций.

в) Интервал интегрирования
разобьем на четное число 2n равных частей, при этом шаг интегрирования h будет равен. На интервале
длиной 2h подынтегральную функцию заменим интерполяционным многочленом второй сте­пени, т. е. удержим в многочлене три первых члена:

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

(3)

Формулы (1), (2) и (3) имеют простой геометрический смысл. В формуле прямоугольников подынтегральная функция f(х) на интервале
заменяется отрезком прямой у = ук, параллельной оси абсцисс, а в формуле трапеций - отрезком прямой
и вычисляется соответственно площадь прямо­угольника и прямолинейной трапеции, которые затем сумми­руются. В формуле Симпсона функция f(х) на интервале
длиной 2h заменяется квадратным трехчленом - параболой
вычисляется площадь криволинейной параболической трапеции, затем площади суммируются.

ЗАКЛЮЧЕНИЕ

В завершении работы, хочется отметить ряд особенностей применения рассмотренных выше методов. Каждый способ приближённого решения определённого интеграла имеет свои преимущества и недостатки, в зависимости от поставленной задачи следует использовать конкретные методы.

Метод замены переменных является одним из основных методов вычисления неопределенных интегралов. Даже в тех случаях, когда мы интегрируем каким-либо другим методом, нам часто приходится в промежуточных вычислениях прибегать к замене переменных. Успех интегрирования зависит в значительной степени от того, сумеем ли мы подобрать такую удачную замену переменных, которая упростила бы данный интеграл.

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

Таким образом, интегрирование всякой рациональной дроби сводится к интегрированию многочлена и нескольких простейших дробей.

Интеграл от любой рациональной функции может быть выражен через элементарные функции в конечном виде, а именно:

    через логарифмы- в случаях простейших дробей 1 типа;

    через рациональные функции- в случае простейших дробей 2 типа

    через логарифмы и арктангенсы- в случае простейших дробей 3 типа

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

Формула Ньютона – Лейбница представляет собой общий подход к нахождению определенных интегралов.

Что касается приемов вычисления определенных интегралов, то они практически ничем не отличаются от всех тех приемов и методов.

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

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

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

Если необходимо получить наиболее точный результат, идеально подходит метод Симпсона .

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

Понравилась статья? Поделитесь с друзьями!
Была ли эта статья полезной?
Да
Нет
Спасибо, за Ваш отзыв!
Что-то пошло не так и Ваш голос не был учтен.
Спасибо. Ваше сообщение отправлено
Нашли в тексте ошибку?
Выделите её, нажмите Ctrl + Enter и мы всё исправим!