МІНІСТЕРСТВО ОСВІТИ І НАУКИ УКРАЇНИ

Національний університет “Львівська політехніка”

Інститут Комп’ютерних наук та інформаційних технологій

Кафедра автоматизованих систем управління

Лабораторна робота № 5-6

з дисципліни

“Математичні методи представлення знань”

на тему:

«Обчислення означених інтегралів»

Виконав:

студент групи КН-29

Коцуба О.

Прийняв:

Биненко Б. І.

Львів – 2011

Обчислення означених інтегралів

Мета роботи: вивчити методи наближених обчислень і запрограмувати алгоритми обчислення означених інтегралів .

Порядок роботи:

1. Попереднє опрацювання теоретичного матеріалу.

2. Отримання допуску до виконання лабораторної роботи.

3. Опрацювання типового навчального завдання (прикладів).

4. Створення проекту для виконання індивідуального завдання.

5. Оформити звіт для захисту лабораторної роботи за зразком:

· назва роботи;

· мета роботи;

· порядок роботи;

· короткі теоретичні відомості;

· алгоритм розв’язування задачі;

· тексти відповідних модулів проекту;

· аналіз отриманих результатів та висновки.

6. Захист лабораторної роботи.

Короткі теоретичні відомості

1. Формули прямокутників.

Нехай на відрізку задана неперервна функція . Потрібно обчислити інтеграл

Розіб’ємо відрізок  на n рівних частин точками , i=0,1,…n-1, довжина кожної з яких дорівнює $IMAGE6$ . Через $IMAGE7$ позначимо значення функції $IMAGE8$ в точках $IMAGE9$ і складемо суми

$IMAGE10$ або $IMAGE6$ $IMAGE12$

Кожна з цих сум є інтегральною сумою для $IMAGE8$ на відрізку і тому наближено виражають означений інтеграл:

$IMAGE15$ (1)

$IMAGE16$ (1/)

Ці формули називаються формулами прямокутників. Із рис. 1 видно, що якщо $IMAGE17$додатна і зростаюча функція, то формула (1) відображає площу ступінчатої фігури, що складена із “ внутрішніх” прямокутників, а формула (1/) – площу фігури, що складена із “зовнішніх” прямокутників.

Похибка методу прямокутників дається формулою (2):

$IMAGE18$ (2)

формула прямокутник лагранж функція

Похибка при цьому буде тим меншою, чим більше число n (тобто чим менший крок поділу $IMAGE6$). Зауважимо, що формули прямокутників дають точні результати для багаточленів першого степеня.

2. Формула трапецій.

Очевидно, що можна отримати більш точне значення інтеграла, якщо дану криву  замінити не ступінчатою лінією, як це мало місце у формулі прямокутників, а вписаною ламаною (рис.2). Тоді площа криволінійної трапеції, обмеженої лініями $IMAGE21$ і  заміниться площами трапецій, обмежених зверху хордами Оскільки площа

  $IMAGE22$

Рис.1 Рис.2

першої трапеції дорівнює $IMAGE23$другої – $IMAGE24$ і т.д.,

то $IMAGE25$

або

$IMAGE26$ (3)

Формула (3) називається формулою трапецій. Число n вибирається довільним, але чим більшим це число буде, а значить, крок $IMAGE6$ меншим, тим з більшою точністю сума в правій частині наближеної рівності (3) буде давати значення інтегралу.

3. Формула парабол (Сімпсона).

Метод Сімпсона найпоширеніший і простіше застосовний для програмування. Його суть полягає в наближенні підінтегральної функції відрізками парабол.

Отже, розглянемо спочатку інтеграл $IMAGE28$, де $IMAGE29$ – парабола; $IMAGE30$, $IMAGE31$, $IMAGE32$ – деякі параметри (або числа).

Тоді

$IMAGE33$

$IMAGE34$

Нехай тепер маємо інтеграл $IMAGE35$, де $IMAGE36$ - неперервна на інтервалі функція. Якщо інтервал розбити на п рівних частинок , i=0,1,…n-1,, то заданий інтеграл І можна записати так:

$IMAGE39$

Якщо на кожному з інтегралів для проміжків $IMAGE40$ функцію $IMAGE36$замінимо параболами $IMAGE42$, що проходять через точки $IMAGE43$ ,то одержимо

$IMAGE44$

Через те, що, формула матиме вигляд:

$IMAGE45$ або

$IMAGE46$

$IMAGE47$ (4)

Формула (4) називається формулою парабол або Сімпсона. Доведено, що похибка обчислень $IMAGE48$ за формулою Сімпсона є такою:

$IMAGE49$ (5)

Проте, цією оцінкою похибки можна користуватись, якщо $IMAGE36$ є хоча б чотири рази диференційовною. Але, якщо $IMAGE36$ навіть чотири рази диференційовна, то часто оцінка четвертої похідної $IMAGE52$ може виявитись досить складною. Тому на практиці переважно користуються таким методом: обчислюють інтеграл, розділяючи інтервал, заданий границями інтегрування, один раз на n рівних частин, а другий раз на т частин. Якщо одержані двоє значень інтеграла мало відрізняються, то результат можна вважати прийнятним. Порівнюючи їх можна оцінити і точність обчислень.

Приклад. Обчислити з точністю до 0,001 інтеграл

$IMAGE53$

Р о з в ’ я з у в а н н я. За формулою (4) маємо:

при  при

	$IMAGE54$	$IMAGE55$			$IMAGE54$	$IMAGE55$		$IMAGE54$	$IMAGE55$
$IMAGE60$	-0,5	0,0000		$IMAGE60$	-0,5	0,00000	$IMAGE62$	0,05	0,0371
$IMAGE63$	-0,4	-0,1203		$IMAGE63$	-0,45	-0,0946	$IMAGE65$	0,10	0,0772
$IMAGE66$	-0,3	-0,1303		$IMAGE66$	-0,40	-0,1203	$IMAGE48$	0,15	0,1200
$IMAGE60$	-0,2	-0,1081		$IMAGE60$	-0,35	-0,1304	$IMAGE65$	0,20	0,1652
$IMAGE66$	-0,1	-0,630		$IMAGE66$	-0,30	-0,1303	$IMAGE48$	0,25	0,2122
$IMAGE60$	0	0,0000		$IMAGE60$	-0,25	-0,1204	$IMAGE48$	0,30	0,2607
$IMAGE60$	0,1	0,0772		$IMAGE60$	-0,20	-0,1081	$IMAGE48$	0,35	0,3103