分岔图感觉不太对，有没有帮忙看看。

Overdosed

最近在练画分岔图，逛了一段时间的论坛找到了几种方法，自己画了一个练练手，感觉这个图有点怪啊。在【1.5~2】这个区间内由单周期直接进入了混沌区域，进入混沌的途径不是应该经过倍周期分岔再进入混沌么？
各位大佬们求解答。必要的话可以把程序贴上来。

kuzb

你是用网上或者论坛里面画分岔图的程序做的吗？

jiangwj

你这图明显不对啊

Overdosed

jiangwj 发表于 2017-1-16 08:26
你这图明显不对啊

我也觉得这个怪怪的，也不知道到底是哪里不对。试了论坛里的两个方法，分岔图的差别太大了。

jiangwj

1. function AA
   
2. 
   
3. global mu w F;
   
4. mu=0.3;w=1.2;
   
5. cF=;
   
6. k=0;
   
7. for F=cF
   
8.     k=k+1;
   
9. =ode45(@duffing,,);
   
10. x1=x(:,1);
    
11. count=find(t>1800);
    
12. Y=x1(count,:);
    
13. t1=t(count,:);
    
14. subplot(5,2,k)
    
15. plot(t1,Y)
    
16. xlabel('t');
    
17. ylabel('x ');
    
18. s=;
    
19. title(s);
    
20. end
    
21. 
    
22. function dx=duffing(t,x)
    
23. global mu w F;
    
24. dx=;
    
25. 杜芬(Duffing)方程 Poincare 截面图
    
26. function aa
    
27. global mu w F;
    
28. mu=0.05;w=1;
    
29. F=7.5;
    
30. x0=;
    
31. tspan = 0:0.1*2*pi:50000;
    
32. =ode45(@duffing,tspan,x0);
    
33. m=zeros(7000,1);
    
34. n=zeros(7000,1);
    
35. for i=1:7000
    
36.     m(i,1)=y(10*i,1);
    
37.     n(i,1)=y(10*i,2);
    
38. end
    
39. plot(m,n,'k.','markersize',1);
    
40. function dx=duffing(t,x)
    
41. global mu w F;
    
42. dx=;

复制代码

我这有一个  你可以试试

Overdosed

jiangwj 发表于 2017-2-20 09:05
我这有一个  你可以试试

这两个图的结果是差不多的。难道是我初值取的不好？