
close all
clear
clc
%构建ODE方程
% DA1=k2*C-k1*A1*R+n*(A2-A1)-r1*A1; %处理细胞内部AHL的改变速率
% DR=k2*C-k1*A*R-b*R+R0;            %处理细胞内部LuxR的改变速率
% DC=k1*A*R-k2*C;                   %处理细胞内部复合物的改变速率
% DCr=a*(C^B/(km^B+C^B))+C0-kd*Cr;  %Crpp酶的改变速率
% DA2=-P*n*(A2-A1)-r2*A2;           %外部AHL的改变速率
%定义初值
% A2=1e-6;
% A1=0;
% R=0;
% C=0;
% Cr=0;
F=[];
G=0:1*10^7:5*10^10;
for i=0:1*10^7:5*10^10
eval(sprintf('[t,y]=ode15s(@eq5,[0,4000],[0;0;0;0;%d]);',i))
y=double(y);
F=[F,max(y(:,4))];
end
plot(G,F,'r') 
xlabel('Initial concentration of external AHL nmol/L')
ylabel('Maximum concentration of CrpP nmol/L')
function dy=eq5(t,y)
k1=0.1; % A+R-C过程速率常数
k2=10;  % C-A+R过程速率常数
n=0.3;  % AHL扩散系数
b=0.02; % LuxR降解速率常数
r1=0.00282; % 内部AHL降解速率常数
r2=0.0004;  % 外部AHL降解速率常数
R0=1;       % 常启动子产生LuxR的基础速率
C0=0.02;    %Crpp产生的基础速率
km=10;      %复合物对启动子的亲和力
kd=0.005;   %Crpp降解速率常数
a=5;        %结合复合物时表达效率
B=1.5;      %
P=4*10^6;      %细胞数量
% y(1)=A1;
% y(2)=R;
% y(3)=C;
% y(4)=Cr;
% y(5)=A2;
dy=zeros(5,1);
dy(1)=k2*y(3)-k1*y(1)*y(2)+n*(y(5)-y(1))-r1*y(1);
dy(2)=k2*y(3)-k1*y(1)*y(2)-b*y(2)+R0;
dy(3)=k1*y(1)*y(2)-k2*y(3);
dy(4)=a*(y(3)^B/(km^B+y(3)^B))-kd*y(4);
dy(5)=-P*n*(y(5)-y(1))-r2*y(5);
end