Date: Jul 12, 2013 11:50 AM
Subject: Re: Incomplete diagram!
"Sanaa" wrote in message <firstname.lastname@example.org>...
> I have a problem that I wish to plot what is called a bifurcation diagram.
> The map under study is given by
> x_((n+1)r)= 1 - ru *x^2_(nr),
> where r is a positive constant may take 0.1,0.2,...
> ru is a positive parameter, and the time t runs from nr to (n+1)r, n=0,1,2,...
> I have the matlab code
> % define the vector of values that r takes
> r_vec = [0.2];
> % define the vector of colors, to plot the dat aof reach value of r in
> % different color
> color_vec = ['b'];
> for rho = 0: 0.005: 2
> % do a loop of values of r
> for j = 1 : size(r_vec,2)
> % define delay
> r = r_vec(j);
> % define coefficient rho
> %rho = 0.5;
> % define the number of deiscrete times in interval [k*r, (k+1)*r]
> n = 10;
> % define the nimber of iterations k = 1,.....,Nit
> Nit = 1000;
> % define the initial vector which is a vector of size n
> x0 = 0.3*ones(1,n);
> % define vectors x_next and x_previous
> x_next = zeros(1,n);
> x_previous = zeros(1,n);
> time = zeros(1,n);
> % initialize x_previous
> x_previous = x0;
> x_next = x_previous;
> time = linspace(0,r,n);
> % set the figure details
> hold on
> for i = 1:Nit
> x_next = 1 - rho*x_previous.^2;
> x_previous = x_next;
> time = linspace(i*r,(i+1)*r,n);
> plot(rho,x_next(1),'.', 'Color', char(color_vec(j)),'linewidth',10)
> hold off
> % save the figure
> % print -depsc map.eps
> The code works but I am wondering why I don't get a complete picture of the diagram, it seems that it is divided into 2 parts!!!!
> Any help please.
I don't have MATLAB installed on this computer, but from the above my guess is to move the statements:
% set the figure details
before your outer loop and see what happens. I suspect the statement "figure(1) is executed twice and that is why your diagram is divided into 2 parts (as you state). But this is just a guess!