Разработка и анализ математических моделей с использованием MATLAB и MAPLE. Коробейников А.Г. - 47 стр.

UptoLike

ВУЗ: 

Университет ИТМО | Санкт-Петербург

Составители: 

Рубрика: 

46
в момент времени t, считая известной начальную численность населения Р
0
= Р(0).
На рис. 6.3 показано множество интегральных кривых,
соответствующих различным значениям начальной численности населения
в пределах от Р
0
= 20 до Р
0
= 300. Обратите внимание, что интегральные
кривые асимптотически приближаются к горизонтальной линии Р = 150.
Действительно, из найденного решения видно, что
(
)
lim 150
t
Pt
→∞
=
независимо
от начального значения Р
0
> 0.
Текст на MATLAB
clc
clear
k=-0.06;
t=0:0.2:150;
length_t = length(t);
P=t;
P0=[ 20; 50; 80; 100; 120; 150; 200; 240; 270; 300];
length_P0 = length(P0)
for i_P0=1:length_P0
for i_t=1:length_t
P(i_t)=150*P0(i_P0)/(P0(i_P0)+(150-P0(i_P0))*exp(k*t(i_t)));
end
plot(t,P);
hold on
end
title('P(t)=150*P0/(P0+(150-P0)*e^-^0^.^0^6^*^t');% \lambda
xlabel('t');
ylabel('P0');
hText=text(10,30,'\leftarrow P0=20');
hText=text(12,210,'\leftarrow P0=300');
grid on

Страницы