Практична і лабораторна робота 3 з дисципліни „Моделювання систем”.
Дослідження логістичного рівняння.
Поведінку популяції, чисельність якої стабілізується на деякому заданому рівні, часто описують за допомогою логістичного рівняння
dx / dt = rx (1- X / k),(1)
де параметр k задає максимальну густину популяції.
Рівняння (1) є нелінійним звичайним диференціальним рівнянням першого порядку. Це рівняння завдяки простоті має аналітичний розв’язок
, (2)
де X 0 – початкове значення X (t).
Неважко помітити, що функція (2) за малих t поводить себе як X 0 ert, тобто є розв’язком лінійного диференціального рівняння dx / dt = rx, і лише із зростанням t виявляється нелінійність рівняння (1).
Можливий розв’язок логістичного рівняння (1) показано на рис.1.
Широке використання логістичного рівняння для моделювання розвитку популяцій пояснюється його простотою і такими властивостями:
– за малих X рівняння (1) зводиться до лінійного рівняння dx / dt = rx, тобто ріст іде експоненціально для малих X;
– для більших X величина X з часом прямує до постійного значення k.
Текст програми на мові MATLAB-6 наступний.
%********* ELABORATION of LOGHISTIC’s EQUATION *************
%********* Основна програма ********************************
сlear;
tlimit=[0, 10]; X0=[0.01];
[t,X]=ode45(@logeq,tlimit,X0);
plot(t,X)
%******** Вміст файлу logeq.m ********************************
function dXdt=logeq(t,X)
dXdt=[ 1.5*X*(1-X/2) ];
В основній програмі після очищення робочої області пам’яті оператором сlear; задається часовий інтервал розв’язування логістичного рівняння tlimit і початкова умова X0.
Далі знаходимо два вектори розв’язку логістичного рівняння [t,X] за допомогою функції ode45, що розв’язує за методом Рунге-Кутта четвертого порядку. Аргументами функції ode45 є ім’я файла з описом диференціального рівняння, двоелементний вектор початку і кінця інтервалу інтегрування та початкова умова.
Останній оператор основної програми будує графік розв’язку у графічному вікні.
Логістичне диференціальне рівняння (1) описане у файлі з іменем logeq.m. Останній оператор цього файлу містить числові значення двох параметрів логістичного рівняння – r =1.5 і k =2. Змінюючи їх значення, а також значення початкової умови треба дослідити поведінку логістичної кривої (2).
Хід лабораторної роботи.
1. Вивчити теоретичну частину.
2. Набрати та відлагодити програму.
3. Отримати логістичну криву для трьох значень параметра r > 0, трьох значень параметра k > 0 та трьох значень початкової умови X (0) > 0, зокрема X (0) < k та X (0) > k.
4. Написати та захистити звіт.