你在function下面，segma1上面增加一行 dx = zeros(2,1); 就可以了
可这样写：
function mydae
[t,y]=ode45(@rigid,[0,8],[0.1,0.1]);
plot(t,y)
xlabel('t')
ylabel('y')
function dy=rigid(t,x)
segmental1=.........%我省略
dx=[
r1*x(1).*(1-x(1)/N1-segmental1*x(2)/N2);
r2*x(2).*(1-segmental2*x(1)/N1-x(2)/N2)];
上面语句保存为函数后试试
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024