Вычисление определенных интегралов по формуле Ньютона-Лейбница не всегда возможно. Многие подынтегральные функции не имеют первообразных в виде элементарных функций, поэтому мы во многих случаях не можем найти точное значение определенного интеграла по формуле Ньютона-Лейбница. С другой стороны, точное значение не всегда и нужно. На практике нам часто достаточно знать приближенное значение определенного интеграла с некоторой заданной степенью точности (например, с точностью до одной тысячной). В этих случаях нам на помощь приходят методы численного интегрирования, такие как метод прямоугольников, метод трапеций, метод Симпсона (парабол) и т.п.
В этой статье подробно разберем метод прямоугольников для приближенного вычисления определенного интеграла.
Сначала остановимся на сути этого метода численного интегрирования, выведем формулу прямоугольников и получим формулу для оценки абсолютной погрешности метода. Далее по такой же схеме рассмотрим модификации метода прямоугольников, такие как метод правых прямоугольников и метод левых прямоугольников. В заключении рассмотрим подробное решение характерных примеров и задач с необходимыми пояснениями.
Навигация по странице.
- Суть метода прямоугольников.
- Метод средних прямоугольников.
- Формула метода средних прямоугольников.
- Оценка абсолютной погрешности метода средних прямоугольников.
- Метод левых прямоугольников и метод правых прямоугольников.
- Примеры применения метода прямоугольников при приближенном вычислении определенных интегралов.
- Метод прямоугольников
- Суть метода прямоугольников
- Метод средних прямоугольников
- Оценка абсолютной погрешности метода средних прямоугольников
- Метод левых прямоугольников и метод правых прямоугольников
- Примеры применения метода прямоугольников при приближенном вычислении определенных интегралов
- Замечание
- Итоги
- Метод прямоугольников
- Сущность метода прямоугольников
- Готовые работы на аналогичную тему
- Методы левых и правых прямоугольников
- Погрешность метода прямоугольников
- 💥 Видео
Видео:Численное интегрирование: Методы Левых Правых прямоугольников, Трапеций, Симпсона c++Скачать
Суть метода прямоугольников.
Пусть функция y = f(x) непрерывна на отрезке [a; b] . Нам требуется вычислить определенный интеграл .
Обратимся к понятию определенного интеграла. Разобьем отрезок [a;b] на n частей точками . Внутри каждого отрезка выберем точку . Так как по определению определенный интеграл есть предел интегральных сумм при бесконечном уменьшении длины элементарного отрезка разбиения , то любая из интегральных сумм является приближенным значением интеграла .
Суть метода прямоугольников заключается в том, что в качестве приближенного значения определенного интеграла берут интегральную сумму (далее мы покажем, какую именно интегральную сумму берут в методе прямоугольников).
Видео:Метод прямоугольников для нахождения определенного интегралаСкачать
Метод средних прямоугольников.
Формула метода средних прямоугольников.
Если отрезок интегрирования [a;b] разбить на РАВНЫЕ части длины h точками (то есть ) и в качестве точек выбрать СЕРЕДИНЫ элементарных отрезков (то есть ), то приближенное равенство можно записать в виде . Это и есть формула метода прямоугольников. Ее еще называют формулой средних прямоугольников из-за способа выбора точек .
называют шагом разбиения отрезка [a;b] .
Приведем графическую иллюстрацию метода средних прямоугольников.
Из чертежа видно, что подынтегральная функция y=f(x) приближается кусочной ступенчатой функцией на отрезке интегрирования.
С геометрической точки зрения для неотрицательной функции y=f(x) на отрезке [a;b] точное значение определенного интеграла представляет собой площадь криволинейной трапеции, а приближенное значение по методу прямоугольников – площадь ступенчатой фигуры.
Оценка абсолютной погрешности метода средних прямоугольников.
Перейдем к оценке абсолютной погрешности метода прямоугольников. Сначала оценим погрешность на элементарном интервале. Погрешность метода прямоугольников в целом будет равна сумме абсолютных погрешностей на каждом элементарном интервале.
На каждом отрезке имеем приближенное равенство . Абсолютную погрешность метода прямоугольников на i -ом отрезке вычисляем как разность между точным и приближенным значением определенного интеграла: . Так как есть некоторое число и , то выражение в силу четвертого свойства определенного интеграла можно записать как . Тогда абсолютная погрешность формулы прямоугольников на i -ом элементарном отрезке будет иметь следующий вид
Если считать, что функция y = f(x) имеет в точке и некоторой ее окрестности производные до второго порядка включительно, то функцию y = f(x) можно разложить в ряд Тейлора по степеням с остаточным членом в форме Лагранжа:
По свойствам определенного интеграла равенства можно интегрировать почленно:
где .
Таким образом, и .
Абсолютная погрешность формулы прямоугольников на отрезке [a; b] равна сумме погрешностей на каждом элементарном интервале, поэтому
и .
Полученное неравенство представляет собой оценку абсолютной погрешности метода прямоугольников.
Видео:Метод левых, правых и средних прямоугольниковСкачать
Метод левых прямоугольников и метод правых прямоугольников.
Перейдем к модификациям метода прямоугольников.
— это формула метода левых прямоугольников.
— это формула метода правых прямоугольников.
Отличие от метода средних прямоугольников заключается в выборе точек не в середине, а на левой и правой границах элементарных отрезков соответственно.
Абсолютная погрешность методов левых и правых прямоугольников оценивается как .
Видео:3.2 Численное интегрирование (лекция)Скачать
Примеры применения метода прямоугольников при приближенном вычислении определенных интегралов.
Перейдем к решению примеров, в которых требуется вычислить приближенное значение определенного интеграла методом прямоугольников.
В основном, встречаются два типа задач. В первом случае задается количество интервалов, на которые разбивается отрезок интегрирования. Во втором случае задается допустимая абсолютная погрешность.
Формулировки задач примерно следующие:
- вычислить приближенно определенный интеграл методом прямоугольников, разбив отрезок интегрирования на n частей;
- Методом прямоугольников найти приближенное значение определенного интеграла с точностью до одной сотой (одной тысячной и т.п.).
Разберем каждый случай.
Сразу оговоримся, что в примерах подынтегральные функции будем брать такие, чтобы можно было найти их первообразные. В этом случае мы сможем вычислить точное значение определенного интеграла и сравнить его с приближенным значением, полученным по методу прямоугольников.
Вычислить определенный интеграл методом прямоугольников, разбив отрезок интегрирования на 10 частей.
В нашем примере a = 4, b = 9, n = 10 , .
Внимательно посмотрим на формулу прямоугольников .
Чтобы ее применить, нам нужно вычислить шаг h и значения функции в точках .
Вычислим шаг: .
Так как , то .
Для i = 1 имеем . Находим соответствующее значение функции .
Для i = 2 имеем . Находим соответствующее значение функции .
И так продолжаем вычисления до i = 10 .
Для удобства представим результаты в виде таблицы.
Подставляем полученные значения в формулу прямоугольников:
Значение исходного определенного интеграла можно вычислить по формуле Ньютона-Лейбница:
.
Первообразная подынтегральной функции была найдена интегрированием по частям.
Как видите, точное значение определенного интеграла отличается от значения, полученного по методу прямоугольников для n = 10 , менее чем на шесть сотых долей единицы.
Вычислите приближенное значение определенного интеграла методами левых и правых прямоугольников с точностью до одной сотой.
По условию имеем a = 1, b = 2 , .
Чтобы применить формулы правых и левых прямоугольников нам необходимо знать шаг h , а чтобы вычислить шаг h необходимо знать на какое число отрезков n разбивать отрезок интегрирования. Так как в условии задачи нам указана точность вычисления 0.01 , то число n мы можем найти из оценки абсолютной погрешности методов левых и правых прямоугольников.
Нам известно, что . Следовательно, если найти n , для которого будет выполняться неравенство , то будет достигнута требуемая степень точности.
Найдем — наибольшее значение модуля первой производной подынтегральной функции на отрезке [1; 2] . В нашем примере это сделать достаточно просто.
Графиком функции производной подынтегральной функции является парабола, ветви которой направлены вниз, на отрезке [1; 2] ее график монотонно убывает. Поэтому достаточно вычислить модули значения производной на концах отрезка и выбрать наибольшее:
В примерах со сложными подынтегральными функциями Вам может потребоваться теория раздела наибольшее и наименьшее значение функции.
Таким образом:
Число n не может быть дробным (так как n – натуральное число – количество отрезков разбиения интервала интегрирования). Поэтому, для достижения точности 0.01 по методу правых или левых прямоугольников, мы можем брать любое n = 9, 10, 11, … Для удобства расчетов возьмем n = 10 .
Формула левых прямоугольников имеет вид , а правых прямоугольников . Для их применения нам требуется найти h и для n = 10 .
Итак,
Точки разбиения отрезка [a; b] определяются как .
Для i = 0 имеем и .
Для i = 1 имеем и .
И так далее до i = 10 .
Полученные результаты удобно представлять в виде таблицы:
Подставляем в формулу левых прямоугольников:
Подставляем в формулу правых прямоугольников:
Вычислим точное значение определенного интеграла по формуле Ньютона-Лейбница:
Очевидно, точность в одну сотую соблюдена.
Во многих случаях нахождение наибольшего значения модуля первой производной (или второй производной для метода средних прямоугольников) подынтегральной функции на отрезке интегрирования является очень трудоемкой процедурой.
Поэтому можно действовать без использования неравенства для оценки абсолютной погрешности методов численного интегрирования. Хотя оценки предпочтительнее.
Для методов правых и левых прямоугольников можно использовать следующую схему.
Берем произвольное n (например, n = 5 ) и вычисляем приближенное значение интеграла. Далее удваиваем количество отрезков разбиения интервала интегрирования, то есть, берем n = 10 , и вновь вычисляем приближенное значение определенного интеграла. Находим разность полученных приближенных значений для n = 5 и n = 10 . Если абсолютная величина этой разности не превышает требуемой точности, то в качестве приближенного значения определенного интеграла берем значение при n = 10 , предварительно округлив его до порядка точности. Если же абсолютная величина разности превышает требуемую точность, то вновь удваиваем n и сравниваем приближенные значения интегралов для n = 10 и n = 20 . И так продолжаем до достижения требуемой точности.
Для метода средних прямоугольников действуем аналогично, но на каждом шаге вычисляем треть модуля разности полученных приближенных значений интеграла для n и 2n . Этот способ называют правилом Рунге.
Вычислим определенный интеграл из предыдущего примера с точностью до одной тысячной по методу левых прямоугольников.
Не будем подробно останавливаться на вычислениях.
Для n = 5 имеем , для n = 10 имеем .
Так как , тогда берем n = 20 . В этом случае .
Так как , тогда берем n = 40 . В этом случае .
Так как , то, округлив 0.01686093 до тысячных, утверждаем, что значение определенного интеграла равно 0.017 с абсолютной погрешностью 0.001 .
В заключении остановимся на погрешности методов левых, правых и средних прямоугольников более детально.
Из оценок абсолютных погрешностей видно, что метод средних прямоугольников даст большую точность, чем методы левых и правых прямоугольников для заданного n . В то же время, объем вычислений одинаков, так что использование метода средних прямоугольников предпочтительнее.
Если говорить о непрерывных подынтегральных функциях, то при бесконечном увеличении числа точек разбиения отрезка интегрирования приближенное значение определенного интеграла теоретически стремиться к точному. Использование методов численного интегрирования подразумевает использование вычислительной техники. Поэтому следует иметь в виду, что при больших n начинает накапливаться вычислительная погрешность.
Еще заметим, если Вам требуется вычислить определенный интеграл с некоторой точностью, то промежуточные вычисления проводите с более высокой точностью. Например, Вам требуется вычислить определенный интеграл с точностью до одной сотой, тогда промежуточные вычисления проводите с точностью как минимум до 0.0001 .
При вычислении определенного интеграла методом прямоугольников (методом средних прямоугольников) пользуемся формулой и оцениваем абсолютную погрешность как .
Для метода левых и правых прямоугольников пользуемся формулами и соответственно. Абсолютную погрешность оцениваем как .
Видео:Метод левых прямоугольников, разбор задачиСкачать
Метод прямоугольников
Не всегда имеется возможность вычисления интегралов по формуле Ньютона-Лейбница. Не все подынтегральные функции имеют первообразные элементарных функций, поэтому нахождение точного числа становится нереальным. При решении таких задач не всегда необходимо получать на выходе точные ответы. Существует понятие приближенного значения интеграла, которое задается методом числового интегрирования типа метода прямоугольников, трапеций, Симпсона и другие.
Данная статья посвящена именно этому разделу с получением приближенных значений.
Будет определена суть метода Симпсона, получим формулу прямоугольников и оценки абсолютной погрешности, метод правых и левых треугольников. На заключительном этапе закрепим знания при помощи решения задач с подробным объяснением.
Видео:Метод средних прямоугольниковСкачать
Суть метода прямоугольников
Если функция y = f ( x ) имеет непрерывность на отрезке [ a ; b ] и необходимо вычислить значение интеграла ∫ a b f ( x ) d x .
Необходимо воспользоваться понятием неопределенного интеграла. Тогда следует разбить отрезок [ a ; b ] на количество n частей x i — 1 ; x i , i = 1 , 2 , . . . . , n , где a = x 0 x 1 x 2 . . . x n — 1 x n = b . В промежутке отрезка x i — 1 ; x i , i = 1 , 2 , . . . , n выберем точку со значением ζ i . Из определения имеем, что существует определенный тип интегральных сумм при бесконечном уменьшении длины элементарного отрезка, который уже разбили. Это выражается формулой λ = m a x i = 1 , 2 , . . . , n ( x i — x i — 1 ) → 0 , тогда получаем, что любая из таких интегральных сумм – приближенное значение интеграла ∫ a b f ( x ) d x ≈ ∑ i = 1 n f ( ζ i ) · ( x i — x i — 1 ) .
Суть метода прямоугольников выражается в том, что приближенное значение считается интегральной суммой.
Видео:Метод левых и метод правых прямоугольниковСкачать
Метод средних прямоугольников
Если разбить интегрируемый отрезок [ a ; b ] на одинаковые части точкой h , то получим a = x 0 , x 1 = x 0 + h , x 2 = x 0 + 2 h , . . . , x — 1 = x 0 + ( n — 1 ) h , x n = x 0 + n h = b , то есть h = x i — x i — 1 = b — a n , i = 1 , 2 , . . . , n . Серединами точек ζ i выбираются элементарные отрезки x i — 1 ; x i , i = 1 , 2 , . . . , n , значит ζ i = x i — 1 + h 2 , i = 1 , 2 , . . . , n .
Тогда приближенное значение ∫ a b f ( x ) d x ≈ ∑ i = 1 n f ( ζ i ) · ( x i — x i — 1 ) записывается таким образом ∫ a b f ( x ) d x ≈ h · ∑ i = 1 n f ( ζ i ) x i — 1 + h 2 . Данная формула называется формулой метода прямоугольников.
Такое название метод получает из-за характера выбора точек ζ i , где шаг разбиения отрезка берется за h = b — a n .
Рассмотрим на приведенном ниже рисунке данный метод.
Чертеж явно показывает, что приближение к кусочной ступенчатой функции
y = f x 0 + h 2 , x ∈ [ x 0 ; x 1 ) f x 1 + h 2 , x ∈ [ x 1 ; x 2 ) . . . f x n — 1 + h 2 , x ∈ [ x n — 1 ; x n ] происходит на всем пределе интегрирования.
С геометрической стороны мы имеем, что неотрицательная функция y = f ( x ) на имеющемся отрезке [ a ; b ] имеет точное значение определенного интеграла и выглядит как криволинейная трапеция, площадь которой необходимо найти. Рассмотрим на рисунке, приведенном ниже.
Видео:метод прямоугольниковСкачать
Оценка абсолютной погрешности метода средних прямоугольников
Для оценки абсолютной погрешности необходимо выполнить ее оценку на заданном интервале. То есть следует найти сумму абсолютных погрешностей каждого интервала. Каждый отрезок x i — 1 ; x i , i = 1 , 2 , . . . , n имеет приближенное равенство ∫ x i — 1 x i f ( x ) d x ≈ f x i — 1 + h 2 · h = f x i — 1 + h 2 · ( x i — x i — 1 ) . Абсолютная погрешность данного метода треугольников δ i , принадлежащей отрезку i , вычисляется как разность точного и приближенного определения интеграла . Имеем, что δ i = ∫ x i — 1 x i f ( x ) d x — f x i — 1 + h 2 · x i — x i — 1 . Получаем, что f x i — 1 + h 2 является некоторым числом, а x i — x i — 1 = ∫ x i — 1 x i d x , тогда выражение f x i — 1 + h 2 · x i — x i — 1 по 4 свойству определения интегралов записывается в форме f x i — 1 + h 2 · x i — x i — 1 = ∫ x — 1 x f x i — 1 + h 2 d x . Отсюда получаем, что отрезок i имеет абсолютную погрешность вида
δ i = ∫ x i — 1 x i f ( x ) d x — f x i — 1 + h 2 · x i — x i — 1 = = ∫ x i — 1 x i f ( x ) d x — ∫ x i — 1 x i x i — 1 + h 2 d x = ∫ x i — 1 x i f ( x ) = — f x i — 1 + h 2 d x
Если взять, что функция y = f ( x ) имеет производные второго порядка в точке x i — 1 + h 2 и ее окрестностях, тогда y = f ( x ) раскладывается в ряд Тейлора по степеням x — x i — 1 + h 2 с остаточным членом в форме разложения по Лагранжу. Получаем, что
f ( x ) = f x i — 1 + h 2 + f ‘ x i — 1 + h 2 · x — x i — 1 + h 2 + + f » ( ε i ) x — x i — 1 + h 2 2 2 ⇔ ⇔ f ( x ) = f ( x i — 1 + h 2 ) = f ‘ x i — 1 + h 2 · x — x i — 1 + h 2 + + f » ( ε i ) x — x i — 1 + h 2 2 2
Исходя из свойства определенного интеграла, равенство может интегрироваться почленно. Тогда получим, что
∫ x i — 1 x i f ( x ) — f x i — 1 + h 2 d x = ∫ x i — 1 x i f ‘ x i — 1 + h 2 · x — x i — 1 + h 2 d x + + ∫ x i — 1 x i f » ε i · x — x i — 1 + h 2 2 2 d x = = f ‘ x i — 1 + h 2 · x — x i — 1 + h 2 2 2 x i — 1 x i + f » ε i · x — x i — 1 + h 2 3 6 x i — 1 x i = = f ‘ x i — 1 + h 2 · x i — h 2 2 2 — x i — 1 — x i — 1 + h 2 2 2 + + f » ε i · x i — h 2 3 6 — x i — 1 — x i — 1 + h 2 3 6 = = f ‘ x i — 1 + h 2 · h 2 8 — h 2 8 + f » ( ε i ) · h 3 48 + h 3 48 = f » ε i · h 3 24
где имеем ε i ∈ x i — 1 ; x i .
Отсюда получаем, что δ i = ∫ x i — 1 x i f ( x ) — f x i — 1 + h 2 d x = f » ε i · h 3 24 .
Абсолютная погрешность формулы прямоугольников отрезка [ a ; b ] равняется сумме погрешностей каждого элементарного интервала. Имеем, что
δ n = ∑ i = 1 n ∫ x i — 1 x i f ( x ) — f x i — 1 + h 2 d x и δ n ≤ m a x x ∈ [ a ; b ] f » ( x ) · n · h 3 24 = m a x x ∈ [ a ; b ] f » ( x ) = b — a 3 24 n 2 .
Неравенство является оценкой абсолютной погрешности метода прямоугольников.
Видео:Метод СимпсонаСкачать
Метод левых прямоугольников и метод правых прямоугольников
Для модификации метода рассмотрим формулы.
∫ a b f ( x ) d x ≈ h · ∑ i = 0 n — 1 f ( x i ) является формулой левых треугольников.
∫ a b f ( x ) d x ≈ h · ∑ i = 1 n f ( x i ) является формулой правых треугольников.
Рассмотрим на примере рисунка, приведенного ниже.
Отличием метода средних прямоугольников считается выбор точек не по центру, а на левой и правой границах данных элементарных отрезков.
Такая абсолютная погрешность методов левых и правых треугольников можно записать в виде
δ n ≤ m a x x ∈ [ a ; b ] f ‘ ( x ) · h 2 · n 2 = m a x x ∈ [ a ; b ] f ‘ ( x ) · ( b — a ) 2 2 n
Видео:3. Численные методы расчета определенного интеграла: прямоугольников, трапеции, парабол (Симпсона)Скачать
Примеры применения метода прямоугольников при приближенном вычислении определенных интегралов
Необходимо рассмотреть решение примеров, где нужно вычислять примерное значение имеющегося определенного интеграла при помощи метода прямоугольников. Рассматривают два типа решения заданий. Суть первого случая – задание количества интервалов для разбивания отрезка интегрирования. Суть второго заключается в наличии допустимой абсолютной погрешности.
Формулировки задач выглядят следующим образом:
- произвести приближенное вычисление определенного интеграла при помощи метода прямоугольников, разбивая на nколичество отрезков интегрирования;
- найти приближенное значение определенного интеграла методом прямоугольников с точностью до одной сотой.
Рассмотрим решения в обоих случаях.
В качестве примера выбрали задания, которые поддаются преобразованию для нахождения их первообразных. Тогда появляется возможность вычисления точного значения определенного интеграла и сравнения с приближенным значением при помощи метода прямоугольников.
Произвести вычисление определенного интеграла ∫ 4 9 x 2 sin x 10 d x при помощи метода прямоугольников, разбивая отрезок интегрирования на 10 частей.
Из условия имеем, что a = 4 , b = 9 , n = 10 , f ( x ) = x 2 sin x 10 . Для применения ∫ a b f ( x ) d x ≈ h · ∑ i = 1 n f x i — 1 + h 2 необходимо вычислить размерность шага h и значение функции f ( x ) = x 2 sin x 10 в точках x i — 1 + h 2 , i = 1 , 2 , . . . , 10 .
Вычисляем значение шага и получаем, что
h = b — a n = 9 — 4 10 = 0 . 5 .
Потому как x i — 1 = a + ( i — 1 ) · h , i = 1 , . . . , 10 , тогда x i — 1 + h 2 = a + ( i — 1 ) · h + h 2 = a + i — 0 . 5 · h , i = 1 , . . . , 10 .
Так как i = 1 , то получаем x i — 1 + h 2 = x 0 + h 2 = a + ( i — 0 . 5 ) · h = 4 + ( 1 — 0 . 5 ) · 0 . 5 = 4 . 25 .
После чего необходимо найти значение функции
f x i — 1 + h 2 = f x 0 + h 2 = f ( 4 . 25 ) = 4 . 25 2 sin ( 4 . 25 ) 10 ≈ — 1 . 616574
При i = 2 получаем x i — 1 + h 2 = x 1 + h 2 = a + i — 0 . 5 · h = 4 + ( 2 — 0 . 5 ) · 0 . 5 = 4 . 75 .
Нахождение соответствующего значения функции получает вид
f x i — 1 + h 2 = f x 1 + h 2 = f ( 4 . 75 ) = 4 . 75 2 sin ( 4 . 75 ) 10 ≈ — 2 . 254654
Вычисления производятся до i = 10 .
Представим эти данные в таблице, приведенной ниже.
i | 1 | 2 | 3 | 4 | 5 |
x i — 1 + h 2 | 4 . 25 | 4 . 75 | 5 . 25 | 5 . 75 | 6 . 25 |
f x i — 1 + h 2 | — 1 . 616574 | — 2 . 254654 | — 2 . 367438 | — 1 . 680497 | — 0 . 129606 |
i | 6 | 7 | 8 | 9 | 10 |
x i — 1 + h 2 | 6 . 75 | 7 . 25 | 7 . 75 | 8 . 25 | 8 . 75 |
f x i — 1 + h 2 | 2 . 050513 | 4 . 326318 | 5 . 973808 | 6 . 279474 | 4 . 783042 |
Значения функции необходимо подставить в формулу прямоугольников. Тогда получаем, что
∫ 4 9 x 2 sin x 10 d x ≈ h · ∑ i = 1 n f x i — 1 + h 2 = = 0 . 5 · — 1 . 616574 — 2 . 25654 — 2 . 367438 — 1 . 680497 — 0 . 129606 + + 2 . 050513 + 4 . 326318 + 5 . 973808 + 6 . 279474 + 4 . 783042 = = 7 . 682193
Исходный интеграл можно вычислить при помощи формулы Ньютона-Лейбница. Получаем, что
∫ 4 9 x 2 · sin x 10 d x = — 1 10 x 2 · cos x + 1 5 x · sin x + 1 5 cos x 4 9 = = 7 5 cos 4 — 4 5 sin 4 — 79 10 cos 9 + 9 5 sin 9 ≈ 7 . 630083
Находим первообразную выражения — 1 10 x 2 · cos x + 1 5 x · sin x + 1 5 cos x соответствующую функции f ( x ) = x 2 sin x 10 . Нахождение производится методом интегрирования по частям.
Отсюда видно, что определенный интеграл отличается от значения, полученном при решении методом прямоугольников, где n = 10 , на 6 долей единицы. Рассмотрим на рисунке, приведенном ниже.
Вычислить приближенного значение определенного интеграла ∫ 1 2 ( — 0 . 03 x 3 + 0 . 26 x — 0 . 26 ) d x при помощи метода левых и правых прямоугольников с точностью до одной сотой.
Из условия мы имеем, что a = 1 , b = 2 и f ( x ) = — 0 . 03 x 3 + 0 . 26 x — 0 . 26 .
Для применения формулы правых и левых прямоугольников нужно знать размерность шага h , а для его вычисления разбиваем отрезок интегрирования на n отрезков. По условию имеем, что точность должна быть до 0 , 01 , тогда нахождение n возможно при помощи оценки абсолютной погрешности методов левых и правых прямоугольников.
Известно, что δ n ≤ m a x x ∈ [ a ; b ] f ‘ ( x ) · ( b — a ) 2 2 n . Для достижения необходимой степени точности необходимо найти такое значение n , для которого неравенство m a x x ∈ [ a ; b ] f ‘ ( x ) · ( b — a ) 2 2 n ≤ 0 . 01 будет выполнено.
Найдем наибольшее значение модуля первой производной, то есть значение m a x x ∈ [ a ; b ] f ‘ ( x ) подынтегральной функции f ( x ) = — 0 . 03 x 3 + 0 . 26 x — 0 . 26 , определенной на отрезке [ 1 ; 2 ] . В нашем случае необходимо выполнить вычисления:
f ‘ ( x ) = — 0 . 03 x 3 + 0 . 26 x — 0 . 26 ‘ = — 0 . 09 x 2 + 0 . 26
Парабола является графиком подынтегральной функции с ветвями, направленными вниз, определенная на отрезке [ 1 ; 2 ] , причем с монотонно убывающим графиком. Необходимо произвести вычисление модулей значений производных на концах отрезков, а из них выбрать наибольшее значение. Получаем, что
f ‘ ( 1 ) = — 0 . 09 · 1 2 + 0 . 26 = 0 . 17 f ‘ ( 2 ) = — 0 . 09 · 2 2 + 0 . 26 = 0 . 1 → m a x x ∈ [ 1 ; 2 ] f ‘ ( x ) = 0 . 17
Решение сложных подынтегральных функций подразумевает обращение к разделу наибольше и наименьшее значение функции.
Тогда получаем, что наибольшее значение функции имеет вид:
m a x x ∈ [ a ; b ] f ‘ ( x ) · ( b — a ) 2 2 n ≤ 0 . 01 ⇔ ⇔ 0 . 17 · ( 2 — 1 ) 2 2 n ≤ 0 . 01 ⇔ 0 . 085 n ≤ 0 . 01 ⇔ n ≥ 8 . 5
Дробность числа n исключается, так как n является натуральным числом. Чтобы прийти к точности 0 . 01 , используя метод правых и левых прямоугольников, не обходимо выбирать любое значение n . Для четкости расчетов возьмем n = 10 .
Тогда формула левых прямоугольников примет вид ∫ a b f ( x ) d x ≈ h · ∑ i = 0 n — 1 f ( x i ) , а правых — ∫ a b f ( x ) d x ≈ h · ∑ i = 1 n f ( x i ) . Для применения их на практике необходимо найти значение размерности шага h и f ( x i ) , i = 0 , 1 , . . . , n , где n = 10 .
h = b — a n = 2 — 1 10 = 0 . 1
Определение точек отрезка [ a ; b ] производится с помощью x i = a + i · h , i = 0 , 1 , . . . , n .
При i = 0 , получаем x i = x 0 = a + i · h = 1 + 0 · 0 . 1 = 1 и f ( x i ) = f ( x 0 ) = f ( 1 ) = — 0 . 03 · 1 3 + 0 . 26 · 1 — 0 . 26 = — 0 . 03 .
При i = 1 , получаем x i = x 1 = a + i · h = 1 + 1 · 0 . 1 = 1 . 1 и f ( x i ) = f ( x 1 ) = f ( 1 . 1 ) = — 0 . 03 · ( 1 . 1 ) 3 + 0 . 26 · ( 1 . 1 ) — 0 . 26 = — 0 . 01393 .
Вычисления производятся до i = 10 .
Вычисления необходимо представить в таблице, приведенной ниже.
i | 0 | 1 | 2 | 3 | 4 | 5 |
x i | 1 | 1 . 1 | 1 . 2 | 1 . 3 | 1 . 4 | 1 . 5 |
f ( x i ) | — 0 . 03 | — 0 . 01393 | 0 . 00016 | 0 . 01209 | 0 . 02168 | 0 . 02875 |
i | 6 | 7 | 8 | 9 | 10 |
x i | 1 . 6 | 1 . 7 | 1 . 8 | 1 . 9 | 2 |
f ( x i ) | 0 . 03312 | 0 . 03461 | 0 . 03304 | 0 . 02823 | 0 . 02 |
Подставим формулу левых треугольников
∫ 1 2 ( — 0 . 03 x 3 + 0 . 26 x — 0 . 26 ) d x ≈ h · ∑ i = 0 n — 1 f ( x i ) = = 0 . 1 · — 0 . 03 — 0 . 01393 + 0 . 00016 + 0 . 01209 + 0 . 02168 + + 0 . 02875 + 0 . 03312 + 0 . 03461 + 0 . 03304 + 0 . 02823 = = 0 . 014775
Подставляем в формулу правых треугольников
∫ 1 2 ( — 0 . 03 x 3 + 0 . 26 x — 0 . 26 ) d x ≈ h · ∑ i = 1 n f ( x i ) = = 0 . 1 · — 0 . 01393 + 0 . 00016 + 0 . 01209 + 0 . 02168 + 0 . 02875 + + 0 . 03312 + 0 . 03461 + 0 . 03304 + 0 . 02823 + 0 . 02 = 0 . 019775
Произведем вычисление по формуле Ньютона-Лейбница:
∫ 1 2 ( — 0 . 03 x 3 + 0 . 26 x — 0 . 26 ) d x = = — 0 . 03 x 4 4 + 0 . 13 x 2 — 0 . 26 x 1 2 = 0 . 0175
Рассмотрим рисунок, приведенный ниже.
Видео:6 Теория: Определенный интеграл Метод левых, правых, центральных прямоугольников, трапеций, СимпсонаСкачать
Замечание
Нахождение наибольшего значения модуля первой производной является трудоемкой работой, поэтому можно исключить использование неравенства для оценивания абсолютной погрешности и методов численного интегрирования. Разрешено применять схему.
Берем значение n = 5 для вычисления приближенного значения интеграла. Необходимо удвоить количество отрезков интегрирования, тогда n = 10 , после чего производится вычисление примерного значения. необходимо найти разность этих значений при n = 5 и n = 10 . Когда разность не соответствует требуемой точности, то приближенным значением считается n = 10 с округлением до десятка.
Когда погрешность превышает необходимую точность, то производится удваивание n и сравнивание приближенных значений. Вычисления производятся до тех пор, пока необходимая точность не будет достигнута.
Для средних прямоугольников выполняются аналогичные действия, но вычисления на каждом шаге требуют разности полученных приближенных значений интеграла для n и 2 n . Такой способ вычисления называется правилом Рунге.
Произведем вычисление интегралов с точностью до одной тысячной при помощи метода левых прямоугольников.
При n = 5 получаем, что ∫ 1 2 ( — 0 . 03 x 3 + 0 . 26 x — 0 . 26 ) d x ≈ 0 . 0116 , а при n = 10 — ∫ 1 2 ( — 0 . 03 x 3 + 0 . 26 x — 0 . 26 ) d x ≈ 0 . 014775 . Так как имеем, что 0 . 0116 — 0 . 014775 = 0 . 003175 > 0 . 001 , возьмем n = 20 . Получаем, что ∫ 1 2 ( — 0 . 03 x 3 + 0 . 26 x — 0 . 26 ) d x ≈ 0 . 01619375 . Имеем 0 . 014775 — 0 . 01619375 = 0 . 00141875 > 0 . 001 , возьмем значение n = 40 , тогда получим ∫ 1 2 ( — 0 . 03 x 3 + 0 . 26 x — 0 . 26 ) d x ≈ 0 . 01686093 . Имеем, что 0 . 1619375 — 0 . 01686093 = 0 . 00066718 0 . 001 , тогда после округления значения проверим, что ∫ 1 2 ( — 0 . 03 x 3 + 0 . 26 x — 0 . 26 ) d x равняется значению 0 , 017 с погрешностью 0 , 001 . Из оценок абсолютных погрешностей видно, что данный метод дает максимальную точность в отличие от метода левых и правых координат для заданного n . Отдается предпочтение методу средних прямоугольников.
Непрерывные подынтегральные функции при бесконечном разделении на отрезки данное приближенно число стремится к точному. Чаще всего такой метод выполняется при помощи специальных программ на компьютере. Поэтому чем больше значение n , тем больше вычислительная погрешность.
Для наиболее точного вычисления необходимо выполнять точные промежуточные действия, желательно с точностью до 0 , 0001 .
Видео:Формула СимпсонаСкачать
Итоги
Для вычисления неопределенного интеграла методом прямоугольников следует применять формулу такого вида ∫ a b f ( x ) d x ≈ h · ∑ i = 1 n f ( ζ i ) x i — 1 + h 2 и оценивается абсолютная погрешность с помощью δ n ≤ m a x x ∈ [ a ; b ] f ‘ ‘ ( x ) · n · h 3 24 = m a x x ∈ [ a ; b ] f ‘ ‘ ( x ) · b — a 3 24 n 2 .
Для решения с помощью методов правых и левых прямоугольников применяют формулы, имеющие вид, ∫ a b f ( x ) d x ≈ h · ∑ i = 0 n — 1 f ( x i ) и ∫ a b f ( x ) d x ≈ h · ∑ i = 1 n f ( x i ) . Абсолютная погрешность оценивается при помощи формулы вида δ n ≤ m a x x ∈ [ a ; b ] f ‘ ( x ) · h 2 · n 2 = m a x x ∈ [ a ; b ] f ‘ ( x ) · b — a 2 2 n .
Видео:Метод трапецийСкачать
Метод прямоугольников
Вы будете перенаправлены на Автор24
Рассмотрим задачу, в которой требуется вычислить определённый интеграл $int^b_a f(x)dx$, при этом функция $f(x)$ является непрерывной на промежутке $left[a;bright]$. Обычно, если существует возможность, интегралы вычисляются через нахождение первообразной, но так как это не всегда возможно, прибегают к использованию приближённых методов.
К наиболее часто используемым приближённым методам относят:
- Метод прямоугольников;
- Метод трапеций;
- Метод Симпсона или иначе метод парабол.
В данной статье мы подробно расcмотрим метод прямоугольников.
Видео:Численное интегрированиеСкачать
Сущность метода прямоугольников
Рассмотрим нахождение определённого интеграла от функции $f(x)$ с точки зрения геометрии. Интеграл $int^b_a f(x)dx$ в данном случае есть не что иное, как площадь фигуры, ограниченной сверху графиком $f(x)$, по бокам прямыми $x=a$ и $x=b$, а снизу осью абсцисс.
Рисунок 1. Метод средних прямоугольников
Для того чтобы найти площадь всей фигуры, можно воспользоваться определением интеграла и разбить всю фигуру на равные сегменты одной и той же длины. Точки на оси абсцисс, которые будут разбивать фигуру, обозначим как $x_i$. Нулевая точка при разбиении $x_0=a$, а конечная точка $x_n=b$. Для того чтобы вычислить длину одного сегмента, воспользуемся формулой:
В методе средних прямоугольников каждый сегмент заменяется на прямоугольник, за высоту которого принимается ордината середины отрезка. Получается, что площадь одного такого прямоугольника равна $S_i= frac cdot f(ξ_i)$, а площадь всей фигуры будет равна:
$int^b_a f(x)dx=fraccdot (f( ξ_0)+f( ξ_1)+. +f( ξ_)$, где $x_i≤ ξ_i≤x_$
Готовые работы на аналогичную тему
Эта формула позволяет не вычислять напрямую площадь искомой фигуры, ограниченной кривой линией, а заменить её приблизительной площадью ступенчатой фигуры, состоящей из прямоугольников.
При использовании метода средних прямоугольников так как $ξ_i=frac<x_i+x_>=x_<i+frac>$, тогда $f( ξ_i)=f(x_<i+frac>)$ обозначим как $y_<i+ frac>$,
и формула примет вид:
Эта формула называется формулой средних прямоугольников.
Видео:Метод трапеций при вычислении определенного интегралаСкачать
Методы левых и правых прямоугольников
Данные методы отличаются от метода средних прямоугольников тем, что здесь в качестве ординаты для элементарного прямоугольника выбирается либо крайнее левое значение функции $f(x)$ (и тогда метод называется методом левых прямоугольников), либо крайнее правое, и тогда метод носит название метода правых прямоугольников.
Формула для применения метода левых прямоугольников выглядит так:
$int^b_a f(x)dx=fraccdot (y_0 + y_1 + y_)left(2right)$
Формула для метода правых прямоугольников:
$int^b_a f(x)dx=fraccdot (y_1 + y_2 + y_n)left(3right)$
Формулы $(1), (2), (3)$ иначе также называются квадратурными составными формулами.
Видео:Формула СимпсонаСкачать
Погрешность метода прямоугольников
Для того чтобы оценить общую погрешность метода прямоугольников, необходимо рассмотреть каждый из элементарных сегментов кривой по отдельности. Общая погрешность в таком случае представляет собой сумму погрешностей всех погрешностей сегментов.
Итак, рассмотрим, чему равна погрешность на одном сегменте.
Площадь одного сегмента вычисляется по приближённой формуле:
Погрешность будем определять по разнице со значением первообразной, вычисленной с помощью формулы Ньютона-Лейбница: $δ_i= int^_<x_> f(x)dx — f(x_+frac<x_-x_>) cdot (x_i-x_)left(5right)$
Так как в левой части равенства $x_-x_$ есть не что иное как $int^x_<x_>dx$ — длина элементарного отрезка, его можно заменить на $dx$. Перепишем правую часть равенства $(4)$, используя это:
Допуская, что фунцкия $f(x)$ дважды дифференцируема в точке $x=x_-x_$ и вокруг неё, разложим её в бесконечную сумму степенных функций, используя ряды Тейлора и формулу Лагранжа:
Применим полученное для подстановки: $f(x)-f(x_+ frac)=f’(x_+frac) cdot(x-(x_ + frac))+f’’(ε_i) cdot frac<(x-(x_+frac))^2>left(7right)$
В конечном итоге для элементарного сегмента $left[x_;x_iright]$ имеем:
Для всей же фигуры погрешность полученной площади составит:
и в конечном виде:
Данная формула используется для получения погрешности при использовании формулы для средних прямоугольников.
Формула для погрешности методов правых и левых прямоугольников выводится аналогичным способом и имеет следующий вид:
Погрешность, полученная с использованием метода правых или левых прямоугольников для вычисления интегралов больше, чем погрешность при использовании метода средних прямоугольников. Поэтому более предпочтительным для приближённого интегрирования является именно метод средних прямоугольников.
Вычислить интеграл $int_1^2 frac=ln2$ с точностью до $0, 001$ используя формулу средних прямоугольников.
Разобьём нашу функцию на 10 равных сегментов.
В начале оценим погрешность вычисления:
В данном случае погрешность меньше либо равна:
$|δ_n|≤0.000042$, следовательно, в данном случае для разбиения можно использовать 10 сегментов.
Разобьём подынтегральную функцию на 10 отрезков, длина каждого из которых $Δx=frac=0,1$ и вычислим значение подынтегральной функции $y(x)=frac$ в середине каждого отрезка:
Сумма всех вычисленных значений функции $f(x)$ составит $6.9284$, а само значение составит:
$int_1^2 frac=frac=0.69284$ — что отвечает требуемому условию о погрешности.
Получи деньги за свои студенческие работы
Курсовые, рефераты или другие работы
Автор этой статьи Дата последнего обновления статьи: 25 02 2021
💥 Видео
Численные методы решения интегралов в MS ExcelСкачать
Как приближённо вычислить интеграл? Формула прямоугольников, трапеций и парабол (Симпсона)Скачать
Формулы для медианы треугольникаСкачать
Алгоритмы. Численное интегрированиеСкачать