siyaoming 发表于 2009-8-18 16:00

求一个延迟的分岔程序

做了好久，没做出来。不知院长们能不能抽出时间来帮帮忙。是这样的。
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 发表于 2009-8-18 22:00

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

siyaoming 发表于 2009-8-19 09:04

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=;T=2*pi/omega;
y1=[];G=[];y2=[];
for h=1:length(f)
sol=dde23(@yanchi,lags,,,f(h));
y1=;y2=;
m=length(y1);
for j=1:m/h
G=;   
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=

无水1324 发表于 2009-8-19 09:40

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

siyaoming 发表于 2009-8-19 09:54

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

无水1324 发表于 2009-8-20 08:03

回复 5楼 siyaoming 的帖子

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

yzsldj 发表于 2009-8-20 08:38

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

siyaoming 发表于 2009-8-20 16:08

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

无水1324 发表于 2009-8-20 16:34

回复 8楼 siyaoming 的帖子

确实不好意思啊，低估了这个问题，你真的可以考虑一下用 7# 的工具箱吧
页: [1]
查看完整版本: 求一个延迟的分岔程序