How to change function to get correct bifurcation behavior?
조회 수: 1 (최근 30일)
이전 댓글 표시
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
답변 (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!