how can I draw bifurcation diagram for given code

95 Ansichten (letzte 30 Tage)
mallela ankamma rao
mallela ankamma rao am 31 Aug. 2022
Beantwortet: Abderrahim. B am 31 Aug. 2022
Good afternoon sir
can you tell me anyone how to draw bifurcation diagram like this
N=1390000000;
zetaa = 0.4;
zetai = 0.3;
zetag = 0.3;
theta = 0.2;
omega = 0.5;
beta = 0.3311;
lamdai = 0.0046;
lamdau = 0.5129;
gammai = 0.0046;
gammau = 0.5129;
gammag = 0.1318;
mui = 0.0046;
muu = 0.4852;
mug = 0.0128;
mu = 0.0014;
pi=2000
Ki = lamdai+gammai+mui+mu; Ku = lamdau+gammau+muu+mu; Kg = gammag+mug+mu; Km = omega+mu;
R0 = (beta*zetag/Km)+(beta*zetai*theta*omega/(Ki*Km))+...
((beta*zetag*lamdai*theta*omega*Ku+beta*zetag*lamdau*(1-theta)*omega*Ki)/(Ki*Kg*Ku*Km));
A1 =(pi-mu*N/R0)/Km;
I1 =(theta*omega)/Ki;
U1=((1-theta)*omega*A1)/Ku;
G1 =(lamdai*theta*omega*Ku+lamdau*(1-theta)*omega*Ki)*A1/(Ki*Kg*Ku);
R1 =(((lamdai*theta*omega)/mu*Ki)+((1-theta)*omega/Ku*mu)+((gammag*lamdai*theta*omega)/(mu*Kg*Ki))...
+(gammag*lamdau*(1-theta)*omega)/(mu*Kg*Ku))*A1;
V1 = beta*(zetaa*A1+zetai*I1+zetag*G1)*(1/N);
plot(R0,V1)
xlabel('R0');
ylabel('V1');
box on
hold on

Antworten (1)

Abderrahim. B
Abderrahim. B am 31 Aug. 2022
Hi!
MATLAB does not have a function to plot this kind of diagrams, but it has capabilities that you can use to write your own specific codes.
Consider using function the below file echanges that you may find useful .
  1. Plotting bifurcation diagram.
  2. Chaotic generators demo
Hope this helps

Kategorien

Find more on Nonlinear Analysis in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by