
- •Тема 5. Наближення та інтерполяція функцій
- •2. Інтерполяційний многочлен Лагранжа
- •3. Похибка інтерполяції многочленом Лагранжа
- •Інтерполяція многочленом Лагранжа з рівновіддаленими вузлами
- •Інтерполяційний многочлен Ньютона
- •5. Наближення функцій, заданих таблично, методом найменших квадратів
- •Алгоритм наближення функції, заданої таблично, методом найменших квадратів
- •Тема 6 Чисельне диференціювання й інтегрування
- •1. Наближення похідних
- •2. Похибка чисельного диференціювання
- •3. Чисельне інтегрування. Поняття про квадратурні формули.
- •4. Найпростіші квадратурні формули
- •1) Формули прямокутників
- •2) Формула трапецій
- •3) Формула Сімпсона (формула парабол)
- •5. Похибки квадратурних формул
- •6. Чисельне розв’язування задачі Коші для звичайних диференціальних рівнянь
- •1) Метод Ейлера
- •2) Метод Рунге-Кутти
1) Формули прямокутників
Вважатимемо,
що вузли
розміщено так, що
,
,
,
.
Крок
в даному випадку дорівнює
.
Позначимо
,
,
– значення даної функції у вузлах.
Замінимо
задану криволінійну трапецію ступінчастою
фігурою, що складається з
прямокутників.
Основи цих прямокутників однакові і
дорівнюють
,
а висоти збігаються із значеннями
в лівих або правих точках частинних
інтервалів.
Для цих двох випадків отримуємо наступні формули прямокутників:
;
.
які називаються формулами лівих і правих прямокутників відповідно.
Широко
розповсюдженим і більш точним є вигляд
формули прямокутників, який використовує
значення функції в середніх точках
елементарних відрізків
.
Тоді отримуємо формулу
прямокутників:
. (10)
яка називається формулою середніх прямокутників.
2) Формула трапецій
Замінимо
криву
не ступінчастою лінією, як у попередньому
випадку, а ламаною, сполучивши сусідні
точки
.
Тоді площа криволінійної трапеції
наближено дорівнює сумі площ прямокутних
трапецій, обмежених зверху відрізками
цієї ламаної. Площа
-ї
частинної прямокутної трапеції дорівнює
,
де
і
– основи
-ї
трапеції, а
– її висота. Тому:
. (11)
Формула (11) називається формулою трапецій.
3) Формула Сімпсона (формула парабол)
Розіб'ємо
відрізок інтегрування
на парне число
(
)
рівних частин з кроком
.
На кожному відрізку
,
,
…,
,…,
підінтегральну функцію
замінимо інтерполяційним многочленом
другого степеня:
,
.
За
можна взяти інтерполяційний многочлен
Лагранжа другого степеня, який проходить
через точки
,
,
:
.
Площа
криволінійної трапеції, обмеженої
параболою, яка проходить через точки
,
,
,
може бути обчислена за допомогою
визначеного інтегралу:
.
Провівши
такі обчислення для кожного елементарного
відрізку
,
візьмемо суму отриманих виразів:
.
Даний
вираз для
і приймається як значення визначеного
інтеграла:
. (12)
Формула (12) називається формулою Сімпсона або формулою парабол.
5. Похибки квадратурних формул
Різницю між лівою і правою частиною квадратурної формули
, (8)
називають
її залишковим членом і позначають
:
.
Абсолютна
похибка
квадратурної формули, очевидно, залежить
від числа
частинних відрізків, на які розбивається
відрізок інтегрування
.
Наступні формули дозволяють оцінювати
абсолютні похибки квадратурних формул
при відомому числі
.
Якщо
функція
має на відрізку
неперервну похідну
і для будь-якого
,
то абсолютна похибка
наближених формул (10), (11), (12) оцінюється
формулою:
.
Якщо
функція
має на відрізку
неперервну другу похідну
і для будь-якого
,
то абсолютна похибка
формули прямокутників оцінюється
формулою:
;
формули прямокутників трапецій:
Абсолютна
похибка
квадратурної формули Сімпсона оцінюється
формулою:
.
Якщо
функція
має на відрізку
неперервну четверту похідну
і для будь-якого
,
то абсолютна похибка
формули Сімпсона оцінюється формулою:
.
Можна
вказати ще один практично зручний спосіб
підрахунку помилки квадратурної
формули
Сімпосона. Якщо
вважати, що
міняється поволі, то отримуємо наближений
вираз для шуканої помилки :
Хай
і
наближені значення інтеграла з кроком
h і кроком 2h. Тоді:
,
тобто
(число
частинних відрізків
кратно 4).
Приклад
3. Обчислити
інтеграл
за
квадратурною формулою:
а) прямокутників;
б) трапецій;
в) Сімпсона.
Оцінити похибку результатів.
Розв’язання.
Точне значення інтеграла
.
Розіб’ємо
відрізок інтегрування
на десять рівних частин:
,
.
Обчислимо значення підінтегральної
функції
в точках розбиття
,
а також в середніх точках
:
-
0
0,0
1,00000
1
0,1
0,90909
0,05
0,95238
2
0,2
0,83333
0,15
0,86957
3
0,3
0,76923
0,25
0,80000
4
0,4
0,71429
0,35
0,74074
5
0,5
0,66667
0,45
0,68966
6
0,6
0,62500
0,55
0,64516
7
0,7
0,58824
0,65
0,60606
8
0,8
0,55556
0,75
0,57143
9
0,9
0,52632
0,85
0,54054
10
1,0
0,50000
0,95
0,51282
а) Обчислимо інтеграл за формулою прямокутників (10):
Оцінимо
похибку результату. У даному випадку
,
.
.
Отже,
добутий результат потрібно записати
так: (залишаємо дві правильні цифри і
одну сумнівну):
.
а) Обчислимо інтеграл за формулою трапецій (11):
.
Оцінимо
похибку результату. У даному випадку
,
.
.
Отже,
добутий результат потрібно записати
так:
.
а) Обчислимо інтеграл за формулою Сімпсона (12):
.
.
Оцінимо
похибку результату. У даному випадку
,
.
.
Отже,
в отриманому результаті не менш як
чотири значущі цифри правильні. Порівнявши
цей результат з точним значенням
інтеграла, переконаємося, що в ньому
правильні всі п’ять значущих цифр.
Отже,
добутий результат потрібно записати
так:
.
Відзначимо, що формула Сімпсона значно точніша формули прямокутників і формули трапецій.