- •Нехай задане рівняння
- •Методи уточнення коренів нелінійних рівнянь Метод поділу проміжку навпіл
- •Значення е задається в межах 10 –410 –6.
- •Метод хорд
- •Комбінований метод хорд та дотичних
- •Метод простої ітерації
- •Метод Ейткена – Стефенсона
- •Метод Стефенсона
- •Методи розв’язування систем лінійних алгебраїчних рівнянь
- •Прямі методи розв’язування систем лінійних алгебраїчних рівнянь
- •Перехід від початкової системи
- •Можна записати, що для всіх
- •Метод Гаусса (рядкова і стовпцева форми розкладу матриці а до трикутного вигляду). Розглянемо систему лінійних рівнянь
- •Метод Гаусса з вибором головного елемента.
- •Ітераційні методи розв’язування систем лінійних алгебраїчних рівнянь
- •Нехай дано лінійну систему
- •Умови збіжності ітераційного процесу Нехай задано зведена до нормального вигляду система лінійних рівнянь
- •Для системи
- •Метод Зейделя
Комбінований метод хорд та дотичних
Метод хорд та дотичних дають наближення кореня з різних сторін (менше і більше від істинного значення). Тому доцільно використовувати обидва способи одночасно, завдяки чому уточнене значення кореня одержується швидше.
Нехай
–
початкове наближення кореня за методом
хорд, а
– за методом дотичних (див.рис.4).
Тоді провівши хорду та дотичну, одержимо відповідні наближення за методом хорд
і за методом дотичних
.
Або в загальному випадку
(11)
(12)
Рис. 4
Якщо припустима абсолютна похибка ε заздалегідь задана, то процес наближення припиняється, доки не буде виявлено, що
Після закінчення процесу за значення кореня х* краще взяти середнє арифметичне одержаних останніх значень
Кращий результат дає наступний порядок обчислень:
Знаходиться наближене значення кореня за методом Ньютона;
Знаходиться наближене значення кореня за методом хорд, використовуючи замість
значення
,
знайдене за методом Ньютона, і процес
повторюється до одержання бажаної
похибки обчислень.
;
.
Рис.5.
Алгоритм комбінованого методу метод хорд та дотичних
Метод простої ітерації
Рис.6.
У цьому методі рівняння заміняється еквівалентним йому рівнянням
(13)
Наприклад, рівняння
зводимо до виду
.
Виберемо за
початкове наближення кореня значення
і підставимо
в праву частину рівняння (1). Одержимо деяке число
(14)
Підставляючи в праву частину рівності (2) замість значення одержимо нове число
Повторюючи процес, будемо мати наступну послідовність
(15)
Якщо ця послідовність
збіжна, то границя цієї послідовності
– корінь рівняння
і може бути обчислений з будь-якою
точністю.
Достатня умова
збіжності методу простої ітерації
формулюється наступним чином: якщо для
всіх
виконується нерівність
(16)
то на проміжку
рівняння
має єдиний корінь і процес ітерації
збігається до цього кореня незалежно
від вибору початкового наближення
.
Таким чином при
практичному знаходженні кореня за
методом ітерації при переході від
рівняння
до (13)
необхідно зобразити
так, щоб похідна
за
абсолютним значенням була якомога менша
одиниці.
Для зведення рівняння до вигляду (13) може бути застосований загальний метод, котрий забезпечує виконання нерівності (16).
Нехай
,
при
,
де m1
– найменше значення похідної
,
;
М1
– найбільше значення похідної на
відрізку [a,
b],
Якщо похідна – від’ємна, то замість рівняння розглянемо рівняння – .
Замінимо це рівняння
еквівалентним йому рівнянням
і виберемо сталу λ
так, щоб забезпечити виконання умови
(16)
.
1)
Розкриваємо
нерівність
Візьмемо праву
нерівність
,
з неї випливає, що
тобто
оскільки
З лівої нерівності
випливає, що
Отже, значення коефіцієнта λ
знаходиться в межах
.
Як правило за λ
приймають значення
де М1
– максимальне значення похідної
на
проміжку
.
Відповідно, ітераційна формула буде мати вигляд
2) Якщо
то можна довести, що
(17)
І відповідний ітераційний процес має вигляд
(18)
Алгоритм методу простої ітерації
