antictl.m

数据挖掘聚类算法

function dx=antictl(t,x)
global f1 f2 d1 b1 k1 k2 b3 d3 k3 g4 k4 d4 g5 k5 d5; 
dx=zeros(5,1);
dx(1)= [f1*f2/(f2^2+(x(1)+x(2))^2)]*x(1)-d1*x(1)-b1*x(1)*x(3)+k1*x(5)*x(2)%健康细胞;
dx(2)= [f1*f2/(f2^2+(x(1)+x(2))^2)]*x(2)+ b1*x(1)*x(3)-d1*x(2)-(k1+k2)*x(2)*x(5)%感染细胞; 
dx(3)=b3*x(2)-d3*x(3)-k3*x(3)*x(4)%病毒;
dx(4)=g4+k4*x(3)*x(4)-k3*x(3)*x(4)-d4*x(4)%抗体响应响应
dx(5)=g5+k5*x(2)*x(5)-d5*x(5)%CTL
%dx(1)=F(N)x(1)-d1x(1)-b1x(1)x(3)+k1x(5)x(2)%健康细胞;
%dx(2)= F(N)x(2)+ b1x(1)x(3)-d1x(2)-(k1+k2)x(2)x(5)%感染细胞;
%dx(3)=b3x(2)–k3x(3)x(4)%病毒;
%dx(4)=g4+k4x(3)x(4)–d4x(4)%抗体响应响应
%dx(5)=g5+k5x(2)x(5)–d5x(5)%CTL