
- •Лабораторна робота №1 Основи програмування в інтегрованій системі MatLab. Розв’язання функціональних рівнянь з однією змінною
- •1.1. Короткі теоретичні відомості по програмуванню в інтегрованій системі автоматизації математичних розрахунків MatLab.
- •2.1. Розв‘язання функціональних рівнянь з однією змінною.
- •2.1.1. Постановка задачі
- •2.1.2. Етап 1: відокремлення коренів
- •2.1.3. Етап 2: уточнення коренів.
- •2.1.3.1. Метод половинного ділення.
- •2.1.3.2. Метод хорд.
- •2.1.3.3. Метод січних
- •2.1.3.4. Метод дотичних (Ньютона)
- •2.1.3.5 Метод простої ітерації
- •1.3. Завдання для виконання лабораторної роботи:
- •Варіанти завдань до лабораторної роботи №1.
- •Міністерство освіти та науки україни
- •Перелік основних програмних продуктів класу matlab.
- •Програмування в системі matlab.
- •Основні типи даних.
- •1 2 3
- •Оператори відношення та їх функції
- •Логічні оператори
- •Рекомендована література до вивчення курсу „Чисельні методи”.
2.1.2. Етап 1: відокремлення коренів
Відокремлення коренів – вста-новлення „тісних” проміжків, кожен з яких вміщує тільки один корінь. Най-більш наочно про демонструвати цей етап можливо графічно. Враховуючи, що дійсні корені рівняння (1.1) – це точки перетину графіка функції F(х) з віссю абсцис, досить побудувати графік F(х) і відмітити на осі OX відрізки, кожен з яких вміщує по одному кореню.
Рис. 1.1. Графічний спосіб відокремлення коренів.
Побудову графіка вдається значно спростити, якщо замінити рівняння (1.1) рівносильним йому рівнянням F1(х)= F2(х). В цьому випадку будують графіки функцій F1(х) та F2(х), потім на вісі OX відмічають відрізки, що вміщують абсциси точок перетину цих графіків.
Рис.
1.2. Варіанти розв’язання рівняння при
заміні функції рівносильними функціями
F1(х)
та F2(х)
на проміжку [ 1, 3
].
Зустрічаються випадки, коли не тільки на першому етапі – при виділенні проміжків, що вміщують корені рівняння, а і при подальшому розв’язанні рівняння розрахунки проводять не за основною функцією F1(x), а за наближеною до неї функцією F2(x) (рис. 1.3).
Рис. 1.3. Реальна функція F1(x) та її наближений аналог-функція F2(x).
При написанні програми, яка виконує відокремлення коренів графічний спосіб не є зручним, тому застосовують найпростіший алгоритм відокремлення коренів: проміжок, на якому відшукуються розв’язки, ділять рівномірно на невеликі відрізки, на кінцях кожного відрізку розраховують значення функції, якщо функція на кінцях відрізку змінює знак – відрізок містить корінь. Величина відрізків поділу залежить від коефіцієнта гладкості функції.
Для відділення коренів можемо скористатися теоремою:
Якщо функція у = F(х) неперервна на інтервалі [a,b] і якщо F(а) та F(b) мають протилежні знаки, тобто F(а) F(b) < 0, то F(х) має хоча б один дійсний корінь на інтервалі [a,b]. Якщо при цьому F(х) має першу похідну, що не змінює знак на інтервалі [a,b], то корінь єдиний.
2.1.3. Етап 2: уточнення коренів.
2.1.3.1. Метод половинного ділення.
Припустимо, що
рівняння
(1.1)
має на відрізку [a, b]
єдиний корінь. Функція F(x)
на цьому відрізку неперервна. Поділимо
відрізок [a,
b]
точкою
навпіл (рис. 1.4). Якщо F(c) 0,
то вибираємо відрізок [a, с]
чи [с, b],
на якому функція змінює знак, тобто
F(а)
F(с)
< 0 чи
F(с)
F(b)
< 0. Для
подальших поділів використовуємо цей
відрізок. На кожному кроці процесу
половинного ділення відрізок зменшується
вдвоє.
Рис. 1.4. Знаходження кореня рівняння методом половинного ділення.
Ітераційний процес продовжуємо, поки довжина відрізка не стане менше заданої точності :
|a – b| < (1.6)
2.1.3.2. Метод хорд.
Даний метод ґрунтується на лінійній інтерполяції функції F(x)=0 по двом значенням, що мають протилежні по знаку значення функції F(a) та F(b). Метод хорд швидше збігається до розв‘язку навіть при досить малих значеннях . Потрібно знайти корінь рівняння F(x)=0 на проміжку [a; b], і якщо відомо, що F(x) неперервна на [a; b] і F(a)∙F(b) < 0. Крім того, перша F'(x) і друга F''(x) похідні функції F(x) зберігають на проміжку [a; b] свій знак. Замінимо функцію F(x) лінійною функцією, яка проходить через вузлові точки (a , F(a)) і (b , F(b)) :
(1.7)
Лінійна функція P(x) на кінцях відрізку [a; b] приймає такі ж значення, як і функція F(x)=0.
Рис. 1.5. Знаходження кореня за методом хорд.
В якості першого наближення при знаходженні кореня функції F(x)=0 візьмемо точне значення кореня функції P(x)=0, тобто х1, яке розрахуємо з рівняння:
(1.8)
При подальшому
дослідженні відрізків [a; х1]
та [х1; b],
виберемо той, на якому функція змінює
знак. З рис. 1.5 бачимо, що таким відрізком
є [a; х1]
. Для вибраного відрізка побудуємо
лінійне наближення функції за формулою
1.7, виконаємо розрахунки кореня для
лінійного наближення за формулою 1.8.
Розрахунки припинимо, коли
.