
- •«Чисельні методи комп’ютерного аналізу»
- •6.050201 «Системна інженерія», спеціалізація:
- •6.050903 «Телекомунікації», спеціалізація:
- •Лабораторна робота №1
- •1.1 Теоретичні відомості
- •1.1.1 Абсолютна і відносна похибки
- •1.1.2 Десятковий запис наближених чисел. Значуща цифра числа. Кількість вірних знаків
- •1.1.3 Зв'язок між кількістю вірних значущих цифр і похибкою числа
- •1.1.4 Пряма і зворотна задачі теорії похибок
- •1.1.4.1 Пряма задача теорії похибок
- •1.1.4.2 Зворотна задача теорії похибок
- •1.2 Завдання на виконання лабораторної роботи
- •1.4 Контрольні питання.
- •Література
- •Лабораторна робота №2
- •2.1 Теоретичні відомості
- •2.1.1 Кінцеві різниці n-х порядків
- •2.1.1.1 Таблиці кінцевих різностей
- •2.1.2 Постановка задачі інтерполяції
- •2.1.2.1 Інтерполяційні формули Ньютона
- •2.1.2.2 Формула Лагранжа
- •2.2 Завдання на проведення лабораторної роботи
- •2.4 Контрольні питання
- •Л абораторна робота №3
- •3.1 Загальні відомості
- •3.1.1 Постановка задачі
- •3.1.2 Методи розв’язку задачі
- •3.1.3 Формули наближеного інтегрування
- •3.1.3.1 Формула прямокутників
- •3.1.3.2 Формула трапецій
- •3.1.3.3 Формула Сімпсона (формула парабол)
- •3.1.3.4 Формула Гауса
- •3.2 Завдання на проведення лабораторної роботи
- •3.4 Контрольні питання
- •Література
- •Лабораторна робота №4
- •4.1 Теоретичні відомості
- •4.1.1 Постановка задачі
- •4.1.2 Відділення коренів. Теорема про оцінку помилки наближеного значення кореня
- •4.1.3 Уточнення кореня методом розподілу відрізка навпіл
- •4.1.4 Метод ітерації
- •4.1.5 Метод Ньютона і його модифікації
- •4.1.6 Метод хорд
- •4.1.7 Комбінований метод дотичних і хорд
- •4.2 Завдання на проведення лабораторної роботи
- •4.4 Контрольні питання
- •Лабораторна робота №5
- •5.1 Теоретичні відомості
- •5.1.1 Постановка задачі
- •5.1.2. Методи розв’язку
- •5.1.2.1. Метод Ейлера-Коші
- •5.1.2.2. Метод Ейлера-Коші з ітераціями
- •5.1.2.3 Модифікований метод Ейлера
- •5.1.2.4. Метод Рунге-Кута
- •5.1.2.5. Явні методи Адамса
- •5.2 Завдання на проведення лабораторної роботи
- •5.4 Контрольні питання
4.1.3 Уточнення кореня методом розподілу відрізка навпіл
Метод розподілу відрізка навпіл (або метод дихотомії) застосовується для уточнення кореня рівняння з наперед заданою точністю, якщо функція задовольняє умові теореми 2.
За початкове наближення вибираємо середину відрізка .
(4.2)
Проводимо
дослідження значення функції на кінцях
відрізків
і
.
Корінь, який шукаємо, знаходиться в тім
відрізку, на кінцях якого функція
здобуває значення протилежних знаків.
За нове наближення вибираємо середину
нового відрізка.
Ітераційний процес продовжується поки не буде досягнуте виконання умов:
(4.3)
Приклад 4.1 Уточнити корінь рівняння
,
який належить відрізку
.
Рішення.
Послідовно маємо:
Можна прийняти
4.1.4 Метод ітерації
Нехай задане рівняння , де – неперервна функція. Щоб знайти дійсні корені цього рівняння, замінимо його рівнозначним:
.
(4.4)
Тепер, щоб рішити рівняння (4.4) застосовують метод послідовних наближень (метод ітерації). Вибирають деяке початкове наближення і послідовно обчислюють наступні наближення:
,
(4.5)
Збіжність
послідовності
забезпечується відповідним вибором
функції
і початкового наближення
.
Вибираючи по-різному функцію
,
можна скласти різні ітераційні методи
рішення рівняння (4.4).
Рисунок 4.1 – Геометрична інтерпретація методу ітерації
а)
;
б)
;
в)
;
г)
.
Метод
ітерації має простий геометричний
зміст. Побудуємо графіки функцій
та
.
Абсциса точки перетину графіків цих
функцій є коренем рівняння (4.4) (рис.
4.1).
На
відрізку
довільно вибираємо точку
і проводимо через неї пряму, рівнобіжну
осі ординат до перетинання з кривою
у точці
.
З точки
проведемо пряму, рівнобіжну вісі абсцис,
до перетинання з прямою
.
У результаті одержимо точку
з ординатою
.
Спроектувавши точку
на вісь
,
знаходимо абсцису
.
Аналогічно через
проводимо пряму, рівнобіжну осі ординат,
до перетинання з кривою
у точці
.
З точки
проводимо пряму, рівнобіжну вісі абсцис,
до перетинання з прямою
у точці
,
абсциса якої
і т.д. У результаті спільні абсциси
точок
і
,
і
є послідовними наближеннями
до
кореня
.
Якщо
,
то послідовні наближення
сходяться до кореня
,
якщо
,
те послідовні наближення віддаляються
від нього.
Швидкість
збіжності залежить також від вибору
початкового наближення
.
Чим ближче до кореня
обране
,
тим швидше буде знайдений результат.
Процес ітерацій можна закінчити, якщо буде виконана наступна нерівність
.
4.1.5 Метод Ньютона і його модифікації
Будемо
вважати, що рівняння
на відрізку
має ізольований корінь
,
тобто
,
а функції
і
неперервні і зберігають свій знак на
.
Нехай
–
-ті
наближення кореня. Розкладемо
у ряд Тейлора в околиці точки
.
Замість рівняння розглянемо рівняння
,
яке
враховує тільки лінійну відносно
частину ряду Тейлора. Вирішивши його
відносно
,
одержимо
.
Узявши знайдене значення за наступне наближення, одержимо
,
(4.6)
Формула
(4.6) визначає метод Ньютона. Він має
просту геометричну інтерпретацію.
Значення,
є абсцисою точки перетинання дотичної
з кривою
у точці
(рис. 4.2). Тому метод Ньютона називають
ще методом дотичних. З рис. 4.2 видно, що
послідовні наближення сходяться до
кореня
монотонно.
Рисунок 4.2 – Геометрична інтерпретація методу Ньютона
За початкове наближення в методі Ньютона слід брати точку
,
у якій
.
Достатні умови збіжності методу Ньютона дає така теорема.
Теорема
3.
Нехай на відрізку
функція
має неперервні зі стійкими знаками
похідні
,
і
.
Тоді існує таке коло
кореня
рівняння
,
що для будь-якого
послідовність
,
обчислена по формулі (5.6), сходиться до
кореня
.
Перевага методу Ньютона перед методом ітерації в тому, що він має більш високу швидкість збіжності.
Недоліком методу Ньютона є те, що на кожній ітерації потрібно обчислювати не тільки значення функції , але і значення її похідної .
Приклад
4.2
Обчислити методом Ньютона від’ємний
корінь рівняння
з п’ятьма вірними знаками.
Рішення.
Вважаючи
в лівій частині рівняння
,
отримаємо
,
,
.
Отже,
корінь
,
який ми шукаємо, знаходиться в інтервалі
.
Звузимо знайдений інтервал. Оскільки
,
то
.
В цьому останньому інтервалі
і
.
Оскільки
і
,
то можемо прийняти за початкове
наближення
.
Послідовні наближення
обчислимо за наступною схемою:
|
|
|
|
|
0 1 2 3 |
-11 -10,3 -10,27 -10,261 |
3453 134,3 37,8 0,2 |
-5183 -4234 -4196 - |
0,7 0,03 0,009 - |
Зупиняючись
на
,
перевіряємо знак значення
.
Оскільки
,
то
,
і будь-яке з цих чисел дає необхідне
наближення.