Page 31 - 6760
P. 31

ПРАКТИЧНЕ ЗАНЯТТЯ № 7


                         Тема: Метод Рунге-Кутта для систем диференціальних
                         рівнянь рівнянь.


                         Технічне забезпечення: ПЕОМ середовище програмування


                                              Короткі теоретичні відомості
                         Для  чисельного  інтегрування  систем  диференціальних
                  рівнянь        та     диференціальних             рівнянь        вищих         порядків

                  застосовуються  ті  ж  методи  і  схеми,  які  наведені  в  роботі  №7,
                  єдиною  відмінністю  є  те,  що  вони  записуються  у  векторній
                  формі; наприклад, для системи:
                                               
                                               x   f   ;,tx  x  0                                (7.1)
                                                                  x
                                                                   0
                  схему четвертого порядку точності записується у вигляді:
                                                         k   f x , ,
                                                                     t
                                                          1      n  n
                                                  k   f   x  5 . 0   ; tk    5 . 0  ,
                                                   2      n       1  n
                                                     k   f   x  5 . 0   ; tk    5 . 0  ,         (7.2)
                                                       3      n       2  n
                                                 k   f  x   5 . 0   ; tk    5 . 0  ;
                                                  4      n        3  n
                                                                
                                                      x     x   k   2k   2k   k  .              (7.3)
                                                       n 1  n      1     2     3   4
                                                                 6
                         Кожен  з  виразів  k   буде  складатись  з  такої  кількості
                                                      i
                  розрахункових  формул,  яка  дорівнює  кількості  невідомих  в
                  системі  рівнянь.  Таким  чином,  для  інтегрування  систем  из
                  кроком  по часу одержується явна обчичлювальна схема.

                         Оскільки  будь-яке  рівняння  вищого  порядку  може  бути
                  зведене до системи диференціальних рівнянь з  n невідомими (n –
                  порядок рівняння), то підхід (7.2 – 7.3) може бути застосований

                  до інтегрування таких рівняннь.
                         Приклад  1.  Записати  рівняння  x   2              x  x  2    x    t f ,    0 xx  ,
                                                                                                           0
                     0
                  x      x  у вигляді системи диференційних рівнянь. Вводимо нову
                           1
                  змінну  x  .  В  такому  випадку  рівняння  записується  у  вигляді
                              
                                  y
                  системи:
                                                           yx 
                                                                          2
                                                          y       xt f  2xy





                                                              31
   26   27   28   29   30   31   32   33   34   35   36