How to change function to get correct bifurcation behavior?

1 visualización (últimos 30 días)

Mohannad Abboushi el 6 de Jun. de 2017

0
Enlazar

Enlace directo a esta pregunta

https://la.mathworks.com/matlabcentral/answers/343591-how-to-change-function-to-get-correct-bifurcation-behavior

Editada: Mohannad Abboushi el 6 de Jun. de 2017

I have a system of ODE's that was modelled and now I am trying to vary the parameters for an ODE of interest (dDSB). The equation is: dDSB:=alpha/H1*dD-k1*c0*(1-h1) where: dD=.001, H1=3000, k1=100, c0=25/3000, h1=(25+2)/3000. Solving for alpha I get 0.0000003333333333*alpha - 991/1200. I then tried to graph the plot:

clf; hold on;
for alpha=0.1:.1:25
Yfix=roots([3.3*10^7*alpha -991/1200]);
I=imag(Yfix)==0;
Yfix=Yfix(I);
Lfix=2*Yfix;
plot(repmat(alpha,size(Yfix)), Yfix, 'o', 'LineWidth',5);
end
hold off