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

lnghng

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

lnghng

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

咕噜噜

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

lnghng

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

咕噜噜

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

无水1324

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

lnghng

是用ode45作的。是相图

咕噜噜

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

无水1324

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

lnghng

帮我看看哪里有错。谢谢了
function f=pendulum3(t,y)
global  A;
f=[y(2);A*cos(2/3*t)-sin(y(1))-y(2)/2];
主程序：
global A
Aa=input('输入A=');
t=input('输入求解区间t=[t0 tf]');
y0=input('输入初值y0=[y01 y02]');
for n=1:length(Aa)
    A=Aa(n);
    [x,y]=ode45('pendulum3',t,y0);
    figure;
    plot(y(500:end,1),y(500:end,2 )
   end

咕噜噜

[x,y]应该是[t,y]吧
另外，你这么循环有什么意义吗

lnghng

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

咕噜噜

可你的循环并没有实现你说的功能