帮忙鉴定下这个分岔图以及程序

zebtra_stripe

帮忙鉴定下这个分岔图以及程序
我的疑问是这个程序作出的分岔图纵坐标是什么含义？（当然是程序如果正确的前提下）
提前谢谢各位高手指教！不胜感激！

clc;
clear;
Z=[];
index2=1000;%已经有500步运行了。

SLRate=.84;%供货调整率
for InventoryRate=linspace(0.5,1,1000);
Y=supplychain(InventoryRate,SLRate,index2);%调用supplychain函数
  %plot3(Y(1000:end,1),Y(1000:end,2),Y(1000:end,3),'-'),grid on;
%归一化
for m=1:3
Max(1,m)=max(Y(:,m));
Y(:,m)=Y(:,m)/Max(1,m);
end
% 舍弃前面迭带的结果，用后面的结果画图
  %取一个截面，用直线穿过截面，获取前后点，取出穿过截面的点
Y(:,1)=Y(:,2)-Y(:,1);
% 对计算结果进行判断，如果点满足x=y，则取点
  for k=2:length(Y)
    f=k-1;
    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 InventoryRate+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 InventoryRate+abs(y)*i];
       end
     end
   end
end
plot(Z,'.','markersize',1)
title('SupplyChain分岔图')
xlabel('InventoryRate'),ylabel('|y| where x=y')

zebtra_stripe

好像几位版主和高手都很忙啊
小弟真的很急需知道这个纵坐标的含义，还请高人提点。

平平

这个我也不明白，我是新手，请多指教

yiquncw

个人感觉这个绘图方法有问题，有些系统不适用，我也用这个方法做过分岔图，效果不好，

zebtra_stripe

高手，那你还有其他的分岔图方法的程序么？
给小弟一份咯。问了这么长时间，也没有回复，更没人知道这纵坐标是什么意思，:'(
正发愁呢。。。:handshake

无水1324

纵坐标的含义：其实就是Poincare截面点

无水1324

to：zebtra_stripe，yiquncw， 平平
作分岔图，第一步你必须把Poincare图搞清楚，然后随参数变化的时候把 Poincare图画出来就是分岔了

zebtra_stripe

受教了。:@)
难怪我觉着庞加莱的程序和分岔图极其相似。想着有点什么关系，就是不大明白具体是什么关系。
但是我还有点疑问就是，庞加莱截面只是去截，也没有规定什么x=y,为什么这里要说明"如果点满足x=y，则取点"
也许问题愚昧，但请不吝赐教！谢谢

无水1324

这里面应该是人为设置了x=y这个截面。所以才要满足这样子的要求，你改为x=2*y+0.1，这种形式也可以，这里面的参数是任意的

laozhou

我画的poincare截面都是对微分方程的数值解以一个周期为步长进行插值截点，效果不错！

htwei

我差不多 抛弃开头一部分瞬态值后 每隔一个周期取一个点  “频闪” 我有个疑问 你的插值是什么意思

[ 本帖最后由 htwei 于 2010-3-3 10:32 编辑 ]

ypp207

能否给个例子。如果用poincare截面来做，截面应该多试几次吧，搞不好与相轨迹相切，你这个系统也是周期激励的？有些疑问，希望详细介绍一下！

无水1324

插值 ,应该是隔一个周期取点，不一定刚好在一个整周期点上。对于选取任意截面的更加有这个可能，
另外ypp207说的是选取截面的一个基本条件。就是截面不能与轨迹相切

purple_paradise

也问一下，对于 多维系统离散，比如4维的，那相图怎么画出来，正常看到的2维的可以（x,y），三维的可以（x,y,z），那更高维的呢？
一个是相图，还有就是分岔图和庞加莱图

无水1324

作二维的图，你就选择两个量化图就好了，不管系统是多少维德