[zhna123]

Dx/dt=a*(1-x/b)*x %菌体浓度模型
function dx = jun(bx,xs)
dx=bx(1)*(1-xs(:,1)/bx(2)).*xs(:,1);

Dp/dt=c*dx/dt+d*x %产物浓度模型
function dp = chan(bp,dxsx)
dp=bp(1)*dxsx(:,1)+bp(2)*dxsx(:,2);

-ds/dt=e*ds/dt-f*dp/dt %底物浓度模型
function ds = di(bs,dxxdp)
ds=-dxxdp(:,1)*bs(1)-dxxdp(:,2)*bs(2);


a b c d e f为待估参数 x p s为已知变量
:cry: