bifurcation2dmap.m

利用Matlab求二维离散系统的分岔图程序。 在非线性动力学及复杂经济学中经常用到类似的问题。也可以用来画其他更高维数的离散映射系统的分岔图。很实用。

%这个小程序虽然很简单，但可以求非线性动力学或者复杂经济学中经常遇到的类似问题。
%也可以求更高维数的离散系统的bifurcation diagram. 可以根据研究需要，选择分岔图中点的颜色和大小。
clear;
   clf;
   r=linspace(0,0.5,4000);
   t=200;
   j=1;
   while j<=4000
       p(1)=2;q(1)=1;i=1;
       while i<=(t-1)      
       p(i+1)=p(i)+0.2*p(i)*[10-2*(0.5+4/3)*p(i)-0.5*q(i)];
       q(i+1)=q(i)+r(j)*q(i)*[10-2*(0.5+1)*q(i)-0.5*p(i)];
       i=i+1;
       end;
      a(j)=p(t);
      b(j)=p(t-1);
      c1(j)=p(t-2);
      c2(j)=p(t-3);
      c(j)=q(t);
      d(j)=q(t-1);
      d1(j)=q(t-2);
      d2(j)=q(t-3);
      j=j+1;
  end;
  plot(r,a,'k.',markersize,1);hold on;
  plot(r,b,'k.',markersize,1);hold on;
  plot(r,c1,'k.',markersize,1);hold on;
  plot(r,c2,'k.',markersize,1);hold on;
  plot(r,d1,'k.',markersize,1);hold on;
  plot(r,d2,'k.',markersize,1);hold on;
  plot(r,c,'k.',markersize,1);hold on;
  plot(r,d,'k.',markersize,1)