先写函数文件如下
function ani=animal(t,x)
ani=[14*x(1)*(1-1/28*x(1)-1/14*x(2));
16*x(2)*(1-1/32*x(2)-1/16*x(1))];
以下是函数文件调用代码
////////////////////////////////////////
[t,x]=ode45('animal',[0,8],[1,1]); %画相轨线
plot(x(:,1),x(:,2),'b')
hold on
[t,x]=ode45('animal',[0,8],[14,16]); %画相轨线
plot(x(:,1),x(:,2),'r')
grid on
xlabel('甲种群密度x1');
ylabel('乙种群密度x2');
title('甲乙种群相轨线');
////////////////////////////////////////////////////
以下是运行结果:
没方程没数据你叫我怎么画,你去百度百科相轨线,把dx/dt = P (t,x,y) dy/dt = Q(t,x,y)给我先