Page 17 - 4299
P. 17

К6. Якщо      ) <0, то a       x  і перейти до К1.
                              f x
                                     f b
                                  j     j           j    j











                                                  а) укорочення інтервалу зліва;















                                               б) укорочення інтервалу справа

                                Рисунок 4.1 – Розв’язок рівняння  (f      ) x =0 методом дихотомії

                       Умова на другому кроці (К2), як правило, не виконується, оскільки, ітераційний процес
               дає  можливість  знайти  лише  наближене  значення  кореня.  Тому  її  слід  замінити  іншою
               умовою  f ( x )  , де >0 – досить мати число, яке визначає точність розв’язку задачі.
                                  
                             j
                                                                              
                     Оскільки на кожному кроці ітерації вихідний відрізок  a ;x  ділиться наполовину, то
                                                                              
                                                                                    j 
                                                                                 j
               він породжує таку послідовність точок, що
                                                          b a
                                                 c   x        , k   0 1 2, , ,....
                                                  j   k     k 1
                                                           2

                     Це означає, що послідовність значень  x     c j   при k   .
                                                              k
                   Для заданої точності розв'язку задачі  (f   ) x =0 -  c   x    із останнього співвідношення
                                                                       j   k
               визначимо  максимальну  кількість  ітерацій  N,  які  необхідно  виконати  для  досягнення
                               
                               6
               точності   10 .
                                                             b a         
                                            N = max : k=     log    / log(  2 )  ,
                                                                                1
                                                                             
                                                                         
               де [...]- ціла частина числа.

                     Розв’язок задачі.
                     Для  розв’язання  нелінійного  алгебраїчного  рівняння      xln 1 x  2  e   x    0  методом
               дихотомії скористаємося машинною програмою, яка працює у середовищі MatLab.
               %================================================
               %Розвязок рівняння f(x)=0 методом діхотомії
               %Вхід: Функція f(x)
               %     a і b - початок і кінець інтервалу [a;b]
               %     epselon - точність розв. задачі f(x)=0
               %Вихід:с - корінь рівяння f(x)=0
               %      Yc - значення функції f(x) при x=c
               %      err - похибка у знаходженні кореня
               %================================================
                                                              16
   12   13   14   15   16   17   18   19   20   21   22