How to change function to get correct bifurcation behavior?

1 次查看（过去 30 天）

Mohannad Abboushi 2017-6-6

0
链接

此问题的直接链接

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

编辑： Mohannad Abboushi 2017-6-6

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