lnghng 发表于 2007-11-27 17:44

求教高手：关于单摆的周期解与混沌解的问题

我不太了解混沌理论中的周期解与混沌解到底是哪各种图。附件是我用matlab作出的单摆的相图，其中阻尼系数是0.25。w/wo是2/3,f=F/mg等于各张图的文件名。请问各位大侠，每张图对应哪种解？？谢谢了。:@L

lnghng 发表于 2007-11-27 17:46

麻烦说明一下周期解是几倍周期解，还有原因？？谢谢了，好人啊。

咕噜噜 发表于 2007-11-27 17:53

相图看不太出来吧
要是一定要说，可能1.2有点像混沌
你还是画poincare映射看看吧

lnghng 发表于 2007-11-27 21:02

我查过一篇文献，上面说1.078是四倍周期解但没有给出图形。poincare映射是什么东西？我不清楚，我现在才开始接触这些东西，麻烦讲详细一点。

咕噜噜 发表于 2007-11-28 09:06

poincare映射你可以搜索论坛找详细解释和计算方法
一般判断是否混沌都是poincare映射，相图和李亚诺普指数结合求解得

无水1324 发表于 2007-11-28 09:57

1、你是用ode45求解的吗？
2、1.45时怎么图像会自相交，这是相图还是？

lnghng 发表于 2007-11-28 13:10

是用ode45作的。是相图

咕噜噜 发表于 2007-11-28 13:49

1.45是相图？好奇怪的相图
能不能提供你的程序上来

无水1324 发表于 2007-11-28 22:17

回复 #7 lnghng 的帖子

你的相图好像一直沿着横坐标变化的哈

lnghng 发表于 2007-11-29 18:37

帮我看看哪里有错。谢谢了
function f=pendulum3(t,y)
globalA;
f=;
主程序：
global A
Aa=input('输入A=');
t=input('输入求解区间t=');
y0=input('输入初值y0=');
for n=1:length(Aa)
    A=Aa(n);
    =ode45('pendulum3',t,y0);
    figure;
    plot(y(500:end,1),y(500:end,2 )
   end

咕噜噜 发表于 2007-11-29 19:17

应该是吧
另外，你这么循环有什么意义吗

lnghng 发表于 2007-11-29 21:09

,还不是一样的吗？循环是为了连续画一系列A的值的图形。

咕噜噜 发表于 2007-11-29 21:26

可你的循环并没有实现你说的功能
页: [1]
查看完整版本: 求教高手：关于单摆的周期解与混沌解的问题