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

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

img

img

输入变量必须是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

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

答：可能是你索引的数组元素超过了你定义的数组大小，请检查一下你的数组大小和你索引的数组元素是否存在冲突。另外，请检查一下你的微分方程模型，是否存在语法错误，或是变量声明错误导致的bug。

方案来自 梦想橡皮擦 狂飙组基于 GPT 编写的 “程秘”

索引超出数组元素的数目通常表示你正在访问一个不存在的数组元素。这通常发生在你尝试访问数组中的第n + 1个元素，但实际上只有n个元素。

这个问题可能有多个原因，包括输入数据格式的问题、数组大小定义的问题等。因此，我需要看到你的代码才能确定问题所在。

该回答引用ChatGPT
尝试下面代码

function dy=eg3(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