马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
本帖最后由 牛小贱 于 2014-7-1 13:21 编辑
分岔图程序如下:
- clear;clc;clf
- global f
- range=linspace(5,20,100);
- for f=range
- y0=[0.001 0.01 0.001 0.010]';
- period=(2*pi)/f;
- step=period/100;
- tspan=[0:step:100*period];
- [t,x]=ode45('Driver_vehicle1',tspan,y0);
- plot(f,x(5000:100:end,1),'k.','markersize',10);
- hold on
- end
复制代码 跑出来的图根本就不是分岔图啊,我测试了几个f=5,8,10三种情况,庞加莱映射图都是几个点,但我跑的分岔图每个f对应的只是一个单一点,跪求大神解决啊。。。
|