I want to draw a bifurcation diagram for system of four delay differential equations
조회 수: 3 (최근 30일)
이전 댓글 표시
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
댓글 수: 0
답변 (0개)
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!