Тема: чисельне інтегрування, квадратурні формули

Лабораторна робота 15.

(2г.)

Мета: Отримати відомості проквадратурні формули, чисельні методи інтегрування

та навчитися застосовувати ці методи до конкретних задач.

Теоретичні відомості.

Якщо функція f (х) неперервна на відрізку [a; b] і відома її первісна F, то справедлива формула Ньютона – Лейбниця . Проте цією формулою неможливо скористатися, якщо первісну F не можна виразити у відомих (традиційно в елементарних) функціях, або якщо функцію f задано таблично або графічно. У цих випадках необхідно будувати методи для наближеного обчислення визначених інтегралів. Найчастіше застосовують квадратурні формули.

Означення. Квадратурні формули – це формули вигляду ≈ . Суму в правій частині формули називають квадратурною сумою, числа х_k і А_k називають відповідно вузлами і коефіцієнтами квадратурної формули.

Найпростіші квадратурні формули можна отримати з наочних міркувань. Наприклад,

нехай функція f (х) близька до лінійної. Тоді природно замінити інтеграл площею трапеції з висотою b – a та основами f (a) і f (b).

f (a)

f (b)

a b

Рисунок 1

В результаті отримаємо формулу трапецій:

≈ (b – a) .

Більш складні квадратурні формули отримують за допомогою апарату інтерполювання. Наприклад, саме так отримана квадратурна формула

≈ (f (a) + 4 f () + f (b)),

яку називають квадратурною формулою Сімпсона.

Якщо функцію f (х) задано на великому проміжку [a; b], то точність розглянутих квадратурних формул стає неприйнятною. Тому для обчислення застосовують відповідну узагальнену квадратурну формулу. Це означає, що відрізок [a; b] ділять на рівні відрізки і на кожному з них застосовують дану квадратурну формулу. Наприклад, узагальнена формула трапецій виглядає так:

≈ (f (a) + 2+ f (b)). (1)

А узагальнена формула Сімпсона має вигляд

≈ (f (a) + 4+ 2+ f (b)). (2)

Хід роботи.

Задача 1. Обчислити наближені значення інтегралу функції f (x) = e^{sin x} cos2x на відрізку [0;1] за узагальненою формулою трапецій з кроками інтегрування h = 0,2 0,1 0,05.

Розв’язання. Спочатку побудуємо електронну таблицю значень функції f (x) у вузлах інтегрування. Надамо чарункам електронної таблиці таких значень:

Символ ↓ означає копіювання попередніх чарунок. Спочатку h = 0,2 і треба копіювати у стовпці А до значення 1, тобто до чарунки А7. В результаті отримаємо таку таблицю:

У стовпці С до отриманої таблиці додамо відповідні коефіцієнти узагальненої формули трапецій:

Аналогічно дістанемо такі таблиці для h = 0,1 і h = 0,05:

На основі цих обчислень побудуємо таблицю для отримання відповідних наближень інтегралу функції f (x) = e^{sin x} cos2x.

Тут у стовпці А крок інтегрування h. У стовпці В інтегральна сума, у стовпці С – відповідне наближене значення інтеграла згідно з узагальненою формулою трапецій (1).

В результаті отримаємо таку таблицю:

Задача 2. Обчислити наближені значення функції f (x) = e^xsinx на відрізку [0;1] за узагальненою формулою Сімпсона з кроками інтегрування h = 0,5 0,25 0,125.

Розв’язання. Зауважимо, що кроки інтегрування h = 0,5 0,25 0,125 є найменшими можливими оскільки за узагальненою формулою Сімпсона (2) кількість відрізків m = , на які вузли інтегрування ділять [0;1], обов’язково є парним. Надамо таких значень чарункам електронної таблиці:

В результаті, наприклад, при h = 0,125 отримаємо таку таблицю:

У стовпці С до отриманої таблиці додамо відповідні коефіцієнти узагальненої формули Сімпсона (2):

Аналогічно з двома іншими таблицями:

На основі цих обчислень побудуємо таблицю для отримання відповідних наближень інтегралу функції f (x) = e^{sin x} cos2x, ідентичну таблиці задачі 1:

Тут у стовпці С – наближення інтегралу згідно з узагальненою формулою Сімпсона (2).

В результаті отримаємо таку таблицю:

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