How to plot parameter vs Lyapunov coefficient
Numerical Bifurcation Analysis Toolbox in Matlab
Brought to you by:
hilmeijer,
willy_govaerts
I am studing a 2D system having bautin bifurcation. I want to plot a parameter vs 1st Lyapunov coefficient in Hopf bifurcation. How to do that? Please provide me where I can get the data of 1st Lyapunov coefficient when varying a system parameter?
Last edit: pallav 2023-12-10
You have not given all context of what you're doing, but I suppose you have computed a Hopf curve and the test function of a degenerate Hopf (Bautin) was active too. From the MatCont folder, first load the data into Workspace:
load Systems/(model name)/diagram/(name of Hopf curve eg H_H(1).mat))
Next plot the test function (which equals the first Lyapunov coefficient) with
plot(x(end,:), h(ii,end)) %typically the index ii=4, but it depends on the precise settings of what you did and did not compute.
Excuse me, but what do you mean by the Hopf curve? I can currently find the Hopf point, but I don't know what to do next.
Thank you very much sir.