Наближене обчислення дійсних коренів многочлена
Аналіз чисельних методів розв'язування рівнянь з однією змінною. Зміст теореми про оцінку похибки наближеного значення кореня. Уточнення степеню концепцією поділу відрізка пополам. Характеристика комбінованого способу дотичних і хорд та простої ітерації.
Рубрика | Математика |
Вид | курсовая работа |
Язык | украинский |
Дата добавления | 26.07.2015 |
Размер файла | 570,9 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
0
Курсова робота
на тему Наближене обчислення дійсних коренів многочленна
ЗМІСТ
ВСТУП
РОЗДІЛ 1. ЧИСЕЛЬНІ МЕТОДИ РОЗВ'ЯЗУВАННЯ РІВНЯНЬ 3 ОДНІЄЮ ЗМІННОЮ
1.1 Постановка задачі
1.2 Відокремлення коренів. Теорема про оцінку похибки наближеного значення кореня
1.3 Уточнення кореня методом поділу відрізка пополам
1.4 Метод дотичних
1.5 Метод хорд
1.6 Комбінований метод дотичних і хорд
1.7 Метод простої ітерації
РОЗДІЛ 2. ПРАКТИЧНА ЧАСТИНА
ВИСНОВКИ
СПИСОК ВИКОРИСТАНИХ ДЖЕРЕЛ
ВСТУП
Центральним в алгебрі многочленів виявляється питання не про практичний пошук коренів многочлена, а питання про їх існування. Відомо, що існують квадратні рівняння з дійсними коефіцієнтами, які не мають дійсних коренів. А чи не знайдеться таке рівняння п'ятого чи більш високого степеня, яке не має жодного кореня навіть серед комплексних чисел, і чи не доведеться для пошуку коренів даних рівнянь переходити від комплексних чисел до більш широкого запасу чисел? Відповідь на це питання дає важлива теорема (основна теорема алгебри), яка стверджує, що будь-яке рівняння з будь-якими числовими коефіцієнтами, не лише з дійсними, але і з комплексними, має хоча б один комплексний корінь.
Твердження, близьке до основної теореми алгебри висловили ще в 17 столітті Жерар (1629) і Декарт (1637). Воно полягало в тому, що кожне алгебраїчне рівняння -го степеня має коренів; якщо ж дійсних коренів менше , при цьому решту коренів слід вважати «уявними». При цьому термін «уявний» не збігався з сучасним поняттям комплексного числа, а означав просто, що потрібну кількість коренів можна собі уявити існуючою.
Більш обережно формулював основну теорему алгебри Ньютон (1707), але у 18 столітті Ейлер чітко сформулював основну теорему алгебри.
Прийнято вважати, що перше строге доведення основної теореми дав Гаусс у 1799. Однак з точки сучасного вчення про неперервність це доведення вимагає деяких доповнень і, крім того, воно стосується лише многочленів з дійсними коефіцієнтами. Але основний хід цього доведення цілком правильний.
Метою курсової роботи є розширення уявлення про наближене обчислення дійсних коренів многочлена.
РОЗДІЛ 1. ЧИСЕЛЬНІ МЕТОДИ РОЗВ'ЯЗУВАННЯ РІВНЯНЬ 3 ОДНІЄЮ ЗМІННОЮ
1.1 Постановка задачі
Нехай задано рівняння з однією змінною
F(x)=0,
де функція F(x) визначена і неперервна на деякому проміжку [а; b].
Розв'язати рівняння означає знайти множину його коренів, тобто таких значень xє[а; b], при яких рівняння (1.1) перетвориться в тотожність. Корінь рівняння (1.1) називають ще нулем функції F(x). Якщо функція - алгебраїчний многочлен, то рівняння (1.1) називається алгебраїчним. Якщо функція F(x) містить тригонометричні, показникові або логарифмічні функції, тоді рівняння (1) називають трансцендентним.
Знайти точні значення коренів заданого рівняння можна лише для найпростіших функцій F(x): алгебраїчних многочленів не вище четвертого степеня, деяких многочленів степеня n ? 5 і деяких трансцендентних функцій.
Універсальних методів для знаходження точних значень коренів алгебраїчних рівнянь степеня п = 5 і трансцендентних рівнянь не існує. Крім того, розв'язуючи практичні задачі, часто дістають рівняння з коефіцієнтами, які є наближеними числами. Тоді постановка задачі знаходження точних коренів не має смислу. Тому важливого значення набувають наближені методи знаходження коренів рівняння з достатньою для практики точністю. Задача знаходження коренів рівняння (1.1) вважається розв'язаною, якщо корені обчислені із наперед заданою точністю.
Нехай х* - точний корінь, а - його наближене значення. Кажуть, що корінь обчислено з наперед заданою точністю є, якщо | х*-| < е. Нехай, наприклад, х*є[а; b] і b-a < е, тоді числа а і b - наближені значення кореня х* відповідно з недостачею і надлишком з точністю є. У цьому випадку за наближене значення з точністю є можна взяти будь-яке число з відрізка [а; b].
Знаходження наближених коренів рівняння (1) складається з двох етапів:
1) відокремлення коренів, тобто знаходження досить малих відрізків, на кожному з яких міститься один і тільки один корінь рівняння;
2) обчислення коренів з наперед заданою точністю. Перший етап називають ще задачею визначення відрізків ізоляції коренів, а другий - уточненням наближених коренів. Перший етап складніший за другий, оскільки для загального випадку немає досить ефективних методів відокремлення коренів. Для знаходження коренів з наперед заданою точністю застосовують методи, які дають можливість уточнювати знайдені наближення коренів.
Корені рівняння (1) можуть бути дійсними і комплексними. Далі розглянуто наближені методи обчислення тільки дійсних коренів рівняння (1.1).
1.2 Відокремлення коренів. Теорема про оцінку похибки наближеного значення кореня
Корінь x* рівняння (1.1) вважається відокремленим на відрізку [а; b], якщо х*є[а;b] і на цьому відрізку дане рівняння не має інших коренів. Щоб відокремити корені рівняння (1.1), треба розбити область визначення даного рівняння на проміжки, на кожному з яких міститься один і тільки один корінь або немає жодного кореня. Відокремлюють корені графічним і аналітичним методами, а також методом послідовного перебору.
Для відокремлення коренів графічним методом будують графік функції у=F(х) і знаходять точки перетину графіка з віссю абсцис та кінці відрізків ізоляції коренів.
Аналітичний метод відокремлення коренів ґрунтується на теоремах з курсу математичного аналізу. Сформулюємо їх.
Теорема 1 (теорема існування кореня). Якщо функція неперервна на [а;b] і набуває на кінцях цього відрізка значень протилежних знаків, тобто F(a)F(b)<0, то всередині відрізка [а;b] існує хоча б один корінь рівняння F(х)=0.
Зазначимо, що теорема не дає відповіді на питання про кількість коренів рівняння (1.1), які належать [a; b]. При виконанні умов теореми рівняння може мати й кілька коренів. На рис. 1.1 зображено графік функції у = F(х), яка задовольняє усі вимоги теореми 1 і має на [a;b] чотири нулі. У досить малому околі точки xз теорему існування кореня застосувати не можна, бо при переході зліва направо через точку x3 знак функції F(х) не змінюється. Точка хз - кратний корінь рівняння (1.1) і його не можна відокремити, користуючись теоремою 1. Тому далі вважатимемо, що F'(х)? 0 для всіх х є [а; b].
Размещено на http://www.allbest.ru/
0
Теорема 2 (теорема існування і єдиності кореня). Якщо функція F(x), неперервна і диференційована на [а; b], набуває на кінцях цього відрізка значень різних знаків, а похідна F'(x) зберігає сталий знак всередині відрізка [а; b], то рівняння F(х) = 0 на цьому відрізку має корінь, причому єдиний.
У відповідності з теоремами 1 і 2 алгоритм відокремлення коренів рівняння (1) можна сформулювати так:
1. Знайти область визначення рівняння.
2. Знайти критичні точки функції F(х).
3. Записати інтервали монотонності функції F(х).
4. Визначити знак функції F(х) на кінцях інтервалів монотонності.
5. Визначити відрізки, на кінцях яких функція F(х) набуває значень протилежних знаків.
6. Знайдені відрізки ізоляції коренів при необхідності звузити.
1.3 Уточнення кореня методом поділу відрізка пополам
Метод поділу відрізка пополам (або метод дихотомії) застосовний для уточнення кореня рівняння F(х)=0 з наперед заданою точністю, допустимою для даної ЕОМ, якщо функція F(х) задовольняє умови теореми 2.
Позначимо через х - точне значення кореня рівняння (1.1) на відрізку [a;b], а е - його граничну абсолютну похибку. Суть методу в тому, що відрізок [а;b] ділять пополам точкою с0=(a+b)/2 і обчислюють F(c0). Якщо F(с0) = 0, то х=с0 є точним значенням кореня. Якщо F(с0) ? 0, але b - а < 2е, то | x*- с0| <е і значення х=с0 буде шуканим наближеним коренем. Якщо F(с0) ? 0 і b-а >2е , тоді розглядають той з двох відрізків [а;с0] і [с0;b]. На кінцях [с0;b] функція F(х) набуває значень протилежних знаків, позначимо цей відрізок [a;b] (рис. 1.2). На відрізку [a;b] функція F(х) задовольняє умови теореми 2. Далі відрізок [a;b] точкою c1=(a+b)/2 ділять пополам і міркують так само, як і раніше.
Таким чином, після кожної ітерація відрізок, на якому розташований корінь, зменшується вдвічі, тобто після n ітерацій він скорочується в 2n раз. Ітераційний процес продовжується доти, поки значення функції після n ітерації не стане меншим по модулі деякого заданого числа е. При цьому абсолютна похибка знайденого кореня не перевищує е. . Методи дихотомії легко реалізується на ЕОМ, але потребує значного обсягу обчислень, щоб досягти високої точності наближеного кореня.
Програма методу має циклічний характер. У програмі при кожному проходженні циклу виконується серія команд
с=(а+b)/2
якщоF(a)F(b)<0 то b:=c інакше a:=с.
Повторення команд циклу продовжують доти, поки не виконається одна з умов F(с) = 0 або b- а < 2є.
1.4 Метод дотичних
Ідея методу полягає в тому, що проводиться дотична до кривої у = F (х) при х=с і шукається точка перетину дотичної з віссю абсцис. При цьому не обов'язково задавати відрізок [а,b], що містить корінь рівняння (1.1), а досить лише знайти деяке початкове наближення кореня х = с0, в якій F(c0)F"(c0)>0
Рівняння дотичної, проведеної до кривої у = F(x) у точці М0 з координатами с0 і F(co), має вид :
у-F(co)=F' (с0) .
Звідси знайдемо наступне наближення кореня c1 як абсцису точку перетину дотичної з віссю 0х (у =0):
с1=с0-F(c0)/F?(c0).
Аналогічно можуть бути знайдені і наступні наближення як точки перетину дотичних з віссю абсцис , проведених у точках M1, M2, і т.д. (рис. 1.3). Формула для n+1-го наближення має вид
cn+1=cn-F(cn)/F'(cn).
Достатні умови збіжності методу Ньютона дає така теорема.
Теорема . Нехай на відрізку [а;b ]функція F(x) має неперервні із сталими знаками похідні F'(x) ? 0,F''(x)? 0 і F(a)F(b) < 0. Тоді існує такий окіл Rа;b] кореня х* рівняння F(x) = 0, що для будь-якого хоє R послідовність ck}, обчислена за формулою (1.3), збігається до кореня х*.
Оцінимо швидкість збіжності методу Ньютона.
Можна довести, що для методу Ньютона справедлива оцінка
/ck - x*/ ,
де M2=, m1=.
З цієї оцінки видно, що для досягнення заданої точності ітераційний процес треба продовжувати доти, поки для двох послідовних наближень ck і c k-1 не виконуватиметься нерівність
/ck - ck-1/ .
Якщо на відрізку а;b] справедлива нерівність М2<2т1 то ітераційний процес можна закінчити, коли виконується умова /ck - ck-1/ .
Перевага методу Ньютона перед методом ітерацій у тому, що він має вищу швидкість збіжності. Так, корінь x* є [2;3 ] рівняння х3 - 2х -5 =0 з точністю е=10-6 і е=10-9 методом Ньютона був обчислений за п'ять і шість ітерацій відповідно, тоді як методом ітерації він був обчислений не менш ніж за шість і десять ітерацій відповідно.
З формули (1.3) видно, що чим більше значення f(x) в околі кореня, тим менша поправка додається до попереднього наближення. Тому метод Ньютона зручно застосовувати тоді, коли в околі кореня графік функції у = f(x) має значну крутість. Крім того, методом Ньютона можна знаходити не тільки дійсні корені рівнянь, а й комплексні. Метод Ньютона легко поширюється і на розв'язування систем нелінійних рівнянь з багатьма невідомими.
Недоліком методу Ньютона є те, що на кожній ітерації треба обчислювати не тільки значення функції f(x), а й значення її похідної. Обчислення похідної f' (x) може бути значно складнішим від обчислення f(x).
1.5 Метод хорд
Метод хорд - один з поширених ітераційних методів. Його ще називають методом лінійного інтерполювання, методом пропорційних частин. Нехай задано рівняння f(х) = 0, де f(х) на відрізку [а;b] має неперервні похідні першого й другого порядків, які зберігають сталі знаки на цьому відрізку, і f(a) f(b) < 0, тобто корінь х* рівняння відокремлений на [а;b].
Ідея методу хорд в тому, що на досить малому відрізку дуга кривої у=f(x) замінюється хордою і абсциса точки перетину хорди з віссю Ох є наближеним значенням кореня.
Нехай для визначеності f(x) > 0. f'(x) > 0, f(a) < 0. f(b) > 0 (рис. 1.4, а). Візьмемо за початкове наближення шуканого кореня х* значення x0= а. Через точки A0 і В проведемо хорду
і за перше наближення кореня х* візьмемо абсцису x1 точки перетину хорди з віссю Ох. Поклавши y=0, х=x1одержуємо
або
Размещено на http://www.allbest.ru/
0
Новий відрізок, що відокремлює корінь, можна визначити, порівнюючи знаки f (а), f (x1) і f (b). Очевидно, що точка x1 ближче до точки х*, чим а, якщо у'у" > 0 (див. рис. 7.4,а), і відрізком, що відокремлює корінь, буде [х1,b] ,у противному випадку, якщо у'у" < 0 (див. рис.1.4, в), відрізком, що відокремлює корінь, буде [а,х1].
Абсциса х2 точки перетину хорди А1В буде другим наближенням кореня. Продовжуючи цей процес необмежено, дістанемо послідовність хо, х1, х2, ..., хn, ...наближених значень кореня х* даного рівняння.
Для хк+1, якщо у'(х)у"(х) > 0 маємо:
к=0,1,2,…
Якщо у'(х)у"(х) < 0 , то для хк+1 можна записати формулу
, к=0,1,2,…
Достатні умови збіжності методу хорд дає така теорема.
Теорема. Нехай на відрізку [а;b] функція f(x) неперервна разом із своїми похідними до другого порядку включно, причому f(a)f(b) < о, а похідні f'(x) і f''(х) зберігають сталі знаки на [а;b], тоді існує такий окіл кореня х* рівняння f(x)= 0, що для будь-якого початкового наближення хо з цього околу послідовність [хk], обчислена за формулою (1.5) або (1.6), збігатиметься до кореня х*.
Для оцінки похибки можна скористатися формулою
/хk - x*/ ,
де M2=,
m1=.
Отже, корінь х* рівняння f(x) =0 буде знайдено методом хорд із наперед заданою точністю є , якщо для двох послідовних наближень хk.і хk-1 справджуватиметься нерівність
.
1.6 Комбінований метод дотичних і хорд
Характерна особливість методів дотичних і хорд та, що послідовності їх наближень монотонні. Причому, якщо для даного рівняння послідовність наближень методу хорд монотонно спадна, то послідовність наближень методу дотичних - монотонне зростаюча, і навпаки. Одночасне застосування цих методів дає змогу наближатися до кореня. рівняння з двох боків, дістаючи наближення з недостачею і надлишком.
Розглянемо рівняння f(x) =0, корінь якого х* є [а;b]. Нехай, наприклад, f(x) > 0, f'(x) > 0, f(а) < 0, f(b)>0.
У даному випадку за початкове наближення в методі хорд вибирають точку х = а, а в методі дотичних - точку b. На відрізку [а;b] застосовують метод дотичних і хорд (спочатку дотичних, а потім хорд). У результаті дістають нові наближення а1, b1 і початковий відрізок ізоляції кореня звузився. Для знаходження нових наближень застосовують метод дотичних і хорд уже на відрізку [а1;b1]. У результаті дістають наближення а2 і b2 відповідно, причому [а2;b2] [a1; b1] [а;b]. Такий процес продовжують доти, поки довжина відрізка [ak;bk] стане меншою або дорівнюватиме величині 2е, де е - наперед задана точність кореня.
За шукане значення кореня беруть півсуму наближень ak і b k, тобто = 0,5(ak + bk), а модуль їх піврізниці дасть граничну абсолютну похибку наближеного кореня, тобто
.
Зазначимо, що на кожному кроці комбінованого методу за нерухомий кінець с у формулі методу хорд треба брати наближення, обчислене на цьому самому кроці за формулою, дотичних.
Формули комбінованого методу дотичних і хорд мають вигляд
bk+1=bk-F(bk)/F'(bk), k=0,1,2…
, k=0,1,2…
За початкове наближення bo у формулі (1.7) методу дотичних беруть той з кінців відрізка [а;b]. в якому значення функції f(x) і її другої похідної мають однакові знаки, тоді протилежний кінець відрізка [а;b] беруть за початкове наближення ао у формулі (1.8) методу хорд.
Завдяки своєрідній комбінації методів дотичних і хорд комбінований метод має вищу швидкість збіжності, ніж методи хорд і дотичних окремо взяті. Так, для рівняння х3-2х-5= 0 його корінь x*є [2;3 ] з точністю до 10-6 і 10-9 комбінованим методом обчислено за дві і три ітерації відповідно, тоді як для методу дотичних потрібно було п'ять і шість ітерацій відповідно, а для методу хорд - дванадцять і дев'ятнадцять ітерацій відповідно.
1.7 Метод простої ітерації
Замінимо рівняння (1.1) еквівалентним йому рівнянням
x = (x).
Це можна зробити різними способами, наприклад
x = x + cf(x), c0.
Припустимо, що вибрано деяке початкове наближення кореня рівняння (1.9). Метод простої ітерації полягає у послідовному знаходженні послідовності {xn} за формулою
xn+1 = (xn), n=0, 1, 2,
Розв'язком рівняння буде границя послідовності, визначеної співвідношенням. Якщо функція (x) має неперервну похідну (x) та
|(x)| q < 1,
то така границя існує і швидкість збіжності визначається нерівністю
| - xn||x1 - x0|
Як видно з нерівності (1.13), швидкість збіжності метода простої ітерації залежить від величини q: чим менше q, тим скоріша збіжність. Тому на практиці важливо представити рівняння (1.9) таким чином, щоб похідна (x) в околі кореня за абсолютною величиною була якомога меншою. Цього дуже часто досягають вибором параметру с у рівнянні (1.10).
РОЗДІЛ 2. ПРАКТИЧНА ЧАСТИНА
Приклади на відокремити корінь рівняння
1. Відокремити корінь рівняння
exsіnx - 1 =0.
Рішення. Побудувати графік функції у =exsіnx - 1 важко. Перепишемо вихідне рівняння у такий спосіб: рівняння похибка корінь ітерація
sin х = е-х .
Побудувавши графіки функцій y=sinx і у= = е-х, бачимо (рис. 2.1), що вони перетинаються в нескінченному числі точок. Отже, наше рівняння має нескінченну множину дійсних корінь, причому всі вони додатні. Можна вказати відрізки, на кожному з яких лежить один і тільки один корінь вихідного рівняння:
[2k; 2k + /2 ] і [2k + /2; (2k + 1)], k = 0, 1, 2, ...
Так, наприклад, найменший додатній корінь розташований на відрізку [0; /2].
Рис.2.1
2. Відокремити корені рівняння x3-x-1=0.
Рішення.
Функція f(x)=x3-х-1 і її похідна , f '(x)=Зх2-1,безупинні на всій числовій осі. Визначимо інтервали монотонності функції f(x). Для цієї мети вирішимо рівняння f '(x)=Зх2-1=0.
Його корені x1 = і x2 = . Отже, інтервалами монотонності функції f(x) є інтервали
(- ; ) ( ; ) ( ; )
При цьому
= x3-х-1 = -
f( )<0, f( )<0
= x3-х-1 = +
Отже, інтервали монотонності
(- ; ) ( ; )
не містять коренів даного рівняння, тому що функція f(х)=х3-х-1 не змінює знак на цих інтервалах, а на інтервалі ( ; ) розташований його єдиний дійсний корінь.
Зі знайденого нескінченного інтервалу виділимо тепер відрізок, що містить єдиний корінь рівняння. У даному випадку для цієї мети досить безпосередньо перевірити знак функції f(x)=x3-х-1 у приналежних ( ; ) целочисленних точках .
Маємо f( )<0, f(1)= 1-1-1 <0, f(2)=8-2-1>0
Отже, єдиний дійсний корінь даного рівняння лежить на відрізку [1, 2].
Поставлена задача вирішена . Подальше звуження границь знайденого інтервалу - це вже - задача про уточнення кореня.
Приклади на застосування методу поділу відрізка пополам
1. Відділити корені аналітично та уточнити один з них методом половинного ділення з точністю до 0,01:
x4 - x3 - 2x2 + 3х - 3 = 0.
- Розв'язання
- Позначимо f(x)= x4 - x3 - 2x2 + 3х - 3.
- Знайдемо похідну = 4x3 - 3x2 - 4х + 3. Обчислимо корені похідної:
- 4x3 - 3x2 - 4х + 3=0; 4x(x2-1)-3(x2-1)=0;
- (x2-1)(4x-3)=0; x1=-1; x2=1; x3=3/4.
- Складемо таблицю знаків функції f(x):
- З таблиці видно, що дане рівняння має два дійсних корені:
- x1(-, -1; x21, +).
- Зменшимо проміжки, в яких знаходяться корені:
- Отже, x1-2, -1; x21, 2.
- Уточнимо один з коренів, наприклад x1-2, -1, методом половинного ділення з точністю до 0,01. Усі обчислення зручно проводити, використовуючи наступну таблицю
- Відповідь: x1-1,73.
- Перепишемо рівняння у вигляді
- .
- Позначимо
- ,
- З графіка видно, що рівняння має два корені: , .
- 4. Запишемо рівняння у вигляді: . Позначимо , , побудуємо графік цих функцій (рис.2. 3). З графіка бачимо, що рівняння має один корінь .
- Для уточнення цього кореня методом спроб виберемо проміжок, на кінцях якого функція має різні знаки. Складемо таблицю:
- Для зручності обрахунків перейдемо до десяткових логарифмів:
- Рис. 2.2
- Рис. 2.3
- Подальші розрахунки проведемо в таблиці:
- Відповідь: .
- Приклади на застосування метод хорд
- 1. Відділити корені графічно та уточнити один з них методом хорд з точністю до 0,001:
- tg(0,55x + 0,1)=x2.
- Розв'язання
- Відділимо корені графічно.
- Побудуємо графіки функцій у1= tg(0,55x + 0,1); у2=х2.
- Отже, додатній корінь рівняння знаходиться у проміжку 0,6; 0,8.
- Для того, щоб уточнити корінь методом хорд, визначимо знаки функції f(x) = tg(0,55x + 0,1)-x2 на кінцях проміжку 0,6; 0,8 та знак другої похідної на цьому проміжку:
- f(0,6)=0,0986; f(0,8)=-0,0406;
- =0,605sin(0,55x+0,1)/(cos3(0,55x+0,1))-2<0 при x0,7; 0,8.
- Тому для обчислень використаємо формулу
- xn+1=xn-, x0=0,6.
- За цією ітераційною схемою матимемо x1=0,742; x2=0,75; x3=0,7502.
- Відповідь: x10,75.
- 2. Методом хорд уточнити корінь рівняння x3-x-1=0,відділений на відрізку [1,2].
- Розв'язання
- f(x)=x3-x-1,
- f'(x)=Зх2-1,
- f”(x)=6.
- Тому що на відрізку [1, 2] маємо f”(x)>0, f(1)<0, f(2)>0, то крапка х1, обумовлена по формулі (5.4)
- =1- =1.1,
- буде лівим кінцем нового відрізка [a1, b1] (х1<x*<b, a1=х1, b1=b.
- Помітимо, що наближене значення с узяте з недоліком, тому що х1<х*, і при округленні з надлишком є небезпека «переступити» через корінь. Як відрізок [a1, b1] варто тепер узяти відрізок [1,1; 2].
- Приклади на застосування методу дотичних
- 1. Методом дотичних уточнити корінь рівняння x3-x-1=0,відділений на відрізку [1, 2].
- Розв'язання
- Маємо f(x)=x3-x-1, f'(x)=Зх2-1, f”(x)=6.
- Тому що на відрізку [1, 2] f”(x)>0, то дотичну треба взяти в правому кінці відрізка (тобто x=b=2, де f(2)=5>0. Ще не обчислюючи c, можна сказати, що точка c буде правим кінцем нового відрізка [a1, b1] (а<xо<c, a1=a, b 1=с). Думаючи c=2, одержимо f'(2)=11 і далі по формулі cn+1=cn-f(cn)/f'(cn): c=2-=1.6
- Наближене значення c узято з надлишком, тому що c>x0. Новий відрізок повинний бути [1; 1,6] (a1=a=1, b 1=c=l,6).
- Ми уточнювали корінь рівняння, заміняв один з кінців вихідного відрізка [а, b] точкою c, більш близької до кореня. Очевидно, більший ефект може бути досягнуть при наближенні до кореня одночасно c двох сторін. Це досягається застосуванням комбінованих методів.
- Приклади на застосування комбінованого методу дотичних і хорд
- Комбінованим методом уточнити корінь рівняння exsіnx - 1 =0 на відрізку
- [0; /2].
- Розв'язання
- За умовою задачі необхідно визначити таке x = , що f()=0, де
- f(x)= exsіnx - 1.
- Функція f(x) на відрізку [0; /2] має неперервні похідні:
- f(x)= ex(sіnx + cosx),
- f(x)= 2 ex cosx.
- Визначимо більш вузький відрізок, на якому функція f(x) змінювала б свій знак, а її похідні не змінювали б свого знаку. Таким відрізком може бути відрізок [0; /3]. Дійсно
- f(0)<0, f(/3)>0; f(x) та f(x) - додатні на [0; /3].
- Оскільки f(0) f(x) < 0, то відповідно до методу хорд, наближення до кореня зліва буде визначатись за формулами
- ,
- де наближення з парним індексом визначаються методом хорд, а наближення з непарним індексом - методом дотичних. Останні визначатимуться за формулами
- x1=/3 - f(/3)/ f(/3),
- x2n+1= x2n-1 - f(x2n-1)/ f(x2n-1),
- Співвідношення (2.1) та (2.3) дають нам початкові наближення:
- x1=0,67; x2=0,424;
- далі використовуємо по черзі формули (2.4) та (2.2), починаючи з (2.4) для n=1, 2, 3,…:
- x3= x1 - f(x1)/ f(x1) = 0,592;
- ; ....
- Як бачимо, ітераційний процес за методом Ньютона (2.4) дає наближення справа до кореня (x1, x3), а ітераційний процес за методом хорд (2.2) - наближення зліва до кореня (x2, x4). Отже, з точністю =0,005 можна вважати коренем даного рівняння значення =0,59.
- 2. Методом хорд і дотичних уточнити корінь рівняння x3-x-1=0, відділений на відрізку [1,2].
- На відрізку [1, 2] маємо f(x)=x3-x-1, f(1)=-1<0, f(2)=5>0;
- f'(x)=Зх2-1, f”(x)=6>0.
- Отже, а1 треба обчислювати по формулі
- - метод хорд, а b1-по формулі b1=b-F(b)/F'(b)- метод дотичних. Запозичаючи результати обчислень по формулам із прикладів 3 і 4, одержимо, що відрізок [a1, b1] повинен бути відрізок [1,1; 1,6].
- Розв'язання
- Позначимо f(x)= x3 - 2x2 + 7х + 3.
- Знайдемо похідну = 3x2 - 4х + 7. Корені похідної відсутні: при будь-яких значеннях х. Отже функція f(x) є зростаючою.
- Складемо таблицю знаків функції f(x):
- Дане рівняння має єдиний корінь, що лежить у проміжку -1, 0.
- Приведемо рівняння до вигляду х=(х) так, щоб виконувалась умова |(x)| q < 1, при х-1, 0. Оскільки max=14 при х-1, 0, то (х)=x- f(x)/14. Тоді (x)=-3x2/14+2x/7+1/2>0 на (-1, 0, а, отже, максимальне значення на -1, 0 досягається на кінцях відрізку. Порівнюючи та матимемо, що q=max. Отже, ітераційний процес матиме вигляд
- xn+1=xn-0,1(xn)3+0,2(xn)2-0,7xn-0,3.
- 2. Уточнити корінь рівняння exsіnx - 1 =0 на відрізку [0; /2] методом простої ітерації.
- Розв'язання
- Запишемо рівняння у вигляді x = x + cf(x), c0:
- x = x + c( exsіnx - 1),
- тобто у нашому випадку (x) = x + c( exsіnx - 1). Виберемо постійну с таким чином, щоб |(x)| =|1+cex(sinx + cosx)| < 1 на відрізку [0; /2]. Для цього дослідимо поведінку функції (x) на [0; /2]. Оскільки (x)=2 cexsinx, на (0; /2)
- (x) < 0 , якщо с < 0;
- (x) > 0 , якщо с > 0.
- Нам краще вибрати випадок, коли (x) < 0 (с < 0), що відповідає монотонно спадаючій функції (x), так як (0)=1 + с > 0 при c > -1. Отже, необхідно вибрати с з інтервалу (-1; 0). Для того, щоб задовольнити умову x = x + cf(x), c0, необхідно визначити найменше та найбільше значення спадаючої функції (x) на відрізку [0; /2]. Найбільше значення: (0)=1 + с > 0, а найменше: (/2) 1 + 4,8с. Бажано, щоб вони обоє були якомога ближче до нуля. Тому слід вибрати значення постійної с таким чином, щоб виконувалась рівність (0) + (/2)=0. З цієї рівності одержимо найоптимальніше значення с=-0,345. Тоді нерівність x = x + cf(x), c0 буде виконуватись при q=1-0,345=0,655.
X |
- |
-1 |
ѕ |
1 |
+ |
|
Sign f(x) |
+ |
- |
- |
- |
+ |
X |
-2 |
-1 |
1 |
2 |
|
sign f(x) |
+ |
- |
- |
+ |
N |
an+ |
bn- |
xn=( an+ bn)/2 |
f(xn) |
|
0 |
-2 |
-1 |
-1,5 |
-3,5625 |
|
1 |
-2 |
-1,5 |
-1,75 |
3,3633 |
|
2 |
-1,75 |
-1,5 |
-1,63 |
-1,814 |
|
3 |
-1,75 |
-1,63 |
-1,69 |
-0,7981 |
|
4 |
-1,75 |
-1,69 |
-1,72 |
-0,2363 |
|
5 |
-1,75 |
-1,72 |
-1,73 |
-0,0406 |
|
6 |
-1,75 |
-1,73 |
-1,74 |
0,1592 |
|
7 |
-1,74 |
-1,73 |
x |
|||
+ |
Зауваження. У даному випадку комбінація методу хорд із методом поділу відрізка пополам дала б більш вузький відрізок [1,1; 1,5]. Однако у загальному випадку комбінування цих двох методів не можна вказати умови, що забезпечують наближення до кореня з двох сторін. Це є наслідком того, що, як уже вказуюся, метод розподілу відрізка навпіл не зв'язаний з поводженням функції f(x)-лівої частини заданого рівняння.
Приклади на застосування методу простої ітерації
1. Відділити корені аналітично та уточнити один з них методом ітерацій:
x3 - 2x2 + 7х + 3 = 0.
X |
- |
-1 |
0 |
+ |
|
Sign f(x) |
- |
- |
+ |
+ |
Подальше розв'язання задачі полягає у виконанні ітераційного процесу
xn+1 = xn - 0,345(sіnxn - 1),
взявши за початкове наближення середину відрізка [0; /2], оскільки корінь ближче до середини, аніж до кінців відрізка.
ВИСНОВКИ
Розглянувши дану тему, стає очевидно, що питання розглянуті в курсовій роботі були актуальні ще в 16 столітті, і не втрачають важливості до нашого часу.
В курсовій роботі було детально розглянуто такі питання: відокремлення коренів; теорема про оцінку похибки наближеного значення кореня; уточнення кореня методом поділу відрізка пополам; метод дотичних; метод хорд; комбінований метод дотичних і хор; метод простої ітерації (доведення та наслідки.
Також було наведено приклади розв'язання рівнянь з допомогою розглянутих методів наближеного обчислення дійсних коренів многочлена.
СПИСОК ВИКОРИСТАНИХ ДЖЕРЕЛ
1. Барковський В.В., Барковська Н.В. Математика для економістів: Вища математика. - К.: Національна академія управління, 1997.
2. Завало С.Т., Костарчук В.М., Хацет Б.І. Алгебра і теорія чисел. 2 том. - К.: Вища школа, 1976. - 384 с.
3. Каплан И.А. Практические занятия по высшей математике. - Харьков: ХГУ, 1967.
4. Карасев А.И. и др. Курс высшей математики для экономических вузов. - М.: Высш. шк., 1982. - Ч. 1,2.
5. Корн Г., Корн Т. Справочник по математике. - М.: Наука, 1980.
6. Костарчук В.М., Хацет Б.І. Курс вищої алгебри. - К.: Рад. Школа, 1964. - 512 с.
7. Кудрявцев В. А., Демидович В. П. Краткий курс высшей математики. -М.: Физматгиз, 1962,1975.
8. Кузнецов Л.А. Сборник заданий по высшей математике (типовые расчеты).- М.: Высш. шк., 1982.
9. Курош А.Г. Курс высшей алгебры. - М.:Наука, 1968. - 432 с.
10. Лященко М.Я., Головань М.С. Чисельні методи : Підручник. - К.: Либідь, 1996. - 288 с.
11. Минорский В. П. Сборник задач по высшей математике. - М.: Наука, 1969.
Размещено на Allbest.ru
Подобные документы
Метод простої ітерації Якобі і метод Зейделя. Необхідна і достатня умова збіжності методу простої ітерації для розв’язання системи лінейних рівнянь. Оцінка похибки. Діагональне домінування матриці як умова збіжності ітерації. Основні переваги цих методів.
презентация [79,9 K], добавлен 06.02.2014Схема класифікації та методи розв'язування рівнянь. Метод половинного ділення. Алгоритм. Метод хорд, Ньютона, їх проблеми. Граф-схема алгоритму Ньютона. Метод простої ітерації. Питання збіжності методу простої ітерації. Теорема про стискаючі відображення.
презентация [310,1 K], добавлен 06.02.2014Методика визначення всіх коренів нелінійного рівняння різними способами: відрізка пополам, хорд, дотичних та ітерацій. Особливості та принципи застосування комп’ютерних технологій в даному процесі. Аналіз отаманих результатів і їх інтерпретація.
лабораторная работа [263,9 K], добавлен 15.12.2015Вивчення рівняння з однією невідомою довільного степеня та способів знаходження коренів таких рівнянь. Доведення основної теореми алгебри. Огляд способу Ньютона встановлення меж дійсних коренів алгебраїчних рівнянь. Відокремлення коренів методом Штурма.
курсовая работа [1,1 M], добавлен 06.10.2012Графічний спосіб розв'язку рівнянь. Комбінований метод пошуку та відокремлення коренів. Метод Ньютона (метод дотичних або лінеаризації). Процедура Ейткена прискорення збіжності. Метод половинного поділу та простих ітерацій уточнення коренів рівняння.
лекция [1,9 M], добавлен 27.07.2013Межі дійсних коренів. Опис та текст програми. Методи наближеного пошуку меж та самих коренів многочлена з дійсними коренями. Метод пошуку точних значень многочленів з числовими коефіцієнтами. Контрольний приклад находження відрізків додатних коренів.
курсовая работа [49,5 K], добавлен 28.03.2009Вивчення теорії наближених обчислень і чисельних методів лінійної алгебри. Опис прямих і ітераційних методів вирішення систем лінійних рівнянь, алгоритмізація і точність наближених обчислень функції. Чисельна інтеграція звичайних диференціальних рівнянь.
лекция [103,6 K], добавлен 06.02.2014Класифікація та типи чисельних методів розв’язування систем лінійних рівнянь і обернення звернення матриць точні, ітераційні та комбіновані. Їх порівняльна характеристика та умови використання в окремих випадках. Вектори та операції над ними, норми.
презентация [85,6 K], добавлен 06.02.2014Аналіз найвідоміших методів розв’язування звичайних диференціальних рівнянь і їх систем, користуючись рекомендованою літературою. Розробка відповідної схеми алгоритму. Розв’язання системи звичайних диференціальних рівнянь в за допомогою MathCAD.
лабораторная работа [412,4 K], добавлен 21.10.2014Основні етапи розв'язування алгебраїчних рівнянь: аналіз задачі, пошук плану розв'язування та його здійснення; перевірка та розгляд інших способів виконання. Раціоналізація розв'язування алгебраїчних рівнянь вищих степенів методом заміни змінних.
курсовая работа [229,8 K], добавлен 13.05.2013