Matlab。试图访问y(4); index超出界限因为numel(y)= 3

时间:2013-05-30 22:28:47

标签: matlab equation-solving matlab-compiler

我的matlab ode45有问题。当我写系统y(1)和y(2)时,只有它可以工作,但是当我使用更多时,它会说 试图访问y(4); index超出界限因为numel(y)= 3 这是我的程序的代码,请检查它 `

function ode_epidemic(  )
global alfa beta hama delta lambda myu
global nfunc
nfunc=0;
alfa=1; beta=1; hama=1; delta=1; lambda=1; myu=1;
options = odeset('RelTol',1e-5);
sol = ode45(@epidemic,[0 20],[5; 0.1; 0],options);
figure
hold on;
xlabel('x','FontSize',11)
ylabel('y','FontSize',11)
plot(sol.x,sol.y(1,:),'-',...
     sol.x,sol.y(2,:),'-.',...
     sol.x,sol.y(3,:),'.','LineWidth',2)
hold off;
grid 'on'
disp('number of iteration')
iter=length(sol.x)
disp('number calculation of F')
nfunc
iter1=iter-1;
h=zeros(iter1,1);
x=zeros(iter1,1);
for i=1:iter1
 h(i)=sol.x(i+1)-sol.x(i); 
 x(i)=sol.x(i);
end
figure
hold on;
xlabel('x');ylabel('y');
plot(x,h,'LineWidth',2);
hold off;
grid 'on'
end

function dy = epidemic(t,y)
global alfa beta hama delta lambda myu
global nfunc
nfunc=nfunc+1;
dy=zeros(4,1);
dy(1)=2*t*y(1)*y(4);
dy(2)=10*t*y(1)^5*y(4);
dy(3)=2*t*y(4);
dy(4)=-2*t*(y(3)-1);
end

`

1 个答案:

答案 0 :(得分:1)

您提供的y矢量只有3个元素。

替换:

ode45(@epidemic,[0 20],[5; 0.1; 0;],options);

ode45(@epidemic,[0 20],[5; 0.1; 0; 0],options);

并且代码运行没有错误。