声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2766|回复: 9

[分形与混沌] 哪位大神帮我看看我的分岔图的程序哪出错了啊

[复制链接]
发表于 2014-2-21 20:05 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
本帖最后由 牛小贱 于 2014-7-1 13:21 编辑

分岔图程序如下:
  1. clear;clc;clf
  2. global f
  3. range=linspace(5,20,100);
  4. for f=range
  5. y0=[0.001 0.01 0.001 0.010]';
  6. period=(2*pi)/f;
  7. step=period/100;  
  8.    tspan=[0:step:100*period];
  9.   [t,x]=ode45('Driver_vehicle1',tspan,y0);
  10.   plot(f,x(5000:100:end,1),'k.','markersize',10);
  11. hold on
  12. end
复制代码
跑出来的图根本就不是分岔图啊,我测试了几个f=5,8,10三种情况,庞加莱映射图都是几个点,但我跑的分岔图每个f对应的只是一个单一点,跪求大神解决啊。。。
QQ图片20140221200245.jpg
回复
分享到:

使用道具 举报

 楼主| 发表于 2014-2-21 22:48 | 显示全部楼层
求助啊。。。我用这段程序跑其他典型混沌能跑出来分岔图,但是我的自己的微分方程就是出不来啊。。。。。
发表于 2014-3-6 09:43 | 显示全部楼层
没头没尾的一个问题,没办法回答

系统没给,庞加莱映射图如何获取没给,庞加莱映射图结果什么样没有
发表于 2014-5-13 20:09 | 显示全部楼层
根据,以前的帖子,你好像,没取到周期点
发表于 2014-5-17 06:37 | 显示全部楼层
kuri2012 发表于 2014-5-13 20:09
根据,以前的帖子,你好像,没取到周期点
  1. plot(f,x(5000:100:end,1),'k.','markersize',10);
复制代码


是取得周期点
 楼主| 发表于 2014-6-30 20:34 | 显示全部楼层

我取的是周期点啊,100个点是一个周期。。
发表于 2014-7-2 19:56 | 显示全部楼层
我恨李雅普诺夫 发表于 2014-6-30 20:34
我取的是周期点啊,100个点是一个周期。。

那有可能你的系统没有出现分叉,只有稳定的单周期
发表于 2014-7-4 19:35 | 显示全部楼层
试着把初始条件改变一下  试试看   
发表于 2014-7-5 18:44 | 显示全部楼层
我恨李雅普诺夫 发表于 2014-2-21 22:48
求助啊。。。我用这段程序跑其他典型混沌能跑出来分岔图,但是我的自己的微分方程就是出不来啊。。。。。

如果用你自己的方程能画出来Poincare映射图可出不来分岔图那就改变初值和参数试试
发表于 2014-7-8 09:23 | 显示全部楼层
{:{44}:}
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-23 15:45 , Processed in 0.061414 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表