用matlab画分叉图

本例还是以达芬振子为例，具体实现过程如下：

第一，建立一个达芬振子议程的M文件，其中force为参数，本例题是求出参数递增时的分叉图型

function df=dafen(t,x,flag,force)
df=[x(2);force*cos(1.2*t)-x(1)^3+x(1)-0.3*x(2)];

第二，画分叉图M文件

clear
ff=0.23:0.002:0.40;
axis([0.22 0.4 0.2 0.6])
options=odeset('RelTol',1e-7);
tt=2*pi/1.2;
hold on
for j=1:length(ff)

    [t,x]=ode45(@dafen,[0:tt/100:80*tt],[0,0],options,[],ff(j));
    i=2000:100:3000
    plot(ff(j),x(i,2),'*')
end