matlab求解分岔图：Error using odearguments

王瑛

根据论坛前辈关于duffing方程分岔图的程序，自己变了一个，但是出现一些错误。错误和源程序如下、希望可以得到大家的帮助。谢谢！

Error using odearguments (line 36)
No default parameters in MASS_SPRING. Use ode45(MASS_SPRING,tspan,y0,...) instead.


Error in ode45 (line 114)
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...


Error in mass_spring_bifur_k (line 20)
[t,Y]=ode45('mass_spring',tspan ,y0);



相关程序如下：

1. %微分方程的程序%
   
2. function dy=mass_spring(t,Y,flag,b,D,g,o,rum,cd,pf,vf)
   
3. global k v p rw I;
   
4. b=0.0028;g=9.807;o=5*pi/6;um=1.256*10^(-6);cd=0.44;pf=0.06;vf=100;rr=1.22;
   
5. D=1/4+1/(1-cos(o))-2/(1-cos(o))^2*log(2/(1+cos(o)));
   
6. dy=[y(2);(y(3)*g+um*(I^2)/4/pi*(log((3*y(3)/4/pi/p)^(1/3)*sin(o)/rw)-D)+cd*pi*((3*y(3)/4/pi/p)^(2/3)-rw^2)*pf*(vf^2)/2-b*y(2)-k*y(1))/y(3);v];
   
7. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   
8. 
   
9. %求解分岔图的程序%
   
10. clear;
    
11. global k v  p rw I;
    
12. T1=0.1347*10^(-3);
    
13. T2=1.413*10^(-5);
    
14. le=0.026;
    
15. rw=0.8*10^(-3);
    
16. VM=p*pi*rw^2*(T1*I+T2*le*I^2);
    
17. v=VM;
    
18. range=[2.7: 0.001: 3.1];
    
19. period=2*pi./v;
    
20. kk=0;
    
21. YY1=[];
    
22. step=2*pi/500;
    
23. for a=range
    
24.   k;
    
25. kk=kk+1;
    
26. %出去前面60个周期的数据，并将最后的结果作为下一次积分的初值。
    
27. tspan=[0:step:60*period];
    
28. y0=[0 0 3e-5];
    
29. [t,Y]=ode45('mass_spring',tspan ,y0);
    
30. y0=Y(end,:);
    
31. j=1;
    
32. for i=60:200
    
33. tspan=[i*period:step:(i+1)*period];
    
34. [t,Y]=ode45(@mass_spring,tspan ,y0,options,b,D,g,o,rum,cd,pf,vf);
    
35. YY1(kk,j)=Y(end,1);
    
36. j=j+1;
    
37. %取出每一个周期内的第一个解的最后一个值。
    
38. y0=Y(end,:);
    
39. end
    
40. end
    
41. bifdata=YY1(:,end-51:end);
    
42. plot(range,bifdata,'k.','markersize',1);

复制代码

ChaChing

LZ才两天在线时间已经2hr了, 表示LZ已先搜/看过! 值得嘉许
个人水平/时间有限, 也非此专业又很少使用ode经验, 但LZ资料齐全
就试试设断点帮看看
不难发现ode45的输入引数tspan为空矩阵
逐步回溯检查(tspan > period > v > VM > I & p)
不难发现可能I及p未给定!?

王瑛

ChaChing 发表于 2013-6-19 21:56

                               
登录/注册后可看大图

LZ才两天在线时间已经2hr了, 表示LZ已先搜/看过! 值得嘉许
个人水平/时间有限, 也非此专业又很少使用ode经 ...

非常感谢！因为也是刚解除matlab没多久，确实是有的参数没有输入。这个问题解决了，其他问题又出现了。我再好好检查检查。如有需要解决的问题，还希望能和大家多交流。谢谢chaching

张某某zyl

最近正在学习

张某某zyl

这是用的庞加莱截面法？