分岔图是一条直线?
离散系统由于模型太大,程序不贴了,就说明一下思路:
第一层for 循环
分岔参数;
第二层for循环
离散映射求解,大概求解2k个点;
然后记录最后的200个点;
end
end
画图部分:
用一层for 循环
分岔参数;
plot(分岔参数i,分别对应记录的200个点)
end
另外说明:我也试过画计算2w个点,取最后的2k个点,结果和上述情况下是差不多的,也是直线的形状,对局部放大部分,和图2一样。
补充编辑:
我把纵坐标 取对数之后,(有负数的,先转正数,取完对数再转负数),如图三所示。
这个能说明什么?
[ 本帖最后由 purple_paradise 于 2010-4-27 14:47 编辑 ] 另外的问题:
画出分岔图和对应的LE图,如下。
关于LE图画法:用定义。由于四维系统雅克比矩阵只含有未知数“分岔参数”,因此,每个分岔参数下,只求出4个特征值。然后把四个特征值取模,lyapunov1= sum(log(abs(eig(J))))/4(对应图2);或者是:lyapunov2=max(log(abs(eig(J))))(对应图3);
此方法求解是否正确?
[ 本帖最后由 purple_paradise 于 2010-4-28 10:14 编辑 ]
页:
[1]