Метод центральных прямоугольников. Метод левых и правых прямоугольников

1. Введение. Постановка задачи……..…………………………2стр.

2. Вывод формулы……………………………………………….3стр.

3. Дополнительный член в формуле прямоугольников……….5стр.

4. Примеры………………………………………………………..7стр.

5. Заключение……………………………………………………..9стр.

6. Список литературы…………………………………………...10стр.

Постановка задачи.

Задача вычисления интегралов возникает во многих областях прикладной математики. В большинстве случаев встречаются определённые интегралы от функций, первообразные которых не выражаются через элементарные функции. Кроме того, в приложениях приходится иметь дело с определёнными интегралами, сами подынтегральные функции не являются элементарными. Распространенными являются также случаи, когда подынтегральная функция задается графиком или таблицей экспериментально полученных значений. В таких ситуациях используют различные методы численного интегрирования, которые основаны на том, что интеграл представляется в виде предела интегральной суммы (суммы площадей), и позволяют определить эту сумму с приемлемой точностью. Пусть требуется вычислить интеграл при условии, что a и b конечны и f(x) является непрерывной функцией на всем интервале (a, b). Значение интеграла I представляет собой площадь, ограниченную кривой f(x),осью x и прямыми x=a, x=b. Вычисление I проводится путем разбиения интервала от a до b на множество меньших интервалов, приближенным нахождением площади каждой полоски, получающейся при таком разбиении, и дальнейшем суммировании площадей этих полосок.

Вывод формулы прямоугольников.

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

З а м е ч а н и е. Пусть функция f(x) непрерывна на сегменте , а

Некоторые точки сегмента . Тогда на этом сегменте найдётся точка такая, что среднее арифметическое .

В самом деле, обозначим через m и M точные грани функции f(x) на сегменте . Тогда для любого номера k справедливы неравенства . Просуммировав эти неравенства по всем номерам и поделив результат на n, получим

Так как непрерывная функция принимает любое промежуточное значение, заключённое между m и M, то на сегменте найдётся точка такая, что

.

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

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

где , а R – дополнительный член. Здесь искомая площадь криволинейной фигуры заменяется площадью некоторой состоящей из прямоугольников ступенчатой фигуры (или – если угодно – определенный интеграл заменяется интегральной суммой). Эта формула и называется формулой прямоугольников.

На практике обычно берут ; если соответствующую среднюю ординату обозначить через , то формула перепишется в виде

.

Дополнительный член в формуле прямоугольников.

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

Справедливо следующее утверждение:

У т в е р ж д е н и е. Если функция f(x) имеет на сегменте непрерывную вторую производную, то на этом сегменте найдётся такая точка

Что дополнительный член R в формуле (1) равен

(2)

Доказательство.

Оценим , считая, что функция f(x) имеет на сегменте [-h, h] непрерывную вторую производную Для этого подвергнем двукратному интегрированию по частям каждый из следующих двух интегралов:

Для первого из этих интегралов получим

Для второго из интегралов аналогично получим

Полусумма полученных для и выражений приводит к следующей формуле:

(3)

Оценим величину , применяя к интегралам и формулу среднего значения и учитывая неотрицательность функций и . Мы получим, что найдутся точка на сегменте [-h, 0] и точка на сегменте

Такие, что

В силу доказанного замечания на сегменте [-h, h] найдётся точка такая, что

Поэтому для полусуммы мы получим следующее выражение:

Подставляя это выражение в равенство (3), получим, что

(4)

. (5)

Так как величина представляет собой площадь некоторого прямоугольника с основанием (рис.1), то формулы (4) и (5) доказывают, что ошибка, совершаемая при замене указанной площадью, имеет порядок

Таким образом, формула тем точнее, чем меньше h. Поэтому для вычисления интеграла естественно представить это интеграл в виде суммы достаточно большого числа n интегралов

И к каждому из указанных интегралов применить формулу (4). Учитывая при этом, что длина сегмента равна , мы получим формулу прямоугольников (1), в которой

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

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

по формуле прямоугольников.

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

П р и м е р 1. Пусть требуется вычислить интеграл .

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

Теперь применим формулу прямоугольников

Таким образом, .

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

П р и м е р 2. Вычислим интеграл с точностью до 0,001.

Применяя формулу Ньютона-Лейбница, получим .

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

Так как для имеем (если ), то

Если взять n=10, то дополнительный член нашей формулы будет Нам придётся внести ещё погрешность, округляя значения функции; постараемся, чтобы границы этой новой погрешности разнились меньше чем на С этой целью достаточно вычислять значение функции с четырьмя знаками, с точностью до 0,00005. Имеем:


Сумма 6,9284.

.

Учитывая, что поправка к каждой ординате (а следовательно и к их среднему арифметическому) содержится между , а также принимая во внимание оценку дополнительного члена , найдём, что содержится между границами и , а следовательно, и подавно между 0,692 и 0,694. Таким образом, .

Заключение.

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

Выше для приближенного вычисления интеграла от функции f(x)

мы исходили из разбиения основного сегмента на достаточно большое число n равных частичных сегментов одинаковой длины h и из последующей замены функции f(x) на каждом частичном сегменте многочленом соответственно нулевого, первого или второго порядка.

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

Список литературы.

1. Фихтенгольц Г.М. Курс дифференциального и интегрального исчисления в 3-х томах, том II. (§§ 332, 335).

2. Ильин В.А., Позняк Э.Г. Основы математического анализа, часть I. Москва «Наука», 1982г. (Глава 12, пп.1, 2, 5).


Учебно-воспитательные задачи:

  • Дидактическая цель. Познакомить учащихся с методами приближённого вычисления определённого интеграла.
  • Воспитательная цель. Тема данного занятия имеет большое практическое и воспитательное значение. Наиболее просто к идее численного интегрирования можно подойти, опираясь на определение определённого интеграла как предела интегральных сумм. Например, если взять какое-либо достаточно мелкое разбиение отрезка [a ; b ] и построить для него интегральную сумму, то её значение можно приближённо принять за значение соответствующего интеграла. При этом важно быстро и правильно производить вычисления с привлечением вычислительной техники.

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

Обеспечение занятия

  • Раздаточный материал. Карточки-задания для самостоятельной работы.
  • ТСО. Мультипроектор, ПК, ноутбуки.
  • Оснащение ТСО. Презентации: “Геометрический смысл производной”, “Метод прямоугольников”, “Метод трапеций”. (Презентации можно взять у автора).
  • Вычислительные средства: ПК, микрокалькуляторы.
  • Методические рекомендации

Вид занятия. Интегрированное практическое.

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

Последовательность занятия.

  1. Формула прямоугольников.
  2. Формула трапеций.
  3. Решение упражнений.

План занятия

  1. Повторение опорных знаний учащихся.

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

  1. Выполнение практической работы.

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

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

Простейшим приближённым методом является метод прямоугольников. Геометрически идея способа вычисления определённого интеграла по формуле прямоугольников состоит в том, что площадь криволинейной трапеции АВСD заменяется суммой площадей прямоугольников, одна сторона которых равна , а друга - .

Если суммировать площади прямоугольников, которые показывают площадь криволинейной трапеции с недостатком [Рисунок1], то получим формулу:

[Рисунок1]

то получим формулу:

Если с избытком

[Рисунок2],

то

Значения у 0 , у 1 ,..., у n находят из равенств , к = 0, 1..., n .Эти формулы называются формулами прямоугольников и дают приближённый результат. С увеличением n результат становится более точным.

Итак, чтобы найти приближённое значение интеграла , нужно:

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


Пример 1. Вычислить по формуле прямоугольников . Найти абсолютную и относительную погрешности вычислений.

Разобьём отрезок [a, b ] на несколько (например, на 6) равных частей. Тогда а = 0, b = 3 ,

х k = a + k х
х
0 = 2 + 0 = 2
х 1 = 2 + 1 = 2,5
х 2 = 2 + 2 =3
х 3 = 2 + 3 = 3
х 4 = 2 + 4 = 4
х 5 = 2 + 5 = 4,5

f (x 0) = 2 2 = 4
f (x 1) = 2 ,5 2 = 6,25
f (x 2) = 3 2 = 9
f (x 3) = 3,5 2 = 12,25
f (x 4) = 4 2 = 16
f (x 5) = 4,5 2 = 20,25.

х 2 2,5 3 3,5 4 4,5
у 4 6,25 9 12,25 16 20,25

По формуле (1):

Для того, чтобы вычислить относительную погрешность вычислений, надо найти точное значение интеграла:



Вычисления проходили долго и мы получили довольно-таки грубое округление. Чтобы вычислить этот интеграл с меньшим приближением, можно воспользоваться техническими возможностями компьютера.

Для нахождения определённого интеграла методом прямоугольников необходимо ввести значения подынтегральной функции f(x) в рабочую таблицу Excel в диапазоне х с заданным шагом х = 0,1.

  1. Составляем таблицу данных и f(x)). х f(x). Аргумент , а в ячейку В1 – слово Функция 2 2,1 ). Затем, выделив блок ячеек А2:А3, автозаполнением получаем все значения аргумента (за правый нижний угол блока протягиваем до ячейки А32, до значения х=5 ).
  2. Далее вводим значения подынтегральной функции. В ячейку В2 необходимо записать её уравнение. Для этого табличный курсор необходимо установить в ячейку В2 и с клавиатуры ввести формулу =А2^2 (при английской раскладке клавиатуры). Нажимаем клавишу Enter . В ячейке В2 появляется 4 . Теперь необходимо скопировать функцию из ячейки В2. Автозаполнением копируем эту формулу в диапазон В2:В32.
    В результате должна быть получена таблица данных для нахождения интеграла.
  3. Теперь в ячейке В33 может быть найдено приближённое значение интеграла. Для этого в ячейку В33 вводим формулу = 0,1*, затем вызываем Мастер функций (нажатием на панели инструментов кнопки Вставка функции (f(x)) . В появившемся диалоговом окне Мастер функции-шаг 1 из 2 слева в поле Категория выбираем Математические. Справа в поле Функция - функцию Сумм. Нажимаем кнопку ОК. Появляется диалоговое окно Сумм. В рабочее поле мышью вводим диапазон суммирования В2:В31. Нажимаем кнопку ОК. В ячейке В33 появляется приближённое значение искомого интеграла с недостатком (37,955 ) .

Сравнивая полученное приближённое значение с истинным значением интеграла (39 ), можно видеть, что ошибка приближения метода прямоугольников в данном случае равна

= |39 - 37 , 955| = 1 ,045

Пример 2. Используя метод прямоугольников, вычислить с заданным шагом х = 0,05.

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

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

[Рисунок3]

Пример 3. Методом трапеций найти с шагом х = 0,1.

  1. Открываем чистый рабочий лист.
  2. Составляем таблицу данных и f(x)). Пусть первый столбец будет значениями х , а второй соответствующими показателями f(x). Для этого в ячейку А1 вводим слово Аргумент , а в ячейку В1 – слово Функция . В ячейку А2 вводится первое значение аргумента – левая граница диапазона (0 ). В ячейку А3 вводится второе значение аргумента – левая граница диапазона плюс шаг построения (0,1 ). Затем, выделив блок ячеек А2:А3, автозаполнением получаем все значения аргумента (за правый нижний угол блока протягиваем до ячейки А33, до значения х=3,1 ).
  3. Далее вводим значения подынтегральной функции. В ячейку В2 необходимо записать её уравнение (в примере синуса). Для этого табличный курсор необходимо установить в ячейку В2. Здесь должно оказаться значение синуса, соответствующее значению аргумента в ячейке А2. Для получения значения синуса воспользуемся специальной функцией: нажимаем на панели инструментов кнопку Вставка функции f(x) . В появившемся диалоговом окне Мастер функции-шаг 1 из 2 слева в поле Категория выбираем Математические. Справа в поле Функция - функцию SIN . Нажимаем кнопку ОК. Появляется диалоговое окно SIN . Наведя указатель мыши на серое поле окна, при нажатой левой кнопке сдвигаем поле вправо, чтобы открыть столбец данных (А ). Указываем значение аргумента синуса щелчком мыши на ячейке А2. Нажимаем кнопку ОК. В ячейке В2 появляется 0. Теперь необходимо скопировать функцию из ячейки В2. Автозаполнением копируем эту формулу в диапазон В2:В33. В результате должна быть получена таблица данных для нахождения интеграла.
  4. Теперь в ячейке В34 может быть найдено приближённое значение интеграла по методу трапеций. Для этого в ячейку В34 вводим формулу = 0,1*((В2+В33)/2+, затем вызываем Мастер функций (нажатием на панели инструментов кнопки Вставка функции (f(x)) . В появившемся диалоговом окне Мастер функции-шаг 1 из 2 слева в поле Категория выбираем Математические. Справа в поле Функция - функцию Сумм. Нажимаем кнопку ОК. Появляется диалоговое окно Сумм. В рабочее поле мышью вводим диапазон суммирования В3:В32. Нажимаем кнопку ОК и ещё раз ОК. В ячейке В34 появляется приближённое значение искомого интеграла с недостатком (1,997 ) .

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

  1. Решение упражнений.

И парадокс состоит в том, что по этой причине (видимо) он довольно редко встречается на практике. Неудивительно, что данная статья появилась на свет через несколько лет после того, как я рассказал о более распространённых методах трапеции и Симпсона , где упомянул о прямоугольниках лишь вскользь. Однако на сегодняшний день раздел об интегралах практически завершён и поэтому настало время закрыть этот маленький пробел. Читаем, вникаем и смотрим видео! ….о чём? Об интегралах, конечно =)

Постановка задачи уже была озвучена на указанном выше уроке, и сейчас мы быстренько актуализируем материал:

Рассмотрим интеграл . Он неберущийся. Но с другой стороны, подынтегральная функция непрерывна на отрезке , а значит, конечная площадь существует. Как её вычислить? Приближённо. И сегодня, как вы догадываетесь – методом прямоугольников.

Разбиваем промежуток интегрирования на 5, 10, 20 или бОльшее количество равных (хотя это не обязательно) отрезков, чем больше – тем точнее будет приближение. На каждом отрезке строим прямоугольник, одна из сторон которого лежит на оси , а противоположная – пересекает график подынтегральной функции. Вычисляем площадь полученной ступенчатой фигуры, которая и будет приближённой оценкой площади криволинейной трапеции (заштрихована на 1-м рисунке) .

Очевидно, что прямоугольники можно построить многими способами, но стандартно рассматривают 3 модификации:

1) метод левых прямоугольников;
2) метод правых прямоугольников;
3) метод средних прямоугольников.

Оформим дальнейшие выкладки в рамках «полноценного» задания:

Пример 1

Вычислить определённый интеграл приближённо:
а) методом левых прямоугольников;
б) методом правых прямоугольников.

Промежуток интегрирования разделить на равных отрезков, результаты вычислений округлять до 0,001

Решение : признАюсь сразу, я специально выбрал такое малое значение – из тех соображений, чтобы всё было видно на чертеже – за что пришлось поплатиться точностью приближений.

Вычислим шаг разбиения (длину каждого промежуточного отрезка) :

Метод левых прямоугольников получил своё называние из-за того,

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

Ни в коем случае не забываем, что округление следует проводить до трёх знаков после запятой – это существенное требование условия , и «самодеятельность» здесь чревата пометкой «оформите задачу, как следует».

Вычислим площадь ступенчатой фигуры, которая равна сумме площадей прямоугольников:


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

При использовании «правого» метода высОты прямоугольников равны значениям функции в правых концах промежуточных отрезков:

Вычислим недостающее значение и площадь ступенчатой фигуры:


– тут, что и следовало ожидать, приближение сильно занижено:

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

В чём их формальное различие? В первой формуле нет слагаемого , а во второй -

На практике рассчитываемые значения удобно заносить в таблицу:


а сами вычисления проводить в Экселе. И быстро, и без ошибок:

Ответ :

Наверное, вы уже поняли, в чём состоит метод средних прямоугольников:

Пример 2

Вычислить приближенно определенный интеграл методом прямоугольников с точностью до 0,01. Разбиение промежутка интегрирования начать с отрезков.

Решение : во-первых, обращаем внимание, что интеграл нужно вычислить с точностью до 0,01 . Что подразумевает такая формулировка?

Если в предыдущей задаче требовалось прОсто округлить результаты до 3 знаков после запятой (а уж насколько они будут правдивы – не важно) , то здесь найденное приближённое значение площади должно отличаться от истины не более чем на .

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

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

Почему? А он при прочих равных условиях (том же самом разбиении) даёт гораздо более точное приближение. Это строго обосновано в теории, и это очень хорошо видно на чертеже:

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

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

Вычисления, как и в предыдущем примере, удобно свести в таблицу. Длина промежуточных отрезков, понятно, та же самая: – и очевидно, что расстояние между серединами отрезков равно этому же числу. Поскольку требуемая точность вычислений составляет , то значения нужно округлять «с запасом» – 4-5 знаками после запятой:


Вычислим площадь ступенчатой фигуры:

Давайте посмотрим, как автоматизировать этот процесс:

Таким образом, по формуле средних прямоугольников:

Как оценить точность приближения? Иными словами, насколько далёк результат от истины (площади криволинейно трапеции) ? Для оценки погрешности существует специальная формула, однако, на практике её применение зачастую затруднено, и поэтому мы будем использовать «прикладной» способ:

Вычислим более точное приближение – с удвоенным количеством отрезков разбиения: . Алгоритм решения точно такой же: .

Найдём середину первого промежуточного отрезка и далее приплюсовываем к полученному значению по 0,3. Таблицу можно оформить «эконом-классом», но комментарий о том, что изменяется от 0 до 10 – всё же лучше не пропускать:


В Экселе вычисления проводятся «в один ряд» (кстати, потренируйтесь) , а вот в тетради таблицу, скорее всего, придётся сделать двухэтажной (если у вас, конечно, не сверхмелкий почерк).

Вычислим суммарную площадь десяти прямоугольников:

Таким образом, более точное приближение:

Которые я и предлагаю вам изучить!

Пример 3: Решение : вычислим шаг разбиения:
Заполним расчётную таблицу:


Вычислим интеграл приближённо методом:
1) левых прямоугольников:
;
2) правых прямоугольников:
;
3) средних прямоугольников:
.

Вычислим интеграл более точно по формуле Ньютона-Лейбница:

и соответствующие абсолютные погрешности вычислений:

Ответ :


Оценка остаточного члена формулы: , или .

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

Инструкция . Введите подынтегральную функцию f(x) , нажмите Решить. Полученное решение сохраняется в файле Word . Также создается шаблон решения в Excel . Ниже представлена видеоинструкция.

Подынтегральная функция f(x)

Пределы интегрирования до .
Точность округления 1 2 3 4 5 6
Количество интервалов разбиения n = или Шаг h =
Метод численного интегрирования функций Формула левых прямоугольников Формула правых прямоугольников Формула средних прямоугольников Формула трапеций Элементарная формула Симпсона Формула Симпсона

Правила ввода функции

Примеры
≡ x^2/(x+2)
cos 2 (2x+π) ≡ (cos(2*x+pi))^2
≡ x+(x-1)^(2/3) Это самая простая квадратурная формула вычисления интеграла, в которой используется одно значение функции
(8.5.1)
где ; h=x 1 -x 0 .
Формула (8.5.1) представляет собой центральную формулу прямоугольников. Вычислим остаточный член. Разложим в ряд Тейлора функцию y=f(x) в точке ε 0:
(8.5.2)
где ; . Проинтегрируем (8.5.2):
(8.5.3)

Во втором слагаемом подынтегральная функция нечетная, а пределы интегрирования симметричны относительно точки ε 0 . Поэтому второй интеграл равен нулю. Таким образом, из (8.5.3) следует .
Т. к. второй множитель подынтегрального выражения не меняет знак, то по теореме о среднем получим , где . После интегрирования получим . (8.5.4)
Сравнивая с остаточным членом формулы трапеций, мы видим, что погрешность формулы прямоугольников в два раза меньше, чем погрешность формулы трапеций. Этот результат верен, если в формуле прямоугольников мы берём значение функции в средней точке.
Получим формулу прямоугольников и остаточный член для интервала . Пусть задана сетка x i =a+ih, i=0,1,...,n, . Рассмотрим сетку ε i =ε 0 +ih, i=1,2,..,n, ε 0 =a-h/2. Тогда . (8.5.5)
Остаточный член .
Геометрически формула прямоугольников может быть представлена следующим рисунком:

Если функция f(x) задана таблично, то используют либо левостороннюю формулу прямоугольников (для равномерной сетки)

либо правостороннюю формулу прямоугольников

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

; .
После интегрирования получим .

Пример . Вычислить интеграл при n=5:
а) по формуле трапеций;
б) по формуле прямоугольников;
в) по формуле Симпсона;
г) по формуле Гаусса;
д) по формуле Чебышева.
Рассчитать погрешность.
Решение. Для 5-ти узлов интегрирования шаг сетки составит 0.125.
При решении будем пользоваться таблицей значений функции. Здесь f(x)=1/x.

x f(x)
x0 0.5 y0 2
x1 0.625 y1 1.6
x2 0.750 y2 1.33
x3 0.875 y3 1.14
x4 1.0 y4 1
a) формула трапеций:
I=h/2×;
I=(0.125/2)×=0.696;
R= [-(b-a)/12]×h×y¢¢(x);
f¢¢(x)=2/(x 3).
Максимальное значение второй производной функции на интервале равно 16: max {f¢¢(x)}, xÎ=2/(0.5 3)=16, поэтому
R=[-(1-0.5)/12]×0.125×16=-0.0833;
б) формула прямоугольников:
для левосторонней формулы I=h×(y0+y1+y2+y3);
I=0.125×(2+1.6+1.33+1.14)=0.759;
R=[(b-a)/6]×h 2 ×y¢¢(x);
R=[(1-0.5)/6]×0.125 2 ×16=0.02;
в) формула Симпсона:
I=(2h/6)×{y0+y4+4×(y1+y3)+2×y2};
I=(2×0.125)/6×{2+1+4×(1.6+1.14)+2×1.33}=0.693;
R=[-(b-a)/180]×h 4 ×y (4) (x);
f (4) (x)=24/(x 5)=768;
R=[-(1-0.5)/180]×(0.125) 4 ×768= - 5.2 e -4;
г) формула Гаусса:
I=(b-a)/2×;
x i =(b+a)/2+t i (b-a)/2
(A i , t i - табличные значения).
t (n=5) A (n=5)
x1 0.9765 y1 1.02 t 1 0.90617985 A 1 0.23692688
x2 0.8846 y2 1.13 t 2 0.53846931 A 2 0.47862868
x3 0.75 y3 1.33 t 3 0 A 3 0.56888889
x4 0.61 y4 1.625 t 4 -0.53846931 A 4 0.47862868
x5 0.52 y5 1.91 t 5 -0.90617985 A 5 0.23692688
I=(1-0.5)/2×(0.2416+0.5408+0.7566+0.7777+0.4525)=0.6923;
д) формула Чебышева:
I=[(b-a)/n] ×S f(x i), i=1..n,
x i =(b+a)/2+[ t i (b-a)]/2 - необходимое приведение интервала интегрирования к интервалу [­­-1;1].
Для n=5
t1 0.832498
t2 0.374541
t3 0
t4 -0.374541
t5 -0.832498
Найдем значения x и значения функции в этих точках:
x1 0,958 f(x1) 1,043
x2 0,844 f(x2) 1,185
x3 0,75 f(x3) 1,333
x4 0,656 f(x4) 1,524
x5 0,542 f(x5) 1,845
Сумма значений функции равна 6,927.
I=(1-0,5)/5×6,927=0,6927.

Формула левых прямоугольников:

Метод средних прямоугольников

Разделим отрезок на n равных частей, т.е. на n элементарных отрезков. Длина каждого элементарного отрезка. Точки деления будут: x 0 =a; x 1 =a+h; x 2 =a+2Ч h,., x n-1 =a+ (n-1) Ч h; x n =b. Эти числа будем называть узлами. Вычислим значения функции f (x) в узлах, обозначим их y 0 , y 1 ,y 2 ,., y n . Cталобыть, y 0 =f (a), y 1 =f (x 1),y 2 =f (x 2),., y n =f (b). Числа y 0 , y 1 ,y 2 ,., y n являются ординатами точек графика функции, соответствующих абсциссам x 0 , x 1 ,x 2 ,., x n. Площадь криволинейной трапеции приближенно заменяется площадью многоугольника, составленного из n прямоугольников. Таким образом, вычисление определенного интеграла сводится к нахождению суммы n элементарных прямоугольников.

Формула средних прямоугольников

Метод правых прямоугольников

Разделим отрезок на n равных частей, т.е. на n элементарных отрезков. Длина каждого элементарного отрезка. Точки деления будут: x 0 =a; x 1 =a+h; x 2 =a+2Ч h,., x n-1 =a+ (n-1) Ч h; x n =b. Эти числа будем называть узлами. Вычислим значения функции f (x) в узлах, обозначим их y 0 , y 1 ,y 2 ,., y n . Cталобыть, y 0 =f (a), y 1 =f (x 1),y 2 =f (x 2),., y n =f (b). Числа y 0 , y 1 ,y 2 ,., y n являются ординатами точек графика функции, соответствующих абсциссам x 0 , x 1 ,x 2 ,., x n. Площадь криволинейной трапеции приближенно заменяется площадью многоугольника, составленного из n прямоугольников. Таким образом, вычисление определенного интеграла сводится к нахождению суммы n элементарных прямоугольников.

Формула правых прямоугольников

Метод Симпсона

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

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


Метод трапеций

Разделим отрезок на n равных частей, т.е. на n элементарных отрезков. Длина каждого элементарного отрезка. Точки деления будут: x 0 =a; x 1 =a+h; x 2 =a+2Ч h,., x n-1 =a+ (n-1) Ч h; x n =b. Эти числа будем называть узлами. Вычислим значения функции f (x) в узлах, обозначим их y 0 , y 1 ,y 2 ,., y n . Cталобыть, y 0 =f (a), y 1 =f (x 1),y 2 =f (x 2),., y n =f (b). Числа y 0 , y 1 ,y 2 ,., y n являются ординатами точек графика функции, соответствующих абсциссам x 0 , x 1 ,x 2 ,., x n

Формула трапеций:

Формула означает, что площадь криволинейной трапеции заменяется площадью многоугольника, составленного из n трапеций (рис.5); при этом кривая заменяется вписанной в нее ломаной.

Последние материалы раздела:

Чудеса Космоса: интересные факты о планетах Солнечной системы
Чудеса Космоса: интересные факты о планетах Солнечной системы

ПЛАНЕТЫ В древние времена люди знали только пять планет: Меркурий, Венера, Марс, Юпитер и Сатурн, только их можно увидеть невооруженным глазом....

Реферат: Школьный тур олимпиады по литературе Задания
Реферат: Школьный тур олимпиады по литературе Задания

Посвящается Я. П. Полонскому У широкой степной дороги, называемой большим шляхом, ночевала отара овец. Стерегли ее два пастуха. Один, старик лет...

Самые длинные романы в истории литературы Самое длинное литературное произведение в мире
Самые длинные романы в истории литературы Самое длинное литературное произведение в мире

Книга длинной в 1856 метровЗадаваясь вопросом, какая книга самая длинная, мы подразумеваем в первую очередь длину слова, а не физическую длину....