- •1 Основні поняття теорії похибок
- •1.1 Поняття похибки
- •1.2 Дії над наближеними числами
- •1.3 Пряма та обернена задачі теорії похибок
- •1.3 Джерела похибок обчислень
- •2 Обчислення значень функцій
- •2.1 Обчислення значень полінома. Схема Горнера
- •2.2 Наближене знаходження сум числових рядів
- •2.3 Обчислення значень аналітичної функції
- •2.4 Обчислення значень показової функції
- •2.5 Обчислення значень логарифмічної функції
- •2.6 Обчислення значень тригонометричних функцій
- •3 Розв’язування систем лінійних алгебраїчних рівнянь (слар)
- •3.1 Концепція методів
- •3.2 Метод простої ітерації
- •3.3 Метод Гаусса-Зейделя
- •4 Розв’язування нелінійних алгебраїчних та трансцендентних рівнянь
- •1. Етап відділення кореня: на цьому етапі відділяється корінь, тобто знаходиться такий відрізок, усередині якого міститься точно один корінь і з цього відрізка береться початкове наближення кореня.
- •2. Етап уточнення кореня: на цьому етапі послідовно уточнюють корінь, тобто знаходять значення із заданою точністю .
- •4.1 Відділення коренів
- •4.2 Метод половинного поділу (метод дихотомії)
- •4.3 Метод хорд (спосіб пропорційних частин)
- •4.4 Метод Ньютона (метод дотичних)
- •4.5 Комбінований метод
- •4.6 Метод ітерації
- •4.7 Метод ітерації для системи двох рівнянь
- •5 Обробка емпіричних даних
- •5.1 Інтерполяція та екстраполяція
- •5.2 Концепція інтерполяції та екстраполяції
- •5.3 Лінійна і квадратична локальні інтерполяції
- •5.4 Глобальна інтерполяція. Многочлен Лагранжа
- •5.5 Глобальна інтерполяція. Многочлен Ньютона
- •5.6 Апроксимація
- •При цьому вимагається, щоб
- •Побудуємо емпіричний точковий графік.
- •Візуальний аналіз побудови дозволяє обрати на роль апроксимуючої функції квадратичну параболу
- •6 Наближене обчислення визначених інтегралів
- •6.1 Концепція чисельного інтегрування
- •6.2 Методи прямокутників та трапецій
- •6.3 Метод Симпсона
- •6.4 Метод Монте-Карло
4.7 Метод ітерації для системи двох рівнянь
Нехай задано систему двох рівнянь з двома невідомими
(4.25)
дійсні корені якої потрібно знайти з заданою точністю.
Припустимо, що система (4.25) має лише ізольовані корені. Кількість цих коренів та їх грубі наближені значення можна встановити, побудувавши криві та , та визначивши наближено координати їх точок перетину.
Нехай - наближені значення коренів системи (4.25), знайдені графічно чи яким-небудь іншим способом (наприклад, грубою прикидкою).
Задамо ітераційний процес, який дозволяє при відомих умовах уточнити дані наближені значення коренів. Для цього представимо систему (4.25) у вигляді
(4.26)
та побудуємо послідовні наближення за такими формулами:
(4.27)
Якщо ітераційний процес збігається, тобто існують границі та , то ці границі та є коренями системи (4.25). Отже, взявши достатньо велике число ітерацій (4.27), ми отримаємо числа та , які будуть відрізнятися від точних коренів та системи (4.25) як завгодно мало. Таким чином, задача буде розв’язаною. Якщо ітераційний процес (4.27) розбіжний, то його не можна використовувати.
Умови збіжності ітераційного процесу (4.27) задає наступна теорема.
Теорема 1 (без доведення). Нехай у деякій замкнутій області (рис. 4.11) існує одна й лише одна пара коренів та системи (4.26).
Якщо:
функції , визначені, неперервно диференційовані у ;
початкові наближення та всі наступні наближення належать ;
у виконані нерівності
то процес послідовних наближень (4.27) збігається до коренів та системи (4.26), тобто та .
Зауваження. Теорема залишається вірною, якщо умову (4.27) замінити умовою
Рис. 4.11 Рис. 4.12
Приклад. Для системи
знайти додатні корені з чотирма значущими цифрами.
Розв’язок. Будуємо графіки функцій та (рис. 4.12).
Наближені значення шуканих коренів .
Для застосування методу ітерації запишемо задану систему рівнянь у наступному вигляді:
Знаходимо частинні похідні:
Тут .
Обмежуючись областю , будемо мати:
Звідси
Отже, якщо послідовні наближення не вийдуть за межі області (що легко встановити у процесі обчислень), то ітераційний процес буде збіжним.
Обчислюємо послідовні наближення за формулами:
Результати обчислень відображені у таблиці:
|
|
|
0 |
3,5 |
2,2 |
1 |
3,479 |
2,259 |
2 |
3,481 |
2,260 |
3 |
3,484 |
2,261 |
4 |
3,486 |
2,261 |
5 |
3,487 |
2,262 |
6 |
3,487 |
2,262 |
Таким чином, можна прийняти
Зауваження. Замість розглянутого процесу послідовних наближень (4.27) інколи зручно користуватися «процесом Зейделя»:
5 Обробка емпіричних даних
В цьому розділі розглянемо тільки основні напрями обробки даних: інтерполяцію, екстраполяцію, апроксимацію, які є базою для вирішення інших задач обробки табличних даних, наприклад, задач згладжування, чисельних методів інтеграції і т.д.