Студопедия

КАТЕГОРИИ:


Архитектура-(3434)Астрономия-(809)Биология-(7483)Биотехнологии-(1457)Военное дело-(14632)Высокие технологии-(1363)География-(913)Геология-(1438)Государство-(451)Демография-(1065)Дом-(47672)Журналистика и СМИ-(912)Изобретательство-(14524)Иностранные языки-(4268)Информатика-(17799)Искусство-(1338)История-(13644)Компьютеры-(11121)Косметика-(55)Кулинария-(373)Культура-(8427)Лингвистика-(374)Литература-(1642)Маркетинг-(23702)Математика-(16968)Машиностроение-(1700)Медицина-(12668)Менеджмент-(24684)Механика-(15423)Науковедение-(506)Образование-(11852)Охрана труда-(3308)Педагогика-(5571)Полиграфия-(1312)Политика-(7869)Право-(5454)Приборостроение-(1369)Программирование-(2801)Производство-(97182)Промышленность-(8706)Психология-(18388)Религия-(3217)Связь-(10668)Сельское хозяйство-(299)Социология-(6455)Спорт-(42831)Строительство-(4793)Торговля-(5050)Транспорт-(2929)Туризм-(1568)Физика-(3942)Философия-(17015)Финансы-(26596)Химия-(22929)Экология-(12095)Экономика-(9961)Электроника-(8441)Электротехника-(4623)Энергетика-(12629)Юриспруденция-(1492)Ядерная техника-(1748)

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




Задача численного интегрирования

ОПРЕДЕЛЕННЫХ ИНТЕГРАЛОВ

Лекция 47. ПРИБЛИЖЕННЫЕ МЕТОДЫ ВЫЧИСЛЕНИЯ

План:

1. Задача численного интегрирования.

2. Формулы прямоугольников.

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

4. Формула парабол (Симпсона).

 

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

Рис. 47.1.
Рассмотрим функцию , непрерывную на отрезке и для наглядности . Тогда в силу геометрического смысла определенного интеграла, он численно равен площади криволинейной трапеции, ограниченной графиком функции y=f(x), прямыми х=а и х = b, отрезком [ a;b ] оси Ох, т.е. = S (рис. 47.1).

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

 

2. Формулы прямоугольников

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

1. С помощью точек хо, х1, х2,…, хп=b разобьём отрезок [ a;b ] на п равных частей. Тогда длина каждого отрезка будет равна . Её называют шагом разбиения и обозначают h.

2. На каждом отрезке [ хо1 ], [ х12 ] … [ хп-1п ] построим прямоугольники высотой f (х 0), f (х 1), f (х 2),…, f (хп-1).

3. Найдем площадь каждого прямоугольника как произведение его длины на ширину:

; .

Тогда сумма площадей всех прямоугольников S будет равна

.

Вынесем за скобки: .

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

(1) – формула прямоугольников.

Если в качестве высоты прямоугольников брать значения функции не в левом, а в правом конце отрезка (рис. 47.2), то формула прямоугольников будет иметь вид:

(2) – формула прямоугольников.

Пример 47.1. Вычислите определенный интеграл :

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

б) по формулам прямоугольников (число точек деления п = 4).

Решение. а) = - точное значение .

б) Выпишем подынтегральную функцию и рассмотрим ее на отрезке [0; 2] (, ). Поскольку число точек деления отрезка равно 4, найдем ширину каждого отрезка (шаг) по формуле : .

Для удобства вычислений все расчеты будем выполнять в электронных таблицах Microsoft Excel. В качестве шапки таблицы можно предложить следующий вариант:

В столбце А будет указываться номер выполняемого шага: 0, 1, 2… п.

В столбце В будут располагаться значения х 0, х 1, х 2 …. хп. Поскольку хо, то в ячейку В2 занесем значение 0. Чтобы найти значение х 1, которое будет находиться в ячейке В3, достаточно к началу промежутка а прибавить ширину шага h. В ячейке В3 будет находиться число 0 + 0,5 = 0,5. Для нахождения каждого последующего значения х к предыдущему необходимо прибавлять ширину шага до тех пор, пока хп не будет равно b.

В столбце С будут содержаться значения функции в точках х 0, х 1, х 2хп, необходимые для расчета значения определенного интеграла по формулам (1) и (2). Чтобы их получить достаточно ввести формулу в ячейку С2. В нашем примере она будет иметь вид: =B2^3. Для заполнения столбца оставшихся значений можно воспользоваться возможностями автозаполнения. Тогда таблица будет иметь вид:

Рассчитаем приближенное значение определенного интеграла по формуле прямоугольников (1) в ячейке D2, по формуле прямоугольников (2) в ячейке Е2. Множитель нами уже вычислен, он равен 0,5, поэтому формулы в ячейках D2 и Е2 будут иметь вид:

D2: =0.5*СУММ(C2:C5) (поскольку суммируются значения функции, начиная с f (х 0)и заканчивая f (хп-1));

Е2: =0.5*СУММ(C3:C6)(суммируются значения функции, начиная с f(х 1) и заканчивая f(хп)).

Результирующая таблица будет следующей:

Видим, что полученные по формулам прямоугольников значения определённого интеграла (2,25 и 6,25) достаточно сильно отличаются от его реального значения (4). Очевидно, что с увеличением числа точек деления приближенное значение определенного интеграла будет ближе к реальному. Но все равно метод прямоугольников относится к числу наименее точных. Рассмотрим другие методы численного интегрирования.

 

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

1. С помощью точек хо, х1, х2,…, хп=b разобьём отрезок [ a;b ] на п равных частей длиной .

2. На каждом отрезке [ хо1 ], [ х12 ] … [ хп-1п ] построим трапеции, соединив отрезками точки (х 0; f (х 0)) и (х 1; f (х 1)), (х 1; f (х 1)) и (х 2; f (х 2)),…, (хп -1; f (х п -1)) и (хп; f (х п))

3. Найдем площадь каждой трапеции как произведение полусуммы её оснований на высоту. Длины оснований первой трапеции будут равны f (х 0) и f (х 1), а высота . Тогда

.

Для второй трапеции длины оснований равны f (х 1) и f (х 2), а высота та же: . Тогда .

Аналогично

Найдем сумму площадей всех трапеций S:

+ +…+ .

Вынесем за скобки: =

= =

= или

.

Поскольку сумма площадей всех трапеций S приближенно равна площади криволинейной трапеции, то можно считать, что

(3)– формула трапеций.

Пример 47.2. Вычислите приближенное значение определенного интеграла по формуле трапеций (число точек деления п = 4).

Решение. Воспользуемся решением примера 47.1. Рассмотрим функцию на отрезке [0; 2], который разбит на четыре части шириной .

В уже созданной в Microsoft Excel таблице в ячейку D4 запишем формулу для расчета приближенного значения определенного интеграла по формуле трапеций (3).

Поскольку шаг равен =0,5, то его нужно умножить на скобку, содержащую полусумму первого и последнего значения функции и сумму всех остальных значений функции из столбца С. Тогда формула в ячейке D4 будет иметь вид: =0.5*((C2+C6)/2+СУММ(C3:C5)).

Расчетная таблица будет следующей:

Полученное по формуле трапеций значение определённого интеграла (4,25) ближе к реальному значению (4), чем значения, вычисленные по формулам прямоугольников. Рассмотрим последний метод – метод парабол - и оценим его точность.

 




Поделиться с друзьями:


Дата добавления: 2014-12-27; Просмотров: 924; Нарушение авторских прав?; Мы поможем в написании вашей работы!


Нам важно ваше мнение! Был ли полезен опубликованный материал? Да | Нет



studopedia.su - Студопедия (2013 - 2024) год. Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав! Последнее добавление




Генерация страницы за: 0.025 сек.