ÎÒ»·Ö֧ͼ³£³£ÏȲî·Ö£¬È»ºóÔËÐкܳ¤Ê±¼äÈ¡×îºóµÄ¼¸°Ù²½¾Í¿ÉÒÔÁË£¬ÀýÈçÒÔÏÂ
clear;
T=1000;dt=0.05;a1=5.0;a2=0.1;b2=2.0;d1=0.4;d2=0.01;e=1.0;r=0;
X=zeros(1,T);Y=zeros(1,T);Z=zeros(1,T);
X(1)=0.4;Y(1)=0.6;Z(1)=5;
hold on
for b1=2:0.02:9
for t=2:T
X(t)=X(t-1)+(X(t-1)*(1-X(t-1))-a1*X(t-1)*Y(t-1)/(1+b1*X(t-1)))*dt;
Y(t)=Y(t-1)+(a1*X(t-1)*Y(t-1)/(1+b1*X(t-1))-a2*Y(t-1)^(e)*Z(t-1)/(1+b2*Y(t-1)^(e)+r*Z(t-1))-d1*Y(t-1))*dt;
Z(t)=Z(t-1)+(a2*Y(t-1)^(e)*Z(t-1)/(1+b2*Y(t-1)^(e)+r*Z(t-1))-d2*Z(t-1))*dt;
end
plot(b1,Z(T-500:T),'k.','MarkerSize',4)
end
µ«ÕâÑù×öµÄ½á¹ûÊÇÖÜÆÚºÍ»ìãçÎÞ·¨·Ö±æ£¬ÓÐÈËÓü«ÖµÃüÁî×ö£¬ºÃÏñ¿ÉÒÔ½â¾öÕâ·½ÃæµÄÎÊÌâ¡£
Çë½ÌһϸßÊÖ£¬°ïÎÒ°ÑÉÏÃæ³ÌÐò¡®plot(b1,Z(T-500:T),'k.','MarkerSize',4) ¡¯Õⲿ·ÖÓü«´óÖµÃüÁî´úÌæÒ»Ï£¬Íò·Ö¸Ðл£¡ |