- •В. Н. Бєловодський Вісім лекцій по чисельних методах
- •7.080407 «Комп'ютерний еколого-економічний моніторинг»)
- •Лекція 1.Елементи теорії погрішностей
- •.Типи й джерела погрішностей
- •1.2.Абсолютні й відносні погрішності наближених чисел
- •1.3. Погрішності виконання арифметичних операцій
- •.Погрішність обчислення функції
- •1.5.Запис наближених чисел
- •1.6.Правила дій над наближеними числами
- •1.7.Погрішності при машинному поданні чисел
- •1.8. Варіанти індивідуальних завдань
- •Лекція 2. Інтерполяція функцій
- •2.1. Постановка задачі
- •2.2. Алгебраїчна інтерполяція, існування й одиничність інтерполяційного багаточлена
- •2.3. Інтерполяційний багаточлен Лагранжа
- •2.4. Кінцеві й розділені різниці
- •2.5. Інтерполяційний багаточлен Ньютона
- •2.6. Порівняльний аналіз інтерполяційних багаточленів
- •2.7. Погрішності інтерполяційних формул
- •2.8. Інтерполяційні формули для рівновіддалених вузлів
- •2.9. Сплайн - інтерполяція
- •2.10. Варіанти індивідуальних завдань
- •Лекція 3. Методи рішення систем лінійних алгебраїчних рівнянь
- •3.1. Попередні зауваження
- •3.2. Точні методи рішення
- •3.3. Наближені методи рішення
- •3.4. Збіжність і погрішність наближених методів
- •3.6. Варіанти індивідуальних завдань
- •Лекція 4. Вирішення нелінійних рівнянь
- •4.1. Попередні зауваження
- •4.2. Методи, засновані на алгебраїчній інтерполяції
- •4.3. Метод послідовних наближень
- •4.4. Завдання
- •Лекція 5. Рішення систем нелінійних рівнянь
- •Метод ітерацій
- •Метод Ньютона
- •5.3. Порівняльний аналіз методів
- •5.4. Завдання
- •Вимоги до програми
- •Вимоги до звіту
- •6.1. Вступні зауваження
- •6.2. Формули Ньютона-Котеса
- •6.3. Найпростіші квадратурні правила
- •6.4. Погрішності квадратурних формул
- •6.5. Поняття про методи Монте-Карло
- •6.6. Завдання
- •Варіанти завдань
- •Лекція 7. Методи рішення звичайних диференціальних рівнянь
- •Вступні зауваження
- •Аналітичні методи
- •Чисельні методи рішення. Правило Рунге
- •Завдання
- •Лекція 8. Основи спектрального аналізу
- •Елементи загальної теорії
- •Дискретна форма рядів Фур'є
- •8.4. Завдання
2.5. Інтерполяційний багаточлен Ньютона
За аналогією з (2.6) розглянемо багаточлен ступеня n .
(2.7) |
Визначимо його значення у вузлових точках. Так, безпосередньо з (2.7) маємо
Далі,
,
але, з урахуванням (2.6) при це рівняється.
Т.о.
Аналогічним образом одержуємо, що
Таким чином, багаточлен (2.7) приймає задані значення в заданих точках і, отже, є інтерполяційним. Він і називається інтерполяційним багаточленом Ньютона. Помітимо, що в силу одиничності (п. 2.2) розходження між інтерполяційними багаточленами Лагранжа й Ньютона чисто зовнішні.
2.6. Порівняльний аналіз інтерполяційних багаточленів
Зовнішні розходження інтерполяційних багаточленів спричиняють ряд особливостей, що впливають на зручність їхнього застосування. Відзначимо деякі з них. Насамперед укажемо, що загальним є те, що нумерація точок абсолютно довільна й ніяк не пов'язана з їхніми числовими значеннями.
До числа достоїнств багаточленів Лагранжа можна віднести ту обставину, що при фіксованих xi, x і зміні yi коефіцієнти не перераховуються. До числа особливостей, – роль коефіцієнтів. Їх можна розглядати, як свого роду, вагові коефіцієнти приyi і по величині можна судити про значимість можливих змінyi. Ця властивість виявляється корисним при теоретичному аналізі.
Як достоїнства багаточлена Ньютона можна відзначити, що при додаванні додаткових вузлових точок відсутня необхідність перерахування коефіцієнтів. У силу відзначеної довільності нумерації вузлових точок до вже наявного багаточлена просто додаються додаткові доданки. Багаточлен Лагранжа такоi особливістi не має.
2.7. Погрішності інтерполяційних формул
Нижче зупинимося лише на основних елементах методики оцінки погрішності інтерполяційних формул, постараючись зберегти її логічну схему.
Нехай y=f(x), - дана функція, а Pn(x)- її інтерполяційний багаточлен. Тоді погрішність Rn(x) інтерполяції функції в точцi х дорівнює
Rn(x)=f(x) – Pn(x),
звідки
f(x)=Pn(x)+Rn(x) (2.8)
Додамо тепер точку x до заданих вузлових точок і розглянемо інтерполяційний багаточлен Pn+1(z) у формі Ньютона, побудований уже по (n+2) точкам х0, х1…хn , x . Його значення в тоці x дорівнює:
Pn+1(x)=f(x)+( x-x0)f(x0, х1)+...+( х-х0)...( х-хn)f(x0,...хn, х)
Але, по побудові, , а першідоданків правої частини являють собою. Таким чином,
.
Порівнюючи тепер це співвідношення з (2.8) ,одержуємо
(2.9) |
Це і є одна з форм подання погрішності апроксимації. Її недоліком є та обставина, що для обчислення необхідне значення, що невідомо. Як вихід з такого положення залишається взяти лише його наближене значення, тобто, однак у цьому випадку співвідношення (2.9) стає вже наближеним.
Якщо ж припустити, що функціядосить гладка й має безперервні похідні до-го порядку включно, то з формули
,
справедливість якої треба по індукції, по теоремі про середній випливає співвідношення
,
де Тоді (2.9) приймає вид
(2.10)
|
Недоліком цього співвідношення є те, що невідомо значення . Однак, якщо відомо вид функції,те корисної може виявитися оцінка
У тому випадку, коли вузлові точки рівновіддалені й ,, те по індукції можна показати, що
(2.11)
Якщо до них додати точку не змінює характер розташування вузлів, то формула залишається справедливої й для розділеної різниці наступного порядку. Т.е.
.
У противному випадку, остання рівність стає наближеним, і зробивши в заміну, одержимо
(2.12)
|