I'm trying to plot the orbit diagram of the logistic map, but the plot doesn't show anything. There are no error messages, so I'm stuck...
close;
clear;
orbit=zeros(1,300);
j=0;
% parameter range
for(r=2.8:0.001:4)
j=j+1;
% random nitiation of iteration
xn1=rand(1);
for(i=1:600)
% calculate logistic map
xn=xn1;
xn1=r*xn*(1-xn);
% wait for transients
if(i>300)
% store the orbit points
orbit(i-300)=xn1;
end
end
plot(r,orbit);
if(j==1)
axis([2.8 4 0 1]);
hold;
end
end

 채택된 답변

Star Strider
Star Strider 2019년 12월 8일

0 개 추천

Try this slight edit of your code:
orbit=zeros(1,300);
j=0;
% parameter range
rv=2.8:0.001:4;
for k = 1:numel(rv)
r = rv(k);
j=j+1;
% random nitiation of iteration
xn1=rand(1);
for i=1:600
% calculate logistic map
xn=xn1;
xn1=r*xn*(1-xn);
% wait for transients
if(i>300)
% store the orbit points
orbit(i-300,k)=xn1;
end
end
% plot(r,orbit,'.');
if(j==1)
axis([2.8 4 0 1]);
hold;
end
end
figure
plot(rv, orbit, '.');
axis([2.8 4 0 1]);

추가 답변 (0개)

카테고리

도움말 센터File Exchange에서 MATLAB에 대해 자세히 알아보기

제품

릴리스

R2019b

태그

질문:

2019년 12월 8일

답변:

2019년 12월 8일

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by