分数阶系统的分岔图怎么画

anqikeli-001

我看了“分岔图绘制不同方法的总结、比较”的帖子，里面的Poincare截面法的程序没有看懂，
[T,Y]=ode45('Lorenz',[0,1],[1;1;1;16;r;4]);  
    [T,Y]=ode45('Lorenz',[0,50],Y(length(Y),:));
上面的第二句是为了得到什么？
而且if Y(k,1)<0   
            if Y(f,1)>0
                y=Y(k,2)-Y(k,1)*(Y(f,2)-Y(k,2))/(Y(f,1)-Y(k,1));
                Z=[Z r+abs(y)*i];
            end  
        else     
            if Y(f,1)<0  
                y=Y(k,2)-Y(k,1)*(Y(f,2)-Y(k,2))/(Y(f,1)-Y(k,1));   
                Z=[Z r+abs(y)*i];   
            end  
好像有点问题吧？？
我现在要做的是分数阶的系统的分岔图，可以用这个方法吗，参数是阶数，希望各位师兄师姐指教！

无水1324

其实第二个就是取第一次求解的值为初值
这样可以有效地消除瞬态响应，不知道这样说你理解不

octopussheng

分数阶系统如何求解啊？没有接触过！

无水1324

只要解可以求出来了，使可以按照这个方法做分岔分析的

xinwilliam

无水将其说到根上了。实际上，画分岔图就是将解描在图上，只不过通过图可以发现分岔现象罢了。

21172485

分数阶的有个台湾人在搞