Курсовая работа: Основные межвидовые взаимодействия, эволюция

[t,NN]=ode45('functhion',[0 3], N0, [], alfa, beta);

plot(t,NN,'g')

end

Очевидно, что неограниченный рост популяции невозможен, начиная с некоторого момента, он будет тормозиться в результате внутривидовой конкуренции.

Нейтрализм:

Численность популяции первого вида не зависит от численности популяции второго вида, поэтому не имеет смысла строить график M(N), построим лишь графики M(t) и N(t).

Дифференциальные уравнения изменения численностей популяций имеют следующий вид:

Равновесие наступает при

Содержимое sistd.m:

function sist=func(t,p)

global a1 a2 b1 b2;

sist=[(a1-b1*p(1))*p(1); (a2-b2*p(2))*p(2)];

Содержимое work1.m

global a1 a2 b1 b2;

a1=2.5;

a2=3.1;

b1=0.002;

b2=0.0023;

N0=1400;

M0=1200;

[t X]=ode45('sistd',[0 5],[N0 M0]);

N=X(:,1);

M=X(:,2);

figure

hold on

grid on

plot(t,N,'b')