# plot the bifurcation into multiple colors on mathematical modelling case

8 vues (au cours des 30 derniers jours)
Latifah Hanum le 7 Fév 2023
I want to plot bifurcation, this is the result I want. Parameter space Phi × u for Pr = 0.002.
I don't have any idea how to make this. I only have this script (attach)
##### 1 commentaireAfficher -1 commentaires plus anciensMasquer -1 commentaires plus anciens
Latifah Hanum le 7 Fév 2023
can you help me to change this script?

Connectez-vous pour commenter.

### Réponse acceptée

Sargondjani le 7 Fév 2023
I dont have time to look into your code, but I usually plot bifurcations using "contour". I compute the eigenvalues, and then define the fieds in the graphs by setting the boundaries of the countour to the bifurcation values. An example of my code:
contour(mm_vec,z1_vec,l1,[1 1],'-k','LineWidth',LW);
contour(mm_vec,z1_vec,l2,[1 1],'-','LineWidth',LW,'Color',[0 0.447 0.741],'HandleVisibility','off');
contour(mm_vec,z1_vec,l3,[1 1],'-b','LineWidth',LW,'HandleVisibility','off');%'Color',[0.85 0.325 0.098]
contour(mm_vec,z1_vec,l1,[-1 -1],'-.k','LineWidth',LW);
contour(mm_vec,z1_vec,l2,[-1 -1],'-.','LineWidth',LW,'Color',[0 0.447 0.741]);
contour(mm_vec,z1_vec,l3,[-1 -1],'-.b','LineWidth',LW,'HandleVisibility','off');%,'Color',[0.85 0.325 0.098]
contour(mm_vec,z1_vec,imag_size,[1e-8 1e-8],':','LineWidth',LW2,'Color',[0.85 0.325 0.098]);%,'fill','on');
In this code mm and z1 are the axis (bifurcation parameters), and l1, l2 and l3 are the eigenvalues, and imag_size is the imaginary part when there exists complex eigenvalues. (Note that these are discrete time bifurcations).
To get smooth plots in case of non-linear graphs your vectors mm and z1 have to be large.
##### 0 commentairesAfficher -2 commentaires plus anciensMasquer -2 commentaires plus anciens

Connectez-vous pour commenter.

### Catégories

En savoir plus sur Nonlinear Dynamics dans Help Center et File Exchange

R2021a

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by