倍周期到混沌的分岔图

kuzb

我还是想问下，图示的倍周期到混沌的图用matlab怎么做出来的，能给个matlab程序链接吗？我之前发过一个帖子，关于这个的。但是用网上的常见的matlab程序做不出来我之前的那个系统的图，不知道什么原因。如果大家做过这个图，能分享一下自己的系统和对应的程序吗？重要的是，系统不是映射，而是关于连续时间的系统，谢谢

detectivegin

我都是自己用origin画的，操作起来很费劲啊

detectivegin

同求程序

kuzb

detectivegin 发表于 2017-3-20 22:22
我都是自己用origin画的，操作起来很费劲啊

你用origin也需要数据啊，你数据不是编程算的？

Eminem

matlab画的比较好

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=;

复制代码

kuzb

Eminem 发表于 2017-3-21 08:43
matlab画的比较好

谢谢，我先试试

detectivegin

kuzb 发表于 2017-3-20 22:42
你用origin也需要数据啊，你数据不是编程算的？

自己编程算的轴心轨迹，然后做这种图

ljq519781936

同求程序

失心控

楼上不是已经给程序了吗  还需要自己动手编的  不能都捡现成的