分岔现象matlab,Matlab用来计算分岔图的程序

【实例简介】

【实例截图】

【核心代码】

function mian6 % 分岔图程序

tic

N=256;

M=200;

% M=100;

w=[150:1:2500];

L=length(w);

%matlabpool(6)

RelTol = 1e-6;%相对误差Relative tolerance

AbsTol = 1e-6;%绝对误差Absolute tolerance

options = odeset('RelTol',RelTol,'AbsTol',ones(1,24)*AbsTol);

for j=1:L

j

% [t,Y] = ode45( @gongshi,[0:(2*pi/M):(288*2*pi)],[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0],w(j));

[t,Y] = ode45( @gongshi,[0:(2*pi/w(j)/M):288*2*pi/w(j)],[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0],options ,w(j));

y1=Y(:,1);

y11=y1((end-N*M):end);

y3=Y(:,3);

y33=y3((end-N*M):end);

for i=1:N

Point1(i,j)=y11(1 (i-1)*M);

Point2(i,j)=y33(1 (i-1)*M);

end

end

w = w';

plot(w,Point1,'.k','MarkerSize',2)

figure

plot(w,Point2,'.k','MarkerSize',2)

save Point1

save Point2

%matlabpool close

toc