how can I draw bifurcation diagram for given code

56 Ansichten (letzte 30 Tage)
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
0 Kommentare-2 ältere Kommentare anzeigen-2 ältere Kommentare ausblenden

Melden Sie sich an, um zu kommentieren.

Antworten (1)

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 .
Hope this helps
0 Kommentare-2 ältere Kommentare anzeigen-2 ältere Kommentare ausblenden

Melden Sie sich an, um zu kommentieren.

Kategorien

Mehr zu MATLAB finden Sie in Help Center und File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by