how to plot bifurcation diagram in matlab

a=7.2*10^-6;

Gm=2.5*10^-2

Sp=2.75*10^-5;

yon=6*10^-2;

b=4.7;

son=4.5*10^-1;

B=10^-4;

vm=0:0.5:2.5;

gp=sinh(vm);

A=10^-10;

soff=1.3*10^-2;

beta=500;

yoff=1.3*10^-2

v=sawtooth(2*pi*1*t)

i=v.*(y*Gm+(1-y)*a.*exp(b*sqrt(abs(v))))

p=i.*v

dydt = B*sinh(v/son).*exp(-y.^2/yon^2).*exp(p/Sp)*heaviside(v)+...

A*sinh(v/soff).*exp(-yoff^2./y.^2).*exp(1/(1+beta*p))*heaviside(-v)

i want to bifurcation diagram for i,v

[Lazaros Moysis]

The best way to plot the bifurcation diagram is to use the local maxima method. Basically, you run a for loop, where each time you change the values of the system's parameter that you want to study. They you simulate the system, find the local maxima for a chosen state, and plot then in a graph. Once the for loop finishes, you can observe how many local maxima the system has for each parameter value. 1 local maxima means period-1 behavior, 2 local maxima means period-2 behavior and so on. Innumerable number of maxima means chaos

https://www.youtube.com/watch?v=jkrmZB9OM2E

https://www.youtube.com/watch?v=p9mg6TOBMY0

Answer 1

[Pranav Verma]

Hi Nune,

Please see the below answers on the same lines which discuss the steps to plot the bifurcation diagrams:

https://www.mathworks.com/matlabcentral/answers/319767-how-to-plot-a-bifurcation-diagram-for-differential-equation-is-there-any-formula-to-plot-the-bifurc

https://www.mathworks.com/matlabcentral/answers/20382-bifurcation-diagram

Thanks

[Atish kumar Sethy]

Hopf bifurcation for holling type problem matlab code

if avilable share the matlab code