求一个延迟的分岔程序

siyaoming

做了好久，没做出来。不知院长们能不能抽出时间来帮帮忙。是这样的。
x1'=x2
x2'=-wo^2*x1-a*x2-a*b*x1^2*x2-w0^2*r*x1^3+f*cos(omega*t)+g1*x1(t-tao1)+g2*x2(t-tao2)
其中wo,a,b,r,f,omega,g2,g1,tao1,tao2都是参数，设置f是变化的其他都是固定的。做关于f变化x2,或者是x1的分岔图。即f为横坐标，x2或者是x1为纵坐标。望院长给出一合适程序。
是用dde23做吧？
在论坛上搜了好久也没有找到类似的程序，要是原创自己也不会，只能希望跟各位院长以及高手学习学习。

无水1324

那我问你一下，对于这个系统你会写出求解的程序吗？
只是要求解方程的。你贴出来，我帮你修改程序

siyaoming

clf
w0=1;a=-1;b=-0.3;r=0.2;omega=1.1;g1=1;g2=1;f=0.1:0.1:5;
lags=[pi/2,2*pi/3];T=2*pi/omega;
y1=[];G=[];y2=[];
for h=1:length(f)
sol=dde23(@yanchi,lags,[0.1 0],[0:T/100:40*T],f(h));
y1=[y1,sol.y(1,2000:100:end)];y2=[y2,sol.y(2,2000:100:end)];
m=length(y1);
for j=1:m/h
G=[G,f(h)];   
end;p=length(G);
N=length(f);
Over_to_step=N-i;
end
figure(1)
plot(G,y1,'k.','markersize',3)
hold on


function dy=yanchi(t,y,Z,f)
w0=1;
a=-1;
b=-0.3;
r=0.2;
omega=1.1;
g1=1;
g2=1;
f=0.1:0.1:5;
dy=[y(2);-w0^2*y(1)-a*y(2)-a*b*y(1)^2*y(2)-w0^2*r*y(1)^3+f*cos(omega*t)+g1*Z(1,1)+g2*Z(2,2)]

无水1324

Ｚ是多少？是时间延迟吗？

siyaoming

这是延迟的表达形式，Z(1,1)表示第一行的延迟是主程序中lags=[pi/2,2*pi/3];中的pi/2,Z(2,2)以此类推。

无水1324

真的很抱歉，我试了好长时间，总是不对，最后我去借了本书，看看，才看到这是个中性方程，它是不能用
dde23的，建议使用simulink去做！

yzsldj

是研究分岔问题吧，应该用DDEBIFTOOL啊

siyaoming

感谢大家的努力，我再看看吧

无水1324

确实不好意思啊，低估了这个问题，你真的可以考虑一下用 7# 的工具箱吧