分岔图是一条直线？

purple_paradise

离散系统
由于模型太大，程序不贴了，就说明一下思路：

第一层for 循环
     分岔参数；
    第二层for循环
        离散映射求解，大概求解2k个点；
        然后记录最后的200个点；
   end
end

画图部分：

用一层for 循环
     分岔参数；
    plot(分岔参数i，分别对应记录的200个点)
end

另外说明：我也试过画计算2w个点，取最后的2k个点，结果和上述情况下是差不多的，也是直线的形状，对局部放大部分，和图2一样。


补充编辑：
我把纵坐标 取对数之后，（有负数的，先转正数，取完对数再转负数），如图三所示。
这个能说明什么？

[ 本帖最后由 purple_paradise 于 2010-4-27 14:47 编辑 ]

purple_paradise

另外的问题：
画出分岔图和对应的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 编辑 ]