Question

Matlab求解微分方程数值解出现错误

Answer 1

你在function下面，segma1上面增加一行 dx = zeros(2,1); 就可以了

追问

为啥啊？什么时候要加那一行啊

追答

ode问题求解时候输出要写成列向量，zeros(2,1)就是为了把dx变成2 x 1的列向量。这样就行了

Answer 2

可这样写：
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)];
上面语句保存为函数后试试