フィルターのクリア
フィルターのクリア

How to change function to get correct bifurcation behavior?

1 回表示 (過去 30 日間)
Mohannad Abboushi
Mohannad Abboushi 2017 年 6 月 6 日
編集済み: 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
How can I get the graph to bifurcate?

サインインしてこの質問に回答する。

回答 (0 件)

カテゴリ

Help Center および File Exchange で Ordinary Differential Equations についてさらに検索

タグ

Community Treasure Hunt

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

Start Hunting!

Translated by