关于#matlab#的问题：matlab微分方程模型，说我索引超出数组元素的数目(2)，程序无法执行

Question

matlab微分方程模型，说我索引超出数组元素的数目(2)。求大家帮我看看我哪里写错了。

Answer 1

输入变量必须是3维的，但是你的是二维的，维数不对
t与y都必须是三维的才行

clc;

R=[0.1,0.2,0.3];
A=[0,0.05,0.08; 0.1,0,0.18;0.06,0.12,0];
B=[2,1.8,1.4];
P0=0.8;
D=[0.02,0.01,0.05];
[t, y]=ode45(@(t,y)eq3(t,y,R,A,B,P0,D),[0,12,15],[100,100,100]);

function dy=eq3(t,y,R,A,B,P0,D)%说明微分变里是3维的
    dy=zeros(3,1);
    %微分方程组
    dy(1)=R(1,1)*y(1)*(1-A(2,1)*y(2)-A(3,1)*y(3))-B(1,1)*P0*(1+sin(pi/6*t))-D(1,1)*y(1);
    dy(2)=R(1,2)*y(2)*(1-A(1,2)*y(1)-A(3,2)*y(3))-B(1,2)*P0*(1+sin(pi/6*t))-D(1,2)*y(2);
    dy(3)=R(1,3)*y(3)*(1-A(1,3)*y(1)-A(2,3)*y(2))-B(1,3)*P0*(1+sin(pi/6*t))-D(1,3)*y(3);
end

如果问题得到解决的话请点 采纳~~