混沌动力学分岔图，代码有误，求教~
function s=der2(t,y,Z,flag,f);ylag1=Z(:,1);
s=[493.8*y(2)*(y(1)+1)-1.2*y(1)*(y(1)+1);
1.458*4*f/3*(1+ylag1(1))-1.458*y(2)-0.486*y(1)-y(2)*(y(1)+1)+0.00243*y(1)*(y(1)+1)];
f=0:0.001:0.2;
hold on
for i=1:length(f)
sol=dde23(@der2,,[-0.1,-0.1],,[],f(i));
y1=sol.y(1,:);
plot(f(i),y1(1:10:end),'k.')
end
我想画的是关于反馈系数的分岔图，模拟了一个例子，但是报错，有没有人知道怎么改？
页:
[1]