Plot ODEs or Bifurcation diagram

Question

0 个投票

Good evening, I am Mawan.

I am a bit troubled with matlab and I hope someone can help.

So I have a model in the form of seven-dimensional system of differential equations. I want to plot a graph between one coefficient of ODEs is the value created varies with one ODE in the system.

I include the following system:

% Parameter Values: Lh=3.285e-2; Psh=7.666e-5; Psv=0.4; bvh=0.8333; bhv=2e-2; Bvh=8.333e-2; sgh=18; sgv=0.6; vh=8.333e-2; vv=0.1; gh=3.704e-3; dh=3.454e-4; roh=1.46e-2; m1h=4.212e-5; m2h=1e-7; m1v=0.1429; m2v=2.279e-4;

%ODEs dy = zeros(7,1); % a column vector

dy(1) = (zt*bhv*y(6)*y(7)/y(4)) * (1-y(1)-y(2)-y(3)) ... - (vh + Psh + (Lh/y(4)))*y(1) ... + dh * y(2) * y(1);

dy(2) = vh*y(1) - (gh+dh+Psh+(Lh/y(4)))*y(2) ... + dh*(y(2))^2;

dy(3) = gh*y(2) - (roh+Psh+(Lh/y(4)))*y(3) ... + dh*y(2)*y(3);

dy(4) = Lh + Psh * y(4) - (m1h+m2h*y(4))*y(4) ... - dh*y(2)*y(4);

dy(5) = zt * (bvh*y(2) + Bvh*y(3)) * (1-y(5)-y(6)) ... - (vv+Psv)*y(5);

dy(6) = vv*y(5) - Psv*y(6);

dy(7) = Psv*y(7) - (m1v + m2v * y(7)) * y(7);

I want to make the plot between the value of the variable zt (7e-4: 8e-4) with y (1). how to make a plot? Thanks for the help.

Yours sincerely

Mawan