Finding equilibrium points for an ODE system
15 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
Hi, I have two functions named cdot and ctdot. I want to find the eqilibrium points which means cdot=ctdot=0. Could you please tell me how can I find points (c,ct) which satisfied in cdot=ctdot=0.
c and ct should be positive between [0,2].
Thanks in advance for any help.
vplc=0.16;
delta=0.1;
Ktau=0.045;
Kc=0.1;
K=0.0075;
Kp=0.15;
gamma=5.5;
kb=0.4;
vss=0.044;
alpha0=delta*6.81e-6/(0.002);
alpha1=delta*2.27e-5/(0.002);
Ke=7;
Vs=0.002;
ks=0.1;
Kf=0.18;
kplc=0.055;
ki=2;
A=(-(vss.*c.^2)./(ks.^2))+((Vs.*K.*gamma.^2.*ct.^2)./(ks.^2))+alpha0+alpha1.*((Ke.^4)./(Ke.^4+(gamma.*ct).^4));
h=(-(0.4.*A.*((Kc.^4).*(Kp.^2))./((p.^2.*c.^2.*gamma.*ct.*Kf))));
jin2=alpha1.*Kce.^4./((gamma.*ct).^4+Kce.^4);
p=(vplc.*c.^2/(c.^2+kplc.^2))./ki;
G1=alpha0+jin2;
G2=((1-h)./tau_max).*c.^4;
Fc=(4.*gamma.*Kf).*((c.^3.*p.^2.*h.*ct)./(Kb.*Kp.^2.*Ktau.^4))-(2.*Vss.*c./Ks.^2);
Fct=((gamma.*Kf.*(c.^4).*(p.^2).*h)./(Kb.*Kp.^2.*Ktau.^4))+((Vs.*K.*gamma.^2)./(Ks.^2))-((4.*gamma.^4.*ct.^3.*alpha1.*Kce.^4)./(Kce.^4+(gamma.*ct).^4).^2);
Fh=(gamma.*Kf.*c.^4.*p.^2.*ct)./(Kb.*Kp.^2.*Ktau.^4);
cdot=(Fct).*(G1)+(Fh).*(G2);
ctdot=-G1.*Fc;
3 Kommentare
Star Strider
am 21 Sep. 2022
‘... find points (c,ct) which satisfied in cdot=ctdot=0’
They don’t intersect, at least in the regions described, as we discussed in how can I find the intersection of two surface. The value of ‘C’ may need to be negative for an intersection to exist.
Siehe auch
Kategorien
Mehr zu Denoising and Compression finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!