
- •Міністерство освіти і науки україни
- •Розділ 1. Інтерполяція функцій
- •1.1. Постановка задачі інтерполяції
- •1.2. Наближене відновлення функції за допомогою інтерполяційного многочлена Лагранжа
- •Питання для самоперевірки
- •Завдання до лабораторної роботи № 1
- •Таблиця варіантів
- •Розділ 2. Чисельне диференціювання функцій за допомогою інтерполяції кубічним сплайном
- •Питання для самоперевірки
- •Індивідуальне завдання № 1
- •Розділ 3. Чисельне інтегрування
- •3.1. Загальні положення
- •3.2. Квадратурна формула Сімпсона
- •Точки поділу
- •Питання для самоперевірки
- •Завдання до лабораторної роботи № 2
- •Варіанти завдань
- •Розділ 4. Наближене розв’язання задачі коші для звичайних диференціальних рівнянь
- •4.1. Метод Ейлера
- •Питання для самоперевірки
- •8.3. Метод Рунге-Кутта четвертого порядку точності
- •Питання для самоперевірки
- •Завдання до лабораторної роботи № 3
- •Індивідуальне завдання № 2
Розділ 1. Інтерполяція функцій
1.1. Постановка задачі інтерполяції
При розв’язанні багатьох практичних задач, що виникають у різних областях, виникає необхідність у використанні теорії наближення функцій або теорії апроксимації функцій.
Суть
задачі апроксимації така. Нехай дана
деяка невідома в аналітичному сенсі
функція
і відома лише її поведінка (наприклад,
дискретні значення в точках
)
на певному відрізку
.
Таку функцію надалі будемо називатиапроксимованою
функцією.
Треба побудувати іншу функцію
,апроксимуючу
функцію,
яка б була близька до функції
з певною похибкою. При цьому потрібне
виконання таких вимог: наявність
дискретних значень функції
;
визначення класу апроксимуючих функцій,
з яких конструюється функція
;
вид критерію згоди між функціями
і
;
оцінка похибки апроксимації.
Вид
функції
залежить від класу розв'язуваних задач.
Наприклад, при дослідженні
напружено-деформованого стану методом
скінченних елементів у задачах механіки
деформованого тіла апроксимуючі функції
представляються комбінацією функцій
,
або сплайнами. Експонентні функції
мають широке застосування у фізиці при
вивченні явищ типу розпаду і нагромадження,
а тригонометричні функції – у механіці
при коливальних процесах.
Критерій згоди або близькості апроксимованої й апроксимуючої функцій визначається з умови мінімуму відстаней між ними. Наприклад, найпоширенішим критерієм є критерій Чебишева:
.
Інший критерій може бути записаний у вигляді:
.
Метод апроксимації, заснований на другому критерії, має назву методу найменших квадратів.
Якщо
при використанні критерію Чебишева
прийняти
,
то це буде означати, що значення
апроксимованої й апроксимуючої функцій
у вузлових точках відрізка
збігаються. Цей спосіб апроксимації
називаєтьсяінтерполюванням
або інтерполяцією.
Питання оцінки похибки апроксимації безумовно залежать від попередніх трьох вимог і розглядаються окремо для конкретного процесу апроксимації.
1.2. Наближене відновлення функції за допомогою інтерполяційного многочлена Лагранжа
Нехай
відомі значення деякої функції
в
різних точках
.
Введемо позначення:
Наприклад, ці значення отримані з експерименту чи знайдені за допомогою досить складних обчислень.
Виникає
задача наближеного відновлення функції
в довільній точці
.
Часто для розв’язання цієї задачі
будується алгебраїчний многочлен
степені
,
що у точках
приймає задані значення, тобто
,
і
називається інтерполяційним. Точки
називаються вузлами інтерполяції.
Наближене
відновлення функції
за формулою
називається
інтерполяцією функції
.
Якщо
розташований поза мінімальним відрізком,
що містить усі вузли інтерполяції, то
заміну функції
за зазначеною формулою називають також
екстраполяцією.
На рисунку 1 приведена блок-схема програми обчислення значення функції у точці за допомогою інтерполяційного многочлена Лагранжа.
У даній блок-схемі: n – кількість вузлових точок; xi – координати вузлових точок; yi – значення функції у вузлових точках; x0 – координата точки, у якій знаходиться значення функції; y0 – значення шуканої функції.
Приклад.
Нехай задані такі значення функції
,
,
,
.
Для
заданої таблиці значень функції
побудуємо інтерполяційний многочлен
Лагранжа у вигляді
.
Розв’язок.
Тут
.
Ясно,
що
.
Для
нашого прикладу
.
;
;
;
.
Тоді
.