
- •Міністерство освіти і науки України
- •Сумський державний університет
- •Конотопський інститут
- •Методичні вказівки
- •Урахування похибок Основні джерела похибок
- •Основні поняття
- •Правила обчислення похибок
- •Методи розв'язування нелінійних рівнянь
- •Метод половинного поділу (бісекцій або діхотомії)
- •Метод січних (хорд, пропорційних частин)
- •Метод Ньютона (дотичних)
- •Метод хорд і дотичних (комбінований метод)
- •Метод простих ітерацій
- •Методи розв'язування систем нелінійних рівнянь
- •Метод простих ітерацій
- •Метод Зейделя
- •Метод Ньютона
- •Модифікований метод Ньютона
- •Розв’язування систем лінійних рівнянь (слар)
- •Метод ітерації
- •Зведення лінійної системи алгебраїчних рівнянь до вигляду, який є зручним для ітерації
- •Метод Зейделя
- •Метод релаксації
- •Наближення функцій
- •Інтерполяція
- •Інтерполяційна формула Лагранжа
- •Оцінка похибки інтерполяційної формули Лагранжа
- •Збіжність функціонального інтерполяційного процесу для неперервних функцій
- •Методика розв’язування задач лінійної інтерполяції
- •Методика розв’язування задачі параболічної інтерполяції
- •Поліноми Чебишова
- •Інші методи інтерполяції. Інтерполяційний багаточлен Ньютона
- •Методи інтегрально-диференціальної інтерполяції
- •Методи інтегрального згладжування
- •Метод найменших квадратів (мнк)
- •Особливості мнк
- •Метод найкращого інтегрального наближення
- •Методи інтерполяції та згладжування на основі сплайнів
- •Інтерполяційні диференціальні кубічні сплайни
- •Метод прогонки
- •Чисельне диференціювання
- •Формули чисельного диференціювання на основі формули Стірлінга
- •Похибки при чисельному диференціюванні
- •Чисельні методи інтегрування функцій
- •Формули Ньютона-Котеса
- •Метод прямокутників
- •Метод трапецій
- •Метод Сімпсона (парабол)
- •Квадратична формула Гауса
- •Порівняння похибок квадратурних формул
- •Методи розв’язування звичайних диференціальних рівнянь та їх систем. Задача Коші
- •Стійкість (коректність) задачі Коші
- •Контрольна робота
Метод хорд і дотичних (комбінований метод)
Розглянемо
рівняння (1) і нехай у точці
виконується
умова
.
Застосуємо в цій точці метод дотичних,
а в точці
х0
=а –
метод хорд. Ітераційні формули
комбінованого
методу мають
вигляд:
(7)
Геометрична
інтерпретація методу. В точці
проводять дотичну до кривої
та
отримують наближення
,
а через точки
та
проводять хорду й отримують наближення
,
тобто на кожному наступному кроці метод
хорд застосовують до нового проміжку
.
Процес продовжують, доки не виконається
умова
<
.
Коренем
рівняння
(1)
буде
.
Метод простих ітерацій
Нехай
відомо, що корінь рівняння (1) лежить на
відрізку
.
Перетворимо рівняння (1) до вигляду
(8)
Таке перетворення може бути виконано різними способами, але для збіжності треба забезпечити виконання умови
<1 (9)
Метод
простих ітерацій або метод послідовних
наближень полягає у тому, що вибираємо
початкове наближення
кореня рівняння
(8),
де
й обчислимо перше наближення за формулою
,
а далі
.
Наступні наближення описує формула
,
(10)
Якщо
існує границя
,то
є коренем рівняння (8).
Теорема
1.
Нехай функція
визначена та диференційована на відрізку
(в області
),
причому всі її значення
.
Тоді, якщо існує правильний дріб
такий, що виконується нерівність (9) при
,
то:
процес ітерації (10) збіжний незалежно від початкового наближення
;
граничне значення
є коренем рівняння (8) на відрізку
(в
області
).
Зауваження
1.
Теорема залишається вірною, якщо
визначена і диференційована в нескінченому
інтервалі
,
причому
повинна задовольняти (9).
Зауваження
2.
В умовах теореми 1 ітерації збігаються
при будь-якому виборі
.
Окрема похибка в обчисленнях, яка не
виходить за межі проміжку
,
не впливає на кінцевий результат. Зростає
лише об'єм обчислень. Тому це надійний
метод обчислень.
Теорема
2.
Нехай функція
визначена і диференційована на деякому
відрізку
,
причому рівняння (8) має корінь, який
лежить у більш вузькому відрізку
,
де
;
.
Тоді, якщо виконується (9) і початкове
наближення
,
то:
всі послідовні наближення знаходяться в інтервалі
:
процес послідовних наближень збіжний, тобто існує
, причому
- єдиний корінь на відрізку
рівняння (8);
виконується оцінка (11).
(11)
Зауваження.
Нехай в деякому околі
кореня
рівняння (8) похідна
зберігає сталий знак і виконана нерівність
(9). Тоді, якщо похідна
додатна, послідовні наближення (10)
збігаються до кореня
монотонно. Якщо похідна
від’ємна, то послідовні наближення
коливаються біля кореня
.
Геометрична інтерпретація Геометрична інтерпретація
методу простих ітерацій методу простих ітерацій
для
випадку
для випадку
Геометрична інтерпретація Геометрична інтерпретація
методу простих ітерацій методу простих ітерацій
для
випадку >
для
випадку
>1
Приклад
1.
Знайти дійсні корені рівняння
з точністю до трьох значущих цифр.
Запишемо
.
Графічним способом встановлюємо, що
рівняння має на відрізку
один дійсний корінь
.
Дотримуючись визначень теореми 2,
задаємо:
і
.
Звідси
.
Так
як
і
,
то при
маємо:
<1.
Якщо
виберемо
,
то всі умови теореми 2 будуть виконані.
Виберемо
і граничну абсолютну похибку
4 і 5 наближення збігаються з точністю до 4 знаків. Тому:
Так
як гранична абсолютна похибка приблизного
кореня
,
включаючи похибку округлення, не
перевищує: 0,0005+0,0001=0,0006<0,0005, то можна
прийняти
.
Приклад
2.
Методом простих ітерацій з точністю
знайти корінь трансцендентного рівняння
на відрізку
.
Запишемо
рівняння у вигляді:
.
На відрізку
маємо
<1.
За початкове наближення обираємо
.
Ітераційний процес запишемо у вигляді:
Послідовно
знаходимо:
На
5 і 6 ітераціях виконується
<
,
тому
.
Звідси
бачимо, що збіжність двостороння і
лінійна, тому що відношення
приблизно однакові при
та рівні
(знаменник геометричної прогресії).
Зауваження
1.
Вихідне рівняння
можна записати у вигляді рівності
,
вибираючи різним способом
.
У деяких випадках
буде менший, а в інших – більший в околі
шуканого кореня
.
Для метода простої ітерації найкращим
є такий запис, для якого
<1,
причому для менших
швидкість збіжності до кореня
є більшою.
Зауваження
2.
Можна виразити
із рівняння
так, щоб для отриманого рівняння
виконувалась умова збіжності
<1
в околі шуканого кореня.
Зауваження
3.
Загальний прийом зведення (1) до (2), для
якого забезпечено виконання
<1.
Нехай
,
причому 0<
.
Замінюємо (1) еквівалентним
,
де
>0-
константа. Тоді:
.
Вибираємо λ, щоб в околі ξ виконувалась
нерівність:
<1.
Звідси
на основі попереднього:
Якщо
,
,
то
<1.
Приклад
3.
Знайти найбільший додатній корінь
рівняння
(1) з точністю
.
Інтервал знаходження кореня
.
Рівняння (1) можна записати у вигляді:
тощо.
Остання
формула є найбільш вдалою, тому що,
взявши інтервал (9;10) і визначивши
матимемо:
,
тому
.
Обчислюємо послідовні наближення
з одним запасним знаком за формулами:
Звідси:
Приклад
4.
Знайти методом простої ітерації корінь
рівняння
з точністю 0,01.
Зауваження
2.
Виконання умови
не
гарантує наближеності до точного
розв’язку.
Приклад
5.
Для
розв'язування рівняння х2
=а можна
прийняти
або
і, відповідно, записати такі ітераційні
процеси:
або
.
Перший
процес взагалі не збігається, а другий
збігається для будь-якого
.
Другий процес збігається дуже швидко,
бо
.
Метод
простих ітерацій має просту геометричну
інтерпретацію. Побудуємо графіки функцій
і
.
Коренем
рівняння (8) є абсциса точки перетину
графіків. Від початкового наближення
х0
будуємо
ламану, абсциси вершин якої є послідовними
наближеннями
кореня
.
На рис. 5 показано випадок, коли
,
а
на рис. 6 — випадок для
Оцінка похибки методу. Оцінимо похибку n-го наближення:
звідки, звівши подібні члени, отримаємо
Якщо
то
;
тоді оцінка похибки наближенняхп
зводиться
до оцінки модуля різниці двох послідовних
наближень.
Зауваження 1. Особливість методу простих ітерацій ― ненакопичення похибки обчислень. Похибка обчислень може вплинути на кількість ітерацій, однак не на точність кінцевого результату.
Зауваження 2. Метод простих ітерацій має лінійну швидкість збіжності.