Лабораторная работа n o 1 «Развитие биологических популяций»

УО «Полоцкий государственный университет имени Евфросинии Полоцкой»

Кафедра физики

ЛАБОРАТОРНАЯ РАБОТА N^o 1

«Развитие биологических популяций»

Группа 21-КФ

Выполнил (а) ________________ Ф.И.О. студента

Проверил (а) ____________ Ф.И.О. преподавателя

Новополоцк

20__г.

Цель работы: Высчитать популяцию схемы Хищник-Жертва
Листинг программы MatLab:
global a1 a2 b12 b21 u
a1 = 5;
b12 = 0.01;
a2 = 3;
b21 = 0.0001;
x10 = 30000;
x20 = 600;
x0 = [x10; x20];
dt = [0 5];
opt = odeset('RelTol', 1e-6);
[t, x] = ode45(@predator, dt, x0, opt);
subplot(2,1,1);
plot(t,x(:,2));
xlabel('time, year');
ylabel('predator fish');
subplot(2,1,2);
plot(t,x(:,1));
xlabel('time,year');
ylabel('rpey fish');
Функция:
function dt = predator(t,x)
global a1 a2 b12 b21 u
if (t< 2)
u = 0;
elseif 2 <= t <= 2.25
u = 1;
else
u = 0;
end
dt = [((a1 - b12 .* x(2)) .* x(1)) - 200 * u; ((-a2 + b21 .* x(1)) .* x(2)) - 200 * u];


Графики:






Вывод: : Высчитал популяцию схемы Хищник-Жертва