Drexel dragonThe Math ForumDonate to the Math Forum



Search All of the Math Forum:

Views expressed in these public forums are not endorsed by Drexel University or The Math Forum.


Math Forum » Discussions » Software » comp.soft-sys.matlab

Topic: why I am not getting a bifurcation diagram?
Replies: 3   Last Post: Feb 20, 2013 1:20 PM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
Sanaa

Posts: 154
Registered: 3/20/12
why I am not getting a bifurcation diagram?
Posted: Feb 20, 2013 7:20 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

Hi,
I am trying to plot a bifurcation diagram for Rossler discrete-time system.
What's wrong in my code?
I am only getting a line!!!!



itermax=300;min=itermax-9;% That is we plot from 291:300 i.e. 10 values of x only.
r=0.1;b=0.2;c=4;
for a=0:0.001:0.4
x=0.1;
x0=x;
y=0.2;
y0=y;
z=0.3;
z0=z;
for n=1:itermax
xn=x0-r*(y0+z0);
yn=y0+r*(x0+a*y0);
zn=z0+r*(b+x0*z0-c*z0);
x=[x xn];
x0=xn;
y=[y yn];
y0=yn;
z=[z zn];
z0=zn;
end
plot(a*ones(10),y(min:itermax),'b.','linewidth',2)
hold on
end
xminmax = xlim;
line(xminmax, [0 0], 'Linewidth', 1,'Color','k')
fsize=15;
xlabel('a','FontSize',fsize)
ylabel('\ity','FontSize',fsize)
hold off
%print(gcf, '-djpeg', '-zbuffer', 'bif01_standard.png');


Any help will be highly appreciated.
Thanks



Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© Drexel University 1994-2014. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.