- •Методичні вказівки до практичних робіт
- •Передмова
- •Практична робота № 1
- •Теоретичні відомості Похибки величин
- •Визначення 1.2. Відносною похибкою (зазвичай визначається в %) називають величину , таку, що:
- •Похибка функції
- •Завдання
- •Контрольні питання
- •Практична робота № 2
- •Теоретичні відомості
- •Завдання
- •Контрольні питання
- •Практична робота № 3
- •Теоретичні відомості
- •Метод хорд
- •Метод Ньютона (дотичних)
- •Метод простої ітерації
- •Граничні оцінки та область існування коренів алгебраїчних рівнянь
- •Метод Лагранжа
- •Метод Ньютона
- •Метод кільця
- •Завдання
- •Контрольні питання
- •Практична робота № 4
- •Теоретичні відомості
- •Багаточлен Лагранжа
- •Інтерполяція за методом Ньютона
- •Інтерполяція сплайнами
- •Завдання
- •Контрольні питання
- •Практична робота № 5
- •Теоретичні відомості
- •Побудова лінійної емпіричної формули
- •Побудова квадратичної емпіричної формули
- •Побудова емпіричної формули найпростіших нелінійних залежностей
- •Завдання
- •Контрольні питання
- •Практична робота № 6
- •Теоретичні відомості
- •Метод прямокутників
- •Метод трапецій
- •Метод Сімпсона
- •Метод Ньютона-Котеса
- •Методи Чебишева і Гауса
- •Завдання
- •Контрольні питання
- •Практична робота № 7
- •Теоретичні відомості Метод Ейлера і модифікований метод Ейлера
- •Метод Рунге-Кутта
- •Метод Мілна
- •Завдання
- •Контрольні питання
- •Практична робота № 8
- •Теоретичні відомості
- •Метод стрільби
- •Кінцево-різницевий метод
- •Методи вирішення задач на власні значення та власні вектори матриць
- •Метод обертань Якобі чисельного вирішення задач на власні значення та власні вектори матриць
- •Завдання
- •Контрольні питання
- •Практична робота № 9
- •Теоретичні відомості
- •Метод сканування
- •Метод золотого перерізу
- •Метод параболічної апроксимації
- •Завдання
- •Контрольні питання
- •Література
Метод простої ітерації
Метод простої ітерації реалізує такий підхід: попередньо початкове рівняння перетворюють до вигляду , що є окремим випадком більш загальної структури . Потім обирають початкове значення і підставляють його в ліву частину рівняння, але , тому що обрано довільно і не є коренем рівняння. Отримане розглядають як чергове наближення до кореня. Його знову підставляють в ліву частину рівняння і отримують наступне значення і т. д., в загальному випадку . Отримана таким чином послідовність при певних умовах може збігатися до кореня
Існують різні способи перетворень рівняння до вигляду ; одні можуть призвести до виконання умови збіжності завжди, інші – в окремих випадках. Найбільш простий спосіб наступний:
,
але він не завжди приводить до успіху. Існує інший спосіб, у відповідності з яким , причому треба обирати так, щоб , де і знак збігався за знаком на .
Похибку рішення можна оцінити з співвідношення:
, де .
У якості критерію закінчення обчислень у методі ітерацій застосовують співвідношення , де - задана похибка рішення.
Часто використовують спрощену умову закінчення пошуку рішення , не обчислюючи максимального значення похідної. У цьому випадку похибка рішення може не відповідати заданій, тобто бути більшою або меншою.
Приклад 3.5
Маємо рівняння . Уточнити корінь з похибкою .
Запишемо . Провівши процедуру відділення коренів, одержимо, що корінь знаходиться в проміжку [0; 0,05], тобто . Приведемо рівняння до виду, зручного для проведення ітерацій . Функцію знаходимо із співвідношення , вважаючи для підвищення збіжності, що , де ; число має той же знак, що і у проміжку [0; 0,05].
Знаходимо:
при .
Приймемо ,
тоді .
За початкове наближення візьмемо , всі інші наближення визначаємо з рівності . Результати зведемо у таблицю 3.5.
Таблиця 3.5
|
|
|
0 |
0 |
0,2614 |
1 |
0,2614 |
0,2266 |
2 |
0,2266 |
0,2309 |
3 |
0,2309 |
0,2303 |
4 |
0,2303 |
0,2304 |
5 |
0,2304 |
|
Відповідь: .
Граничні оцінки та область існування коренів алгебраїчних рівнянь
У багатьох випадках досить знати граничне значення кореня. Це необхідно у різних додатках (наприклад, при оцінці стійкості систем керування), а також для більш ефективного відділення коренів нелінійного рівняння. Цю проблему можна вирішити, не знаходячи всі корені, а скориставшись спеціальними методами.
Метод Лагранжа
Метод зводиться до визначення верхньої границі додатних коренів за формулою:
, (3.1)
де - номер першого за порядком від’ємного коефіцієнта в поліномі лівої частини рівняння ; - найбільша з абсолютних величин від’ємних коефіцієнтів ; при цьому передбачається, що .
Приклад 3.6
Визначити граничне значення додатних коренів рівняння:
.
У цьому рівнянні , тому спочатку помножимо обидві частини рівняння на -1. Одержимо:
.
Тут ; ; ; ; .
Отже, ; .
Застосовуючи формулу (3.1), одержимо:
.