zebtra_stripe 发表于 2009-12-26 09:21

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

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

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=;
       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=;
       end
   end
   end
end
plot(Z,'.','markersize',1)
title('SupplyChain分岔图')
xlabel('InventoryRate'),ylabel('|y| where x=y')


zebtra_stripe 发表于 2009-12-29 16:04

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

平平 发表于 2009-12-30 22:48

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

yiquncw 发表于 2010-1-4 13:53

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

zebtra_stripe 发表于 2010-1-4 15:32

回复 地板 yiquncw 的帖子

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

无水1324 发表于 2010-1-5 08:14

回复 沙发 zebtra_stripe 的帖子

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

无水1324 发表于 2010-1-5 08:16

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

zebtra_stripe 发表于 2010-1-5 11:07

回复 7楼 无水1324 的帖子

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

无水1324 发表于 2010-1-5 15:04

回复 8楼 zebtra_stripe 的帖子

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

laozhou 发表于 2010-1-13 14:31

回复 9楼 无水1324 的帖子

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

htwei 发表于 2010-3-3 10:30

回复 10楼 laozhou 的帖子

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

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

ypp207 发表于 2010-4-13 12:00

回复 10楼 laozhou 的帖子

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

无水1324 发表于 2010-4-13 16:55

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

purple_paradise 发表于 2010-4-26 12:07

回复 13楼 无水1324 的帖子

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

无水1324 发表于 2010-4-26 16:59

回复 14楼 purple_paradise 的帖子

作二维的图，你就选择两个量化图就好了，不管系统是多少维德
页: [1] 2
查看完整版本: 帮忙鉴定下这个分岔图以及程序