I want to draw a bifurcation diagram for system of four delay differential equations

3 次查看（过去 30 天）

Akhtar Jan 2022-7-8

clc

lags = 0.6;

sol = dde23(@ddex1de,lags,@ddex1hist,[0,8]);

t = sol.x;

y = sol.y;

plot(t,y(1,:));

function dydt = ddex1de(t,y,Z)

ylag1 = Z(:,1);

k1 = 2;

k_1 = 1;

k2 = 1;

k3 = 1;

k_3 = 1;

k4 = 2;

k5 = 1;

E1 = 1;

E2 = 1;

dydt = [-k1*y(1)*E1+k_1*y(2)+k4*y(3)-k5*ylag1(1);

k1*y(1)*E1-(k_1+k2)*y(2);

k3*y(4)*E2-(k_3+k4)*y(3);

k2*y(2)+k_3*y(3)-k3*y(4)*E2+k5*ylag1(1)];

end

function S = ddex1hist(t)

S = [5;0;0;0];

end

MATLAB Programming

在 Help Center 和 File Exchange 中查找有关 Programming 的更多信息

R2022a

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

Translated by