Discover MakerZone

MATLAB and Simulink resources for Arduino, LEGO, and Raspberry Pi

Learn more

Discover what MATLAB® can do for your career.

Opportunities for recent engineering grads.

Apply Today

Thread Subject:
bifurcation diagram

Subject: bifurcation diagram

From: J G

Date: 16 Sep, 2010 11:52:04

Message: 1 of 3

I have a table of data which I'd like to plot in a bifurcation diagram to illustrate equilibrium points (y-axis) and a changing parameter (x-axis). I don't think there is a function for this? Do I just use a general plot function to do this, changing my data into vectors? Thanks if you can help.

Subject: bifurcation diagram

From: Grzegorz Knor

Date: 16 Sep, 2010 12:07:04

Message: 2 of 3

"J G" <jemma.geoghegan@gmail.com> wrote in message <i6t0d4$geq$1@fred.mathworks.com>...
> I have a table of data which I'd like to plot in a bifurcation diagram to illustrate equilibrium points (y-axis) and a changing parameter (x-axis). I don't think there is a function for this? Do I just use a general plot function to do this, changing my data into vectors? Thanks if you can help.


From Wikipedia, the free encyclopedia:

r_values=2.5:.001:4;
iterations_per_value = 500;
output_matrix = zeros(iterations_per_value,length(r_values));
for count = 1 : length(r_values)
        r = r_values(count);
        y = .5;
        for i=1:500;
            y = r*(y*(1-y));
        end
 
        for i=1:iterations_per_value
            y = r*(y*(1-y));
            output_matrix(i,count) = y;
        end
end
plot(r_values,output_matrix,'.k');
grid on;

Subject: bifurcation diagram

From: J G

Date: 16 Sep, 2010 12:27:05

Message: 3 of 3

Thank you for your help! Would you mind writing the pseudocode (I generally follow it but I've just started using matlab). Thanks again.

>
>
> From Wikipedia, the free encyclopedia:
>
> r_values=2.5:.001:4;
> iterations_per_value = 500;
> output_matrix = zeros(iterations_per_value,length(r_values));
> for count = 1 : length(r_values)
> r = r_values(count);
> y = .5;
> for i=1:500;
> y = r*(y*(1-y));
> end
>
> for i=1:iterations_per_value
> y = r*(y*(1-y));
> output_matrix(i,count) = y;
> end
> end
> plot(r_values,output_matrix,'.k');
> grid on;

Tags for this Thread

No tags are associated with this thread.

What are tags?

A tag is like a keyword or category label associated with each thread. Tags make it easier for you to find threads of interest.

Anyone can tag a thread. Tags are public and visible to everyone.

Contact us